Development of the CHARMM Force Field for Lipids

The development of the CHARMM additive all-atom lipid force field (FF) is traced from the early 1990's to the most recent version (C36) published in 2010. Though simulations with early versions yielded useful results, they failed to reproduce two important quantities: a zero surface tension at the experimental bilayer surface area, and the signature splitting of the deuterium order parameters in the glycerol and upper chain carbons. Systematic optimization of parameters based on high level quantum mechanical data and free energy simulations have resolved these issues, and bilayers with a wide range of lipids can be simulated in tensionless ensembles using C36. Issues associated with other all-atom lipid FFs, success and limitations in the C36 FF and ongoing developments are also discussed.     

Giant Residual Dipolar 13C–1H Couplings in High‐Spin Organoiron Complexes: Elucidation of Their Structures in Solution by 13C NMR Spectroscopy

High‐spin Fe^II–alkyl complexes with bis(pyridylimino)isoindolato ligands were synthesized and their paramagnetic ¹H and ¹³C NMR spectra were analyzed comprehensively. The experimental ¹³C—¹H coupling values are temperature (T^?1)‐ as well as magnetic‐field (B²)‐dependent and deviate considerably from typical scalar ¹J_CH couplings constants. This deviation is attributed to residual dipolar couplings (RDCs), which arise from partial alignment of the complexes in the presence of a strong magnetic field. The analysis of the experimental RDCs allows an unambiguous assignment of all ¹³C NMR resonances and, additionally, a structural refinement of the conformation of the complexes in solution. Moreover the RDCs can be used for the analysis of the alignment tensor and hence the tensor of the anisotropy of the magnetic susceptibility.     

Comparison between self-guided Langevin dynamics and molecular dynamics simulations for structure refinement of protein loop conformations

This article presents a comparative analysis of two replica‐exchange simulation methods for the structure refinement of protein loop conformations, starting from low‐resolution predictions. The methods are self‐guided Langevin dynamics (SGLD) and molecular dynamics (MD) with a Nosé–Hoover thermostat. We investigated a small dataset of 8‐ and 12‐residue loops, with the shorter loops placed initially from a coarse‐grained lattice model and the longer loops from an enumeration assembly method (the Loopy program). The CHARMM22 + CMAP force field with a generalized Born implicit solvent model (molecular‐surface parameterized GBSW2) was used to explore conformational space. We also assessed two empirical scoring methods to detect nativelike conformations from decoys: the all‐atom distance‐scaled ideal‐gas reference state (DFIRE‐AA) statistical potential and the Rosetta energy function. Among the eight‐residue loop targets, SGLD out performed MD in all cases, with a median of 0.48 Å reduction in global root‐mean‐square deviation (RMSD) of the loop backbone coordinates from the native structure. Among the more challenging 12‐residue loop targets, SGLD improved the prediction accuracy over MD by a median of 1.31 Å, representing a substantial improvement. The overall median RMSD for SGLD simulations of 12‐residue loops was 0.91 Å, yielding refinement of a median 2.70 Å from initial loop placement. Results from DFIRE‐AA and the Rosetta model applied to rescoring conformations failed to improve the overall detection calculated from the CHARMM force field. We illustrate the advantage of SGLD over the MD simulation model by presenting potential‐energy landscapes for several loop predictions. Our results demonstrate that SGLD significantly outperforms traditional MD in the generation and populating of nativelike loop conformations and that the CHARMM force field performs comparably to other empirical force fields in identifying these conformations from the resulting ensembles. Published 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

RDC‐enhanced structure calculation of a β‐heptapeptide in methanol

Residual dipolar couplings (RDCs) are a rich source of structural information that goes beyond the range covered by the nuclear Overhauser effect or scalar coupling constants. They can only be measured in partially oriented samples. RDC studies of peptides in organic solvents have so far been focused on samples in chloroform or DMSO. Here, we show that stretched poly(vinyl acetate) can be used for the partial alignment of a linear β‐peptide with proteinogenic side chains in methanol. ¹D_CH, ¹D_NH, and ²D_HH RDCs were collected with this sample and included as restraints in a simulated annealing calculation. Incorporation of RDCs in the structure calculation process improves the long‐range definition in the backbone of the resulting 3₁₄‐helix and uncovers side‐chain mobility. Experimental side‐chain RDCs of the central leucine and valine residues are in good agreement with predicted values from a local three‐state model. Copyright © 2016 John Wiley & Sons, Ltd.     

