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.     

Force field dependence of phospholipid headgroup and acyl chain properties: Comparative molecular dynamics simulations of DMPC bilayers

The reliability of molecular simulations largely depends on the quality of the empirical force field parameters. Force fields used in lipid simulations continue to be improved to enhance the agreement with experiments for a number of different properties. In this work, we have carried out molecular dynamics simulations of neat DMPC bilayers using united‐atom Berger force field and three versions of all‐atom CHARMM force fields. Three different systems consisting of 48, 72, and 96 lipids were studied. Both particle mesh Ewald (PME) and spherical cut‐off schemes were used to evaluate the long‐range electrostatic interactions. In total, 21 simulations were carried out and analyzed to find out the dependence of lipid properties on force fields, system size, and schemes to calculate long‐range interactions. The acyl chain order parameters calculated from Berger and the recent versions of CHARMM simulations have shown generally good agreement with the experimental results. However, both sets of force fields deviate significantly from the experimentally observed P‐C dipolar coupling values for the carbon atoms that link the choline and glycerol groups with the phosphate groups. Significant differences are also observed in several headgroup parameters between CHARMM and Berger simulations. Our results demonstrate that when changes were introduced to improve CHARMM force field using PME scheme, all the headgroup parameters have not been reoptimized. The headgroup properties are likely to play a significant role in lipid–lipid, protein–lipid, and ligand–lipid interactions and hence headgroup parameters in phospholipids require refinement for both Berger and CHARMM force fields. © 2009 Wiley Periodicals, Inc.J Comput Chem, 2010     

Evaluations of AMBER force field parameters by MINA approach for copper-based nucleases

We describe the development of the AMBER force field parameters for 46 nucleases involving most kinds of copper nucleases with high DNA affinities and specificities by MINA approach that could evaluate accurate force constants for batch bonds/angles on the basis of energies of three adjacent lengths/angles geometries. The molecular mechanics (MM) and molecular dynamic simulations on adducts of the 21 representative copper-based nucleases with DNA are in excellent agreement with those of experimental results. Furthermore, to validate the evaluated parameters, the studied structures performed frequency analysis together with normal mode calculations in quantum mechanics and MM calculations. The force field parameters evaluated in this work provide an extension of AMBER force field, and the results of molecular dynamics simulations of adduct of copper nuclease and duplex DNA illustrate the potential utility of these parameters.     

A molecular mechanics force field for lignin

A CHARMM molecular mechanics force field for lignin is derived. Parameterization is based on reproducing quantum mechanical data of model compounds. Partial atomic charges are derived using the RESP electrostatic potential fitting method supplemented by the examination of methoxybenzene:water interactions. Dihedral parameters are optimized by fitting to critical rotational potentials and bonded parameters are obtained by optimizing vibrational frequencies and normal modes. Finally, the force field is validated by performing a molecular dynamics simulation of a crystal of a lignin fragment molecule and comparing simulation-derived structural features with experimental results. Together with the existing force field for polysaccharides, this lignin force field will enable full simulations of lignocellulose.     

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.     

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.     

MATCH: an atom-typing toolset for molecular mechanics force fields

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.     

Development of the force field parameters for phosphoimidazole and phosphohistidine

Phosphorylation of histidine-containing proteins is a key step in the mechanism of many phosphate transfer enzymes (kinases, phosphatases) and is the first stage in a wide variety of signal transduction cascades in bacteria, yeast, higher plants, and mammals. Studies of structural and dynamical aspects of such enzymes in the phosphorylated intermediate states are important for understanding the intimate molecular mechanisms of their functioning. Such information may be obtained via molecular dynamics and/or docking simulations, but in this case appropriate force field parameters for phosphohistidine should be explicitly defined. In the present article we describe development of the GROMOS96 force field parameters for phosphoimidazole molecule--a realistic model of the phosphohistidine side chain. The parameterization is based on the results of ab initio quantum chemical calculations with subsequent refinement and testing using molecular mechanics and molecular dynamics simulations. The set of force constants and equilibrium geometry is employed to derive force field for the phosphohistidine moiety. Resulting parameters and topology are incorporated into the molecular modeling package GROMACS and used in molecular dynamics simulations of a phosphohistidine-containing protein in explicit solvent.     

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     

Simulation of multiphase systems utilizing independent force fields to control intraphase and interphase behavior

