Solution conformations of structured peptides: continuum electrostatics versus distance-dependent dielectric functions

To compare different implicit solvent potentials, the folding thermodynamics of the helical peptide RN24 and the β-hairpin peptide BH8 are studied by molecular dynamics simulation with adaptive umbrella sampling. As the potential energy functions, the analytical continuum solvent (ACS) potential and three simplified variants, termed EPSR₁, EPSR₄, and EPSR₁₀, are used. The ACS potential is a combination of the standard CHARMM force field for the internal energy (bonds, angles, dihedrals) and the van der Waals energy with the analytical continuum electrostatic (ACE) potential and a non-polar solvation potential. The EPSR potentials differ from the ACS potential by the use of Coulomb's law with a distance-dependent dielectric function to calculate the electrostatic energy. With the ACS potential, quantitative agreement with experiment is obtained for the helix propensity (RN24: 62% calculated vs 50–60% experiment) and the β-hairpin propensity (BH8: 33% calculated vs 19–37% experiment) of the peptides. During the simulations with the EPSR potentials, no significant formation of secondary structure is observed. It is shown that the preference for coil conformations over conformations with secondary structure by the EPSR potentials is due to an overestimation of the energy of salt bridge formation, independent of the magnitude of the Coulomb energy relative to the other energy terms. Possible improvements of the distance-dependent dielectric functions which may permit their application to the simulation of peptide folding, are discussed. Received: 11 July 1998 / Accepted: 22 September 1998 / Published online: 17 December 1998     

Determination of conformational free energies of peptides by multidimensional adaptive umbrella sampling

We improve the multidimensional adaptive umbrella sampling method for the computation of conformational free energies of biomolecules. The conformational transition between the alpha-helical and beta-hairpin conformational states of an alanine decapeptide is used as an example. Convergence properties of the weighted-histogram-analysis-based adaptive umbrella sampling can be improved by using multiple replicas in each adaptive iteration and by using adaptive updating of the bounds of the umbrella potential. Using positional root-mean-square deviations from structures of the alpha-helical and beta-hairpin reference states as reaction coordinates, we obtained well-converged free energy surfaces of both the in-vacuum and in-solution decapeptide systems. From the free energy surfaces well-converged relative free energies between the two conformational states can be derived. Advantages and disadvantages of different methods for obtaining conformational free energies as well as implications of our results in studying conformational transitions of proteins and in improving force field are discussed.     

Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations

A new adaptive umbrella sampling technique for molecular dynamics simulations is described. The high efficiency of the technique renders multidimensional adaptive umbrella sampling possible and thereby enables uniform sampling of the conformational space spanned by several degrees of freedom. The efficiency is achieved by using the weighted histogram analysis method to combine the results from different simulations, by a suitable extrapolation scheme to define the umbrella potential for regions that have not been sampled, and by a criterion to identify simulations during which the system was not in equilibrium. The technique is applied to two test systems, the alanine dipeptide and the threonine dipeptide, to sample the configurational space spanned by one or two dihedral angles. The umbrella potentials applied at the end of each adaptive umbrella sampling run are equal to the negative of the corresponding potentials of mean force. The trajectories obtained in the simulations can be used to calculate dynamical variables that are of interest. An example is the distribution of the distance between the HN and the Hβ proton that can be important for the interpretation of NMR experiments. Factors influencing the accuracy of the calculated quantities are discussed. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1450–1462, 1997     

Implementation of an adaptive umbrella sampling method for the calculation of multidimensional potential of mean force of chemical reactions in solution

We describe the implementation of an adaptive umbrella sampling method, making use of the weighted histogram analysis method, for computing multidimensional potential of mean force for chemical reaction in solution. The approach is illustrated by investigating the effect of aqueous solution on the free energy surface for the proton transfer reaction of [H(3)N-H-NH(3)](+) using a combined quantum mechanical and molecular mechanical AM1/TIP3P potential.     

Determination of free energy profiles by repository based adaptive umbrella sampling: bridging nonequilibrium and quasiequilibrium simulations

We propose a new adaptive sampling approach to determine free energy profiles with molecular dynamics simulations, which is called as "repository based adaptive umbrella sampling" (RBAUS). Its main idea is that a sampling repository is continuously updated based on the latest simulation data, and the accumulated knowledge and sampling history are then employed to determine whether and how to update the biasing umbrella potential for subsequent simulations. In comparison with other adaptive methods, a unique and attractive feature of the RBAUS approach is that the frequency for updating the biasing potential depends on the sampling history and is adaptively determined on the fly, which makes it possible to smoothly bridge nonequilibrium and quasiequilibrium simulations. The RBAUS method is first tested by simulations on two simple systems: a double well model system with a variety of barriers and the dissociation of a NaCl molecule in water. Its efficiency and applicability are further illustrated in ab initio quantum mechanics/molecular mechanics molecular dynamics simulations of a methyl-transfer reaction in aqueous solution.     

