Influence of variation of a side chain on the folding equilibrium of a β‐peptide: Limitations of one‐step perturbation

In a recent study (Lin et al., Helv. Chim. Acta 2011, 94 , 597), the one‐step perturbation method was applied to tackle a challenging computational problem, that is, the calculation of the folding free enthalpies ΔG_F,U of six hepta‐β‐peptides with different, Ala, Val, Leu, Ile, Ser, or Thr, side chains in the fifth residue. The ΔG_F,U values obtained using one‐step perturbation based on a single molecular dynamics simulation of a judiciously chosen reference state with soft‐core atoms in the side chain of the fifth residue showed an overall accuracy of about k_BT for the four peptides with nonpolar side chains, but twice as large deviations were observed for the peptides with polar side chains. Here, alternative reference‐state Hamiltonians that better cover the conformational space relevant to these peptides are investigated, and post simulation rotational sampling of the χ₁ and χ₂ torsional angles of the fifth residue is carried out to sample different orientations of the side chain. A reference state with rather soft atoms yields accurate ΔG_F,U values for the peptides with the Ser and Thr side chains, but it failed to correctly predict the folding free enthalpy for one peptide with a nonpolar side chain, that is, Leu. Based on the results and those of earlier studies, possible ways to improve the accuracy of the efficient one‐step perturbation technique to compute free enthalpies of folding are discussed. © 2013 Wiley Periodicals, Inc.     

Using one‐step perturbation to predict the effect of changing force‐field parameters on the simulated folding equilibrium of a β‐peptide in solution

Computer simulation using molecular dynamics is increasingly used to simulate the folding equilibria of peptides and small proteins. Yet, the quality of the obtained results depends largely on the quality of the force field used. This comprises the solute as well as the solvent model and their energetic and entropic compatibility. It is, however, computational very expensive to perform test simulations for each combination of force‐field parameters. Here, we use the one‐step perturbation technique to predict the change of the free enthalpy of folding of a β‐peptide in methanol solution due to changing a variety of force‐field parameters. The results show that changing the solute backbone partial charges affects the folding equilibrium, whereas this is relatively insensitive to changes in the force constants of the torsional energy terms of the force field. Extending the cut‐off distance for nonbonded interactions beyond 1.4 nm does not affect the folding equilibrium. The same result is found for a change of the reaction‐field permittivity for methanol from 17.7 to 30. The results are not sensitive to the criterion, e.g., atom‐positional RMSD or number of hydrogen bonds, that is used to distinguish folded and unfolded conformations. Control simulations with perturbed Hamiltonians followed by backward one‐step perturbation indicated that quite large perturbations still yield reliable results. Yet, perturbing all solvent molecules showed where the limitations of the one‐step perturbation technique are met. The evaluated methodology constitutes an efficient tool in force‐field development for molecular simulation by reducing the number of required separate simulations by orders of magnitude. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Free‐energy differences between states with different conformational ensembles

Multiple conformations separated by high‐energy barriers represent a challenging problem in free‐energy calculations due to the difficulties in achieving adequate sampling. We present an application of thermodynamic integration (TI) in conjunction with the local elevation umbrella sampling (LE/US) method to improve convergence in alchemical free‐energy calculations. TI‐LE/US was applied to the guanosine triphosphate (GTP) to 8‐Br‐GTP perturbation, molecules that present high‐energy barriers between the anti and syn states and that have inverted preferences for those states. The convergence and reliability of TI‐LE/US was assessed by comparing with previous results using the enhanced‐sampling one‐step perturbation (OSP) method. A linear interpolation of the end‐state biasing potentials was sufficient to dramatically improve sampling along the chosen reaction coordinate. Conformational free‐energy differences were also computed for the syn and anti states and compared to experimental and theoretical results. Additionally, a coupled OSP with LE/US was carried out, allowing the calculation of conformational and alchemical free energies of GTP and 8‐substituted GTP analogs. © 2013 Wiley Periodicals, Inc.     

Using enveloping distribution sampling to compute the free enthalpy difference between right- and left-handed helices of a β-peptide in solution