Fixed‐charge empirical force fields have been developed and widely used over the past three decades for all‐atom molecular simulations. Most simulation programs providing these methods enable only one set of force field parameters to be used for the entire system. Whereas this is generally suitable for single‐phase systems, the molecular environment at the interface between two phases may be sufficiently different from the individual phases to require a different set of parameters to be used to accurately represent the system. Recently published simulations of peptide adsorption to material surfaces using the CHARMM force field have clearly demonstrated this issue by revealing that calculated values of adsorption free energy substantially differ from experimental results. Whereas nonbonded parameters could be adjusted to correct this problem, this cannot be done without also altering the conformational behavior of the peptide in solution, for which CHARMM has been carefully tuned. We have developed a dual‐force‐field approach (Dual‐FF) to address this problem and implemented it in the CHARMM simulation package. This Dual‐FF method provides the capability to use two separate sets of nonbonded force field parameters within the same simulation: one set to represent intraphase interactions and a separate set to represent interphase interactions. Using this approach, we show that interfacial parameters can be adjusted to correct errors in peptide adsorption free energy without altering peptide conformational behavior in solution. This program thus provides the capability to enable both intraphase and interphase molecular behavior to be accurately and efficiently modeled in the same simulation. © 2012 Wiley Periodicals, Inc.     

The QM‐MM interface for CHARMM‐deMon

We present a new QM/MM interface for fast and efficient simulations of organic and biological molecules. The CHARMM/deMon interface has been developed and tested to perform minimization and atomistic simulations for multi‐particle systems. The current features of this QM/MM interface include readability for molecular dynamics, tested compatibility with Free Energy Perturbation simulations (FEP) using the dual topology/single coordinate method. The current coupling scheme uses link atoms, but further extensions of the code to incorporate other available schemes are planned. We report the performance of different levels of theory for the treatment of the QM region, while the MM region was represented by a classical force‐field (CHARMM27) or a polarizable force‐field based on a simple Drude model. The current QM/MM implementation can be coupled to the dual‐thermostat method and the VV2 integrator to run molecular dynamics simulations. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Force field development for cofactors in the photosystem II

We present a set of force field (FF) parameters compatible with the AMBER03 FF to describe five cofactors in photosystem II (PSII) of oxygenic photosynthetic organisms: plastoquinone‐9 (three redox forms), chlorophyll‐a, pheophytin‐a, heme‐b, and β‐carotene. The development of a reliable FF for these cofactors is an essential step for performing molecular dynamics simulations of PSII. Such simulations are important for the calculation of absorption spectrum and the further investigation of the electron and energy transfer processes. We have derived parameters for partial charges, bonds, angles, and dihedral‐angles from solid theoretical models using systematic quantum mechanics (QM) calculations. We have shown that the developed FF parameters are in good agreement with both ab initio QM and experimental structural data in small molecule crystals as well as protein complexes. © 2012 Wiley Periodicals, Inc.     

Towards in Silico Liquid Crystals. Realistic Transition Temperatures and Physical Properties for n‐Cyanobiphenyls via Molecular Dynamics Simulations

We study the important n‐cyanobiphenyl (with n= 4–8) series of mesogens, using modelling and molecular dynamics simulations. We are able to obtain spontaneously ordered nematics upon cooling isotropic samples of 250 molecules. By using the united‐atom force field developed herein, we show that the experimental isotropic–nematic transition temperatures are reproduced within 4 K, allowing a molecular‐level interpretation of the odd–even effect along the series. Other properties, like densities, orientational order parameters and NMR residual dipolar couplings are also reproduced well, demonstrating the feasibility of predictive in silico modelling of nematics from the molecular structure.     

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.     

A potential with low point charges for pure siliceous zeolites

A modified CHARMM force‐field (ZHB potential) with low point charges for silica was previously proposed by Zimmerman et al. (J. Chem. Theory Comput. 2011, 7, 1695). The ZHB potential is advantageous for quantum mechanics/molecular mechanics simulations as it minimizes the electron spill‐out problems. In the same spirit, here we propose a modified ZHB potential (MZHB) by reformulating its bonding potential, while retaining the nonbonding potential as in the ZHB force‐field. We show that several structural and dynamic properties of silica, like the IR spectrum, distribution functions, mechanical properties, and negative thermal expansion computed using the MZHB potential agree well with experimental data. Further, transferability of MZHB is also tested for reproducing the crystallographic structures of several polymorphs of silica. © 2015 Wiley Periodicals, Inc.     

Force-field development and molecular dynamics simulations of ferrocene-peptide conjugates as a scaffold for hydrogenase mimics

