Molecular Structure Refinement Based on Residual Dipolar Couplings: A Comparison of the Molecular Rotational-Sampling Method with the Alignment-Tensor Approach.

In NMR experiments, residual dipolar couplings (RDCs) in a molecule can be measured by averaging the dipolar couplings (DCs) over the rotational motion of a molecule in an environment that induces a slight anisotropic orientation distribution of the molecule. Since the shape of the anisotropic distribution cannot be measured, it is standard practice to use a particular orientation distribution of the molecule with respect to the magnetic field, in the form of a so-called alignment tensor (AT), to calculate RDC-values for the molecule. Since the same alignment tensor is commonly used to calculate the different RDCs of a molecule, this approach rests on the assumption that the rotational motion of the molecule is decoupled from its internal motions and that the molecule is rigid. The validity of these two assumptions is investigated for a small, simple molecule, using a relatively rigid atomic interaction function or force field and a more flexible one. By simulating the molecule using an orientation-biasing force an anisotropic rotational distribution can be generated, for which RDCs can be obtained. Using these RDCs as target RDCs when applying one of the two approaches of structure refinement based on RDCs, it can be investigated how well the target RDCs are approximated in the RDC restraining and whether the corresponding nonuniform orientation distribution is reproduced. For the relatively rigid version of the molecule, the AT approach reproduces the target RDC-values, although the nonuniform orientation distribution of the angle θab,H between the vector r⃗ab connecting two atoms a and b in the molecule and the vector representing the direction of the magnetic field H⃗ as generated in the orientation-biasing simulation cannot be reproduced in the AT RDC-restraining simulation. For the relatively flexible version of the molecule, the AT approach fails to reproduce both the target RDC values and the nonuniform orientation distribution. For biomolecules with flexible parts, the application of the AT approach is thus not recommended. Instead, a method based on sampling of the rotational and internal degrees of freedom of the molecule should be applied in molecular structure determination or refinement based on measured RDCs.

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. It avoids the use of an alignment tensor and, therefore, the assumptions that the overall rotation of the molecule is decoupled from its internal motions and that the molecule is rigid. Two simple molecules, a relatively rigid and a very flexible cyclo-octane molecule with eight aliphatic side chains containing 24 united atoms, serve as so-called "toy model" test systems. The method demonstrates the influence of molecular flexibility, force-field dominance, and the number of RDC restraints available on the outcome of structure refinement based on RDCs. Magnetic-field rotational sampling is basically equivalent but more efficient than explicitly sampling the rotational degrees of freedom of the molecule. In addition, the performance of the method is less dependent on the number NRDC of measured RDC-values available. The restraining forces bias the overall orientation distribution of the molecule correctly. This study suggests that the information content of RDCs with respect to molecular structure is limited.

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. Applying recently proposed algorithms to constrain bond angles or dihedral angles, it is investigated, using the protein trypsin inhibitor as test molecule, whether bond angles and dihedral angles involving hydrogen atoms or even stiff proper (torsional) dihedral angles as well as improper ones (maintaining particular tetrahedral or planar geometries) may be constrained without generating too many artificial side effects. Constraining the relative positions of the hydrogen atoms in the protein allows for a lengthening of the time step by a factor of two. Additionally constraining the improper dihedral angles and the stiff proper (torsional) dihedral angles in the protein does not allow for an increase of the MD or SD time step.

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 3 J-couplings and 200 S2 order parameters. These data are used to compare the performance of the three approaches. It is found that a molecular dynamics (MD) simulation in explicit water approximates the experimental data much better than stochastic dynamics (SD) simulation in vacuo without or with a solvent-accessible-surface-area (SASA) implicit-solvation term added to the force field. This is due to the missing energetic and entropic contributions and hydrogen-bonding capacities of the water molecules and the missing dielectric screening effect of this high-permittivity solvent. Omission of explicit water molecules leads to compaction of the protein, an increased internal strain, distortion of exposed loop and turn regions and excessive intra-protein hydrogen bonding. As a consequence, the conformation and dynamics of groups on the surface of the protein, which may play a key role in protein-protein interactions or ligand or substrate binding, may be incorrectly modelled. It is thus recommended to include water molecules explicitly in structure refinement of proteins in aqueous solution based on nuclear magnetic resonance (NMR) or other experimentally measured data.

On the use of 3J-coupling NMR data to derive structural information on proteins.