Poisson-Boltzmann calculations versus molecular dynamics simulations for calculating the electrostatic potential of a solvated peptide

We performed a very long molecular dynamics simulation of a peptide in explicit water molecules and ions and averaged the electrostatic potential caused by peptide, water and ions at eight points in the vicinity of the peptide. These electrostatic potential values were directly compared to the potential calculated by solving the non-linear Poisson-Boltzmann equation for the system, which describes the solvent using continuum electrostatics. We analyze the contribution of dielectric constant, conformational flexibility and solvation effects on the electrostatic potential at these eight points. Received: 24 April 1998 / Accepted: 4 August 1998 / Published online: 23 November 1998     

Molecular modelling of amphipathic basic model peptides: Interactions in aqueous solutions and in the presence of lipids

Molecular modelling calculations based on experimental data obtained in solution and in small unilamellar vesicles are used to study interactions between amphiphilic basic peptides and membranes. The behaviour of such peptides during the initial and final stages of the adsorption process is our primary interest. Primary sequences of 20 amino acid residues were designed with equal numbers of basic lysines and hydrophobic leucines in order to get an amphipathic α helix. First, in solution, aggregates with an increasing number (up to nine) of helical monomers were built up and the hydrophobic solvent accessible surface per monomer was analysed on energy minimised structures. This showed that aggregates with 5–8 of monomers should be equally probable, in reasonable accordance with experimental data. In addition, models of membranes with 21 dimyristoyl-phosphatidylcholine lipids were constructed; amphiphilic peptides were merged into these assemblies with their axes parallel to the monolayer surface and the whole lipid/peptide complex was submitted to a few steps of simulated annealing and further energy minimisation techniques in order to equilibrate alkyl chains in the vicinity of the peptide. These simulations yield an estimation of the penetration depth for the peptide in the monolayer of ∼3.2 ?, whereas experimental approaches to this question were not productive. The modification in the peptide net electrical charge by interchanging Leu in Lys residues in such systems is also examined: for low-charged peptides the penetration depth increases. Received: 20 May 1998 / Accepted : 3 September 1998 / Published online: 7 December 1998     

Using the local elevation method to construct optimized umbrella sampling potentials: Calculation of the relative free energies and interconversion barriers of glucopyranose ring conformers in water

A method is proposed to combine the local elevation (LE) conformational searching and the umbrella sampling (US) conformational sampling approaches into a single local elevation umbrella sampling (LEUS) scheme for (explicit‐solvent) molecular dynamics (MD) simulations. In this approach, an initial (relatively short) LE build‐up (searching) phase is used to construct an optimized biasing potential within a subspace of conformationally relevant degrees of freedom, that is then used in a (comparatively longer) US sampling phase. This scheme dramatically enhances (in comparison with plain MD) the sampling power of MD simulations, taking advantage of the fact that the preoptimized biasing potential represents a reasonable approximation to the negative of the free energy surface in the considered conformational subspace. The method is applied to the calculation of the relative free energies of β‐D ‐glucopyranose ring conformers in water (within the GROMOS 45A4 force field). Different schemes to assign sampled conformational regions to distinct states are also compared. This approach, which bears some analogies with adaptive umbrella sampling and metadynamics (but within a very distinct implementation), is shown to be: (i) efficient (nearly all the computational effort is invested in the actual sampling phase rather than in searching and equilibration); (ii) robust (the method is only weakly sensitive to the details of the build‐up protocol, even for relatively short build‐up times); (iii) versatile (a LEUS biasing potential database could easily be preoptimized for small molecules and assembled on a fragment basis for larger ones). © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

The topology of catchment regions of potential energy hypersurfaces

A simple proof is presented for a fundamental topological property of catchment regions of potential energy hypersurfaces: each catchment region C(λ,i), representing a chemical species and its conformational range on the potential energy hypersurface, is simply k-connected for each dimension k=1,2,…3N−6−λ, where λ is the index of the catchment region. The consequences of this property on the structure of the fundamental group of reaction mechanisms (the one-dimensional homotopy group of reaction paths) is discussed. Received: 8 July 1998 / Accepted: 6 October 1998 / Published online: 23 February 1999     

