Helix stability of oligoglycine, oligoalanine, and oligo-β-alanine dodecamers reflected by hydrogen-bond persistence

Accurate Host-Guest Binding Free Energies Using the AMOEBA Polarizable Force Field

A grand challenge of computational biophysics is accurate prediction of interactions between molecules. Molecular dynamics (MD) simulations have recently gained much interest as a tool to directly compute rigorous intermolecular binding affinities. The choice of a fixed point-charge or polarizable multipole force field used in MD is a topic of ongoing discussion. To compare alternative methods, we participated in the SAMPL7 and SAMPL8 Gibb octaacid host-guest challenges to assess the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) polarizable multipole force field. Advantages of AMOEBA over fixed charge models include improved representation of molecular electrostatic potentials and better description of water occupying the unligated host cavity. Prospective predictions for 26 host-guest systems exhibit a mean unsigned error vs experiment of 0.848 kcal/mol across all absolute binding free energies, demonstrating excellent agreement between computational and experimental results. In addition, we explore two topics related to the inclusion of ions in MD simulations: use of a neutral co-alchemical protocol and the effect of salt concentration on binding affinity. Use of the co-alchemical method minimally affects computed energies, but salt concentration significantly perturbs our binding results. Higher salt concentration strengthens binding through classical charge screening. In particular, added Na+ ions screen negatively charged carboxylate groups near the binding cavity, thereby diminishing repulsive coulomb interactions with negatively charged guests. Overall, the AMOEBA results demonstrate the accuracy available through a force field providing a detailed energetic description of the four octaacid hosts and 13 charged organic guests. Use of the AMOEBA polarizable atomic multipole force field in conjunction with an alchemical free energy protocol can achieve chemical accuracy in application to realistic molecular systems.

Polypeptide Chain Growth Mechanisms and Secondary Structure Formation in Glycine Gas-Phase Deposition on Silica Surfaces

Peptide formation by amino acids condensation represents a crucial reaction in the quest of the origins of life as well as in synthetic chemistry. However, it is still poorly understood in terms of efficiency and reaction mechanism. In the present work, peptide formation has been investigated through thermal condensation of gas-phase glycine in fluctuating silica environments as a model of prebiotic environments. In-situ IR spectroscopy measurements under a controlled atmosphere reveal that a humidity fluctuating system subjected to both temperature and water activity variations results in the formation of more abundant peptides compared to a dehydrated system subjected only to temperature fluctuations cycles. A model is proposed in which hydration steps result in the hydrolysis and redistribution of the oligomers formed during previous deposition in dry conditions. This results in the formation of self-assembled aggregates with well-defined secondary structures (especially β-sheets). Upon further monomers feeding, structural elements are conserved in newly growing chains, with indications of templated polymerization. The structural dynamics of peptides were also evaluated. Rigid self-assembled structures with a high resistance to further wetting/drying cycles and inaccessibility to isotopic exchange were present in the humidity fluctuating system compared to more flexible structures in the dehydrated system. The resistance and growth of self-assembled structures were also investigated for an extended duration of Gly deposition using isotope labeling.

Searching for Aquamelt Behavior among Silklike Biomimetics during Fibrillation under Flow