Residual dipolar couplings in short peptides reveal systematic conformational preferences of individual amino acids

Residual dipolar couplings (RDCs) observed by NMR in solution under weak alignment conditions can monitor average net orientations and order parameters of individual bonds. By their simple geometrical dependence, RDCs bear particular promise for the quantitative characterization of conformations in partially folded or unfolded proteins. We have systematically investigated the influence of amino acid substitutions X on the conformation of unfolded model peptides EGAAXAASS as monitored by their (1)H(Nu)-(15)N and (1)H(alpha)-(13)C(alpha) RDCs detected at natural abundance of (15)N and (13)C in strained polyacrylamide gels. In total, 14 single amino acid substitutions were investigated. The RDCs show a specific dependence on the substitution X that correlates to steric or hydrophobic interactions with adjacent amino acids. In particular, the RDCs for the glycine and proline substitutions indicate less or more order, respectively, than the other amino acids. The RDCs for aromatic substitutions tryptophane and tyrosine give evidence of a kink in the peptide backbone. This effect is also observable for orientation by Pf1 phages and corroborated by variations in (13)C(alpha) secondary shifts and (3)J(HNH)(alpha) scalar couplings in isotropic samples. RDCs for a substitution with the beta-turn sequence KNGE differ from single amino acid substitutions. Terminal effects and next neighbor effects could be demonstrated by further specific substitutions. The results were compared to statistical models of unfolded peptide conformations derived from PDB coil subsets, which reproduce overall trends for (1)H(Nu)-(15)N RDCs for most substitutions, but deviate more strongly for (1)H(alpha)-(13)C(alpha) RDCs. The outlined approach opens the possibility to obtain a systematic experimental characterization of the influence of individual amino acid/amino acid interactions on orientational preferences in polypeptides.     

All‐atom polarizable force field for DNA based on the classical drude oscillator model

Presented is a first generation atomistic force field (FF) for DNA in which electronic polarization is modeled based on the classical Drude oscillator formalism. The DNA model is based on parameters for small molecules representative of nucleic acids, including alkanes, ethers, dimethylphosphate, and the nucleic acid bases and empirical adjustment of key dihedral parameters associated with the phosphodiester backbone, glycosidic linkages, and sugar moiety of DNA. Our optimization strategy is based on achieving a compromise between satisfying the properties of the underlying model compounds in the gas phase targeting quantum mechanical (QM) data and reproducing a number of experimental properties of DNA duplexes in the condensed phase. The resulting Drude FF yields stable DNA duplexes on the 100‐ns time scale and satisfactorily reproduce (1) the equilibrium between A and B forms of DNA and (2) transitions between the BI and BII substates of B form DNA. Consistency with the gas phase QM data for the model compounds is significantly better for the Drude model as compared to the CHARMM36 additive FF, which is suggested to be due to the improved response of the model to changes in the environment associated with the explicit inclusion of polarizability. Analysis of dipole moments associated with the nucleic acid bases shows the Drude model to have significantly larger values than those present in CHARMM36, with the dipoles of individual bases undergoing significant variations during the MD simulations. Additionally, the dipole moment of water was observed to be perturbed in the grooves of DNA. © 2014 Wiley Periodicals, Inc.     

J-GFT NMR for precise measurement of mutually correlated nuclear spin-spin couplings

