Additive CHARMM force field for naturally occurring modified ribonucleotides.
Ontology highlight
ABSTRACT:
SUBMITTER: Xu Y
PROVIDER: S-EPMC4801715 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform "all-CHARMM" simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems.
| S-EPMC2888302 | biostudies-literature
Project description:We present an extension of the CHARMM hexopyranose monosaccharide additive all-atom force field to enable modeling of glycosidic-linked hexopyranose polysaccharides. The new force field parameters encompass 1→1, 1→2, 1→3, 1→4, and 1→6 hexopyranose glycosidic linkages, as well as O-methylation at the C(1) anomeric carbon, and are developed to be consistent with the CHARMM all-atom biomolecular force fields for proteins, nucleic acids, and lipids. The parameters are developed in a hierarchical fashion using model compounds containing the key atoms in the full carbohydrates, in particular O-methyl-tetrahydropyran and glycosidic-linked dimers consisting of two molecules of tetrahyropyran or one of tetrahydropyran and one of cyclohexane. Target data for parameter optimization include full two-dimensional energy surfaces defined by the Φ/Ψ glycosidic dihedral angles in the disaccharide analogs as determined by quantum mechanical MP2/cc-pVTZ single point energies on MP2/6-31G(d) optimized structures (MP2/cc-pVTZ//MP2/6-31G(d)). In order to achieve balanced, transferable dihedral parameters for the Φ/Ψ glycosidic dihedral angles, surfaces for all possible chiralities at the ring carbon atoms involved in the glycosidic linkages are considered, resulting in over 5000 MP2/cc-pVTZ//MP2/6-31G(d) conformational energies. Also included as target data are vibrational frequencies, pair interaction energies and distances with water molecules, and intramolecular geometries including distortion of the glycosidic valence angle as a function of the glycosidic dihedral angles. The model-compound optimized force field parameters are validated on full disaccharides through comparison of molecular dynamics results to available experimental data. Good agreement is achieved with experiment for a variety of properties including crystal cell parameters and intramolecular geometries, aqueous densities, and aqueous NMR coupling constants associated with the glycosidic linkage. The newly-developed parameters allow for the modeling of linear, branched, and cyclic hexopyranose glycosides both alone and in heterogenous systems including proteins, nucleic acids and/or lipids when combined with existing CHARMM biomolecular force fields.
| S-EPMC2757763 | biostudies-literature
Project description:An additive all-atom empirical force field for aldopentofuranoses, methyl-aldopentofuranosides (Me-aldopentofuranosides), and fructofuranose carbohydrates, compatible with existing CHARMM carbohydrate parameters, is presented. Building on existing parameters transferred from cyclic ethers and hexopyranoses, parameters were further developed using target data for complete furanose carbohydrates as well as O-methyl tetrahydrofuran. The bond and angle equilibrium parameters were adjusted to reproduce target geometries from a survey of furanose crystal structures, and dihedral parameters were fit to over 1700 quantum mechanical (QM) MP2/cc-pVTZ//MP2/6-31G(d) conformational energies. The conformational energies were for a variety of complete furanose monosaccharides and included two-dimensional ring pucker energy surfaces. Bonded parameter optimization led to the correct description of the ring pucker for a large set of furanose compounds, while furanose-water interaction energies and distances reproduced QM HF/6-31G(d) results for a number of furanose monosaccharides, thereby validating the nonbonded parameters. Crystal lattice unit cell parameters and volumes, aqueous-phase densities, and aqueous NMR ring pucker and exocyclic data were used to validate the parameters in condensed-phase environments. Conformational sampling analysis of the ring pucker and exocyclic group showed excellent agreement with experimental NMR data, demonstrating that the conformational energetics in aqueous solution are accurately described by the optimized force field. Overall, the parameters reproduce available experimental data well and are anticipated to be of utility in future computational studies of carbohydrates, including in the context of proteins, nucleic acids, and/or lipids when combined with existing CHARMM biomolecular force fields.
| S-EPMC2741538 | biostudies-literature
Project description:Presented is an extension of the CHARMM additive carbohydrate all-atom force field to enable modeling of polysaccharides containing furanose sugars. The new force field parameters encompass 1 ? 2, 1 ? 3, 1 ? 4, and 1 ? 6 pyranose-furanose linkages and 2 ? 1 and 2 ? 6 furanose-furanose linkages, building on existing hexopyranose and furanose monosaccharide parameters. The model compounds were chosen to be monomers or glycosidic-linked dimers of tetrahydropyran (THP) and tetrahydrofuran (THF) as to contain the key atoms in full carbohydrates. Target data for optimization included two-dimensional quantum mechanical (QM) potential energy scans of the ?/? glycosidic dihedral angles, with geometry optimization at the MP2/6-31G(d) level followed by MP2/cc-pVTZ single-point energies. All possible chiralities of the model compounds at the linkage carbons were considered, and for each geometry, the THF ring was constrained to the favorable south or north conformations. Target data also included QM vibrational frequencies and pair interaction energies and distances with water molecules. Force field validation included comparison of computed crystal properties, aqueous solution densities, and NMR J-coupling constants to experimental reference values. Simulations of infinite crystals showed good agreement with experimental values for intramolecular geometries as well as for crystal unit cell parameters. Additionally, aqueous solution densities and available NMR data were reproduced to a high degree of accuracy, thus validating the hierarchically optimized parameters in both crystalline and aqueous condensed phases. The newly developed parameters allow for the modeling of linear, branched, and cyclic pyranose/furanose polysaccharides both alone and in heterogeneous systems including proteins, nucleic acids, and/or lipids when combined with existing additive CHARMM biomolecular force fields.
| S-EPMC2958709 | biostudies-literature
Project description:Presented is an extension of the CHARMM additive all-atom carbohydrate force field to enable the modeling of phosphate and sulfate linked to carbohydrates. The parameters are developed in a hierarchical fashion using model compounds containing the key atoms in the full carbohydrates. Target data for parameter optimization included full two-dimensional energy surfaces defined by the glycosidic dihedral angle pairs in the phosphate/sulfate model compound analogs of hexopyranose monosaccharide phosphates and sulfates, as determined by quantum mechanical (QM) MP2/cc-pVTZ single point energies on MP2/6-31+G(d) optimized structures. In order to achieve balanced, transferable dihedral parameters for the dihedral angles, surfaces for all possible anomeric and conformational states were included during the parametrization process. In addition, to model physiologically relevant systems both the mono- and di-anionic charged states were studied for the phosphates. This resulted in over 7000 MP2/cc-pVTZ//MP2/6-31G+(d) model compound conformational energies which, supplemented with QM geometries, were the main target data for the parametrization. Parameters were validated against crystals of relevant monosaccharide derivatives obtained from the Cambridge Structural Database (CSD) and larger systems, namely inositol-(tri/tetra/penta) phosphates non-covalently bound to the pleckstrin homology (PH) domain and oligomeric chondroitin sulfate in solution and in complex with cathepsin K protein.
| S-EPMC3367516 | biostudies-literature
Project description:CHARMM is one of the most widely used biomolecular force fields. Although developed in close connection with a dedicated molecular simulation engine of the same name, it is also usable with other codes. GROMACS is a well-established, highly optimized, and multipurpose software for molecular dynamics, versatile enough to accommodate many different force field potential functions and the associated algorithms. Due to conceptional differences related to software design and the large amount of numeric data inherent to residue topologies and parameter sets, conversion from one software format to another is not straightforward. Here, we present an automated and validated means to port the CHARMM force field to a format read by the GROMACS engine, harmonizing the different capabilities of the two codes in a self-documenting and reproducible way with a bare minimum of user interaction required. Being based entirely on the upstream data files, the presented approach does not involve any hard-coded data, in contrast with previous attempts to solve the same problem. The heuristic approach used for perceiving the local internal geometry is directly applicable for analogous transformations of other force fields.
| S-EPMC10369483 | biostudies-literature
Project description:Parametrization of the additive all-atom CHARMM force field for acyclic polyalcohols, acyclic carbohydrates and inositol is conducted. Initial parameters were transferred from the alkanes and hexopyranose carbohydrates, with subsequent development and optimization of parameters unique to the molecules considered in this study. Using the model compounds acetone and acetaldehyde, nonbonded parameters for carbonyls were optimized targeting quantum mechanical interaction data for solute-water pairs and pure solvent thermodynamic data. Bond and angle parameters were adjusted by comparing optimized geometries to small molecule crystal survey data and by performing vibrational analyses on acetone, acetaldehyde and glycerol. C-C-C-C, C-C-C-O, C-C-OH and O-C-C-O torsional parameters for polyol chains were fit to quantum mechanical dihedral potential energy scans comprising over 1500 RIMP2/cc-pVTZ//MP2/6-31G(d) conformations using an automated Monte Carlo simulated annealing procedure. Comparison of computed condensed-phase data, including crystal lattice parameters and densities, NMR proton-proton couplings, densities and diffusion coefficients of aqueous solutions, to experimental data validated the optimized parameters. Parameter development for these compounds proved particularly challenging because of the flexibility of the acyclic sugars and polyalcohols as well as the intramolecular hydrogen bonding between vicinal hydroxyls for all of the compounds. The newly optimized additive CHARMM force field parameters are anticipated to be of utility for atomic level of detail simulations of acyclic polyalcohols, acyclic carbohydrates and inositol in solution.
| S-EPMC2760998 | biostudies-literature
Project description:Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.
| S-EPMC4712441 | biostudies-literature
Project description:Accurate force-field (FF) parameters are key to reliable prediction of properties obtained from molecular modeling (MM) and molecular dynamics (MD) simulations. With ever-widening applicability of MD simulations, robust parameters need to be generated for a wider range of chemical species. The CHARMM General Force Field program (CGenFF, https://cgenff.umaryland.edu/) is a tool for obtaining initial parameters for a given small molecule based on analogy with the available CGenFF parameters. However, improvement of these parameters is often required and performing their optimization remains tedious and time consuming. In addition, tools for optimization of small molecule parameters in the context of the Drude polarizable FF are not yet available. To overcome these issues, the FFParam package has been designed to facilitate the parametrization process. The package includes a graphical user interface (GUI) created using Qt libraries. FFParam supports Gaussian and Psi4 for performing quantum mechanical calculations and CHARMM and OpenMM for MM calculations. A Monte Carlo simulated annealing (MCSA) algorithm has been implemented for automated fitting of partial atomic charge, atomic polarizabilities and Thole scale parameters. The LSFITPAR program is called for automated fitting of bonded parameters. Accordingly, FFParam provides all the features required for generation and analysis of CHARMM and Drude FF parameters for small molecules. FFParam-GUI includes a text editor, graph plotter, molecular visualization, and text to table converter to meet various requirements of the parametrization process. It is anticipated that FFParam will facilitate wider use of CGenFF as well as promote future use of the Drude polarizable FF.
| S-EPMC7323454 | biostudies-literature
Project description:The C36 CHARMM lipid force field has been extended to include sphingolipids, via a combination of high-level quantum mechanical calculations on small molecule fragments, and validation by extensive molecular dynamics simulations on N-palmitoyl and N-stearoyl sphingomyelin. NMR data on these two molecules from several studies in bilayers and micelles played a strong role in the development and testing of the force field parameters. Most previous force fields for sphingomyelins were developed before the availability of the detailed NMR data and relied on x-ray diffraction of bilayers alone for the validation; these are shown to be too dense in the bilayer plane based on published chain order parameter data from simulations and experiments. The present simulations reveal O-H:::O-P intralipid hydrogen bonding occurs 99% of the time, and interlipid N-H:::O=C (26-29%, depending on the lipid) and N-H:::O-H (17-19%). The interlipid hydrogen bonds are long lived, showing decay times of 50 ns, and forming strings of lipids, and leading to reorientational correlation time of nearly 100 ns. The spontaneous radius of curvature for pure N-palmitoyl sphingomyelin bilayers is estimated to be 43-100 Å, depending on the assumptions made in assigning a bending constant; this unusual positive curvature for a two-tailed neutral lipid is likely associated with hydrogen bond networks involving the NH of the sphingosine group.
| S-EPMC4119286 | biostudies-literature