Recently, the method of enveloping distribution sampling (EDS) to efficiently obtain free enthalpy differences between different molecular systems from a single simulation has been generalized to compute free enthalpy differences between different conformations of a system [Z. X. Lin, H. Y. Liu, S. Riniker, and W. F. van Gunsteren, J. Chem. Theory Comput. 7, 3884 (2011)]. However, the efficiency of EDS in this case is hampered if the parts of the conformational space relevant to the two end states or conformations are far apart and the conformational diffusion from one state to the other is slow. This leads to slow convergence of the EDS parameter values and free enthalpy differences. In the present work, we apply the EDS methodology to a challenging case, i.e., to calculate the free enthalpy difference between a right-handed 2.7(10∕12)-helix and a left-handed 3(14)-helix of a hexa-β-peptide in solution from a single simulation. No transition between the two helices was detected in a standard EDS parameter update simulation, thus enhanced sampling techniques had to be applied, which included adiabatic decoupling (AD) of solute and solvent motions in combination with increasing the solute temperature, and lowering the shear viscosity of the solvent. AD was found to be unsuitable to enhance the sampling of the solute conformations in the EDS parameter update simulations. Lowering the solvent shear viscosity turned out to be useful during EDS parameter update simulations, i.e., it did speed up the conformational diffusion of the solute, more transitions between the two helices were observed. This came at the cost of more CPU time spent due to the shorter time step needed for simulations with the lower solvent shear viscosity. Using an improved EDS parameter update scheme, parameter convergence was five-fold enhanced. The resulting free enthalpy difference between the two helices calculated from EDS agrees well with the result obtained through direct counting from a long MD simulation, while the EDS technique significantly enhances the sampling of both helices over non-helical conformations.     

A new equation of state for polar and nonpolar pure fluids

Chung, T.H., Khan, M.M., Lee, L.L. and Starling, K.E., 1984. A new equation of state for polar and nonpolar pure fluids. Fluid Phase Equilibria, 17: 351–372A new equation of state based on the concept of perturbation theory and the hard-convex-body equation of state has been developed successfully for nonpolar compounds. The equation can predict the thermodynamic properties (density, enthalpy departure and vapor pressure) of a wide range of pure fluids from small, spherical (argon-like) molecules to large, structurally complex molecules. For nonpolar compounds, the equation employs three parameters: the shape, size and energy parameters. For normal paraffins, the size parameter (hard-core volume) is related to the measurable van der Waals volume given by Bondi. For most other compounds, it is related to the critical volume. The shape-parameter values reflect the structure and degree of acentricity of the compound of interest. The equation has been extended to polar and associating compounds by using the mean-potential model. For polar compounds, a fourth parameter is required. The equation has been tested extensively for polar (dipolar and quadrupolar) and hydrogen-bonding compounds. The applicability of this equation for such a wide variety of substances provides an important first step in the development of a composition-dependent equation of state for mixtures.     

A study of molecular vibrational relaxation mechanism in condensed phase based upon mixed quantum-classical molecular dynamics. II. Noncollisional mechanism for the relaxation of a polar solute in supercritical water

Mixed quantum-classical molecular dynamics method has been applied to vibrational relaxation of a hydrophilic model NO in supercritical water at various densities along an isotherm above the critical temperature. The relaxation rate was determined based on Fermi's golden rule at each state point and showed an inverse S-shaped curve as a function of bulk density. The hydration number was also calculated as a function of bulk density based on the calculated radial distribution function, which showed a good correlation with the relaxation rate. Change of the survival probability of the solute vibrational state was analyzed as a function of time together with the trajectory of the solvent water and the interaction with it. We will show that the solvent molecule resides near the solute molecule for a while and the solvent contributes to the relaxation by the random-noiselike Coulombic interaction only when it stays near the solute. After the solvent leaves the solute, it shows no contribution to the relaxation. The relaxation mechanism for this system is significantly different from the collisional one found for a nonpolar solute in nonpolar solvent in Paper I. Then, the relaxation rate is determined, on average, by the hydration number or local density of the solvent. Thus, the density dependence of the relaxation rate for the polar solute in supercritical water is apparently similar to that found for the nonpolar solute in nonpolar solvent, although the molecular process is quite different from each other.     

Self-assembling cyclic peptides: molecular dynamics studies of dimers in polar and nonpolar solvents

