Impact of 2'-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA

Here, we present an update of the CHARMM27 all-atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson-Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2'-hydroxyl moiety to almost exclusively sample the O3' orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3' orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2'-hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3' orientations of the 2'-hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J-coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2'-hydroxyl moiety and support a model whereby the 2'-hydroxyl can enhance the probability of conformational transitions in RNA.     

CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields

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. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

OPLS all-atom force field for carbohydrates

The OPLS all-atom (AA) force field for organic and biomolecular systems has been expanded to include carbohydrates. Starting with reported nonbonded parameters of alcohols, ethers, and diols, torsional parameters were fit to reproduce results from ab initio calculations on the hexopyranoses, α,β-d -glucopyranose, α,β-d -mannopyranose, α,β-d -galactopyranose, methyl α,β-d -glucopyranoside, and methyl α,β-d -mannopyranoside. In all, geometry optimizations were carried out for 144 conformers at the restricted Hartree–Fock (RHF)/6–31G* level. For the conformers with a relative energy within 3 kcal/mol of the global minima, the effects of electron correlation and basis-set extension were considered by performing single-point calculations with density functional theory at the B3LYP/6–311+G** level. The torsional parameters for the OPLS-AA force field were parameterized to reproduce the energies and structures of these 44 conformers. The resultant force field reproduces the ab initio calculated energies with an average unsigned error of 0.41 kcal/mol. The α/β ratios as well as the relative energies between the isomeric hexopyranoses are in good accord with the ab initio results. The predictive abilities of the force field were also tested against RHF/6–31G* results for d -allopyranose with excellent success; a surprising discovery is that the lowest energy conformer of d -allopyranose is a β anomer. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1955–1970, 1997     

An all-atom force field for metallocenes

A new all-atom force field, for the molecular modeling of metallocenes was constructed. Quantum chemical calculations were performed to obtain several force field terms not yet defined in the literature. The remainder were transferred from the OPLS-AA/AMBER framework. The parametrization work included the obtention of geometrical parameters, torsion energy profiles, and distributions of atomic charges that blend smoothly with the OPLS-AA specification for a variety of organic molecular fragments. Validation was carried out by comparing simulated and experimental data for five different ferrocene derivatives in the crystalline phase. The present model can be regarded as a step toward a general force field for metallocenes, built in a systematic way, easily integrated with OPLS-AA, and transferable between different metal-ligand combinations.     

FFParam: Standalone package for CHARMM additive and Drude polarizable force field parametrization of small molecules

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.     

Extension of the CHARMM general force field to sulfonyl‐containing compounds and its utility in biomolecular simulations

Presented is an extension of the CHARMM General Force Field (CGenFF) to enable the modeling of sulfonyl‐containing compounds. Model compounds containing chemical moieties such as sulfone, sulfonamide, sulfonate, and sulfamate were used as the basis for the parameter optimization. Targeting high‐level quantum mechanical and experimental crystal data, the new parameters were optimized in a hierarchical fashion designed to maintain compatibility with the remainder of the CHARMM additive force field. The optimized parameters satisfactorily reproduced equilibrium geometries, vibrational frequencies, interactions with water, gas phase dipole moments, and dihedral potential energy scans. Validation involved both crystalline and liquid phase calculations showing the newly developed parameters to satisfactorily reproduce experimental unit cell geometries, crystal intramolecular geometries, and pure solvent densities. The force field was subsequently applied to study conformational preference of a sulfonamide based peptide system. Good agreement with experimental IR/NMR data further validated the newly developed CGenFF parameters as a tool to investigate the dynamic behavior of sulfonyl groups in a biological environment. CGenFF now covers sulfonyl group containing moieties allowing for modeling and simulation of sulfonyl‐containing compounds in the context of biomolecular systems including compounds of medicinal interest. © 2012 Wiley Periodicals, Inc.     

Additive CHARMM force field for naturally occurring modified ribonucleotides

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Modification of the CHARMM force field for DMPC lipid bilayer