Bifurcation of reaction pathways: the set of valley ridge inflection points of a simple three-dimensional potential energy surface

This paper serves for the better understanding of the branching phenomenon of reaction paths of potential energy hypersurfaces in more than two dimensions. We apply the recently proposed reduced gradient following (RGF) method for the analysis of potential energy hypersurfaces having valley-ridge inflection (VRI) points. VRI points indicate the region of possible reaction path bifurcation. The relation between RGF and the so-called global Newton search for stationary points (Branin method) is shown. Using a 3D polynomial test surface, a whole 1D manifold of VRI points is obtained. Its relation to RGF curves, steepest descent and gradient extremals is discussed as well as the relation of the VRI manifold to bifurcation points of these curves. Received: 8 July 1998 / Accepted: 24 August 1998 / Published online: 23 November 1998     

Molecular dynamics coupled with a virtual system for effective conformational sampling

An enhanced conformational sampling method is proposed: virtual‐system coupled canonical molecular dynamics (VcMD). Although VcMD enhances sampling along a reaction coordinate, this method is free from estimation of a canonical distribution function along the reaction coordinate. This method introduces a virtual system that does not necessarily obey a physical law. To enhance sampling the virtual system couples with a molecular system to be studied. Resultant snapshots produce a canonical ensemble. This method was applied to a system consisting of two short peptides in an explicit solvent. Conventional molecular dynamics simulation, which is ten times longer than VcMD, was performed along with adaptive umbrella sampling. Free‐energy landscapes computed from the three simulations mutually converged well. The VcMD provided quicker association/dissociation motions of peptides than the conventional molecular dynamics did. The VcMD method is applicable to various complicated systems because of its methodological simplicity. © 2018 Wiley Periodicals, Inc.     

Diatomic potential well depths from analyses of high-resolution electron energy spectra for autoionizing collision complexes

High-resolution energy spectra of electrons released in Penning ionization collisions of metastable rare gas atoms Rg*(ns) (Rg = He, Ne, Ar, Kr, Xe) with several open-shell and closed-shell atoms are analyzed to determine the well depth of the potential energy curve which describes the respective autoionizing collision complex. We thereby elucidate trends in the chemical interaction of Rg* with closed-shell target atoms A(ns ²) and establish a basis for detailed comparison with the respective interactions involving open-shell, ground state alkali atoms A(ns). From electron energy spectra due to␣associative ionization (RgH⁺ formation) in Rg* + H(1s) collisions, we determine binding energies for the RgH⁺(¹Σ) ground state potential (Rg = Ne, Ar, Kr, Xe) with uncertainties around 0.03 eV. Received: 30 June 1998 / Accepted: 5 August 1998 / Published online: 28 October 1998     

Rapid free energy calculation of peptide self-assembly by REMD umbrella sampling

We extend umbrella sampling with replica exchange steps to calculate free energies that are important in the self-assembly of peptides. This leads to a more than 10-fold speed up over conventional umbrella sampling, thereby providing a practical method to calculate these free-energy differences. This approach can also observe first-order phase transitions and pinpoint the location of the concomitant boundary. When conformational changes are involved, this method can handle peptides up to approximately 7 residues, providing a rapid and accurate assessment of the thermodynamic properties of model systems, and can thus be used to answer fundamental questions about peptide self-assembly. When no major conformational changes are involved, we expect the size limit to be equal to that of standard molecular dynamics.     

Analysis of the statistical error in umbrella sampling simulations by umbrella integration

Umbrella sampling simulations, or biased molecular dynamics, can be used to calculate the free-energy change of a chemical reaction. We investigate the sources of different sampling errors and derive approximate expressions for the statistical errors when using harmonic restraints and umbrella integration analysis. This leads to generally applicable rules for the choice of the bias potential and the sampling parameters. Numerical results for simulations on an analytical model potential are presented for validation. While the derivations are based on umbrella integration analysis, the final error estimate is evaluated from the raw simulation data, and it may therefore be generally applicable as indicated by tests using the weighted histogram analysis method.     

Conformations of ethyl esters versus thiolesters