Values of 3J-couplings as obtained from NMR experiments on proteins cannot easily be used to determine protein structure due to the difficulty of accounting for the high sensitivity of intermediate 3J-coupling values (4-8 Hz) to the averaging period that must cover the conformational variability of the torsional angle related to the 3J-coupling, and due to the difficulty of handling the multiple-valued character of the inverse Karplus relation between torsional angle and 3J-coupling. Both problems can be solved by using 3J-coupling time-averaging local-elevation restraining MD simulation. Application to the protein hen egg white lysozyme using 213 backbone and side-chain 3J-coupling restraints shows that a conformational ensemble compatible with the experimental data can be obtained using this technique, and that accounting for averaging and the ability of the algorithm to escape from local minima for the torsional angle induced by the Karplus relation, are essential for a comprehensive use of 3J-coupling data in protein structure determination.

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. At comparable computational effort, the use of bond-length constraints in proteins leads to better energy conservation and less distorted properties than the two MTS algorithms investigated.

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. DCs are unsuitable to maintain particular, for example, linear or flat, geometries of molecules. The proposed algorithm can easily handle bond-length, bond-angle, and dihedral-angle constraints simultaneously, as when calculating a potential of mean force along a dihedral-angle degree of freedom.

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 S2CH and S2NH 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 S2CH and S2NH order parameters towards experimentally derived target values S2CH (exp) and S2NH (exp). Application of S2CH and S2NH order-parameter restraining MD simulation to bond vectors in 63 side chains of the protein hen egg white lysozyme using 51 S2CH (exp) target values and 28 S2NH (exp) target values shows that a conformational ensemble compatible with the experimentally derived data can be obtained by using this technique. It is observed that S2CH order-parameter restraining of C-H bonds in methyl groups is less reliable than S2NH order-parameter restraining because of the possibly less valid assumptions and approximations used to derive experimental S2CH (exp) values from NMR relaxation measurements and the necessity to adopt the assumption of uniform rotational motion of methyl C-H bonds around their symmetry axis and of the independence of these motions from each other. The restrained simulations demonstrate that side chains on the protein surface are highly dynamic. Any hydrogen bonds they form and that appear in any of four different crystal structures, are fluctuating with short lifetimes in solution.

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. In this overview, advantages and shortcomings of assumptions and approximations commonly used when simulating bio-molecular systems are considered. What the developers of bio-molecular force fields and simulation software can do to facilitate and broaden research involving bio-molecular simulations is also discussed.

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. It can easily be used to obtain a potential of mean force along a dihedral-angle coordinate.

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. X-ray diffraction of crystals has provided structures of EpoA unbound and bound to αß-tubulin. Since both crystal structures derived from X-ray diffraction intensities do not completely satisfy the three available sets of NOE distance bounds for EpoA, molecular dynamics (MD) simulations have been employed to obtain conformational ensembles in aqueous and in DMSO solution that are compatible with the respective NOE data. It was found that EpoA displays a larger conformational variability in DMSO than in water and the two conformational ensembles show little overlap. Yet, they both provide conformational scaffolds that are energetically accessible at physiological temperature and pressure.

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 Qexp of experimentally observable quantities Q. One must consider 1) the accuracy of Qexp , 2) the accuracy of the function Q(rN ) used to calculate a Q-value based on a molecular configuration rN of N particles, 3) the sensitivity of the function Q(rN ) to the configuration rN , 4) the relative time scales of the simulation and experiment, 5) the degree to which the calculated and experimental properties are equivalent, and 6) the degree to which the system simulated matches the experimental conditions. Experimental data is limited in scope and generally corresponds to averages over both time and space. A critical analysis of the various factors influencing the apparent degree of (dis)agreement between simulations and experiment is presented and illustrated using examples from the literature. What can be done to enhance the validation of molecular simulation is also discussed.

Interpretation of Seemingly Contradictory Data: Low NMR S2 Order Parameters Observed in Helices and High NMR S2 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 S2 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.5) indicating flexibility, whereas residue 143 in the centre of a long flexible loop region has a high order parameter (0.82). Using S2 order parameter restraining MD simulations, this paradox has been resolved. Low S2 values in helices are due to the presence of a mixture of 310 -helical and α-helical hydrogen bonds. High S2 values in relatively disordered parts of a protein may be due to fluctuating networks of hydrogen bonds between the backbone and the side chains, which restrict the motion of N-H bond vectors.

The key to predicting the stability of protein mutants lies in an accurate description and proper configurational sampling of the folded and denatured states.

