A twist in the road less traveled: The AMBER ff15ipq-m force field for protein mimetics.

We present a new force field, AMBER ff15ipq-m, for simulations of protein mimetics in applications from therapeutics to biomaterials. This force field is an expansion of the AMBER ff15ipq force field that was developed for canonical proteins and enables the modeling of four classes of artificial backbone units that are commonly used alongside natural α residues in blended or "heterogeneous" backbones: chirality-reversed D-α-residues, the Cα-methylated α-residue Aib, homologated ß-residues (ß3) bearing proteinogenic side chains, and two cyclic ß residues (ßcyc; APC and ACPC). The ff15ipq-m force field includes 472 unique atomic charges and 148 unique torsion terms. Consistent with the AMBER IPolQ lineage of force fields, the charges were derived using the Implicitly Polarized Charge (IPolQ) scheme in the presence of explicit solvent. To our knowledge, no general force field reported to date models the combination of artificial building blocks examined here. In addition, we have derived Karplus coefficients for the calculation of backbone amide J-coupling constants for ß3Ala and ACPC ß residues. The AMBER ff15ipq-m force field reproduces experimentally observed J-coupling constants in simple tetrapeptides and maintains the expected conformational propensities in reported structures of proteins/peptides containing the artificial building blocks of interest-all on the µs timescale. These encouraging results demonstrate the power and robustness of the IPolQ lineage of force fields in modeling the structure and dynamics of natural proteins as well as mimetics with protein-inspired artificial backbones in atomic detail.

Integrating NMR, SAXS, and Atomistic Simulations: Structure and Dynamics of a Two-Domain Protein.

Multidomain proteins with two or more independently folded functional domains are prevalent in nature. Whereas most multidomain proteins are linked linearly in sequence, roughly one-tenth possess domain insertions where a guest domain is implanted into a loop of a host domain, such that the two domains are connected by a pair of interdomain linkers. Here, we characterized the influence of the interdomain linkers on the structure and dynamics of a domain-insertion protein in which the guest LysM domain is inserted into a central loop of the host CVNH domain. Expanding upon our previous crystallographic and NMR studies, we applied SAXS in combination with NMR paramagnetic relaxation enhancement to construct a structural model of the overall two-domain system. Although the two domains have no fixed relative orientation, certain orientations were found to be preferred over others. We also assessed the accuracies of molecular mechanics force fields in modeling the structure and dynamics of tethered multidomain proteins by integrating our experimental results with microsecond-scale atomistic molecular dynamics simulations. In particular, our evaluation of two different combinations of the latest force fields and water models revealed that both combinations accurately reproduce certain structural and dynamical properties, but are inaccurate for others. Overall, our study illustrates the value of integrating experimental NMR and SAXS studies with long timescale atomistic simulations for characterizing structural ensembles of flexibly linked multidomain systems.

Links between the charge model and bonded parameter force constants in biomolecular force fields.

The ff15ipq protein force field is a fixed charge model built by automated tools based on the two charge sets of the implicitly polarized charge method: one set (appropriate for vacuum) for deriving bonded parameters and the other (appropriate for aqueous solution) for running simulations. The duality is intended to treat water-induced electronic polarization with an understanding that fitting data for bonded parameters will come from quantum mechanical calculations in the gas phase. In this study, we compare ff15ipq to two alternatives produced with the same fitting software and a further expanded data set but following more conventional methods for tailoring bonded parameters (harmonic angle terms and torsion potentials) to the charge model. First, ff15ipq-Qsolv derives bonded parameters in the context of the ff15ipq solution phase charge set. Second, ff15ipq-Vac takes ff15ipq's bonded parameters and runs simulations with the vacuum phase charge set used to derive those parameters. The IPolQ charge model and associated protocol for deriving bonded parameters are shown to be an incremental improvement over protocols that do not account for the material phases of each source of their fitting data. Both force fields incorporating the polarized charge set depict stable globular proteins and have varying degrees of success modeling the metastability of short (5-19 residues) peptides. In this particular case, ff15ipq-Qsolv increases stability in a number of α-helices, correctly obtaining 70% helical character in the K19 system at 275 K and showing appropriately diminishing content up to 325 K, but overestimating the helical fraction of AAQAA3 by 50% or more, forming long-lived α-helices in simulations of a ß-hairpin, and increasing the likelihood that the disordered p53 N-terminal peptide will also form a helix. This may indicate a systematic bias imparted by the ff15ipq-Qsolv parameter development strategy, which has the hallmarks of strategies used to develop other popular force fields, and may explain some of the need for manual corrections in this force fields' evolution. In contrast, ff15ipq-Vac incorrectly depicts globular protein unfolding in numerous systems tested, including Trp cage, villin, lysozyme, and GB3, and does not perform any better than ff15ipq or ff15ipq-Qsolv in tests on short peptides. We analyze the free energy surfaces of individual amino acid dipeptides and the electrostatic potential energy surfaces of each charge model to explain the differences.