In this paper, we elucidate a generic mechanism behind strain-induced phase transition in aqueous solutions of silk-inspired biomimetics by atomistic molecular dynamics simulations. We show the results of modeling of homopeptides polyglycine Gly₃₀ and polyalanine Ala₃₀ and a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅, i.e., the simplest and yet relevant sequences that could mimic the behavior of natural silk under stress conditions. First, we analyze hydrophobicities of the sequences by calculating the Gibbs free energy of hydration and inspecting the interchain hydrogen bonding and hydration by water. Second, the force-extension profiles are scanned and compared with the results for poly(ethylene oxide), the synthetic polymer for which the aquamelt behavior has been proved recently. Additionally, the conformational transitions of oligopeptides from coiled to extended states are characterized by a generalized order parameter and by the dependence of the solvent-accessible surface area of the chains on applied stretching. Fibrillation itself is surveyed using both the two-dimensional interchain pair correlation function and the SAXS/WAXS patterns for the aggregates formed under stress. These are compared with experimental data found in the literature on fibril structure of silk composite materials doped with oligoalanine peptides. Our results show that tensile stress introduced into aqueous oligopeptide solutions facilitates interchain interactions. The oligopeptides display both a greater resistance to extension as compared to poly(ethylene oxide) and a reduced ability for hydrogen bonding of the stretched chains between oligomers and with water. Fiber formation is proved for all simulated objects, but the most structured one is made of a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅: For this sequence, we obtain the highest degree of the secondary structure motifs in the fiber. We conclude that this is the most promising candidate among considered sequences to find the aquamelt behavior in further experimental studies.

Large Increase in the Dielectric Constant and Partial Loss of Coherence Increases Tunneling Rates across Molecular Wires

Although the dielectric behavior of monolayers is important in a large range of applications, its role in charge transport studies involving molecular junctions is largely ignored. This paper describes a large increase in the relative static dielectric constant (ε_r) by simply increasing the thickness of well-organized monolayers of oligoglycine and oligo(ethylene glycol) from 7 up to 14. The resulting large capacitance of 3.5-5.1 μF/cm² is thickness-independent, which is highly attractive for field-effect transistor applications. This increase of ε_r results in a linear increase of the thermal activation energy by a factor of 6, which suggests that the mechanism of charge transport gradually changes from coherent to (partially) incoherent tunneling. The comparisons of oligoglycine (which readily forms hydrogen bonds with neighboring molecules) and methyl terminated oligo(ethylene glycol) (which lacks hydrogen bond donors) monolayers, kinetic isotope effects, and relative humidity-dependent measurements all indicate the importance of strong hydrogen bonds involving ionic species and strongly bonded water in the unusual dielectric behavior and the incoherent tunneling mechanism. This partial loss of coherence of the charge carriers can explain the unusually small tunneling decay coefficients across long molecular wires, and the length-dependent increase of ε_r of monolayers opens up interesting new applications.

Absolute binding free energies for the SAMPL6 cucurbit[8]uril host-guest challenge via the AMOEBA polarizable force field

As part of the SAMPL6 host-guest blind challenge, the AMOEBA force field was applied to calculate the absolute binding free energy for a cucurbit[8]uril host complexed with 14 diverse guests, ranging from small, rigid structures to drug molecules. The AMOEBA results from the initial submission prompted an investigation into aspects of the methodology and parameterization employed. Lessons learned from the blind challenge include: a double annihilation scheme (electrostatics and van der Waals) is needed to obtain proper sampling of guest conformations, annihilation of key torsion parameters of the guest are recommended for flexible guests, and a more thorough analysis of torsion parameters is warranted. When put in to practice with the AMOEBA model, the lessons learned improved the MUE from 2.63 to 1.20 kcal/mol and the RMSE from 3.62 to 1.68 kcal/mol, respectively. Overall, the AMOEBA protocol for determining absolute binding free energies benefitted from participation in the SAMPL6 host-guest blind challenge and the results suggest the implementation of the methodology in future host-guest calculations.

Conformation Analysis of GalNAc-Appended Sugar Amino Acid Foldamers as Glycopeptide Mimics

Sugar amino acid (SAA)-based foldamers with well-defined secondary structures were appended with N-acetylgalactosamine (GalNAc) sugars to access sequence-defined, multidentate glycoconjugates with full control over number, spacing and position. Conformation analysis of these glycopeptides by extensive NMR spectroscopic studies revealed that the appended GalNAc units had a profound influence on the native conformational behaviour of the SAA foldamers. Whereas the 2,5-cis glycoconjugate showed a helical structure in water, comprising of two consecutive 16-membered hydrogen bonds, its 2,5-trans congener displayed an unprecedented 16/10-mixed turn structure not seen before in any glycopeptide foldamer.

