Methods for Classical-Mechanical Molecular Simulation in Chemistry: Achievements, Limitations, Perspectives.

More than a half century ago it became feasible to simulate, using classical-mechanical equations of motion, the dynamics of molecular systems on a computer. Since then classical-physical molecular simulation has become an integral part of chemical research. It is widely applied in a variety of branches of chemistry and has significantly contributed to the development of chemical knowledge.

Molecular structure refinement based on residual dipolar couplings using magnetic-field rotational sampling.

A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using magnetic-field rotational sampling of the rotational degrees of freedom of a molecule in conjunction with molecule-internal configurational sampling. By applying rotational sampling, as is occurring in the experiment, leading to observable RDCs, the method stays close to the experiment.

A Method to Derive Structural Information on Molecules from Residual Dipolar Coupling NMR Data.

A method for structure refinement of molecules based on residual dipolar coupling (RDC) data is proposed. It calculates RDC values using rotational and molecule-internal configurational sampling instead of the common refinement procedure that is based on the approximation of the nonuniform rotational distribution of the molecule by a single alignment tensor representing the average nonuniformity of this distribution. Using rotational sampling, as is occurring in the experiment leading to observable RDCs, the method stays close to the experiment.

Molecular dynamics simulation or structure refinement of proteins: are solvent molecules required? A case study using hen lysozyme.

In protein simulation or structure refinement based on values of observable quantities measured in (aqueous) solution, solvent (water) molecules may be explicitly treated, omitted, or represented by a potential of mean-solvation-force term, depending on protein coordinates only, in the force field used. These three approaches are compared for hen egg white lysozyme (HEWL). This 129-residue non-spherical protein contains a variety of secondary-structure elements, and ample experimental data are available: 1630 atom-atom Nuclear Overhauser Enhancement (NOE) upper distance bounds, 213  J-couplings and 200 S order parameters.

On the use of intra-molecular distance and angle constraints to lengthen the time step in molecular and stochastic dynamics simulations of proteins.

Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. An example of such a time-saving technique is the application of constraints to particular degrees of freedom when integrating Newton's or Langevin's equations of motion in molecular dynamics (MD) or stochastic dynamics (SD) simulations, respectively. The application of bond-length constraints is standard practice in protein simulations and allows for a lengthening of the time step by a factor of three.

On the use of multiple-time-step algorithms to save computing effort in molecular dynamics simulations of proteins.

Computer simulation of proteins in aqueous solution at the atomic level of resolution is still limited in time span and system size due to limited computing power available and thus employs a variety of time-saving techniques that trade some accuracy against computational effort. Examples of such time-saving techniques are the application of constraints to particular degrees of freedom or the use of a multiple-time-step (MTS) algorithm distinguishing between particular forces when integrating Newton's equations of motion. The application of two types of MTS algorithms to bond-stretching forces versus the remaining forces in molecular dynamics (MD) simulations of a protein in aqueous solution or of liquid water is investigated and the results in terms of total energy conservation and the influence on various other properties are compared to those of MD simulations of the same systems using bond-length, and for water bond-angle, constraints.

On the Effect of the Various Assumptions and Approximations used in Molecular Simulations on the Properties of Bio-Molecular Systems: Overview and Perspective on Issues.

Computer simulations of molecular systems enable structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic or supra-molecular level and plays an increasingly important role in chemistry, biology and physics. To interpret the results of such simulations appropriately, the degree of uncertainty and potential errors affecting the calculated properties must be considered. Uncertainty and errors arise from (1) assumptions underlying the molecular model, force field and simulation algorithms, (2) approximations implicit in the interatomic interaction function (force field), or when integrating the equations of motion, (3) the chosen values of the parameters that determine the accuracy of the approximations used, and (4) the nature of the system and the property of interest.

A method to apply bond-angle constraints in molecular dynamics simulations.

An algorithm to apply bond-angle constraints in molecular dynamics simulations of macromolecules or molecular liquids is presented. It uses Cartesian coordinates and determines the Lagrange multipliers required for maintaining the constraints iteratively. It constitutes an alternative to the use of only distance constraints (DCs) between particles to maintain a particular geometry.

