Long timestep dynamics of peptides by the dynamics driver approach

Proteins: Structure, Function and Bioinformatics - Tập 21 Số 4 - Trang 282-302 - 1995
Philippe Derreumaux1,2, Tamar Schlick1,2,3
1Chemistry Department, New York University, New York, New York 10012
2Courant Institute of Mathematical Sciences, New York University, New York, New York, 10012
3The Howard Hughes Medical Institute, New York University, New York, New York 10012

Tóm tắt

AbstractPrevious experience with the Langevin/implicit‐Euler scheme for dynamics (“LI”) on model systems (butane, water) has shown that LI is numerically stable for timesteps in the 5–20 fs range but quenches high‐frequency modes. To explore applications to polypeptides, we apply LI to model systems (several dipeptides, a tetrapeptide, and a 13‐residue oligoalanine) and also develop a new dynamics driver approach (“DA”). The DA scheme, based on LI, addresses the important issue of proper sampling, which is unlikely to be solved by small‐time step integration methods or implicit methods with intrinsic damping at room temperature, such as LI. Equilibrium averages, time‐dependent molecular properties, and sampling trends at room temperature are reported for both LI and DA dynamics simulations, which are then compared to those generated by a standard explicit discretization of the Langevin equation with a 1 fs timestep. We find that LI's quenching effects are severe on both the fast and slow (due to vibrational coupling) frequency modes of all‐atom polypeptides and lead to more restricted dynamics at moderate timesteps (40 fs). The DA approach empirically counteracts these damping effects by adding random atomic perturbations to the coordinates at each step (before the minimization of a dynamics function). By restricting the energetic fluctuations and controlling the kinetic energy, we are able with a 60 fs timestep to generate continuous trajectories that sample more of the relevant conformational space and also reproduce reasonably Boltzmann statistics. Although the timescale for transition may be accelerated by the DA approach, the transitional. information obtained for the alanine dipeptide and the tetrapeptide is consistent with that obtained by several other theoretical approaches that focus specifically on the determination of pathways. While the trajectory for oligoalanine by the explicit scheme over the nanosecond timeframe remains in the vicinity of the full αR‐helix starting structure, and a high‐temperature (6000°K) MD trajectory departs slowly from the a helical structure, the DA‐generated trajectory for the same CPU time exhibits unfolding and refolding and reveals a range of conformations with an intermediate helix content. Significantly, this range of states is more consistent with spectroscopic experiments on small peptides, as well as the cooperative two‐state model for helix–coil transition. The good, near‐Boltzmann statistics reported for the smaller systems above, in combination with the interesting oligoalanine results, suggest that DA is a promising tool for efficiently exploring conformational spaces of biomolecules and exploring folding/unfolding processes of polypeptides. © 1995 Wiley‐Liss, Inc.

Từ khóa


Tài liệu tham khảo

van Gunsteren W. F., 1989, Computer Simulations of Biomolecular Systems

10.1080/07391102.1989.10508503

Kollman P. A., 1989, Decisions in force field development: An alternative to those described by Roterman et al, J. Biomol. Struct. Dyn., 7, 421

10.1080/07391102.1991.10507872

Gibson K. D., 1991, Decisions in force field development: Reply to Kollman and Dill, J. Biomol. Struct. Dyn., 8, 1103

10.1080/07391102.1991.10507873

10.1080/08927029308022162

10.1002/jcc.540150207

10.1016/0166-1280(93)87151-3

10.1017/CBO9781139167864

10.1002/jcc.540141011

10.1016/0021-9991(77)90098-5

10.1063/1.460004

10.1063/1.465633

10.1126/science.220.4598.671

10.1002/bip.360240308

10.1103/PhysRevA.45.8894

Duane S. Kennedy A. D. Pendleto B. J. Roweth D.Hybrid Monte Carlo. Phys. Lett. B195–216 1987.

10.1063/1.465904

10.1016/0301-4622(92)85023-W

10.1073/pnas.83.9.2782

10.1073/pnas.84.19.6611

10.1073/pnas.76.8.3585

10.1002/bip.360330803

10.1002/cpa.3160420706

10.1002/cpa.3160420807

10.1063/1.459935

10.1016/0009-2614(92)85028-9

10.1016/0022-2836(92)90263-J

10.1021/bi00088a014

10.1002/jcc.540040211

10.1021/ja00333a021

10.1021/ja00016a010

10.1021/j100265a018

10.1063/1.458491

10.1021/ja00502a001

10.1021/j100128a027

10.1002/bip.360320508

10.1021/ja00294a005

10.1002/prot.340030408

10.1021/j100188a054

10.1016/0022-2836(92)90264-K

10.1021/ja00077a066

10.1002/bip.360340109

10.1016/0009-2614(84)80098-6

Press W. H., 1986, Numerical Recipes

10.1080/00268978800101881

10.1145/128745.150973

10.1145/128745.150975

10.1002/jcc.540150506

10.1063/1.1699114

10.1021/bi00139a031

10.1021/cr00023a008

10.1021/ma00234a015

10.1021/ja00291a014

10.1021/bi00238a033

10.1002/bip.360310911

10.1016/S0022-2836(05)80172-X

10.1126/science.2237416

10.1021/bi00230a009

10.1002/bip.360311202

10.1063/1.1730390

10.1063/1.1731802

10.1063/1.467422

10.1126/science.2402636

Derreumaux P. Karplus M. Schlick T.Simulations of the loop opening/closing motion of triose phosphate isomerase by the dynamics driver approach. In preparation.

10.1107/S0021889891004399

Thông tin
Thông tin xuất bản

Proteins: Structure, Function and Bioinformatics

Tập 21 Số 4

282-302

Thông tin tác giả