Characterization of biaryl torsional energetics and its treatment in OPLS all-atom force fields.
Ontology highlight
ABSTRACT:
SUBMITTER: Dahlgren MK
PROVIDER: S-EPMC3707137 | biostudies-literature | 2013 May
REPOSITORIES: biostudies-literature
Publications
Characterization of biaryl torsional energetics and its treatment in OPLS all-atom force fields.
Dahlgren Markus K MK Schyman Patric P Tirado-Rives Julian J Jorgensen William L WL
Journal of chemical information and modeling 20130513 5
The frequency of biaryl substructures in a database of approved oral drugs has been analyzed. This led to designation of 20 prototypical biaryls plus 10 arylpyridinones for parametrization in the OPLS all-atom force fields. Bond stretching, angle-bending, and torsional parameters were developed to reproduce the MP2 geometries and torsional energy profiles. The transferability of the new parameters was tested through their application to three additional biaryls. The torsional energetics for the ...[more]
Similar Datasets
Project description:The development and validation of new peptide dihedral parameters are reported for the OPLS-AA force field. High accuracy quantum chemical methods were used to scan φ, ψ, χ1, and χ2 potential energy surfaces for blocked dipeptides. New Fourier coefficients for the dihedral angle terms of the OPLS-AA force field were fit to these surfaces, utilizing a Boltzmann-weighted error function and systematically examining the effects of weighting temperature. To prevent overfitting to the available data, a minimal number of new residue-specific and peptide-specific torsion terms were developed. Extensive experimental solution-phase and quantum chemical gas-phase benchmarks were used to assess the quality of the new parameters, named OPLS-AA/M, demonstrating significant improvement over previous OPLS-AA force fields. A Boltzmann weighting temperature of 2000 K was determined to be optimal for fitting the new Fourier coefficients for dihedral angle parameters. Conclusions are drawn from the results for best practices for developing new torsion parameters for protein force fields.
| S-EPMC4504185 | biostudies-literature
Project description:The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and ?-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.
| S-EPMC5816952 | biostudies-literature
Project description:We introduce a toolset of program libraries collectively titled multipurpose atom-typer for CHARMM (MATCH) for the automated assignment of atom types and force field parameters for molecular mechanics simulation of organic molecules. The toolset includes utilities for the conversion of multiple chemical structure file formats into a molecular graph. A general chemical pattern-matching engine using this graph has been implemented whereby assignment of molecular mechanics atom types, charges, and force field parameters are achieved by comparison against a customizable list of chemical fragments. While initially designed to complement the CHARMM simulation package and force fields by generating the necessary input topology and atom-type data files, MATCH can be expanded to any force field and program, and has core functionality that makes it extendable to other applications such as fragment-based property prediction. In this work, we demonstrate the accurate construction of atomic parameters of molecules within each force field included in CHARMM36 through exhaustive cross validation studies illustrating that bond charge increment rules derived from one force field can be transferred to another. In addition, using leave-one-out substitution it is shown that it is also possible to substitute missing intra and intermolecular parameters with ones included in a force field to complete the parameterization of novel molecules. Finally, to demonstrate the robustness of MATCH and the coverage of chemical space offered by the recent CHARMM general force field (Vanommeslaeghe, et al., J Comput Chem 2010, 31, 671), one million molecules from the PubChem database of small molecules are typed, parameterized, and minimized.
| S-EPMC3228871 | biostudies-literature
Project description:Traditional approaches to specifying a molecular mechanics force field encode all the information needed to assign force field parameters to a given molecule into a discrete set of atom types. This is equivalent to a representation consisting of a molecular graph comprising a set of vertices, which represent atoms labeled by atom type, and unlabeled edges, which represent chemical bonds. Bond stretch, angle bend, and dihedral parameters are then assigned by looking up bonded pairs, triplets, and quartets of atom types in parameter tables to assign valence terms and using the atom types themselves to assign nonbonded parameters. This approach, which we call indirect chemical perception because it operates on the intermediate graph of atom-typed nodes, creates a number of technical problems. For example, atom types must be sufficiently complex to encode all necessary information about the molecular environment, making it difficult to extend force fields encoded this way. Atom typing also results in a proliferation of redundant parameters applied to chemically equivalent classes of valence terms, needlessly increasing force field complexity. Here, we describe a new approach to assigning force field parameters via direct chemical perception. Rather than working through the intermediary of the atom-typed graph, direct chemical perception operates directly on the unmodified chemical graph of the molecule to assign parameters. In particular, parameters are assigned to each type of force field term (e.g., bond stretch, angle bend, torsion, and Lennard-Jones) based on standard chemical substructure queries implemented via the industry-standard SMARTS chemical perception language, using SMIRKS extensions that permit labeling of specific atoms within a chemical pattern. We use this to implement a new force field format, called the SMIRKS Native Open Force Field (SMIRNOFF) format. We demonstrate the power and generality of this approach using examples of specific molecules that pose problems for indirect chemical perception and construct and validate a minimalist yet very general force field, SMIRNOFF99Frosst. We find that a parameter definition file only ∼300 lines long provides coverage of all but <0.02% of a 5 million molecule drug-like test set. Despite its simplicity, the accuracy of SMIRNOFF99Frosst for small molecule hydration free energies and selected properties of pure organic liquids is similar to that of the General Amber Force Field, whose specification requires thousands of parameters. This force field provides a starting point for further optimization and refitting work to follow.
| S-EPMC6245550 | biostudies-literature
Project description:DFT calculations have been used to develop improved descriptions of the torsional energetics for nucleosides and nucleotides in the OPLS-AA force field. Scans of nucleotide dihedral angles (γ, χ, and β) and methyl phosphates provided the bases for the new torsional parameters. In addition, the angle-bending parameters of phosphodiesters and ribose were updated, and adjustments were made to existing carbohydrate torsions to better capture the sugar puckering landscape of ribose. MD simulations of nucleosides with the new parameters demonstrate a significant improvement in the ribose sugar puckering and χ angle distributions. Additionally, energy-minimization of protein-nucleotide crystal structures with the new parameters produced accurate poses.
| S-EPMC5821430 | biostudies-literature
Project description:F(o)F(1)-ATPase is a rotary motor protein synthesizing ATP from ADP driven by a cross-membrane proton gradient. The proton flow through the membrane-embedded F(o) generates the rotary torque that drives the rotation of the asymmetric shaft of F(1). Mechanical energy of the rotating shaft is used by the F(1) catalytic subunit to synthesize ATP. It was suggested that elastic power transmission with transient storage of energy in some compliant part of the shaft is required for the observed high turnover rate. We used atomistic simulations to study the spatial distribution and structural determinants of the F(1) torsional elasticity at the molecular level and to comprehensively characterize the elastic properties of F(1)-ATPase. Our fluctuation analysis revealed an unexpected heterogeneity of the F(1) shaft elasticity. Further, we found that the measured overall torsional moduli of the shaft arise from two distinct contributions, the intrinsic elasticity and the effective potential imposed on the shaft by the catalytic subunit. Separation of these two contributions provided a quantitative description of the coupling between the rotor and the catalytic subunit. This description enabled us to propose a minimal quantitative model of the F(1) energetics along the rotary degrees of freedom near the resting state observed in the crystal structures. As opposed to the usually employed models where the motor mechanical progression is described by a single angular variable, our multidimensional treatment incorporates the spatially inhomogeneous nature of the shaft and its interactions with the stator and offers new insight into the mechanoenzymatics of F(1)-ATPase.
| S-EPMC3088579 | biostudies-literature
Project description:All-atom molecular dynamics simulations have been applied in the recent past to explore the free energetics underlying ion transport processes in biological ion channels. Roux and co-workers, Kuyucak and co-workers, Busath and co-workers, and others have performed rather elegant and extended time scale molecular dynamics simulations using current state-of-the-art fixed-charge (nonpolarizable) force fields to calculate the potential of mean force defining the equilibrium flux of ions through prototypical channels such as gramicidin A. An inescapable conclusion of such studies has been the gross overestimation of the equilibrium free energy barrier, generally predicted to be from 10 to 20 kcal/mol depending on the force field and simulation protocol used in the calculation; this translates to an underestimation of experimentally measurable single channel conductances by several orders of magnitude. Next-generation polarizable force fields have been suggested as possible alternatives for more quantitative predictions of the underlying free energy surface in such systems. (1) Presently, we consider ion permeation energetics in the gramicidin A channel using a novel polarizable force field. Our results predict a peak barrier height of 6 kcal/mol relative to the channel entrance; this is significantly lower than the uncorrected value of 12 kcal/mol for nonpolarizable force fields such as GROMOS and CHARMM27 which do not account for electronic polarization. These results provide promising initial indications substantiating the long-conjectured importance of polarization effects in describing ion-protein interactions in narrow biological channels.
| S-EPMC2818143 | biostudies-literature
Project description:In this manuscript we expand significantly on our earlier communication by investigating the bilayer self-assembly of eight different types of phospholipids in unbiased molecular dynamics (MD) simulations using three widely used all-atom lipid force fields. Irrespective of the underlying force field, the lipids are shown to spontaneously form stable lamellar bilayer structures within 1 microsecond, the majority of which display properties in satisfactory agreement with the experimental data. The lipids self-assemble via the same general mechanism, though at formation rates that differ both between lipid types, force fields and even repeats on the same lipid/force field combination. In addition to zwitterionic phosphatidylcholine (PC) and phosphatidylethanolamine (PE) lipids, anionic phosphatidylserine (PS) and phosphatidylglycerol (PG) lipids are represented. To our knowledge this is the first time bilayer self-assembly of phospholipids with negatively charged head groups is demonstrated in all-atom MD simulations.
| S-EPMC5091658 | biostudies-literature
Project description:As the past century has been characterized by waves of viral pandemics, there is an ever-growing role for molecular simulation-based research. In this study, we utilize all-atom molecular dynamics to simulate an enterovirus-D68 capsid and examine the dependency of viral capsid dynamics and properties on AMBER and CHARMM force fields. Out of the six force fields studied, we note that CHARMM36m and CHARMM36 generate secondary structures that are most consistent with protein structural data and sample the largest conformational space. The ion distribution and radius of gyration of the capsid are similar across all force fields investigated.
| S-EPMC8978650 | biostudies-literature
Project description:Experimental studies on intrinsically disordered and unfolded proteins have shown that in isolation they typically have low populations of secondary structure and exhibit distance scalings suggesting that they are at near-theta-solvent conditions. Until recently, however, all-atom force fields failed to reproduce these fundamental properties of intrinsically disordered proteins (IDPs). Recent improvements by refining against ensemble-averaged experimental observables for polypeptides in aqueous solution have addressed deficiencies including secondary structure bias, global conformational properties, and thermodynamic parameters of biophysical reactions such as folding and collapse. To date, studies utilizing these improved all-atom force fields have mostly been limited to a small set of unfolded or disordered proteins. Here, we present data generated for a diverse library of unfolded or disordered proteins using three progressively improved generations of Amber03 force fields, and we explore how global and local properties are affected by each successive change in the force field. We find that the most recent force field refinements significantly improve the agreement of the global properties such as radii of gyration and end-to-end distances with experimental estimates. However, these global properties are largely independent of the local secondary structure propensity. This result stresses the need to validate force fields with reference to a combination of experimental data providing information about both local and global structure formation.
| S-EPMC7507668 | biostudies-literature