On the Use of Side-Chain NMR Relaxation Data to Derive Structural and Dynamical Information on Proteins: A Case Study Using Hen Lysozyme.

Values of and order parameters derived from NMR relaxation measurements on proteins cannot be used straightforwardly to determine protein structure because they cannot be related to a single protein structure, but are defined in terms of an average over a conformational ensemble. Molecular dynamics simulation can generate a conformational ensemble and thus can be used to restrain and order parameters towards experimentally derived target values (exp) and (exp). Application of and order-parameter restraining MD simulation to bond vectors in 63 side chains of the protein hen egg white lysozyme using 51 (exp) target values and 28 (exp) target values shows that a conformational ensemble compatible with the experimentally derived data can be obtained by using this technique.

Algorithms to apply dihedral-angle constraints in molecular or stochastic dynamics simulations.

Various algorithms to apply dihedral-angle constraints in molecular dynamics or stochastic dynamics simulations of molecular systems are presented, investigated, and tested. They use Cartesian coordinates and determine the Lagrangian multipliers necessary for maintaining the constraints iteratively. The most suitable algorithm to maintain a dihedral-angle constraint is numerically compared to the alternative to use distance constraints to this end.

Conformational Properties of the Chemotherapeutic Drug Analogue Epothilone A: How to Model a Flexible Protein Ligand Using Scarcely Available Experimental Data.

Epothilones are among the most potent chemotherapeutic drugs used for the treatment of cancer. Epothilone A (EpoA), a natural product, is a macrocyclic molecule containing 34 non-hydrogen atoms and a thiazole side chain. NMR studies of EpoA in aqueous solution, unbound as well as bound to αβ-tubulin, and unbound in dimethyl sulfoxide (DMSO) solution have delivered sets of nuclear Overhauser effect (NOE) atom-atom distance bounds, but no structures based on NMR data are present in structural data banks.

Validation of Molecular Simulation: An Overview of Issues.

Computer simulation of molecular systems enables structure-energy-function relationships of molecular processes to be described at the sub-atomic, atomic, supra-atomic, or supra-molecular level. To interpret results of such simulations appropriately, the quality of the calculated properties must be evaluated. This depends on the way the simulations are performed and on the way they are validated by comparison to values Q of experimentally observable quantities Q.

Using Complementary NMR Data Sets To Detect Inconsistencies and Model Flaws in the Structure Determination of Human Interleukin-4.

The derivation of protein structure from values of observable quantities measured in NMR experiments is a rather nontrivial task due to (i) the limited number of data compared to degrees of freedom of a protein, (ii) the uncertainty inherent to the function connecting an observable quantity to molecular structure, (iii) the finite quality of biomolecular models and force fields used in structure refinement, and (iv) the conformational freedom of a protein in aqueous solution, which requires extensive conformational sampling and appropriate conformational averaging when calculating or restraining to sets of NMR data. The protein interleukin-4 (IL-4) has been taken as a test case using NOE distances, S order parameters, and J-couplings as test data and the former two types of data as restraints. It is shown that, by combining sets of different, complementary NMR data as restraints in MD simulations, inconsistencies in the data or flaws in the model and procedures used to derive protein structure from NMR data can be detected.

Interpretation of Seemingly Contradictory Data: Low NMR S Order Parameters Observed in Helices and High NMR S Order Parameters in Disordered Loops of the Protein hGH at Low pH.

At low pH, human growth hormone (hGH) adopts a partially folded state, in which the native helices are maintained, but the long loop regions and side-chain packing become disordered. Some of the S order parameters for backbone N-H vectors derived from NMR relaxation measurements on hGH at low pH initially seem contradictory. Three isolated residues (15, 20, and 171) in helices A and D exhibit low order parameter values (<0.

Deriving Structural Information from Experimentally Measured Data on Biomolecules.

During the past half century, the number and accuracy of experimental techniques that can deliver values of observables for biomolecular systems have been steadily increasing. The conversion of a measured value Q of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: 1) insufficient or inaccurate values of Q , 2) inaccuracies in the function Q(r→) used to relate the quantity Q to structure r→ , 3) how to account for the averaging inherent in the measurement of Q , 4) how to handle the possible multiple-valuedness of the inverse r→(Q) of the function Q(r→) , to mention a few. These apply to a variety of observable quantities Q and measurement techniques such as X-ray and neutron diffraction, small-angle and wide-angle X-ray scattering, free-electron laser imaging, cryo-electron microscopy, nuclear magnetic resonance, electron paramagnetic resonance, infrared and Raman spectroscopy, circular dichroism, Förster resonance energy transfer, atomic force microscopy and ion-mobility mass spectrometry.

