A single number, i.e. RMSD, is no substitute for actually visualizing your
system. Also keep in mind that 30.0 and even 5.0 are very strong restraint
force constants, and the system will hardly move at all during these steps.
Can you see what is happening during the trajectory? Also, did you use the
same exact protocol when you did the constant protonation simulations? Did
you use the same force constants? It seems strange that the RMSD for the
constant protonation was larger in the heating/equilibration for constant
protonation than it was for constant pH (I think you must have used
different restraints).
Finally, what frame are you aligning to? Both should be aligned to the same
structure to get a legitimate comparison.
Good luck!
Jason
On Tue, Jul 27, 2010 at 6:15 AM, Sangita Kachhap <sangita.imtech.res.in>wrote:
> I have done the constant pH 4 simulation of lysozyme.
> 1) I minimized the structure using restraint on backbone atom, restraint_wt
> = 30.0
> for 3000 steps.
> 2) I heat the system gradualy for 1K/1ps using restraint on backbone atom,
> restraint_wt = 5.0
> 3) Equilibration for 200ps using restraint on backbone atom, restraint_wt =
> 1.0
> 4) Production run for 10ns without restraint.
>
> I also did the whole thing using icnstph=0 (constant protonation)to compare
> both
> result
>
> I processed the trajectory using ptraj and got the following result:
>
> Heating
> maximum rmsd in const pH MD - 0.98
> maximum rmsd in const protonation - 2.1
>
> Equilibration
> maximum rmsd in const pH MD - 1.15
> maximum rmsd in const protonation - 3.8
>
> Production
> maximum rmsd in const pH MD - 5.5 (in 5ns)
> maximum rmsd in const pH MD - 10.9 (in next 5ns)
> maximum rmsd in const protonation - 5 (in 5ns)
> maximum rmsd in const protonation - 6 (in next 5ns)
>
> In constant pH simulation rmsd in not stabilizing but in case of constant
> protonation it is slightly stable.
>
> why rmsd is not stabilizing upto 10ns production run can anyone help me.
>
> Thanks With regard
> Sangita Kachhap
> JRF
> BIC,IMTECH
> CHANDIGARH
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
Jason M. Swails
Quantum Theory Project,
University of Florida
Ph.D. Graduate Student
352-392-4032
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Tue Jul 27 2010 - 09:30:04 PDT
