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

Cumulative asparagine to aspartate deamidation fails to perturb γD-crystallin structure and stability

Deamidation frequently is invoked as an important driver of crystallin aggregation and cataract formation. Here, we characterized the structural and biophysical consequences of cumulative Asn to Asp changes in γD-crystallin. Using NMR spectroscopy, we demonstrate that N- or C-terminal domain-confined or fully Asn to Asp changed γD-crystallin exhibits essentially the same 1H-15N HSQC spectrum as the wild-type protein, implying that the overall structure is retained. Only a very small thermodynamic destabilization for the overall Asn to Asp γD-crystallin variants was noted by chaotropic unfolding, and assessment of the colloidal stability, by measuring diffusion interaction parameters, yielded no substantive differences in association propensities. Furthermore, using molecular dynamics simulations, no significant changes in dynamics for proteins with Asn to Asp or iso-Asp changes were detected. Our combined results demonstrate that substitution of all Asn by Asp residues, reflecting an extreme case of deamidation, did not affect the structure and biophysical properties of γD-crystallin. This suggests that these changes alone cannot be the major determinant in driving cataract formation.

Effects of altered backbone composition on the folding kinetics and mechanism of an ultrafast-folding protein

Sequence-encoded protein folding is a ubiquitous biological process that has been successfully engineered in a range of oligomeric molecules with artificial backbone chemical connectivity. A remarkable aspect of protein folding is the contrast between the rapid rates at which most sequences in nature fold and the vast number of conformational states possible in an unfolded chain with hundreds of rotatable bonds. Research efforts spanning several decades have sought to elucidate the fundamental chemical principles that dictate the speed and mechanism of natural protein folding. In contrast, little is known about how protein mimetic entities transition between an unfolded and folded state. Here, we report effects of altered backbone connectivity on the folding kinetics and mechanism of the B domain of Staphylococcal protein A (BdpA), an ultrafast-folding sequence. A combination of experimental biophysical analysis and atomistic molecular dynamics simulations performed on the prototype protein and several heterogeneous-backbone variants reveal the interplay among backbone flexibility, folding rates, and structural details of the transition state ensemble. Collectively, these findings suggest a significant degree of plasticity in the mechanisms that can give rise to ultrafast folding in the BdpA sequence and provide atomic level insights into how protein mimetic chains adopt an ordered folded state.

Accelerated Molecular Dynamics for Peptide Folding: Benchmarking Different Combinations of Force Fields and Explicit Solvent Models

Accelerated molecular dynamics (aMD) protocols were assessed on predicting the secondary structure of eight peptides, of which two are helical, three are β-hairpins, and three are disordered. Protocols consisted of combinations of three force fields (ff99SB, ff14SB, ff19SB) and two explicit solvation models (TIP3P and OPC), and were evaluated in two independent aMD simulations, one starting from an extended conformation, the other starting from a misfolded conformation. The results of these analyses indicate that all three combinations performed well on helical peptides. As for β-hairpins, ff19SB performed well with both solvation methods, with a slight preference for the TIP3P solvation model, even though performance was dependent on both peptide sequence and initial conformation. The ff19SB/OPC combination had the best performance on intrinsically disordered peptides. In general, ff14SB/TIP3P suffered the strongest helical bias.

Conformational landscapes of artificial peptides predicted by various force fields: are we ready to simulate β-amino acids?

With the introduction of artificial peptides as antimicrobial agents and organic catalysts, numerous efforts have been made to design foldamers with desirable structures and functions. Computational tools are a helpful proxy for revealing the dynamic structures at atomic resolution and understanding foldamer's complex structure-function relationships. However, the performance of conventional force fields in predicting the structures of artificial peptides has not been systematically evaluated. In this study, we critically assessed three popular force fields, AMBER ff14SB, CHARMM36m, and OPLS-AA/L, in predicting conformational propensities of a β-peptide foldamer at monomer and hexamer levels. Simulation results were compared to those obtained from quantum chemistry calculations and experimental data. We also utilised replica exchange molecular dynamics simulations to investigate the energy landscape of each force field and assess the similarities and differences between force fields. We compared different solvent systems in the AMBER ff14SB and CHARMM36m frameworks and confirmed the unanimous role of hydrogen bonds in shaping energy landscapes. We anticipate that our data will pave the way for further improvements to force fields and for understanding the role of solvents in peptide folding, crystallisation, and engineering.