G-matrix Fourier transform (GFT) NMR spectroscopy is presented for accurate and precise measurement of chemical shifts and nuclear spin-spin couplings correlated according to spin system. The new approach, named "J-GFT NMR", is based on a largely extended GFT NMR formalism and promises to have a broad impact on projection NMR spectroscopy. Specifically, constant-time J-GFT (6,2)D (HA-CA-CO)-N-HN was implemented for simultaneous measurement of five mutually correlated NMR parameters, that is, 15N backbone chemical shifts and the four one-bond spin-spin couplings 13Calpha-1Halpha, 13Calpha-13C', 15N-13C', and 15N-1HNu. The experiment was applied for measuring residual dipolar couplings (RDCs) in an 8 kDa protein Z-domain aligned with Pf1 phages. Comparison with RDC values extracted from conventional NMR experiments reveals that RDCs are measured with high precision and accuracy, which is attributable to the facts that (i) the use of constant time evolution ensures that signals do not broaden whenever multiple RDCs are jointly measured in a single dimension and (ii) RDCs are multiply encoded in the multiplets arising from the joint sampling. This corresponds to measuring the couplings multiple times in a statistically independent manner. A key feature of J-GFT NMR, i.e., the correlation of couplings according to spin systems without reference to sequential resonance assignments, promises to be particularly valuable for rapid identification of backbone conformation and classification of protein fold families on the basis of statistical analysis of dipolar couplings.     

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.     

Probing Long‐Range Anisotropic Interactions: a General and Sign‐Sensitive Strategy to Measure 1H–1H Residual Dipolar Couplings as a Key Advance for Organic Structure Determination

Residual dipolar couplings (RDCs) are amongst the most powerful NMR parameters for organic structure elucidation. In order to maximize their effectiveness in increasingly complex cases such as flexible compounds, a maximum of RDCs between nuclei sampling a large distribution of orientations is needed, including sign information. For this, the easily accessible one‐bond ¹H–¹³C RDCs alone often fall short. Long‐range ¹H–¹H RDCs are both abundant and typically sample highly complementary orientations, but accessing them in a sign‐sensitive way has been severely obstructed due to the overflow of ¹H–¹H couplings. Here, we present a generally applicable strategy that allows the measurement of a large number of ¹H–¹H RDCs, including their signs, which is based on a combination of an improved PSYCHEDELIC method and a new selective constant‐time β‐COSY experiment. The potential of ¹H–¹H RDCs to better determine molecular alignment and to discriminate between enantiomers and diastereomers is demonstrated.     

Increasing the accuracy of solution NMR structures of membrane proteins by application of residual dipolar couplings. High-resolution structure of outer membrane protein A

The structure determination of membrane proteins is one of the most challenging applications of solution NMR spectroscopy. The paucity of distance information available from the highly deuterated proteins employed requires new approaches in structure determination. Here we demonstrate that significant improvement in the structure accuracy of the membrane protein OmpA can be achieved by refinement with residual dipolar couplings (RDCs). The application of charged polyacrylamide gels allowed us to obtain two alignments and accurately measure numerous heteronuclear dipolar couplings. Furthermore, we have demonstrated that using a large set of RDCs in the refinement can yield a structure with 1 A rms deviation to the backbone of the high-resolution crystal structure. Our simulations with various data sets indicate that dipolar couplings will be critical for obtaining accurate structures of membrane proteins.     

Long‐Range Residual Dipolar Couplings: A Tool for Determining the Configuration of Small Molecules

Together with NOE and J coupling, one‐bond residual dipolar coupling (RDC), which reports on the three‐dimensional orientation of an internuclear vector in the molecular frame, plays an important role in the conformation and configuration analysis of small molecules in solution by NMR spectroscopy. When the molecule has few C? H bonds, or too many bonds are in parallel, the available RDCs may not be sufficient to obtain the alignment tensor used for structure elucidation. Long‐range RDCs that connect nuclei over multiple bonds are normally not parallel to the single bonds and therefore complement one‐bond RDCs. Herein we present a method for extracting the long‐range RDC of a chosen proton or group of protons to all remotely connected carbon atoms, including non‐protonated carbon atoms. Alignment tensors fitted directly to the total long‐range couplings (T=J+D) enabled straightforward analysis of both the long‐range and one‐bond RDCs for strychnine.     

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.     

Observation of CH⋅⋅⋅π Interactions between Methyl and Carbonyl Groups in Proteins

