Calculating solution redox free energies with ab initio quantum mechanical/molecular mechanical minimum free energy path method.
Ontology highlight
ABSTRACT:
SUBMITTER: Zeng X
PROVIDER: S-EPMC2736613 | biostudies-other | 2009 Apr
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Macroscopic pKa values were calculated for all compounds in the SAMPL6 blind prediction challenge, based on quantum chemical calculations with a continuum solvation model and a linear correction derived from a small training set. Microscopic pKa values were derived from the gas-phase free energy difference between protonated and deprotonated forms together with the Conductor-like Polarizable Continuum Solvation Model and the experimental solvation free energy of the proton. pH-dependent microstate free energies were obtained from the microscopic pKas with a maximum likelihood estimator and appropriately summed to yield macroscopic pKa values or microstate populations as function of pH. We assessed the accuracy of three approaches to calculate the microscopic pKas: direct use of the quantum mechanical free energy differences and correction of the direct values for short-comings in the QM solvation model with two different linear models that we independently derived from a small training set of 38 compounds with known pKa. The predictions that were corrected with the linear models had much better accuracy [root-mean-square error (RMSE) 2.04 and 1.95 pKa units] than the direct calculation (RMSE 3.74). Statistical measures indicate that some systematic errors remain, likely due to differences in the SAMPL6 data set and the small training set with respect to their interactions with water. Overall, the current approach provides a viable physics-based route to estimate macroscopic pKa values for novel compounds with reasonable accuracy.
| S-EPMC6240492 | biostudies-literature
Project description:Despite being at the heart of chemical thought, the curly arrow notation of reaction mechanisms has been treated with suspicion-the connection with rigorous molecular quantum mechanics being unclear. The connection requires a view of the wavefunction that goes beyond molecular orbitals and rests on the most fundamental property of electrons. The antisymmetry of electronic wavefunctions requires that an N-electron wavefunction repeat itself in 3N dimensions, thus exhibiting tiles. Inspection of wavefunction tiles permits insight into structure and mechanism. Here, we demonstrate that analysis of the wavefunction tile along a reaction coordinate reveals the electron movements depicted by the curly arrow notation for several reactions. The Diels-Alder reaction is revealed to involve the separation and counter propagation of electron spins. This unprecedented method of extracting the movements of electrons during a chemical reaction is a breakthrough in connecting traditional depictions of chemical mechanism with state-of-the-art quantum chemical calculations.
| S-EPMC5897577 | biostudies-literature
Project description:The majority of macromolecular crystal structures are determined using the method of molecular replacement, in which known related structures are rotated and translated to provide an initial atomic model for the new structure. A theoretical understanding of the signal-to-noise ratio in likelihood-based molecular replacement searches has been developed to account for the influence of model quality and completeness, as well as the resolution of the diffraction data. Here we show that, contrary to current belief, molecular replacement need not be restricted to the use of models comprising a substantial fraction of the unknown structure. Instead, likelihood-based methods allow a continuum of applications depending predictably on the quality of the model and the resolution of the data. Unexpectedly, our understanding of the signal-to-noise ratio in molecular replacement leads to the finding that, with data to sufficiently high resolution, fragments as small as single atoms of elements usually found in proteins can yield ab initio solutions of macromolecular structures, including some that elude traditional direct methods.
| S-EPMC5389281 | biostudies-literature
Project description:Crystal structure determination requires solving the phase problem. This can be accomplished using ab initio direct methods for small molecules and macromolecules at resolutions higher than 1.2 Å, whereas macromolecular structure determination at lower resolution requires either molecular replacement using a homologous structure or experimental phases using a derivative such as covalent labeling (for example selenomethionine or mercury derivatization) or heavy-atom soaking (for example iodide ions). Here, a case is presented in which crystals were obtained from a 30.8 kDa protein sample and yielded a 1.6 Å resolution data set with a unit cell that could accommodate approximately 8 kDa of protein. Thus, it was unclear what had been crystallized. Molecular replacement with pieces of homologous proteins and attempts at iodide ion soaking failed to yield a solution. The crystals could not be reproduced. Sequence-independent molecular replacement using the structures available in the Protein Data Bank also failed to yield a solution. Ultimately, ab initio structure solution proved successful using the program ARCIMBOLDO, which identified two α-helical elements and yielded interpretable maps. The structure was the C-terminal dimerization domain of the intended target from Mycobacterium smegmatis. This structure is presented as a user-friendly test case in which an unknown protein fragment could be determined using ARCIMBOLDO.
| S-EPMC6130419 | biostudies-literature
Project description:Mechanical properties of FeCrMn-based steels are of major importance for practical applications. In this work, we investigate mechanical properties of disordered paramagnetic fcc FeCr 10 ⁻ 16 Mn 12 ⁻ 32 alloys using density functional theory. The effects of composition and temperature changes on the magnetic state, elastic properties and stacking fault energies of the alloys are studied. Calculated dependencies of the lattice and elastic constants are used to evaluate the effect of the solid solution strengthening by Mn and Cr using a modified Labusch-Nabarro model and a model for concentrated alloys. The effect of Cr and Mn alloying on the stacking fault energies is calculated and discussed in connection to possible deformation mechanisms.
| S-EPMC6480156 | biostudies-literature
Project description:Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. The situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.
| S-EPMC8864787 | 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:We use a combination of database analysis and quantum chemical studies to investigate the role of cis and trans Hoogsteen:Hoogsteen (H:H) base pairs and associated higher-order structures in RNA. We add three new examples to the list of previously identified base-pair combinations belonging to these families and, in addition to contextual classification and characterization of their structural and energetic features, we compare their interbase interaction energies and propensities toward participation in triplets and quartets. We find that some base pairs, which are nonplanar in their isolated minimum energy geometries, attain planarity and stability upon triplet formation. A:A H:H trans is the most frequent H:H combination in RNA structures. This base pair occurs at many distinct positions in known rRNA structures, where it helps in the interaction of ribosomal domains in the 50S subunit. It is also present as a part of tertiary interaction in tRNA structures. Although quantum chemical studies suggest an intrinsically nonplanar geometry for this base pair in isolated form, it has the tendency to attain planar geometry in RNA crystal structures by forming higher-order tertiary interactions or in the presence of additional base-phosphate interactions. The tendency of this base pair to form such additional interactions may be helpful in bringing together different segments of RNA, thus making it suitable for the role of facilitator for RNA folding. This also explains the high occurrence frequency of this base pair among all H:H interactions.
| S-EPMC2856888 | biostudies-literature
Project description:In silico screening of drug-target interactions is a key part of the drug discovery process. Changes in the drug scaffold via contraction or expansion of rings, the breaking of rings, and the introduction of cyclic structures from acyclic structures are commonly applied by medicinal chemists to improve binding affinity and enhance favorable properties of candidate compounds. These processes, commonly referred to as scaffold hopping, are challenging to model computationally. Although relative binding free energy (RBFE) calculations have shown success in predicting binding affinity changes caused by perturbing R-groups attached to a common scaffold, applications of RBFE calculations to modeling scaffold hopping are relatively limited. Scaffold hopping inevitably involves breaking and forming bond interactions of quadratic functional forms, which is highly challenging. A novel method for handling ring opening/closure/contraction/expansion and linker contraction/expansion is presented here. To the best of our knowledge, RBFE calculations on linker contraction/expansion have not been previously reported. The method uses auxiliary restraints to hold the atoms at the ends of a bond in place during the breaking and forming of the bonds. The broad applicability of the method was demonstrated by examining perturbations involving small-molecule macrocycles and mutations of proline in proteins. High accuracy was obtained using the method for most of the perturbations studied. The rigor of the method was isolated from the force field by validating the method using relative and absolute hydration free energy calculations compared to standard simulation results. Unlike other methods that rely on λ-dependent functional forms for bond interactions, the method presented here can be employed using modern molecular dynamics software without modification of codes or force field functions.
| S-EPMC8215533 | biostudies-literature
Project description:Recently developed high-entropy alloys (HEAs) consisting of multiple principal elements represent a new field of metallurgy and have demonstrated appealing properties for a wide range of applications. Using ab initio alloy theory, we reveal the alloying effect on the elastic properties and the ideal tensile strength (ITS) in the [001] direction of four body-centered cubic (bcc) refractory HEAs based on Zr, V, Ti, Nb, and Hf. We find that these HEAs show high elastic anisotropy and large positive Cauchy pressure, suggesting good extrinsic ductility. Starting from ZrNbHf, it is found that the ITS decreases with equimolar Ti addition. On the other hand, if both Ti and V are added to ZrNbHf, the ITS is enhanced by about 42%. An even more captivating effect is the ITS increase by about 170%, if Ti and V are substituted for Hf. The alloying effect on the ITS is explained by the d-band filling. An intrinsic brittle-to-ductile transition is found in terms of the failure mode under uniaxial tension. These investigations suggest that intrinsically ductile HEAs with high ideal strength can be achieved by controlling the proportion of group four elements to group five elements.
| S-EPMC4510575 | biostudies-other