BACKGROUND: The contribution of particular hydrogen bonds to the stability of a protein fold can be investigated experimentally as well as computationally by the construction of protein mutants which lack particular hydrogen-bond donors or acceptors with a subsequent determination of their structural stability. However, the comparison of experimental data with computational results is not straightforward. One of the difficulties is related to the representation of the unfolded state conformation. METHODS: A series of molecular dynamics simulations of the 34-residue WW domain of protein Pin1 and 20 amide-to-ester mutants started from the X-ray crystal structure and the NMR solution structure are analysed in terms of backbone-backbone hydrogen bonding and differences in free enthalpies of folding in order to provide a structural interpretation of the experimental data available. RESULTS: The contribution of the different ß-sheet hydrogen bonds to the relative stability of the mutants with respect to wild type cannot be directly inferred from experimental thermal denaturation temperatures or free enthalpies of chaotrope denaturation for the different mutants, because some ß-sheet hydrogen bonds show sizeable variation in occurrence between the different mutants. CONCLUSIONS: A proper representation of unfolded state conformations appears to be essential for an adequate description of relative stabilities of protein mutants. GENERAL SIGNIFICANCE: The simulations may be used to link the structural Boltzmann ensembles to relative free enthalpies of folding between mutants and wild-type protein and show that unfolded conformations have to be treated with a sufficient level of detail in free energy calculations of protein stability. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

On the use of time-averaging restraints when deriving biomolecular structure from ³J -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 (CF3)-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 CF3-substituted derivative, that was designed to stabilize the dominant octreotide conformer responsible for receptor binding, turned out to have a lower affinity. The obtained spectroscopic data were interpreted as to show an increased flexibility of the CF3 derivative compared to the unsubstituted octreotide, which could then explain the lower affinity. In this article, we use MD simulation without and with time-averaged NOE distance and time-averaged local-elevation 3J-coupling restraining representing experimental NMR data to determine the conformational properties of the different peptides in the different solvents for which experimental data are available, that are compatible with the NOE atom-atom distance bounds and the 3JHNHα-couplings as derived from the NMR measurements. The conformational ensembles show that the CF3 substitution in combination with the change of solvent from water to methanol leads to a decrease in flexibility and a shift in the populations of the dominant conformers that are compatible with the experimental data.

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. In order to delineate the influence of the ß-sheet core of the fibril upon the structural variability of the loop, the latter is also simulated without the ß-sheet core, but with its N- and C-terminal residues restrained at their positions in the fibril. The analysis of the trajectories shows that the structural variability of the loop is restricted by the ß-sheet core, least at its N-terminal end and most in the middle of the trimer.

Directed evolution of a model primordial enzyme provides insights into the development of the genetic code.

The contemporary proteinogenic repertoire contains 20 amino acids with diverse functional groups and side chain geometries. Primordial proteins, in contrast, were presumably constructed from a subset of these building blocks. Subsequent expansion of the proteinogenic alphabet would have enhanced their capabilities, fostering the metabolic prowess and organismal fitness of early living systems. While the addition of amino acids bearing innovative functional groups directly enhances the chemical repertoire of proteomes, the inclusion of chemically redundant monomers is difficult to rationalize. Here, we studied how a simplified chorismate mutase evolves upon expanding its amino acid alphabet from nine to potentially 20 letters. Continuous evolution provided an enhanced enzyme variant that has only two point mutations, both of which extend the alphabet and jointly improve protein stability by >4 kcal/mol and catalytic activity tenfold. The same, seemingly innocuous substitutions (Ile→Thr, Leu→Val) occurred in several independent evolutionary trajectories. The increase in fitness they confer indicates that building blocks with very similar side chain structures are highly beneficial for fine-tuning protein structure and function.

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 Qexp of an observable quantity Q into structural information is, however, a task beset with theoretical and practical problems: 1) insufficient or inaccurate values of Qexp , 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 Qexp , 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. The process of deriving structural information from measured data is reviewed with an eye to non-experts and newcomers in the field using examples from the literature of the effect of the various choices and approximations involved in the process. A list of choices to be avoided is provided.

Polarizable coarse-grained models for molecular dynamics simulation of liquid cyclohexane.

Force field parameters for polarizable coarse-grained (CG) supra-atomic models of liquid cyclohexane are proposed. Two different bead sizes were investigated, one representing two fine-grained (FG) CH(2)r united atoms of the cyclohexane ring, and one representing three FG CH(2)r united atoms. Electronic polarizability is represented by a massless charge-on-spring particle connected to each CG bead. The model parameters were calibrated against the experimental density and heat of vaporization of liquid cyclohexane, and the free energy of cyclohexane hydration. Both models show good agreement with thermodynamic properties of cyclohexane, yet overestimate the self-diffusion. The dielectric properties of the polarizable models agree very well with experiment.