The CHARMM force field for DMPC lipids was modified in order to improve agreement with experiment for a number of important properties of hydrated lipid bilayer. The modification consists in introduction of a scaling factor 0.83 for 1-4 electrostatic interactions (between atoms separated by three covalent bonds), which provides correct transgauche ratio in the alkane tails, and recalculation of the headgroup charges on the basis of HF/6-311(d,p) ab-initio computations. Both rigid TIP3P and flexible SPC water models were used with the new lipid model, showing similar results. The new model in a 75 ns simulation has shown a correct value of the area per lipid at zero surface tension, as well as good agreement with the experiment for the electron density, structure factor, and order parameters, including those in the headgroup part of lipids.     

Optimization of CHARMM force field parameters for the chalcone fragment

In this work,we developed the CHARMM all-atom force field parameters for the nonstandard biological residue chalcone,followed by the standard protocol for the CHARMM27 force field development.Target data were generated via ab initio calculations at the MP2/6-31G* and HF/6-31G* levels.The reference data included interaction energies between water and the model compound F(a fragment of chalcone).Bond,angle,and torsion parameters were derived from the ab initio calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters of standard residues.The optimized CHARMM parameters perform well in reproducing the target data.We expect that the extension of the CHARMM27 force field parameters for chalcone will facilitate the molecular simulation studies of the reaction mechanism of intramolecular cyclization of chalcone catalyzed by chalcone isomerase.     

Optimization of the CHARMM additive force field for DNA: Improved treatment of the BI/BII conformational equilibrium

The B-form of DNA can populate two different backbone conformations: BI and BII, defined by the difference between the torsion angles ε and ζ (BI = ε-ζ < 0 and BII = ε-ζ > 0). BI is the most populated state, but the population of the BII state, which is sequence dependent, is significant and accumulating evidence shows that BII affects the overall structure of DNA, and thus influences protein-DNA recognition. This work presents a reparametrization of the CHARMM27 additive nucleic acid force field to increase the sampling of the BII form in MD simulations of DNA. In addition, minor modifications of sugar puckering were introduced to facilitate sampling of the A form of DNA under the appropriate environmental conditions. Parameter optimization was guided by quantum mechanical data on model compounds, followed by calculations on several DNA duplexes in the condensed phase. The selected optimized parameters were then validated against a number of DNA duplexes, with the most extensive tests performed on the EcoRI dodecamer, including comparative calculations using the Amber Parm99bsc0 force field. The new CHARMM model better reproduces experimentally observed sampling of the BII conformation, including sampling as a function of sequence. In addition, the model reproduces the A form of the 1ZF1 duplex in 75 % ethanol, and yields a stable Z-DNA conformation of duplex (GTACGTAC) in its crystal environment. The resulting model, in combination with a recent reoptimization of the CHARMM27 force field for RNA, will be referred to as CHARMM36.     

CHARMM force field and molecular dynamics simulations of protonated polyethylenimine

As a gene delivery vector, polyethylenimine (PEI) shows one of the highest transfection efficiencies, while effectively protecting DNA from enzyme degradation. The distinctive charge pattern of protonated PEI is widely considered responsible for fundamental process such as DNA condensation into PEI/DNA polyplexes (which are able to enter cells via endocytosis), proton sponge effect (which triggers the release of polyplexes from endosome), and release of DNA from polyplexes (to be further processed inside the nucleus). Our investigations are largely motivated by the crucial need for a realistic molecular mechanics force field (FF) for PEI, and, accordingly, we focus on two major issues: (1) development of a new atomistic (CHARMM) FF for PEI in different protonation states, rigorously derived from high‐quality ab initio calculations performed on model polymers, and (2) molecular dynamics investigations of solvated PEI, providing a detailed picture of the dynamic structuring thereof in dependence on their size and protonation state. The modeled PEI chains are essentially described in terms of gyration radius, end‐to‐end distance, persistence length, radial distribution functions, coordination numbers, and diffusion coefficients. They turn out to be more rigid than in other computational studies and we find diffusion coefficients in fair agreement with experimental data. The developed atomistic FF proves adequate for the realistic modeling of the size and protonation behavior of linear PEI, either as individual chains or composing polyplexes. © 2017 Wiley Periodicals, Inc.     

CHARMM36 all‐atom additive protein force field: Validation based on comparison to NMR data

