Vibrational subsystem analysis: A method for probing free energies and correlations in the harmonic limit.
Ontology highlight
ABSTRACT:
SUBMITTER: Woodcock HL
PROVIDER: S-EPMC2676147 | biostudies-other | 2008 Dec
REPOSITORIES: biostudies-other
Similar Datasets
Project description:We present a methodology to compute, at reduced computational cost, Gibbs free energies, enthalpies, and entropies of adsorption from molecular dynamics. We calculate vibrational partition functions from vibrational energies, which we obtain from the vibrational density of states by projection on the normal modes. The use of a set of well-chosen reference structures along the trajectories accounts for the anharmonicities of the modes. For the adsorption of methane, ethane, and propane in the H-CHA zeolite, we limit our treatment to a set of vibrational modes localized at the adsorption site (zeolitic OH group) and the alkane molecule interacting with it. Only two short trajectories (1-20 ps) are required to reach convergence (<1 kJ/mol) for the thermodynamic functions. The mean absolute deviations from the experimentally measured values are 2.6, 2.8, and 4.7 kJ/mol for the Gibbs free energy, the enthalpy, and the entropy term (-TΔS), respectively. In particular, the entropy terms show a major improvement compared to the harmonic approximation and almost reach the accuracy of the previous use of anharmonic frequencies obtained with curvilinear distortions of individual modes. The thermodynamic functions so obtained follow the trend of the experimental values for methane, ethane, and propane, and the Gibbs free energy of adsorption at experimental conditions is correctly predicted to change from positive for methane (5.9 kJ/mol) to negative for ethane (-4.8 kJ/mol) and propane (-7.1 kJ/mol).
| S-EPMC8444336 | biostudies-literature
Project description:Molecular vibrations, especially low frequency motions, may be used as an indication of the rigidity or the flatness of the protein folding energy landscape. We have studied the vibrational properties of native folded as well as random coil structures of more than 60 polypeptides. The picture we obtain allows us to perceive how and why the energy landscape progressively rigidifies while still allowing potential flexibility. Compared with random coil structures, both alpha-helices and beta-hairpins are vibrationally more flexible. The vibrational properties of loop structures are similar to those of the corresponding random coil structures. Inclusion of an alpha-helix tends to rigidify peptides and so-called building blocks of the structure, whereas the addition of a beta-structure has less effect. When small building blocks coalesce to form larger domains, the protein rigidifies. However, some folded native conformations are still found to be vibrationally more flexible than random coil structures, for example, beta(2)-microglobulin and the SH3 domain. Vibrational free energy contributes significantly to the thermodynamics of protein folding and affects the distribution of the conformational substates. We found a weak correlation between the vibrational folding energy and the protein size, consistent with both previous experimental estimates and theoretical partition of the heat capacity change in protein folding.
| S-EPMC1301155 | biostudies-other
Project description:In this work, regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit are demonstrated for the first time. The performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with two Pople basis sets (6-311++G** and 6-311++G(3df,2pd)), the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, 6) was tested.The BLYP-calculated harmonic frequencies were found to be markedly closer than the B3LYP-calculated harmonic frequencies to the experimentally derived values, while the calculated anharmonic frequencies consistently underestimated the observed wavenumbers. The different basis set families gave very similar estimated values for the CBS parameters. The anharmonic frequencies calculated with B3LYP/aug-pc-3 were consistently significantly higher than those obtained with the pc-3 basis set; applying the aug-pcseg-n basis set family alleviated this problem. Utilization of B3LYP/aug-pcseg-n basis sets instead of B3LYP/aug-cc-pVXZ, which is computationally less expensive, is suggested for medium-sized molecules. Harmonic BLYP/pc-2 calculations produced fairly accurate ethylene frequencies. Graphical Abstract In this study, the performance of the VPT2 scheme implemented using density functional theory (DFT-BLYP and DFT-B3LYP) in combination with the polarization-consistent basis sets pc-n, aug-pc-n, and pcseg-n (n = 0, 1, 2, 3, 4), and the correlation-consistent basis sets cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5, and 6) was tested. For the first time, we demonstrated regular convergence patterns of the structural, harmonic, and VPT2-calculated anharmonic vibrational parameters of ethylene towards the Kohn-Sham complete basis set (KS CBS) limit.
| S-EPMC4723629 | biostudies-literature
Project description:The thermodynamic reasons why membrane proteins form stable complexes inside the hydrophobic lipid bilayer remain poorly understood. This is largely because of a lack of membrane-protein systems amenable for equilibrium studies and a limited number of methods for measuring these reactions. Recently, we reported the equilibrium dimerization of the CLC-ec1 Cl-/H+ transporter in lipid bilayers (Chadda et al. 2016. eLife https://doi.org/10.7554/eLife.17438), which provided a new type of model system for studying protein association in membranes. The measurement was conducted using the subunit-capture approach, involving passive dilution of the protein in large multilamellar vesicles, followed by single-molecule photobleaching analysis of the Poisson distribution describing protein encapsulation into extruded liposomes. To estimate the fraction of dimers (FDimer ) as a function of protein density, the photobleaching distributions for the nonreactive, ideal monomer and dimer species must be known so that random co-capture probabilities can be accounted for. Previously, this was done by simulating the Poisson process of protein reconstitution into a known size distribution of liposomes composed of Escherichia coli polar lipids (EPLs). In the present study, we investigate the dependency of FDimer and ?G° on the modeling through a comparison of different liposome size distributions (EPL versus 2:1 POPE/POPG). The results show that the estimated FDimer values are comparable, except at higher densities when liposomes become saturated with protein. We then develop empirical controls to directly measure the photobleaching distributions of the nonreactive monomer (CLC-ec1 I201W/I422W) and ideal dimer (WT CLC-ec1 cross-linked by glutaraldehyde or CLC-ec1 R230C/L249C cross-linked by a disulfide bond). The measured equilibrium constants do not depend on the correction method used, indicating the robustness of the subunit-capture approach. This strategy therefore presents a model-free way to quantify protein dimerization in lipid bilayers, offering a simplified strategy in the ongoing effort to characterize equilibrium membrane-protein reactions in membranes.
| S-EPMC5806674 | biostudies-literature
Project description:A simple and robust formulation of the path-independent confinement method for the calculation of free energies is presented. The simplified confinement method (SCM) does not require matrix diagonalization or switching off the molecular force field, and has a simple convergence criterion. The method can be readily implemented in molecular dynamics programs with minimal or no code modifications. Because the confinement method is a special case of thermodynamic integration, it is trivially parallel over the integration variable. The accuracy of the method is demonstrated using a model diatomic molecule, for which exact results can be computed analytically. The method is then applied to the alanine dipeptide in vacuum, and to the ?-helix ? ?-sheet transition in a 16-residue peptide modeled in implicit solvent. The SCM requires less effort for the calculation of free energy differences than previous formulations because it does not require computing normal modes. The SCM has a diminished advantage for determining absolute free energy values, because it requires decreasing the MD integration step to obtain accurate results. An approximate confinement procedure is introduced, which can be used to estimate directly the configurational entropy difference between two macrostates, without the need for additional computation of the difference in the free energy or enthalpy. The approximation has convergence properties similar to those of the standard confinement method for the calculation of free energies. The use of the approximation requires about 5 times less wall-clock simulation time than that needed to compute enthalpy differences to similar precision from an MD trajectory. For the biomolecular systems considered in this study, the errors in the entropy approximation are under 10%. Practical applications of the methods to proteins are currently limited to implicit solvent simulations.
| S-EPMC3569517 | biostudies-literature
Project description:A minimal observable length is a common feature of theories that aim to merge quantum physics and gravity. Quantum mechanically, this concept is associated with a nonzero minimal uncertainty in position measurements, which is encoded in deformed commutation relations. In spite of increasing theoretical interest, the subject suffers from the complete lack of dedicated experiments and bounds to the deformation parameters have just been extrapolated from indirect measurements. As recently proposed, low-energy mechanical oscillators could allow to reveal the effect of a modified commutator. Here we analyze the free evolution of high-quality factor micro- and nano-oscillators, spanning a wide range of masses around the Planck mass mP (≈ 22 μg). The direct check against a model of deformed dynamics substantially lowers the previous limits on the parameters quantifying the commutator deformation.
| S-EPMC4557370 | biostudies-literature
Project description:CXC and CC chemokines are involved in numerous biological processes, and their function in situ may be significantly influenced by heterodimer formation, as was recently reported, for example, for CXC chemokines CXCL4/PF4 and CXCL8/IL8 that interact to form heterodimers that modulate chemotactic and cell proliferation activities. Here we used molecular dynamics simulations to determine relative association free energies (overall average and per residue) for homo- and heterodimer pairs of CXC (CXCL4/PF4, CXCL8/IL8, CXCL1/Gro-alpha, and CXCL7/NAP-2) and CC (CCL5/RANTES, CCL2/MCP-1, and CCL8/MCP-2) chemokines. Even though structural homology among monomer folds of all CXC and CC chemokines permits heterodimer assembly, our calculated association free energies depend upon the particular pair of chemokines in terms of the net electrostatic and nonelectrostatic forces involved, as well as (for CC/CXC mixed chemokines) the selection of dimer type (CC or CXC). These relative free energies indicate that association of some pairs of chemokines is more favorable than others. Our approach is validated by correlation of calculated and experimentally determined free energies. Results are discussed in terms of CXC and CC chemokine function and have significant biological implications.
| S-EPMC2527121 | biostudies-literature
Project description:Vibrational modes of molecules are fundamental properties determined by intramolecular bonding, atomic masses, and molecular geometry, and often serve as important channels for dissipation in nanoscale processes. Although single-molecule junctions have been used to manipulate electronic structure and related functional properties of molecules, electrical control of vibrational mode energies has remained elusive. Here we use simultaneous transport and surface-enhanced Raman spectroscopy measurements to demonstrate large, reversible, voltage-driven shifts of vibrational mode energies of C60 molecules in gold junctions. C60 mode energies are found to vary approximately quadratically with bias, but in a manner inconsistent with a simple vibrational Stark effect. Our theoretical model instead suggests that the mode shifts are a signature of bias-driven addition of electronic charge to the molecule. These results imply that voltage-controlled tuning of vibrational modes is a general phenomenon at metal-molecule interfaces and is a means of achieving significant shifts in vibrational energies relative to a pure Stark effect.
| S-EPMC3910565 | biostudies-literature
Project description:A quantum mechanical/molecular mechanical minimum free energy path (QM/MM-MFEP) method was developed to calculate the redox free energies of large systems in solution with greatly enhanced efficiency for conformation sampling. The QM/MM-MFEP method describes the thermodynamics of a system on the potential of mean force surface of the solute degrees of freedom. The molecular dynamics (MD) sampling is only carried out with the QM subsystem fixed. It thus avoids "on-the-fly" QM calculations and thus overcomes the high computational cost in the direct QM/MM MD sampling. In the applications to two metal complexes in aqueous solution, the new QM/MM-MFEP method yielded redox free energies in good agreement with those calculated from the direct QM/MM MD method. Two larger biologically important redox molecules, lumichrome and riboflavin, were further investigated to demonstrate the efficiency of the method. The enhanced efficiency and uncompromised accuracy are especially significant for biochemical systems. The QM/MM-MFEP method thus provides an efficient approach to free energy simulation of complex electron transfer reactions.
| S-EPMC2736613 | biostudies-other
Project description:We report a theoretical description and numerical tests of the extended-system adaptive biasing force method (eABF), together with an unbiased estimator of the free energy surface from eABF dynamics. Whereas the original ABF approach uses its running estimate of the free energy gradient as the adaptive biasing force, eABF is built on the idea that the exact free energy gradient is not necessary for efficient exploration, and that it is still possible to recover the exact free energy separately with an appropriate estimator. eABF does not directly bias the collective coordinates of interest, but rather fictitious variables that are harmonically coupled to them; therefore is does not require second derivative estimates, making it easily applicable to a wider range of problems than ABF. Furthermore, the extended variables present a smoother, coarse-grain-like sampling problem on a mollified free energy surface, leading to faster exploration and convergence. We also introduce CZAR, a simple, unbiased free energy estimator from eABF trajectories. eABF/CZAR converges to the physical free energy surface faster than standard ABF for a wide range of parameters.
| S-EPMC5402294 | biostudies-literature