Further along the Road Less Traveled: AMBER ff15ipq, an Original Protein Force Field Built on a Self-Consistent Physical Model.

We present the AMBER ff15ipq force field for proteins, the second-generation force field developed using the Implicitly Polarized Q (IPolQ) scheme for deriving implicitly polarized atomic charges in the presence of explicit solvent. The ff15ipq force field is a complete rederivation including more than 300 unique atomic charges, 900 unique torsion terms, 60 new angle parameters, and new atomic radii for polar hydrogens. The atomic charges were derived in the context of the SPC/Eb water model, which yields more-accurate rotational diffusion of proteins and enables direct calculation of nuclear magnetic resonance (NMR) relaxation parameters from molecular dynamics simulations. The atomic radii improve the accuracy of modeling salt bridge interactions relative to contemporary fixed-charge force fields, rectifying a limitation of ff14ipq that resulted from its use of pair-specific Lennard-Jones radii. In addition, ff15ipq reproduces penta-alanine J-coupling constants exceptionally well, gives reasonable agreement with NMR relaxation rates, and maintains the expected conformational propensities of structured proteins/peptides, as well as disordered peptides-all on the microsecond (µs) time scale, which is a critical regime for drug design applications. These encouraging results demonstrate the power and robustness of our automated methods for deriving new force fields. All parameters described here and the mdgx program used to fit them are included in the AmberTools16 distribution.

Structural Insight into Fungal Cell Wall Recognition by a CVNH Protein with a Single LysM Domain.

MGG_03307 is a lectin isolated from Magnaporte oryzae, a fungus that causes devastating rice blast disease. Its function is associated with protecting M. oryzae from the host immune response in plants. To provide the structural basis of how MGG_03307 protects the fungus, crystal structures of its CVNH-LysM module were determined in the absence and presence of GlcNAc-containing cell wall chitin constituents, which can act as pathogen-associated molecular patterns. Our structures revealed that glycan binding is accompanied by a notable conformational change in the LysM domain and that GlcNAc3 and GlcNAc4 are accommodated similarly. GlcNAc5 and GlcNAc6 interact with the LysM domain in multiple conformations, as evidenced by solution nuclear magnetic resonance studies. No dimerization of MoCVNH3 via its LysM domain was observed upon binding to GlcNAc6, unlike in multiple LysM domain-containing proteins. Importantly, we define a specific consensus binding mode for the recognition of GlcNAc oligomers by single LysM domains.

In vivo structure-activity relationship studies support allosteric targeting of a dual specificity phosphatase.

Dual specificity phosphatase 6 (DUSP6) functions as a feedback attenuator of fibroblast growth factor signaling during development. In vitro high throughput chemical screening attempts to discover DUSP6 inhibitors have yielded limited success. However, in vivo whole-organism screens of zebrafish identified compound 1 (BCI) as an allosteric inhibitor of DUSP6. Here we designed and synthesized a panel of analogues to define the structure-activity relationship (SAR) of DUSP6 inhibition. In vivo high-content analysis in transgenic zebrafish, coupled with cell-based chemical complementation assays, identified structural features of the pharmacophore of 1 that were essential for biological activity. In vitro assays of DUSP hyperactivation corroborated the results from in vivo and cellular SAR. The results reinforce the notion that DUSPs are druggable through allosteric mechanisms and illustrate the utility of zebrafish as a model organism for in vivo SAR analyses.

Evaluating the strength of salt bridges: a comparison of current biomolecular force fields.

Recent advances in computer hardware and software have made rigorous evaluation of current biomolecular force fields using microsecond-scale simulations possible. Force fields differ in their treatment of electrostatic interactions, including the formation of salt bridges in proteins. Here we conducted an extensive evaluation of salt bridge interactions in the latest AMBER, CHARMM, and OPLS force fields, using microsecond-scale molecular dynamics simulations of amino acid analogues in explicit solvent. We focused on salt bridges between three different pairs of oppositely charged amino acids: Arg/Asp, Lys/Asp, and His(+)/Asp. Our results reveal considerable variability in the predicted KA values of the salt bridges for these force fields, as well as differences from experimental data: almost all of the force fields overestimate the strengths of the salt bridges. When amino acids are represented by side-chain analogues, the AMBER ff03 force field overestimates the KA values the least, while for complete amino acids, the AMBER ff13α force field yields the lowest KA value, most likely caused by an altered balance of side-chain/side-chain and side-chain/backbone contacts. These findings confirm the notion that the implicit incorporation of solvent polarization improves the accuracy of modeling salt bridge interactions.