N-terminal diproline and charge group effects on the stabilization of helical conformation in alanine-based short peptides: CD studies with water and methanol as solvent

Protein folding problem remains a formidable challenge as main chain, side chain and solvent interactions remain entangled and have been difficult to resolve. Alanine-based short peptides are promising models to dissect protein folding initiation and propagation structurally as well as energetically. The effect of N-terminal diproline and charged side chains is assessed on the stabilization of helical conformation in alanine-based short peptides using circular dichroism (CD) with water and methanol as solvent. A1 (Ac-Pro-Pro-Ala-Lys-Ala-Lys-Ala-Lys-Ala-NH₂ ) is designed to assess the effect of N-terminal homochiral diproline and lysine side chains to induce helical conformation. A2 (Ac-Pro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) and A3 (Ac-dPro-Pro-Glu-Glu-Ala-Ala-Lys-Lys-Ala-NH₂ ) with N-terminal homochiral and heterochiral diproline, respectively, are designed to assess the effect of Glu...Lys (i, i + 4) salt bridge interactions on the stabilization of helical conformation. The CD spectra of A1, A2 and A3 in water manifest different amplitudes of the observed polyproline II (PPII) signals, which indicate different conformational distributions of the polypeptide structure. The strong effect of solvent substitution from water to methanol is observed for the peptides, and CD spectra in methanol evidence A2 and A3 as helical folds. Temperature-dependent CD spectra of A1 and A2 in water depict an isodichroic point reflecting coexistence of two conformations, PPII and β-strand conformation, which is consistent with the previous studies. The results illuminate the effect of N-terminal diproline and charged side chains in dictating the preferences for extended-β, semi-extended PPII and helical conformation in alanine-based short peptides. The results of the present study will enhance our understanding on stabilization of helical conformation in short peptides and hence aid in the design of novel peptides with helical structures. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.

High-resolution crystal structures of protein helices reconciled with three-centered hydrogen bonds and multipole electrostatics

Theoretical and experimental evidence for non-linear hydrogen bonds in protein helices is ubiquitous. In particular, amide three-centered hydrogen bonds are common features of helices in high-resolution crystal structures of proteins. These high-resolution structures (1.0 to 1.5 Å nominal crystallographic resolution) position backbone atoms without significant bias from modeling constraints and identify Φ = -62°, ψ = -43 as the consensus backbone torsional angles of protein helices. These torsional angles preserve the atomic positions of α-β carbons of the classic Pauling α-helix while allowing the amide carbonyls to form bifurcated hydrogen bonds as first suggested by Némethy et al. in 1967. Molecular dynamics simulations of a capped 12-residue oligoalanine in water with AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications), a second-generation force field that includes multipole electrostatics and polarizability, reproduces the experimentally observed high-resolution helical conformation and correctly reorients the amide-bond carbonyls into bifurcated hydrogen bonds. This simple modification of backbone torsional angles reconciles experimental and theoretical views to provide a unified view of amide three-centered hydrogen bonds as crucial components of protein helices. The reason why they have been overlooked by structural biologists depends on the small crankshaft-like changes in orientation of the amide bond that allows maintenance of the overall helical parameters (helix pitch (p) and residues per turn (n)). The Pauling 3.6₁₃ α-helix fits the high-resolution experimental data with the minor exception of the amide-carbonyl electron density, but the previously associated backbone torsional angles (Φ, Ψ) needed slight modification to be reconciled with three-atom centered H-bonds and multipole electrostatics. Thus, a new standard helix, the 3.6_13/10-, Némethy- or N-helix, is proposed. Due to the use of constraints from monopole force fields and assumed secondary structures used in low-resolution refinement of electron density of proteins, such structures in the PDB often show linear hydrogen bonding.