Ethyl formate and other substituted ethyl formates exist in stable anti and gauche conformations about the COCC dihedral angle, according to microwave spectroscopic studies. Similar studies of ethyl thiolformates characterize stable gauche conformations about the corresponding CSCC dihedral angle in every compound studied, but the anti conformation is found only in ethyl fluorothiolformate and chlorothiolformate. Ab initio calculations that include electron correlation via MP2 or the B3LYP density functional model have been carried out for ethyl and methyl formate and thiolformate and their fluoroformate analogs. These calculations reveal that the potential energy minima at gauche and anti COCC configurations are well developed in every case. However, although the gauche minimum for the CSCC torsional angle is clearly defined, the potential function near the anti CSCC configuration corresponds to a potential energy plateau rather than a minimum. In the case of ethyl fluorothiolformate, a modest well is predicted at the anti CSCC configuration, in agreement with experimental results. Received: 7 July 1998 / Accepted: 21 September 1998 / Published online: 15 February 1999     

Analysis of the RGD sequence in protein structures: comparison to the conformations of the RGDW and DRGDW peptides determined by molecular dynamics simulations

Geometric properties of the RGD sequence in a data set of protein crystal and NMR structures deposited in the Protein Data Bank were examined to identify structural characteristics that are related to cell adhesion activity. Interatomic distances and dihedral angles are examined. These geometric measures are then used in an analysis of the conformations of the RGDW and _DRGDW peptides obtained from molecular dynamics simulations (Stote RH, et al. (2000) J Phys ChemB 104:1624). This analysis leads to the suggestion that differences in the accessible conformations contribute to the difference in biological activity between the RGDW and the _DRGDW peptides. Received: 15 August 2000 / Accepted: 4 October 2000 / Published online: 21 March 2001     

Stationary points on the H2CO potential energy surface: dependence on theoretical level

A total of 36 stationary points have been located on the H₂CO potential energy surface by means of gradient extremal following. These 36 points are believed to represent all the important stationary points on this surface. There is no indication that the structure of the surface becomes less complicated as the size of the basis set is enlarged at the Hartree-Fock level of theory, but many of the second- and third-order saddle points disappear when electron correlation is introduced. Of the ten first-order saddle points (transition structures) located, the majority have reaction paths entering the associated minima in a side-on approach, i.e. these cannot be located by uphill walking from the minimum. Received: 5 February 1998 / Accepted: 21 May 1998 / Published online: 29 July 1998     

Biomolecular free energy profiles by a shooting/umbrella sampling protocol, "BOLAS"

We develop an efficient technique for computing free energies corresponding to conformational transitions in complex systems by combining a Monte Carlo ensemble of trajectories generated by the shooting algorithm with umbrella sampling. Motivated by the transition path sampling method, our scheme "BOLAS" (named after a cowboy's lasso) preserves microscopic reversibility and leads to the correct equilibrium distribution. This makes possible computation of free energy profiles along complex reaction coordinates for biomolecular systems with a lower systematic error compared to traditional, force-biased umbrella sampling protocols. We demonstrate the validity of BOLAS for a bistable potential, and illustrate the method's scope with an application to the sugar repuckering transition in a solvated deoxyadenosine molecule.     

Statistical Analysis of Ion Mobility Spectrometry. I. Unbiased and Guided Replica-Exchange Molecular Dynamics

Achieving (bio)macromolecular structural assignment from the interpretation of ion mobility spectrometry (IMS) experiments requires successful comparison with computer modeling. Replica-exchange molecular dynamics simulations with suitable force fields not only offer a convenient framework to locate relevant conformations, especially in the case of multiple-funnel energy landscapes, but they are also well suited to statistical analyses. In the present paper, we discuss two extensions of the method used to improve its efficiency in the context of IMS. Two doubly-protonated polyalanines [RA₄XA₄K + 2H]²⁺ with X = V and D appear as favorable cases for which the calculated collision cross-section distributions naturally agree with the measurements, providing reliable candidate structures. For these compounds, a careful consideration of other order parameters based on the weighted histogram method resolves several otherwise hidden underlying conformational families. In the case of a much larger peptide exhibiting bistability, assignment is more difficult but could be achieved by guiding the sampling with an umbrella potential using the square gyration radius as the biasing coordinate. Applied to triply protonated bradykinine, the two presented methods indicate that different conformations compatible with the measurements are very close in energy.     

Using simplified protein representation as a reference potential for all-atom calculations of folding free energy

An effective approach for evaluating folding free-energy surfaces of explicit all-atom models is developed and examined. This approach is based on using the potential of a simplified protein model as a reference potential for calculating the free energy of the corresponding explicit model. Preliminary results are presented for the folding free energy of a 12-residue helix. The potential of the method for studies of protein-folding processes is discussed, emphasizing the ability to determine the difference between the results of simplified and explicit models. This can help in establishing the validity of simplified folding models. Received: 30 December 1998 / Accepted: 27 January 1999 / Published online: 5 May 1999