Deterministic global optimization and torsion angle dynamics for molecular structure prediction

Determination, prediction, and understanding of structures, using the energy landscapes of chemical systems – Part II

In the past decade, new theoretical approaches have been developed to determine, predict and understand the struc-ture of chemical compounds. The central element of these methods has been the investigation of the energy landscape of chemical systems. Applications range from extended crystalline and amorphous compounds over clusters and molecular crystals to proteins. In this review, we are going to give an introduction to energy landscapes and methods for their investigation, together with a number of examples. These include structure prediction of extended and mo-lecular crystals, structure prediction and folding of proteins, structure analysis of zeolites, and structure determination of crystals from powder diffraction data.

A comparative study of molecular dynamics in Cartesian and in internal coordinates: dynamical instability in the latter caused by nonlinearity of the equations of motion

The stability of a general molecular dynamics (MD) integration scheme is examined for simulations in generalized (internal plus external) coordinates (GCs). An analytic expression is derived for the local error in energy during each integration time step. This shows that the explicit dependence of the mass-matrix on GCs, which makes the system's Lagrange equations of motion nonlinear, causes MD simulations in GCs to be less stable than those in Cartesian coordinates (CCs). In terms of CCs, the corresponding mass-matrix depends only on atomic masses and thus atomistic motion is subject to the linear Newton equations, which makes the system more stable. Also investigated are two MD methods in GCs that utilize nonzero elements of the vibrational spectroscopic B-matrices. One updates positions and velocities in GCs that are iteratively adjusted so as to conform to the velocity Verlet equivalent in GCs. The other updates positions in GCs and velocities in CCs that are adjusted to satisfy the internal constraints of the new constrained WIGGLE MD scheme. The proposed methods are applied to an isolated n-octane molecule and their performances are compared with those of several CCMD schemes. The simulation results are found to be consistent with the analytic stability analysis. Finally, a method is presented for computing nonzero elements of B-matrices for external rotations without imposing the Casimir-Eckart conditions.

Application of torsion angle molecular dynamics for efficient sampling of protein conformations

We investigate the application of torsion angle molecular dynamics (TAMD) to augment conformational sampling of peptides and proteins. Interesting conformational changes in proteins mainly involve torsional degrees of freedom. Carrying out molecular dynamics in torsion space does not only explicitly sample the most relevant degrees of freedom, but also allows larger integration time steps with elimination of the bond and angle degrees of freedom. However, the covalent geometry needs to be fixed during internal coordinate dynamics, which can introduce severe distortions to the underlying potential surface in the extensively parameterized modern Cartesian-based protein force fields. A "projection" approach (Katritch et al. J Comput Chem 2003, 24, 254-265) is extended to construct an accurate internal coordinate force field (ICFF) from a source Cartesian force field. Torsion crossterm corrections constructed from local molecular fragments, together with softened van der Waals and electrostatic interactions, are used to recover the potential surface and incorporate implicit bond and angle flexibility. MD simulations of dipeptide models demonstrate that full flexibility in both the backbone phi/psi and side chain chi1 angles are virtually restored. The efficacy of TAMD in enhancing conformational sampling is then further examined by folding simulations of small peptides and refinement experiments of protein NMR structures. The results show that an increase of several fold in conformational sampling efficiency can be reliably achieved. The current study also reveals some complicated intrinsic properties of internal coordinate dynamics, beyond energy conservation, that can limit the maximum size of the integration time step and thus the achievable gain in sampling efficiency.