The increasing importance of hydrogenase enzymes in the new energy research field has led us to examine the structure and dynamics of potential hydrogenase mimics, based on a ferrocene-peptide scaffold, using molecular dynamics (MD) simulations. To enable this MD study, a molecular mechanics force field for ferrocene-bearing peptides was developed and implemented in the CHARMM simulation package, thus extending the usefulness of the package into peptide-bioorganometallic chemistry. Using the automated frequency-matching method (AFMM), optimized intramolecular force-field parameters were generated through quantum chemical reference normal modes. The partial charges for ferrocene were derived by fitting point charges to quantum-chemically computed electrostatic potentials. The force field was tested against experimental X-ray crystal structures of dipeptide derivatives of ferrocene-1,1'-dicarboxylic acid. The calculations reproduce accurately the molecular geometries, including the characteristic C2-symmetrical intramolecular hydrogen-bonding pattern, that were stable over 0.1 micros MD simulations. The crystal packing properties of ferrocene-1-(D)alanine-(D)proline-1'-(D)alanine-(D)proline were also accurately reproduced. The lattice parameters of this crystal were conserved during a 0.1 micros MD simulation and match the experimental values almost exactly. Simulations of the peptides in dichloromethane are also in good agreement with experimental NMR and circular dichroism (CD) data in solution. The developed force field was used to perform MD simulations on novel, as yet unsynthesized peptide fragments that surround the active site of [Ni-Fe] hydrogenase. The results of this simulation lead us to propose an improved design for synthetic peptide-based hydrogenase models. The presented MD simulation results of metallocenes thereby provide a convincing validation of our proposal to use ferrocene-peptides as minimal enzyme mimics.     

Polarizable empirical force field for sulfur‐containing compounds based on the classical Drude oscillator model

Condensed‐phase computational studies of molecules using molecular mechanics approaches require the use of force fields to describe the energetics of the systems as a function of structure. The advantage of polarizable force fields over nonpolarizable (or additive) models lies in their ability to vary their electronic distribution as a function of the environment. Toward development of a polarizable force field for biological molecules, parameters for a series of sulfur‐containing molecules are presented. Parameter optimization was performed to reproduce quantum mechanical and experimental data for gas phase properties including geometries, conformational energies, vibrational spectra, and dipole moments as well as for condensed phase properties such as heats of vaporization, molecular volumes, and free energies of hydration. Compounds in the training set include methanethiol, ethanethiol, propanethiol, ethyl methyl sulfide, and dimethyl disulfide. The molecular volumes and heats of vaporization are in good accordance with experimental values, with the polarizable model performing better than the CHARMM22 nonpolarizable force field. Improvements with the polarizable model were also obtained for molecular dipole moments and in the treatment of intermolecular interactions as a function of orientation, in part due to the presence of lone pairs and anisotropic atomic polarizability on the sulfur atoms. Significant advantage of the polarizable model was reflected in calculation of the dielectric constants, a property that CHARMM22 systematically underestimates. The ability of this polarizable model to accurately describe a range of gas and condensed phase properties paves the way for more accurate simulation studies of sulfur‐containing molecules including cysteine and methionine residues in proteins. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Establishing effective simulation protocols for β‐ and α/β‐peptides. III. Molecular mechanical model for acyclic β‐amino acids

All‐atom molecular mechanics (MM) force field parameters are developed for the backbone of acyclic β‐amino acid using an improved version of the multiobjective evolutionary algorithm (MOEA). The MM model is benchmarked using β³‐homo‐Alanine (β³‐hAla) diamide in water with SCC‐DFTB/MM simulations as the reference. Satisfactory agreements are found between the MM and SCC‐DFTB/MM results regarding the distribution of key dihedral angles for the β³‐hAla diamide in water. The MM model is further applied to a β‐hepta‐peptide in methanol solution. The calculated NOE values and ³J coupling constants averaged over different trajectories are consistent with experimental data. By contrast, simulations using parameters directly transferred from the CHARMM22 force field for proteins lead to much worse agreement, which highlights the importance of careful parameterization for non‐natural peptides, for which the improved MOEA is particularly useful. Finally, as an initial application of the new force field parameters, the behaviors of a short random copolymer consisting of β amino acids in bulk solution and membrane/water interface are studied using a generalized Born implicit solvent model (GBSW). Results for four selected sequences show that segregation of hydrophobic and cationic groups occur easily at the membrane/solution interface for all sequences. The sequence that features alternating short blocks exhibits signs of lower stability at the interface compared to other sequences. These results confirm the hypothesis in recent experimental studies that β‐amino‐acid based random copolymers can develop a high degree of amphiphilicity without regular three‐dimensional structure. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Parameterization of the GPR119 Receptor Agonist AR231453

The GPR119 receptor is a class A G protein‐coupled receptor expressed mainly in pancreatic beta cells. Since GPR119 receptor activation ameliorates Type 2 Diabetes through an increase in glucose‐dependent insulin release, the development of new GPR119 receptor agonists would be worthwhile. A better understanding of the way agonists interact with the receptor would help to design better ligands for the receptor. It also would help to better understand the agonist mechanism of action. An understanding of how agonists interact with the receptor can be acquired using molecular dynamics simulations, which cannot be performed without having force field parameters for the ligand molecule. This study presents the development of CHARMM force field parameters for AR231453, the prototypical first potent and orally available GPR119 agonist, using the Force Field Tool Kit. The parameters are validated through Normal Mode Analysis calculations and molecular dynamics simulations in combination with infrared spectroscopy. © 2017 Wiley Periodicals, Inc.