This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.105.071522v1 90/4/1191    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Nutt, D. R. Articles by Meuwly, M. Search for Related Content PubMed PubMed Citation Articles by Nutt, D. R. Articles by Meuwly, M.

Studying Reactive Processes with Classical Dynamics: Rebinding Dynamics in MbNO

Department of Chemistry, University of Basel, 4056 Basel, Switzerland

Correspondence: Address reprint requests to Markus Meuwly, E-mail: m.meuwly{at}unibas.ch.

A new surface-crossing algorithm suitable for describing bond-breaking and bond-forming processes in molecular dynamics simulations is presented. The method is formulated for two intersecting potential energy manifolds which dissociate to different adiabatic states. During simulations, crossings are detected by monitoring an energy criterion. If fulfilled, the two manifolds are mixed over a finite number of time steps, after which the system is propagated on the second adiabat and the crossing is carried out with probability one. The algorithm is extensively tested (almost 0.5 µs of total simulation time) for the rebinding of NO to myoglobin. The unbound surface (Fe···NO) is represented using a standard force field, whereas the bound surface (Fe–NO) is described by an ab initio potential energy surface. The rebinding is found to be nonexponential in time, in agreement with experimental studies, and can be described using two time constants. Depending on the asymptotic energy separation between the manifolds, the short rebinding timescale is between 1 and 9 ps, whereas the longer timescale is about an order of magnitude larger. NO molecules which do not rebind within 1 ns are typically found in the Xenon-4 pocket, indicating the high affinity of NO to this region in the protein.