On 4/11/14, 7:47 AM, neha_bharty wrote:
Hello
I am new user of gromacs. I am using gromacs 4.6.3
I wanted to minimize the energy of a protein. Before energy minimization the
prockeck result for Ramachandran plot was 89% in allowed region but after
energy minimization the result become worse i.e. 87%.
I was performing energy minimization to refine the model. But the model
become worse after that.

That's not a huge change, nor is energy minimization rigorous enough to fully refine a model. EM will probably optimize the side chains pretty well, but the backbone may require actual MD simulation to sample better conformations. The changes that take place during EM are typically very small.

I Model the protein using homology modeling and needed it for docking

I am also doing it with different protein molecules but again the same
problem occurs.
I don't know that its normal or there is problem with my parameter or force
field.
I am also mentioning the detail of my energy minimization.

The details are:
Force Field Used: CHARMM27 all-atom force field
Water Model: TIP3P


CHARMM36 may be a better option. The recalibration of the CMAP terms could be relevant here. See http://mackerell.umaryland.edu/CHARMM_ff_params.html to download the force field files.

-Justin

em.mdp file is:

; VARIOUS PREPROCESSING OPTIONS =
title                    =
cpp                      = /lib/cpp
include                  =
define                   = -DFLEXIBLE

; RUN CONTROL PARAMETERS =
integrator               = steep
; start time and timestep in ps =
tinit                    = 0
dt                       = 0.002
nsteps                   = 15000

; ENERGY MINIMIZATION OPTIONS =
emtol                    = 10
emstep                   = 0.1
nstcgsteep               = 1000

coulombtype     =pme
nstenergy       =10

Please help me out for this.


--
View this message in context: 
http://gromacs.5086.x6.nabble.com/Re-Energy-minimization-result-problem-tp5015793.html
Sent from the GROMACS Users Forum mailing list archive at Nabble.com.


--
==================================================

Justin A. Lemkul, Ph.D.
Ruth L. Kirschstein NRSA Postdoctoral Fellow

Department of Pharmaceutical Sciences
School of Pharmacy
Health Sciences Facility II, Room 601
University of Maryland, Baltimore
20 Penn St.
Baltimore, MD 21201

jalem...@outerbanks.umaryland.edu | (410) 706-7441
http://mackerell.umaryland.edu/~jalemkul

==================================================
--
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.

Reply via email to