[gmx-users] protein stability as a dimer

[Justin A. Lemkul]

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.

-Justin

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

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

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================