[gmx-users] protein stability as a dimer

Justin A. Lemkul jalemkul at vt.edu
Thu Jul 1 13:47:14 CEST 2010

sonali dhindwal wrote:
> Hello All,
> I have done simulation for 1ns on a protein dimer using GROMOS96 43a1 
> force field.
> I want to study if the protein is stable as  a dimer or not, so can you 
> please give me some suggestion as to what analysis  I could do for the same.
> I checked g_rms after the simulation, graph which I am getting is like, 
> first rmsd is increases rapidly to 0.2nm for for 0.1 ns and
> then till 0.8 ns it still increased to 0.3 nm and
> still there are more fluctuations but and it decreses again to 0.28 nm 
> at 0.1 ns.
> So what could i infer about the stability of the protein as a dimer from 
> this data.

If you've only done 100 ps of simulation, not much.  I would also question why 
you consider it a "decrease" from 0.3 nm to 0.28 nm rather than just a 
fluctuation.  Run a longer simulation to gather meaningful RMSD values.

As for other analysis to do, do some homework in the literature, and also 
consider your own goals for running the simulation.  Certainly you have a purpose.


> Thanks in advance.
> Regards
> --
> Sonali Dhindwal


Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080


More information about the gromacs.org_gmx-users mailing list