Re: AMBER: Simulation of small protein

From: Carlos Simmerling <>
Date: Thu, 23 Oct 2003 12:13:09 -0400

it depends on whether you're seeing backbone changes-
we've publsihed that there is significant overstabilization of alpha
helix in parm99+GB. on a long time scale, all of
our sequences increased in helical content and the
energy dropped. Seeing this during 200ps is a bit surprising

without knowing what the topology looks like it is
hard to say whether more compact is reasonable (is it
a rather open structure to start with?)

----- Original Message -----
From: "Endres, Robert G." <>
To: <>
Sent: Thursday, October 23, 2003 11:52 AM
Subject: AMBER: Simulation of small protein

> Dear AMBER users,
> After heating/equilibration of a NMR structure of a small (36 res.)
protein with weak constraints, I was doing a short (0.2 ns) MD simulation
(parm99.dat) at 300K with the GBSA implicit solvent model.
> I was a bit surprised that the total energy dropped further by almost 100
kcal/mol compared to the starting structure (close to NMR structure), and
the ESURF term (proportional to the total surface) decreased from 17 to 12
> So the protein got more compact during the MD simulation. The RMSD started
from zero and increased to 3-5 A.
> It seems that the protein changed quite strongly the structure during the
simulation. Does anyone have experience with this, i.e. is this "normal" or
"expected", is this a problem with the force field/solvent model or NMR
> Many thanks for suggestions,
> Robert
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to
> To unsubscribe, send "unsubscribe amber" to

The AMBER Mail Reflector
To post, send mail to
To unsubscribe, send "unsubscribe amber" to
Received on Thu Oct 23 2003 - 17:53:01 PDT
Custom Search