Protein structure and function is dependent on myriad noncovalent interactions. Direct detection and characterization of these weak interactions in large biomolecules, such as proteins, is experimentally challenging. Herein, we report the first observation and measurement of long-range “through-space” scalar couplings between methyl and backbone carbonyl groups in proteins. These J couplings are indicative of the presence of noncovalent C−H⋅⋅⋅π hydrogen-bond-like interactions involving the amide π network. Experimentally detected scalar couplings were corroborated by a natural bond orbital analysis, which revealed the orbital nature of the interaction and the origins of the through-space J couplings. The experimental observation of this type of CH⋅⋅⋅π interaction adds a new dimension to the study of protein structure, function, and dynamics by NMR spectroscopy.     

Force field‐dependent solution properties of glycine oligomers

Molecular simulations can be used to study disordered polypeptide systems and to generate hypotheses on the underlying structural and thermodynamic mechanisms that govern their function. As the number of disordered protein systems investigated with simulations increase, it is important to understand how particular force fields affect the structural properties of disordered polypeptides in solution. To this end, we performed a comparative structural analysis of Gly₃ and Gly₁₀ in aqueous solution from all atom, microsecond molecular dynamics (MD) simulations using the CHARMM 27 (C27), CHARMM 36 (C36), and Amber ff12SB force fields. For each force field, Gly₃ and Gly₁₀ were simulated for at least 300 ns and 1 μs, respectively. Simulating oligoglycines of two different lengths allows us to evaluate how force field effects depend on polypeptide length. Using a variety of structural metrics (e.g., end‐to‐end distance, radius of gyration, dihedral angle distributions), we characterize the distribution of oligoglycine conformers for each force field and show that each sample conformation space differently, yielding considerably different structural tendencies of the same oligoglycine model in solution. Notably, we find that C36 samples more extended oligoglycine structures than both C27 and ff12SB. © 2015 Wiley Periodicals, Inc.     

Combination of the CHARMM27 force field with united-atom lipid force fields

Computer simulations offer a valuable way to study membrane systems, from simple lipid bilayers to large transmembrane protein complexes and lipid-nucleic acid complexes for drug delivery. Their accuracy depends on the quality of the force field parameters used to describe the components of a particular system. We have implemented the widely used CHARMM22 and CHARMM27 force fields in the GROMACS simulation package to (i) combine the CHARMM22 protein force field with two sets of united-atom lipids parameters; (ii) allow comparisons of the lipid CHARMM27 force field with other lipid force fields or lipid-protein force field combinations. Our tests do not show any particular issue with the combination of the all-atom CHARMM22 force field with united-atoms lipid parameters, although pertinent experimental data are lacking to assess the quality of the lipid-protein interactions. The conversion utilities allow automatic generation of GROMACS simulation files with CHARMM nucleic acids and protein parameters and topologies, starting from pdb files using the standard GROMACS pdb2gmx method. CMAP is currently not implemented.     

Internucleotide J-couplings and chemical shifts of the N-H···N hydrogen-bonds in the radiation-damaged guanine-cytosine base pairs

Internucleotide ^2hJ_NN spin‐spin couplings and chemical shifts (δ(¹H) and Δδ(¹⁵N)) of N? H···N H‐bond units in the natural and radiation‐damaged G‐C base pairs were predicted using the appropriate density functional theory calculations with a large basis set. Four possible series of the damaged G‐C pairs (viz., dehydrogenated and deprotonated G‐C pairs, GC^?? and GC^?+ radicals) were discussed carefully in this work. Computational NMR results show that radicalization and anionization of the base pairs can yield strong effect on their ^2hJ_NN spin scalar coupling constants and the corresponding chemical shifts. Thus, variations of the NMR parameters associated with the N? H···N H‐bonds may be taken as an important criterion for prejudging whether the natural G‐C pair is radiation‐damaged or not. Analysis shows that ^2hJ_NN couplings are strongly interrelated with the energy gaps (ΔE_LP→σ*) and the second‐order interaction energies (E(2)) between the donor N lone‐pair (LP_N) and the acceptor σ^*_{N? H} localized NBO orbitals, and also are sensitive to the electron density distributions over the σ*(N? H) orbital, indicating that ^2hJ_NN couplings across the N? H···N H‐bonds are charge‐transfer‐controlled. This is well supported by variation of the electrostatic potential surfaces and corresponding charge transfer amount between G and C moieties. It should be noted that although the NMR spectra for the damaged G‐C pair radicals are unavailable now and the states of the radicals are usually detected by the electron spin resonance, this study provides a correlation of the properties of the damaged DNA species with some of the electronic parameters associated with the NMR spectra for the understanding of the different state character of the damaged DNA bases. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011.     