The self-assembly of cyclic D,L-alpha-peptides into hollow nanotubes is a crucial mechanistic step in their application as antibacterial and drug-delivery agents. To understand this process, molecular dynamics (MD) simulations were performed on dimers of cyclic peptides formed from cyclo [(-L-Trp-D-N-MeLeu-)4-]2 and cyclo [(-L-Trp-D-Leu-)4-]2 subunits in nonpolar (nonane) and polar (water) solvent. The dimers were observed to be stable only in nonpolar solvent over the full 10 ns length of the MD trajectory. The behavior of the dimers in different solvents is rationalized in terms of the intersubunit hydrogen bonding, hydrogen bonding with the solvent, and planarity of the rings. It is shown that the phi and psi dihedral angles of a single uncapped ring in nonane lie in the beta-sheet region of the Ramachandran plot, and the ring stays in a flat conformation. Steered MD (SMD) simulations based on Jarzynski's equality were performed to obtain the potential of mean force as a function of the distance between the two rings of the capped dimer in nonane. It is also shown that a single peptide subunit prefers to reside close to the nonane/water interface rather than in bulk solvent because of the amphiphilic character of the peptide ring. The present MD results build the foundation for using MD simulations to study the mechanism of the formation of cyclic peptide nanotubes in lipid bilayers.     

Coulomb replica‐exchange method: Handling electrostatic attractive and repulsive forces for biomolecules

We propose a new type of the Hamiltonian replica‐exchange method (REM) for molecular dynamics (MD) and Monte Carlo simulations, which we refer to as the Coulomb REM (CREM). In this method, electrostatic charge parameters in the Coulomb interactions are exchanged among replicas while temperatures are exchanged in the usual REM. By varying the atom charges, the CREM overcomes free‐energy barriers and realizes more efficient sampling in the conformational space than the REM. Furthermore, this method requires only a smaller number of replicas because only the atom charges of solute molecules are used as exchanged parameters. We performed Coulomb replica‐exchange MD simulations of an alanine dipeptide in explicit water solvent and compared the results with those of the conventional canonical, replica exchange, and van der Waals REMs. Two force fields of AMBER parm99 and AMBER parm99SB were used. As a result, the CREM sampled all local‐minimum free‐energy states more frequently than the other methods for both force fields. Moreover, the Coulomb, van der Waals, and usual REMs were applied to a fragment of an amyloid‐β peptide (Aβ) in explicit water solvent to compare the sampling efficiency of these methods for a larger system. The CREM sampled structures of the Aβ fragment more efficiently than the other methods. We obtained β‐helix, α‐helix, 3₁₀‐helix, β‐hairpin, and β‐sheet structures as stable structures and deduced pathways of conformational transitions among these structures from a free‐energy landscape. © 2012 Wiley Periodicals, Inc.     

Comparison of retention on traditional alkyl, polar endcapped, and polar embedded group stationary phases

The development of new RP stationary phases containing polar groups has provided the chromatographer with a variety of stationary phase choices with differing selectivities. Polar endcapped and polar embedded group stationary phases have found use in solving a wide variety of separation problems, especially for the efficient separation of organic bases as well as separations necessitating the use of highly aqueous mobile phases. In this report, the retention thermodynamics of small, nonpolar solutes on traditional alkyl, polar endcapped, and polar embedded group stationary phases are compared. It is found that the nonpolar (methylene) transfer enthalpy is less favorable when polar embedded group phases are used, when compared to traditional or polar endcapped phases. In contrast, the transfer enthalpy of a phenyl group is found to be more favorable when a polar endcapped phase is used. In addition, the retention characteristics of these phases are compared using a set of solutes with differing solvatochromic parameters. Hydrogen-bond acids appear to have enhanced retention on polar embedded group phases, while hydrogen-bond bases have enhanced retention on polar endcapped phases.     

Dispersion corrected interaction of polar and nonpolar fluids confined within carbon nanotubes: Density functional theoretical analysis using Grimme's D3 scheme

