We present an algorithm for simulating the long time scale dynamics of
proteins and other macromolecules, Our method applies the concept of
multiple time step integration to the diffusive Langevin equation, in
which short time scale dynamics are replaced by friction and noise, Th
e macromolecular force held is represented at atomic resolution, Slow
motions are modeled by constrained Langevin dynamics with very large t
ime steps, while faster degrees of freedom are kept in local thermal e
quilibrium, In the limit of a sufficiently large molecule, our algorit
hm is shown to reduce the CPU time required by two orders of magnitude
, We test the algorithm on two systems, alanine dipeptide and bovine p
ancreatic trypsin inhibitor (BPTI), and find that it accurately calcul
ates a variety of equilibrium and dynamical properties, In the case of
BPTI, the CPU time required is reduced by nearly a factor of 60 compa
red to a conventional, unconstrained Langevin simulation using the sam
e force held. (C) 1998 Wiley-Liss, Inc.