Identification of a short ACE2-derived stapled peptide targeting the SARS-CoV-2 spike protein

The design and synthesis of a series of peptide derivatives based on a short ACE2 α-helix 1 epitope and subsequent [i - i+4] stapling of the secondary structure resulted in the identification of a 9-mer peptide capable to compete with recombinant ACE2 towards Spike RBD in the micromolar range. Specifically, SARS-CoV-2 Spike inhibitor screening based on colorimetric ELISA assay and structural studies by circular dichroism showed the ring-closing metathesis cyclization being capable to stabilize the helical structure of the 9-mer ³⁴HEAEDLFYQ⁴² epitope better than the triazole stapling via click chemistry. MD simulations showed the stapled peptide being able not only to bind the Spike RBD, sterically interfering with ACE2, but also showing higher affinity to the target as compared to parent epitope.

Chemical Shifts of Artificial Monomers Used to Construct Heterogeneous-Backbone Protein Mimetics in Random Coil and Folded States

The construction of protein-sized synthetic chains that blend natural amino acids with artificial monomers to create so-called heterogeneous-backbones is a powerful approach to generate complex folds and functions from bio-inspired agents. A variety of techniques from structural biology commonly used to study natural proteins have been adapted to investigate folding in these entities. In NMR characterization of proteins, proton chemical shift is a straightforward to acquire, information-rich metric that bears directly on a variety of properties related to folding. Leveraging chemical shift to gain insight into folding requires a set of reference chemical shift values corresponding to each building block type (i.e., the 20 canonical amino acids in the case of natural proteins) in a random coil state and knowledge of systematic changes in chemical shift associated with particular folded conformations. Although well documented for natural proteins, these issues remain unexplored in the context of protein mimetics. Here, we report random coil chemical shift values for a library of artificial amino acid monomers frequently used to construct heterogeneous-backbone protein analogues as well as a spectroscopic signature associated with one monomer class, β3-residues bearing proteinogenic side chains, adopting a helical folded conformation. Collectively, these results will facilitate the continued utilization of NMR for the study of structure and dynamics in protein-like artificial backbones.

Heterogeneous-Backbone Proteomimetic Analogues of Lasiocepsin, a Disulfide-Rich Antimicrobial Peptide with a Compact Tertiary Fold

The emergence of resistance to clinically used antibiotics by bacteria presents a significant problem in public health. Natural antimicrobial peptides (AMPs) are a valuable source of antibiotics that act by a mechanism less prone to the evolutionary development of resistance. In an effort to realize the promise of AMPs while overcoming limitations such as poor biostability, researchers have sought sequence-defined oligomers with artificial amide-based backbones that show membrane-disrupting functions similar to natural agents. Most of this precedent has focused on short peptidomimetic analogues of unstructured chains or secondary folds; however, the natural antimicrobial arsenal includes a number of small- and medium-sized proteins that act via an ordered tertiary structure. Generating proteomimetic analogues of these scaffolds poses a challenge due to the increased complexity of the target for mimicry. Here, we report the development of heterogeneous-backbone variants of lasiocepsin, a 27-residue disulfide-rich AMP found in bee venom that adopts a compact tertiary fold. Iterative cycles of design, synthesis, and biological evaluation yielded analogues of the natural domain with ∼30 to 40% artificial backbone content, comparable antibacterial activity, reduced host cell toxicity, and improved stability to proteolytic degradation. High-resolution structures determined for several variants by NMR provide insights into the interplay among backbone composition, tertiary fold, and biological properties. Collectively, the results reported here broaden the scope of protein functional mimicry by artificial backbone analogues of tertiary folding patterns and suggest protein backbone engineering as a means to tune protein function by exerting site-specific control over protein folded structure.

Development and Validation of Fluorinated, Aromatic Amino Acid Parameters for Use with the AMBER ff15ipq Protein Force Field