The structural and flow characteristics of fluids within carbon nanotube (CNT) is dictated by the interaction of fluid molecules within the nanocavity of CNT. Therefore, in the present study, dispersion corrected density functional theory has been used to investigate the structure and interaction of polar and nonpolar molecules within CNT. The present study shows that there is profound effect on the interaction due to dispersion. The interaction energy of the confined water was found to be reduced with increasing distance of the water molecule from the wall of the CNT. The water is preferentially adsorbed over methane due to stronger interaction with CNT over methane. Further, water is preferentially adsorbed over methanol molecule when interaction is calculated without dispersion but after inclusion of dispersion interaction, the calculated results show that the methanol–CNT interaction is stronger than that of water molecule and hence preferentially adsorbed within the CNT as revealed from MD simulation. The present calculation reveals that that the effect of CNT confinement on the IR spectra of the single file water is quite considerable compared to the IR spectra of tetrahedral bulk water cluster. Therefore, the present results might be useful for the separation of polar molecule from nonpolar molecule during fabrication of CNT‐based filter and purification system.     

Rapid approximate calculation of water binding free energies in the whole hydration domain of (bio)macromolecules

The evaluation of water binding free energies around solute molecules is important for the thermodynamic characterization of hydration or association processes. Here, a rapid approximate method to estimate water binding free energies around (bio)macromolecules from a single molecular dynamics simulation is presented. The basic idea is that endpoint free‐energy calculation methods are applied and the endpoint quantities are monitored on a three‐dimensional grid around the solute. Thus, a gridded map of water binding free energies around the solute is obtained, that is, from a single short simulation, a map of favorable and unfavorable water binding sites can be constructed. Among the employed free‐energy calculation methods, approaches involving endpoint information pertaining to actual thermodynamic integration calculations or endpoint information as exploited in the linear interaction energy method were examined. The accuracy of the approximate approaches was evaluated on the hydration of a cage‐like molecule representing either a nonpolar, polar, or charged water binding site and on α‐ and β‐cyclodextrin molecules. Among the tested approaches, the linear interaction energy method is considered the most viable approach. Applying the linear interaction energy method on the grid around the solute, a semi‐quantitative thermodynamic characterization of hydration around the whole solute is obtained. Disadvantages are the approximate nature of the method and a limited flexibility of the solute. © 2016 Wiley Periodicals, Inc.     

Conformational equilibria of terminally blocked single amino acids at the water-hexane interface. A molecular dynamics study

The conformational equilibria of the acetyl and methyl amide terminally blocked L-alanine, L-leucine and L-glutamine amino acids are examined in vacuum, in bulk water, and at the water-hexane interface, using multi-nanosecond molecular dynamics simulations. The two-dimensional probability distribution functions of finding the peptides at different dihedral angles of the backbone, phi and psi, are calculated, and free energy differences between different conformational states are determined. All three peptides are interfacially active, i.e. tend to accumulate at the interface even though they are not amphiphilic. Conformational states stable in both gas phase and water are also stable in the interfacial environment. Their populations, however, cannot be simply predicted from the knowledge of conformational equilibria in the bulk phases, indicating that the interface exerts a unique effect on the peptides. Conformational preferences in the interfacial environment arise from the interplay between electrostatic and hydrophobic effects. As in an aqueous solution, electrostatic solute-solvent interactions lead to the stabilization of more polar peptide conformations. The hydrophobic effect is manifested at the interface by a tendency to segregate polar and nonpolar moieties of the solute into the aqueous and the hexane phases, respectively. For the terminally blocked glutamine, this favors conformations for which such a segregation is compatible with the formation of strong, backbone-side chain intramolecular hydrogen bonds on the hexane side of the interface. The influence of the hydrophobic effect can be also noted in the orientational preferences of the peptides at the interface. The terminally blocked leucine is oriented such that its nonpolar side chain is buried in hexane, whereas the polar side chain of glutamine is immersed in water. The free energies of rotating the peptides along the axis parallel to the interface by more than 90 degrees are substantial. This indicates that peptide folding at interfaces is strong by driven by the tendency to adopt amphiphilic structures.     

Influence of Variation of a Side Chain on the Folding Equilibrium of a β‐Peptide