On the use of time-averaging restraints when deriving biomolecular structure from ³ -coupling values obtained from NMR experiments.

Deriving molecular structure from [Formula: see text]-couplings obtained from NMR experiments is a challenge due to (1) the uncertainty in the Karplus relation [Formula: see text] connecting a [Formula: see text]-coupling value to a torsional angle [Formula: see text], (2) the need to account for the averaging inherent to the measurement of [Formula: see text]-couplings, and (3) the sampling road blocks that may emerge due to the multiple-valuedness of the inverse function [Formula: see text] of the function [Formula: see text]. Ways to properly handle these issues in structure refinement of biomolecules are discussed and illustrated using the protein hen egg white lysozyme as example.

A molecular dynamics simulation investigation of the relative stability of the cyclic peptide octreotide and its deprotonated and its (CF)-Trp substituted analogs in different solvents.

The cyclic octa-peptide octreotide and its derivatives are used as diagnostics and therapeutics in relation to particular types of cancers. This led to investigations of their conformational properties using spectroscopic, NMR and CD, methods. A CF-substituted derivative, that was designed to stabilize the dominant octreotide conformer responsible for receptor binding, turned out to have a lower affinity.

Investigation of the structural preference and flexibility of the loop residues in amyloid fibrils of the HET-s prion.

The structural variability of a 16-residue loop (residues 246-261) which is in part disordered and connects two layers of the β-solenoid formed by the prion-form of HET-s and its prion domain HET-s(218-289) is investigated using molecular dynamics computer simulation. A system of three HET-s(218-289) molecules in a β-sheet structure as in the fibril is simulated in aqueous solution. The trajectory structures appear to be consistent with the Cα chemical shift data obtained.

On the pitfalls of peer review.

The review process of academic, scientific research and its basic tenets is considered, thereby distinguishing between (i) reviewing of manuscripts to be published in the scientific literature, (ii) reviewing of research proposals to be financed by funding agencies, (iii) reviewing of educational or research institutions with respect to their proper functioning, and (iv) reviewing of scientists with the aim of appointing or tenuring faculty.

Supra-Atomic Coarse-Grained GROMOS Force Field for Aliphatic Hydrocarbons in the Liquid Phase.

A supra-atomic coarse-grained (CG) force field for liquid n-alkanes is presented. The model was calibrated using experimental thermodynamic data and structural as well as energetic properties for 14 n-alkanes as obtained from atomistic fine-grained (FG) simulations of the corresponding hydrocarbons using the GROMOS 45A3 biomolecular force field. A variation of the nonbonded force-field parameters obtained from mapping the FG interactions onto the CG degrees of freedom to fit the density and heat of vaporization to experimental values turned out to be mandatory for a correct reproduction of these data by the CG model, while the bonded force-field parameters for the CG model could be obtained from a Boltzmann-weighted fit with some variations with respect to the corresponding properties from the FG simulations mapped onto the CG degrees of freedom.

Effects of Polarizable Solvent Models upon the Relative Stability of an α-Helical and a β-Hairpin Structure of an Alanine Decapeptide.

The free enthalpy of changing a nonpolarizable solvent into a polarizable solvent is calculated for an alanine decapeptide solvated in water, methanol, chloroform, or carbon tetrachloride. Introducing polarizability into a water solvent does not change the relative stability between the α-helix and the β-hairpin, while for methanol and chloroform the α-helix is stabilized by about 1 kJ mol(-1) per residue and for carbon tetrachloride by about 2 kJ mol(-1) per residue. These results suggest that the less polar the solvent is, the more the α-helical structure is stabilized with respect to the β-hairpin structure by the use of a polarizable solvent model instead of a nonpolarizable one.