We developed force field parameters for fluorinated, aromatic amino acids enabling molecular dynamics (MD) simulations of fluorinated proteins. These parameters are tailored to the AMBER ff15ipq protein force field and enable the modeling of 4, 5, 6, and 7F-tryptophan, 3F- and 3,5F-tyrosine, and 4F- or 4-CF3-phenylalanine. The parameters include 181 unique atomic charges derived using the implicitly polarized charge (IPolQ) scheme in the presence of SPC/Eb explicit water molecules and 9 unique bond, angle, or torsion terms. Our simulations of benchmark peptides and proteins maintain expected conformational propensities on the μs time scale. In addition, we have developed an open-source Python program to calculate fluorine relaxation rates from MD simulations. The extracted relaxation rates from protein simulations are in good agreement with experimental values determined by 19F NMR. Collectively, our results illustrate the power and robustness of the IPolQ lineage of force fields for modeling the structure and dynamics of fluorine-containing proteins at the atomic level.

Do molecular dynamics force fields accurately model Ramachandran distributions of amino acid residues in water?

Molecular dynamics (MD) is a powerful tool for studying intrinsically disordered proteins, however, its reliability depends on the accuracy of the force field. We assess Amber ff19SB, Amber ff14SB, OPLS-AA/M, and CHARMM36m with respect to their capacity to capture intrinsic conformational dynamics of 14 guest residues x (=G, A, L, V, I, F, Y, D^P, E^P, R, C, N, S, T) in GxG peptides in water. The MD-derived Ramachandran distribution of each guest residue is used to calculate 5 J-coupling constants and amide I' band profiles to facilitate a comparison to spectroscopic data through reduced χ² functions. We show that the Gaussian model, optimized to best fit the experimental data, outperforms all MD force fields by an order of magnitude. The weaknesses of the MD force fields are: (i) insufficient variability of the polyproline II (pPII) population among the guest residues; (ii) oversampling of antiparallel at the expense of transitional β-strand region; (iii) inadequate sampling of turn-forming conformations for ionizable and polar residues; and (iv) insufficient guest residue-specificity of the Ramachandran distributions. Whereas Amber ff19SB performs worse than the other three force fields with respect to χ² values, it accounts for residue-specific pPII content better than the other three force fields. Additional testing of residue-specific RSFF1 and Amber ff14SB combined with TIP4P/2005 on six guest residues x (=A, I, F, D^P, R, S) reveals that residue specificity derived from protein coil libraries or an improved water model alone do not result in significantly lower χ² values.

Implications of the unfolded state in the folding energetics of heterogeneous-backbone protein mimetics

Sequence-encoded folding is the foundation of protein structure and is also possible in synthetic chains of artificial chemical composition. In natural proteins, the characteristics of the unfolded state are as important as those of the folded state in determining folding energetics. While much is known about folded structures adopted by artificial protein-like chains, corresponding information about the unfolded states of these molecules is lacking. Here, we report the consequences of altered backbone composition on the structure, stability, and dynamics of the folded and unfolded states of a compact helix-rich protein. Characterization through a combination of biophysical experiments and atomistic simulation reveals effects of backbone modification that depend on both the type of artificial monomers employed and where they are applied in sequence. In general, introducing artificial connectivity in a way that reinforces characteristics of the unfolded state ensemble of the prototype natural protein minimizes the impact of chemical changes on folded stability. These findings have implications in the design of protein mimetics and provide an atomically detailed picture of the unfolded state of a natural protein and artificial analogues under non-denaturing conditions.

Biophysical experiments and atomistic simulation reveal impacts of protein backbone alteration on the ensemble that defines the unfolded state. These effects have implications on folded stability of protein mimetics.

Analysis of folded structure and folding thermodynamics in heterogeneous-backbone proteomimetics

Recent years have seen a growing number of examples of designed oligomeric molecules with artificial backbone connectivity that are capable of adopting complex folded tertiary structures analogous to those seen in natural proteins. A range of experimental techniques from structural biology and biophysics have been brought to bear in the study of these proteomimetic agents. Here, we discuss some considerations encountered in the characterization of high-resolution folded structure as well as folding thermodynamics of protein-like artificial backbones. We provide an overview of the use of X-ray crystallography and NMR spectroscopy in such systems and review example applications of these methods in the primary literature. Further, we provide detailed protocols for two experiments that have proved useful in our prior and ongoing efforts to compare folding thermodynamics between natural protein domains and heterogeneous-backbone counterparts.