The ability to design well‐folding β‐peptides with a specific biological activity requires detailed insight into the relationship between the β‐amino acid sequence and the three‐dimensional structure of the peptide. Here, we present a molecular‐dynamics (MD) study of the influence of a variation of a side chain on the folding equilibrium of a β‐heptapeptide that folds into a 3₁₄‐helical structure. The side chain of the 5th residue, a valine, was changed into five differently branched side chains of different lengths and polarity, Ala, Leu, Ile, Ser, and Thr. Two computational techniques, long‐time MD simulations and the one‐step perturbation method, were used to obtain free enthalpies of folding. The simulations show that all six peptides exhibit similar folding behavior, and that their dominant fold is the same, i.e., a 3₁₄‐helix. Despite the similarities of their structural properties, a small stabilization effect of ca. 2 kJ mol^?1 on the folding equilibrium of the 3₁₄‐helical structure due to a branching C_γ‐atom in the β³‐side chain is observed. These results confirm those of previous circular dichroism (CD) studies. The length of side chain and its polarity seem to have no apparent (de)stabilization effect. Application of the cost‐effective one‐step perturbation method to predict free‐enthalpy differences appeared to yield an overall accuracy of about k_BT, which is not sufficient to detect the small stabilization effect.     

Solvent effect on the absorption spectra of coumarin 120 in water: A combined quantum mechanical and molecular mechanical study

The solvent effect on the absorption spectra of coumarin 120 (C120) in water was studied utilizing the combined quantum mechanical∕molecular mechanical (QM∕MM) method. In molecular dynamics (MD) simulation, a new sampling scheme was introduced to provide enough samples for both solute and solvent molecules to obtain the average physical properties of the molecules in solution. We sampled the structure of the solute and solvent molecules separately. First, we executed a QM∕MM MD simulation, where we sampled the solute molecule in solution. Next, we chose random solute structures from this simulation and performed classical MD simulation for each chosen solute structure with its geometry fixed. This new scheme allowed us to sample the solute molecule quantum mechanically and sample many solvent structures classically. Excitation energy calculations using the selected samples were carried out by the generalized multiconfigurational perturbation theory. We succeeded in constructing the absorption spectra and realizing the red shift of the absorption spectra found in polar solvents. To understand the motion of C120 in water, we carried out principal component analysis and found that the motion of the methyl group made the largest contribution and the motion of the amino group the second largest. The solvent effect on the absorption spectrum was studied by decomposing it in two components: the effect from the distortion of the solute molecule and the field effect from the solvent molecules. The solvent effect from the solvent molecules shows large contribution to the solvent shift of the peak of the absorption spectrum, while the solvent effect from the solute molecule shows no contribution. The solvent effect from the solute molecule mainly contributes to the broadening of the absorption spectrum. In the solvent effect, the variation in C-C bond length has the largest contribution on the absorption spectrum from the solute molecule. For the solvent effect on the absorption spectrum from the solvent molecules, the solvent structure around the amino group of C120 plays the key role.     

Assessment of enveloping distribution sampling to calculate relative free enthalpies of binding for eight netropsin-DNA duplex complexes in aqueous solution

The performance of enveloping distribution sampling (EDS) simulations to estimate free enthalpy differences associated with seven alchemical transformations of A-T into G-C base pairs at the netropsin binding site in the minor groove of a 13-base pair DNA duplex in aqueous solution is evaluated. It is demonstrated that sufficient sampling can be achieved with a two-state EDS Hamiltonian even for large perturbations such as the simultaneous transformation of up to three A-T into three G-C base pairs. The two parameters required to define the EDS reference state Hamiltonian are obtained automatically using a modified version of a scheme presented in earlier work. The sensitivity of the configurational sampling to a variation of these parameters is investigated in detail. Although for relatively small perturbations, that is, one base pair, the free enthalpy estimate depends only weakly on the EDS parameters, the sensitivity is stronger for the largest perturbation. Yet, EDS offers various convenient measures to evaluate the degree of sampling and thus the reliability of the free enthalpy estimate and appears to be an efficient alternative to the conventional thermodynamic integration methodology to obtain free energy differences for molecular systems.     

An efficient method for computing excess free energy of liquid