Protein structure and dynamics can be characterized on the atomistic level with both nuclear magnetic resonance (NMR) experiments and molecular dynamics (MD) simulations. Here, we quantify the ability of the recently presented CHARMM36 (C36) force field (FF) to reproduce various NMR observables using MD simulations. The studied NMR properties include backbone scalar couplings across hydrogen bonds, residual dipolar couplings (RDCs) and relaxation order parameter, as well as scalar couplings, RDCs, and order parameters for side‐chain amino‐ and methyl‐containing groups. It is shown that the C36 FF leads to better correlation with experimental data compared to the CHARMM22/CMAP FF and suggest using C36 in protein simulations. Although both CHARMM FFs contains the same nonbond parameters, our results show how the changes in the internal parameters associated with the peptide backbone via CMAP and the χ₁ and χ₂ dihedral parameters leads to improved treatment of the analyzed nonbond interactions. This highlights the importance of proper treatment of the internal covalent components in modeling nonbond interactions with molecular mechanics FFs. © 2013 Wiley Periodicals, Inc.     

CHARMM force field parameters for simulation of reactive intermediates in native and thio-substituted ribozymes

Force field parameters specifically optimized for residues important in the study of RNA catalysis are derived from density-functional calculations, in a fashion consistent with the CHARMM27 all-atom empirical force field. Parameters are presented for residues that model reactive RNA intermediates and transition state analogs, thio-substituted phosphates and phosphoranes, and bound Mg(2+) and di-metal bridge complexes. Target data was generated via density-functional calculations at the B3LYP/6-311++G(3df,2p)// B3LYP/6-31++G(d,p) level. Partial atomic charges were initially derived from CHelpG electrostatic potential fitting and subsequently adjusted to be consistent with the CHARMM27 charges. Lennard-Jones parameters were determined to reproduce interaction energies with water molecules. Bond, angle, and torsion parameters were derived from the density-functional calculations and renormalized to maintain compatibility with the existing CHARMM27 parameters for standard residues. The extension of the CHARMM27 force field parameters for the nonstandard biological residues presented here will have considerable use in simulations of ribozymes, including the study of freeze-trapped catalytic intermediates, metal ion binding and occupation, and thio effects.     

Ab initio parameterization of an all-atom polarizable and dissociable force field for water

A novel all-atom, dissociative, and polarizable force field for water is presented. The force field is parameterized based on forces, stresses, and energies obtained form ab initio calculations of liquid water at ambient conditions. The accuracy of the force field is tested by calculating structural and dynamical properties of liquid water and the energetics of small water clusters. The transferability of the force field to dissociated states is studied by considering the solvation of a proton and the ionization of water at extreme conditions of pressure and temperature. In the case of the solvated proton, the force field properly describes the presence of both Eigen and Zundel configurations. In the case of the pressure-induced ice VIII/ice X transition and the temperature-induced transition to a superionic phase, the force field is found to describe accurately the proton symmetrization and the melting of the proton sublattice, respectively.     

Revised CHARMM force field parameters for iron‐containing cofactors of photosystem II

Photosystem II is a complex protein–cofactor machinery that splits water molecules into molecular oxygen, protons, and electrons. All‐atom molecular dynamics simulations have the potential to contribute to our general understanding of how photosystem II works. To perform reliable all‐atom simulations, we need accurate force field parameters for the cofactor molecules. We present here CHARMM bonded and non‐bonded parameters for the iron‐containing cofactors of photosystem II that include a six‐coordinated heme moiety coordinated by two histidine groups, and a non‐heme iron complex coordinated by bicarbonate and four histidines. The force field parameters presented here give water interaction energies and geometries in good agreement with the quantum mechanical target data. © 2017 Wiley Periodicals, Inc.     

CHARMM‐GUI ligand reader and modeler for CHARMM force field generation of small molecules

Reading ligand structures into any simulation program is often nontrivial and time consuming, especially when the force field parameters and/or structure files of the corresponding molecules are not available. To address this problem, we have developed Ligand Reader & Modeler in CHARMM‐GUI. Users can upload ligand structure information in various forms (using PDB ID, ligand ID, SMILES, MOL/MOL2/SDF file, or PDB/mmCIF file), and the uploaded structure is displayed on a sketchpad for verification and further modification. Based on the displayed structure, Ligand Reader & Modeler generates the ligand force field parameters and necessary structure files by searching for the ligand in the CHARMM force field library or using the CHARMM general force field (CGenFF). In addition, users can define chemical substitution sites and draw substituents in each site on the sketchpad to generate a set of combinatorial structure files and corresponding force field parameters for throughput or alchemical free energy simulations. Finally, the output from Ligand Reader & Modeler can be used in other CHARMM‐GUI modules to build a protein‐ligand simulation system for all supported simulation programs, such as CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Ligand Reader & Modeler is available as a functional module of CHARMM‐GUI at http://www.charmm-gui.org/input/ligandrm . © 2017 Wiley Periodicals, Inc.     