Side chain orientation from methyl 1H-1H residual dipolar couplings measured in highly deuterated proteins

High-level deuteration is a prerequisite for the study of high molecular weight systems using liquid-state NMR. Here, we present new experiments for the measurement of proton-proton dipolar couplings in CH(2)D methyl groups of (13)C labeled, highly deuterated (70-80%) proteins. (1)H-(1)H residual dipolar couplings (RDCs) have been measured in two alignment media for 57 out of 70 possible methyl containing residues in the 167-residue flavodoxin-like domain of the E. coli sulfite reductase. These data yield information on the orientation of the methyl symmetry axis with respect to the molecular alignment frame. The alignment tensor characteristics were obtained very accurately from a set of backbone RDCs measured on the same protein sample. To demonstrate that accurate structural information is obtained from these data, the measured methyl RDCs for Valine residues are analyzed in terms of chi(1) torsion angles and stereospecific assignment of the prochiral methyl groups. On the basis of the previously determined backbone solution structure of this protein, the methyl RDC data proved sufficient to determine the chi(1) torsion angles in seven out of nine valines, assuming a single-rotamer model. Methyl RDCs are complementary to other NMR data, for example, methyl-methyl NOE, to determine side chain conformation in high molecular weight systems.     

Discriminating the helical forms of peptides by NMR and molecular dynamics simulation

The HNCO NMR pulse sequence was applied to three selectively labeled (15)N and (13)C isotopic homologues of the peptide Ac-WAAAH(AAARA)(3)A-NH(2) to probe directly for hydrogen bonds between residues 8 and 11 (characteristic of a 3(10)-helix), 8 and 12 (alpha-helix), and 8 and 13 (pi-helix). The experiments demonstrate conclusively, and in agreement with circular dichroism studies, that the center of the peptide is alpha-helical; there is no discernible 3(10)- or pi-helix at these specific positions. Molecular dynamics simulations of the preceding peptide and Ac-(AAAAK)(3)A-NH(2) in water using the potential energy parameter set CHARMM22/CMAP correctly yield an alpha-helix, in contrast to simulations with the set CHARMM22, which result in a pi-helix.     

Correlation of protein structure and dynamics to scalar couplings across hydrogen bonds

NMR-observable scalar couplings across hydrogen bonds in nucleic acids and proteins present a quantitative measure for the geometry and--by the implicit experimental time averaging--dynamics of hydrogen bonds. We have carried out in-depth molecular dynamics (MD) simulations with various force fields on three proteins: ubiquitin, the GB1 domain of protein G, and the SMN Tudor domain, for which experimental h3JNC' scalar couplings of backbone hydrogen bonds and various high-resolution X-ray structures are available. Theoretical average values for h3JNC' were calculated from the snapshots of these MD simulations either by density functional theory or by a geometric parametrization (Barfield, M. J. Am. Chem. Soc. 2002, 124, 4158-4168). No significant difference was found between the two methods. The results indicate that time-averaging using explicit water solvation in the MD simulations improves significantly the agreement between experimental and theoretical values for the lower resolution structures ubiquitin (1.8 A), Tudor domain (1.8 A), and protein G (2.1 A). Only marginal improvement is found for the high-resolution structure (1.1 A) of protein G. Hence, experimental h3JNC' values are compatible with a static, high-resolution structural model. The MD averaging of the low-resolution structures moves the averages of the rHO distance and the H...O=C angle theta closer to their respective values in the high-resolution structures, thereby improving the agreement using experimental h3JNC' data. In contrast, MD averaging with implicit water models deteriorates the agreement with experiment for all proteins. The differing behavior can be explained by an artifactual lengthening of H-bonds caused by the implicit water models.