We present a new theoretical method for efficient calculation of free energy of liquid. This interaction entropy method allows one to compute entropy and free energy of liquid from standard single step MD (molecular dynamics) simulation directly in liquid state without the need to perform MD simulations at many intermediate states as required in thermodynamic integration or free energy perturbation methods. In this new approach, one only needs to evaluate the interaction energy of a single (fixed) liquid molecule with the rest of liquid molecules as a function of time from a standard MD simulation of liquid and the fluctuation of distribution of this interaction energy is then used to calculate the interaction entropy of the liquid. Explicit theoretical derivation of this interaction entropy approach is provided and numerical calculations for the benchmark liquid water system were carried out using three different water models. Numerical analysis of the result was performed and comparison of the computational result with experimental data and other theoretical results were provided. Excellent agreement of calculated free energies with the experimental data using TIP4P model is obtained for liquid water.     

Coupled semiempirical quantum mechanics and molecular mechanics (QM/MM) calculations on the aqueous solvation free energies of ionized molecules

The aqueous solvation free energies of ionized molecules were computed using a coupled quantum mechanical and molecular mechanical (QM/MM) model based on the AM1, MNDO, and PM3 semiempirical molecular orbital methods for the solute molecule and the TIP3P molecular mechanics model for liquid water. The present work is an extension of our model for neutral solutes where we assumed that the total free energy is the sum of components derived from the electrostatic/polarization terms in the Hamiltonian plus an empirical “nonpolar” term. The electrostatic/polarization contributions to the solvation free energies were computed using molecular dynamics (MD) simulation and thermodynamic integration techniques, while the nonpolar contributions were taken from the literature. The contribution to the electrostatic/polarization component of the free energy due to nonbonded interactions outside the cutoff radii used in the MD simulations was approximated by a Born solvation term. The experimental free energies were reproduced satisfactorily using variational parameters from the vdW terms as in the original model, in addition to a parameter from the one-electron integral terms. The new one-electron parameter was required to account for the short-range effects of overlapping atomic charge densities. The radial distribution functions obtained from the MD simulations showed the expected H-bonded structures between the ionized solute molecule and solvent molecules. We also obtained satisfactory results by neglecting both the empirical nonpolar term and the electronic polarization of the solute, i.e., by implementing a nonpolarization model. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1028–1038, 1999     

A solvation‐free‐energy functional: A reference‐modified density functional formulation

The three‐dimensional reference interaction site model (3D‐RISM) theory, which is one of the most applicable integral equation theories for molecular liquids, overestimates the absolute values of solvation‐free‐energy (SFE) for large solute molecules in water. To improve the free‐energy density functional for the SFE of solute molecules, we propose a reference‐modified density functional theory (RMDFT) that is a general theoretical approach to construct the free‐energy density functional systematically. In the RMDFT formulation, hard‐sphere (HS) fluids are introduced as the reference system instead of an ideal polyatomic molecular gas, which has been regarded as the appropriate reference system of the interaction‐site‐model density functional theory for polyatomic molecular fluids. We show that using RMDFT with a reference HS system can significantly improve the absolute values of the SFE for a set of neutral amino acid side‐chain analogues as well as for 504 small organic molecules. © 2015 Wiley Periodicals, Inc.     

Enthalpy distributions and densities of states in linear polymers

Using polypropylene as an example, we applied a method we have recently developed to calculate the probability distribution of enthalpy from the temperature dependence of the heat capacity. The method involves the use of local temperature expansions of the heat capacity to calculate a set of moments of the enthalpy distribution. Using the maximum‐entropy method, one can then construct the enthalpy distribution for the system. The method is completely model free. The enthalpy distribution so obtained is the analogue of the Maxwell–Boltzmann distribution of kinetic energies for a gas, and like that function, tells one the distribution of enthalpies that an average unit in the polymer chain can have, a quantity that is crucial to understanding the chemical and physical properties of a polymer. Given the enthalpy distribution, one can then calculate the Gibbs free energy and the density of states that correspond to a particular value of enthalpy, thus giving one an expanded thermodynamics of the system in addition to the usual average quantities. We illustrate the fact that the Gibbs free‐energy distribution for this system scales as a simple function of temperature and that the density‐of ‐states function yields a simple empirical partition function for the system giving both the average thermodynamics and the distribution functions. © 2001 John Wiley & Sons, Inc. J Polym Sci Part B: Polym Phys 39: 1513–1518, 2001