Empirical force field calculations on methylammonium chloride and its deuterated derivatives

A valence force field for the methylammonium ion is presented. The force field was determined from a modified least squares fit of the observed Raman frequencies of methylammonium chloride and three of its deuterated derivatives. In addition, a Urey—Bradley force field was tested; it was found to be impossible to obtain a good fit to the observed frequencies without the inclusion of a non-bonded hydrogen—hydrogen interaction in the potential energy. The reason is that the unmodified Urey—Bradley force field does not include any bend—bend interactions, which are necessary in order to reproduce the rocking modes.     

CHARMM fluctuating charge force field for proteins: I parameterization and application to bulk organic liquid simulations

A first-generation fluctuating charge (FQ) force field to be ultimately applied for protein simulations is presented. The electrostatic model parameters, the atomic hardnesses, and electronegativities, are parameterized by fitting to DFT-based charge responses of small molecules perturbed by a dipolar probe mimicking a water dipole. The nonbonded parameters for atoms based on the CHARMM atom-typing scheme are determined via simultaneously optimizing vacuum water-solute geometries and energies (for a set of small organic molecules) and condensed phase properties (densities and vaporization enthalpies) for pure bulk liquids. Vacuum solute-water geometries, specifically hydrogen bond distances, are fit to 0.19 A r.m.s. error, while dimerization energies are fit to 0.98 kcal/mol r.m.s. error. Properties of the liquids studied include bulk liquid structure and polarization. The FQ model does indeed show a condensed phase effect in the shifting of molecular dipole moments to higher values relative to the gas phase. The FQ liquids also appear to be more strongly associated, in the case of hydrogen bonding liquids, due to the enhanced dipolar interactions as evidenced by shifts toward lower energies in pair energy distributions. We present results from a short simulation of NMA in bulk TIP4P-FQ water as a step towards simulating solvated peptide/protein systems. As expected, there is a nontrivial dipole moment enhancement of the NMA (although the quantitative accuracy is difficult to assess). Furthermore, the distribution of dipole moments of water molecules in the vicinity of the solutes is shifted towards larger values by 0.1-0.2 Debye in keeping with previously reported work.     

Modeling the effect of ion-induced shock waves and DNA breakage with the reactive CHARMM force field

Ion-induced DNA damage is an important effect underlying ion beam cancer therapy. This article introduces the methodology of modeling DNA damage induced by a shock wave caused by a projectile ion. Specifically it is demonstrated how single- and double strand breaks in a DNA molecule could be described by the reactive CHARMM (rCHARMM) force field implemented in the program MBN Explorer. The entire workflow of performing the shock wave simulations, including obtaining the crucial simulation parameters, is described in seven steps. Two exemplary analyses are provided for a case study simulation serving to: (a) quantify the shock wave propagation and (b) describe the dynamics of formation of DNA breaks. The article concludes by discussing the computational cost of the simulations and revealing the possible maximal computational time for different simulation set-ups.     

Molecular dynamics effects on luminescence properties of oligothiophene derivatives: a molecular mechanics-response theory study based on the CHARMM force field and density functional theory

CHARMM force field parameter values for a class of oligothiophene derivatives have been derived with reference to density functional theory/B3LYP potential energy surfaces. The force field parametrization of these luminescent conjugated polyelectrolytes includes the electronic ground state as well as the strongly light absorbing first excited state. In conjunction with quantum chemical response theory calculations of transition state properties, a molecular dynamical model of the Stokes shift is obtained. The theoretical model is benchmarked against experimental data recorded at room temperature which refer to sodium salts of p-HTAA and p-FTAA with distilled water as a solvent. For p-HTAA the theoretically predicted Stokes shift of 112 nm is in good agreement with the experimental result of 124 nm, given the approximations about exciton localization that were introduced to obtain a force field for the excited state.