Conformations of phenylalanine in the tripeptides AFA and GFG probed by combining MD simulations with NMR, FTIR, polarized Raman, and VCD spectroscopy

Solid-state packing dictates the unexpected solubility of aromatic peptides

The understanding and prediction of the solubility of biomolecules, even of the simplest ones, reflect an open question and unmet need. Short aromatic tripeptides are among the most highly aggregative biomolecules. However, in marked contrast, Ala-Phe-Ala (AFA) was surprisingly found to be non-aggregative and could be solubilized at millimolar concentrations. Here, aiming to uncover the underlying molecular basis of its high solubility, we explore in detail the solubility, aggregation propensity, and atomic-level structure of the tripeptide. We demonstrate an unexpectedly high water solubility of AFA reaching 672 mM, two orders of magnitude higher than reported previously. The single crystal structure reveals an anti-parallel β sheet conformation devoid of any aromatic interactions. This study provides clear mechanistic insight into the structural basis of solubility and suggests a simple and feasible tool for its estimation, bearing implications for design of peptide drugs, peptides materials, and advancement of peptide nanotechnology.

Multiconformational analysis of tripeptides upon consideration of implicit and explicit hydration effects

During the last two decades, numerous observed data obtained by various physical techniques, also supported by molecular modeling approaches, have highlighted the structuring features of tripeptides, as well as their aggregation properties. Herein, we focus on the structural dynamics of four trimers, i.e., Gly-Gly-Gly, Gly-Ala-Gly, Ala-Ala-Ala and Ala-Phe-Ala, in an aqueous environment. Density functional theory calculations (DFT) were carried out to assess the stability of four types of secondary structures, i.e., β-strand, polyproline-II (pP-II), α-helix and γ-turn, of which the formation had been described in these tripeptides. Both implicit and explicit hydration effects were analyzed on the conformational and energetic features of trimers. It has been shown that the use of M062X functional (versus B3LYP) improve the stability of intramolecular H-bonds, especially in inverse γ-turn structures, as well as the energetic and conformational equilibrium in all tripeptides. Explicit hydration reflected by the presence of five water molecules around the backbone polar sites (NH₃⁺, N-H, CO and NH₂) considerably changes the conformational landscapes of the trimers. Characteristic intramolecular and intermolecular interactions evidenced by the calculations, were emphasized.

Conformation Dependence of Diphenylalanine Self-Assembly Structures and Dynamics: Insights from Hybrid-Resolution Simulations

The molecular design of peptide-assembled nanostructures relies on extensive knowledge pertaining to the relationship between conformational features of peptide constituents and their behavior regarding self-assembly, and characterizing the conformational details of peptides during their self-assembly is experimentally challenging. Here, we demonstrate that a hybrid-resolution modeling method can be employed to investigate the role that conformation plays during the assembly of terminally capped diphenylalanines (FF) through microsecond simulations of hundreds or thousands of peptides. Our simulations discovered tubular or vesicular nanostructures that were consistent with experimental observation while reproducing critical self-assembly concentration and secondary structure contents in the assemblies that were measured in our experiments. The atomic details provided by our method allowed us to uncover diverse FF conformations and conformation dependence of assembled nanostructures. We found that the assembled morphologies and the molecular packing of FFs in the observed assemblies are linked closely with side-chain angle and peptide bond orientation, respectively. Of various conformations accessible to soluble FFs, only a select few are compatible with the assembled morphologies in water. A conformation resembling a FF crystal, in particular, became predominant due to its ability to permit highly ordered and energetically favorable FF packing in aqueous assemblies. Strikingly, several conformations incompatible with the assemblies arose transiently as intermediates, facilitating key steps of the assembly process. The molecular rationale behind the role of these intermediate conformations were further explained. Collectively, the structural details reported here advance the understanding of the FF self-assembly mechanism, and our method shows promise for studying peptide-assembled nanostructures and their rational design.

Self-assembly of tripeptides into γ-turn nanostructures

Self-assembling phenylalanine-based peptides have garnered interest owing to their potential for creating new functional materials. Here, we designed four diastereomers, l-Phe-l-Phe-l-Phe (FFF), d-Phe-l-Phe-l-Phe (fFF), l-Phe-d-Phe-l-Phe (FfF) and l-Phe-l-Phe-d-Phe (FFf), to analyze the effect of the d-isomer on the self-assembly. Using SEM, TG, VCD, and solid-state NMR measurements, we found that only FFf forms a γ-turn conformation and self-assembles into a nanoplate with higher thermal stability. The supramolecular structure of FFf consists of intra- and intermolecular hydrogen bonds and π-π stackings. From our results, we have discovered that FFf forms a new type of self-assembling γ-turn conformation, clarifying the structural role of a d-amino acid residue in supramolecular formation.

Amino acid conformations control the morphological and chiral features of the self-assembled peptide nanostructures: Young investigators perspective

HYPOTHESIS

A variety of nanostructures with different chiral features can be self-assembled from short peptides with highly similar sequences. We hypothesize that these supramolecular nanostructures are ruled by the constituent amino acid residues which adopt their conformations under the influence of intra-/inter-molecular interactions during peptide self-assembly.

APPROACH

Through reviewing recent advances in the self-assembly of short peptides and focusing on the relationship between amino acid conformations, peptide secondary structures and intra-/inter-molecular interactions within the supramolecular architectures, we aim to rationalize the complex interactive processes involved in the self-assembly of short, designed peptides.

RESULTS

Given the highly complexing interactive processes, the adoption of amino acid conformations and their control over peptide self-assembly consist of 4 main steps: (1) Each amino acid residue adopts its unique conformation in a specific sequence; (2) The sequence exhibits its own main chain geometry and determines the propensity of the intermolecular alignment within the building block; (3) The structural propensity of the building block and the packing mode between them determine the self-assembled structural features such as twisting, growth and chirality; (4) In addition to intra-/inter-molecular interactions, inter-sheet and inter-building block interactions could also affect the residue conformations and nanostructures, causing structural readjustment.

Optimization of Protein Backbone Dihedral Angles by Means of Hamiltonian Reweighting

Molecular dynamics simulations depend critically on the accuracy of the underlying force fields in properly representing biomolecules. Hence, it is crucial to validate the force-field parameter sets in this respect. In the context of the GROMOS force field, this is usually achieved by comparing simulation data to experimental observables for small molecules. In this study, we develop new amino acid backbone dihedral angle potential energy parameters based on the widely used 54A7 parameter set by matching to experimental J values and secondary structure propensity scales. In order to find the most appropriate backbone parameters, close to 100 000 different combinations of parameters have been screened. However, since the sheer number of combinations considered prohibits actual molecular dynamics simulations for each of them, we instead predicted the values for every combination using Hamiltonian reweighting. While the original 54A7 parameter set fails to reproduce the experimental data, we are able to provide parameters that match significantly better. However, to ensure applicability in the context of larger peptides and full proteins, further studies have to be undertaken.

Realistic sampling of amino acid geometries for a multipolar polarizable force field

The Quantum Chemical Topological Force Field (QCTFF) uses the machine learning method kriging to map atomic multipole moments to the coordinates of all atoms in the molecular system. It is important that kriging operates on relevant and realistic training sets of molecular geometries. Therefore, we sampled single amino acid geometries directly from protein crystal structures stored in the Protein Databank (PDB). This sampling enhances the conformational realism (in terms of dihedral angles) of the training geometries. However, these geometries can be fraught with inaccurate bond lengths and valence angles due to artefacts of the refinement process of the X-ray diffraction patterns, combined with experimentally invisible hydrogen atoms. This is why we developed a hybrid PDB/nonstationary normal modes (NM) sampling approach called PDB/NM. This method is superior over standard NM sampling, which captures only geometries optimized from the stationary points of single amino acids in the gas phase. Indeed, PDB/NM combines the sampling of relevant dihedral angles with chemically correct local geometries. Geometries sampled using PDB/NM were used to build kriging models for alanine and lysine, and their prediction accuracy was compared to models built from geometries sampled from three other sampling approaches. Bond length variation, as opposed to variation in dihedral angles, puts pressure on prediction accuracy, potentially lowering it. Hence, the larger coverage of dihedral angles of the PDB/NM method does not deteriorate the predictive accuracy of kriging models, compared to the NM sampling around local energetic minima used so far in the development of QCTFF.

IMPIPS: the immune protection-inducing protein structure concept in the search for steric-electron and topochemical principles for complete fully-protective chemically synthesised vaccine development

Determining immune protection-inducing protein structures (IMPIPS) involves defining the stereo-electron and topochemical characteristics which are essential in MHC-p-TCR complex formation. Modified high activity binding peptides (mHABP) were thus synthesised to produce a large panel of IMPIPS measuring 26.5 ±3.5Å between the farthest atoms fitting into Pockets 1 to 9 of HLA-DRβ1* structures. They displayed a polyproline II-like (PPII_L) structure with their backbone O and N atoms orientated to establish H-bonds with specific residues from HLA-DRβ1*-peptide binding regions (PBR). Residues having specific charge and gauche⁺ orientation regarding p3χ1, p5χ2, and p7χ1 angles determined appropriate rotamer orientation for perfectly fitting into the TCR to induce an appropriate immune response. Immunological assays in Aotus monkeys involving IMPIPS mixtures led to promising results; taken together with the aforementioned physicochemical principles, non-interfering, long-lasting, protection-inducing, multi-epitope, multistage, minimal subunit-based chemically-synthesised peptides can be designed against diseases scourging humankind.

Local order in the unfolded state: conformational biases and nearest neighbor interactions

The discovery of Intrinsically Disordered Proteins, which contain significant levels of disorder yet perform complex biologically functions, as well as unwanted aggregation, has motivated numerous experimental and theoretical studies aimed at describing residue-level conformational ensembles. Multiple lines of evidence gathered over the last 15 years strongly suggest that amino acids residues display unique and restricted conformational preferences in the unfolded state of peptides and proteins, contrary to one of the basic assumptions of the canonical random coil model. To fully understand residue level order/disorder, however, one has to gain a quantitative, experimentally based picture of conformational distributions and to determine the physical basis underlying residue-level conformational biases. Here, we review the experimental, computational and bioinformatic evidence for conformational preferences of amino acid residues in (mostly short) peptides that can be utilized as suitable model systems for unfolded states of peptides and proteins. In this context particular attention is paid to the alleged high polyproline II preference of alanine. We discuss how these conformational propensities may be modulated by peptide solvent interactions and so called nearest-neighbor interactions. The relevance of conformational propensities for the protein folding problem and the understanding of IDPs is briefly discussed.

pH-Independence of trialanine and the effects of termini blocking in short peptides: a combined vibrational, NMR, UVCD, and molecular dynamics study

Several lines of evidence now well establish that unfolded peptides in general, and alanine in specific, have an intrinsic preference for the polyproline II (pPII) conformation. Investigation of local order in the unfolded state is, however, complicated by experimental limitations and the inherent dynamics of the system, which has in some cases yielded inconsistent results from different types of experiments. One method of studying these systems is the use of short model peptides, and specifically short alanine peptides, known for predominantly sampling pPII structure in aqueous solution. Recently, He et al. ( J. Am. Chem. Soc. 2012 , 134 , 1571 - 1576 ) proposed that unblocked tripeptides may not be suitable models for studying conformational propensities in unfolded peptides due to the presence of end effect, that is, electrostatic interactions between investigated amino acid residues and terminal charges. To determine whether changing the protonation states of the N- and C-termini influence the conformational manifold of the central amino acid residue in tripeptides, we have examined the pH-dependence of unblocked trialanine and the conformational preferences of alanine in the alanine dipeptide. To this end, we measured and globally analyzed amide I' band profiles and NMR J-coupling constants. We described conformational distributions as the superposition of two-dimensional Gaussian distributions assignable to specific subspaces of the Ramachandran plot. Results show that the conformational ensemble of trialanine as a whole, and the pPII content (χpPII = 0.84) in particular, remains practically unaffected by changing the protonation state. We found that compared to trialanine, the alanine dipeptide has slightly lower pPII content (χpPII = 0.74) and an ensemble more reminiscent of the unblocked Gly-Ala-Gly model peptide. In addition, a two-state thermodynamic analysis of the conformational sensitive Δε(T) and (3)J(H(N)H(α))(T) data obtained from electronic circular dichroism and H NMR spectra indicate that the free energy landscape of trialanine is similar in all protonation states. MD simulations for the investigated peptides corroborate this notion and show further that the hydration shell around unblocked trialanine is unaffected by the protonation/deprotonation of the C-terminal group. In contrast, the alanine dipeptide shows a reduced water density around the central residue as well as a less ordered hydration shell, which decreases the pPII propensity and reduces the lifetime of sampled conformations.

Disorder and order in unfolded and disordered peptides and proteins: a view derived from tripeptide conformational analysis. I. Tripeptides with long and predominantly hydrophobic side chains

We performed a conformational analysis of the central residues of three tripeptides glycyl-L-isoleucyl-glycine (GIG), glycyl-L-tyrosyl-glycine (GYG) and glycyl-L-arginyl-glycine (GRG) in aqueous solution, based on a global analysis of amide I' band profiles and NMR J-coupling constants. The results are compared with recently reported distributions of GVG, GFG and GEG. For GIG and GYG, we found that even though the polyproline II (pPII) fraction is below 0.5, it is still the most populated conformation, whereas GVG and GFG show both a larger β-strand fraction. For GRG, we observed a clear dominance of pPII over β-strand, reminiscent of observations for GEG and GKG. This finding indicates that terminal charges on otherwise hydrophobic residue side chains stabilize pPII over β-strand conformations. For all peptides investigated we found that a variety of compact and turn-like conformations constitute nearly 20 percent of their conformational distributions. Attempts to analyze our data with a simple two-state pPII-->/<--β model therefore do not yield any satisfactory reproduction of experimental results. A comparison of the obtained GxG ensembles with conformational distributions of GxG segments in truncated coil libraries (helices and sheets omitted) revealed a much larger fraction of type II β(i+2) and type III β like conformations for the latter. Thus, a comparison of conformational distributions of unfolded peptide segments in solution and in coil libraries reveal interesting information on how the interplay between intrinsic propensities of amino acid residues and non-local interactions in polypeptide chains determine the conformations of loop segments in proteins.

Electrostatics of hydrogen exchange for analyzing protein flexibility

Electrostatic interactions at the protein-aqueous interface modulate the reactivity of solvent-exposed backbone amides by a factor of at least a billion fold. The brief (∼10 ps) lifetime of the peptide anion formed during the hydroxide-catalyzed exchange reaction helps enable the experimental rates to be robustly predictable by continuum dielectric methods. Since this ability to predict the structural dependence of exchange reactivity also applies to the protein amide hydrogens that are only rarely exposed to the bulk solvent phase, electrostatic analysis of the experimental exchange rates provides an effective assessment of whether a given model ensemble is consistent with the properly weighted Boltzmann conformational distribution of the protein native state.

Assessment of the intrinsic conformational preferences of dipeptide amino acids in aqueous solution by combined umbrella sampling/MBAR statistics. A comparison with experimental results

The propensities of 19 amino acid dipeptides have been calculated by a distributed umbrella sampling molecular dynamics simulation procedure using the OPLS-AA force field. The potential of mean force maps was estimated with the multiple Bennett acceptance ratio statistics. The resulting propensities compare satisfactorily well with very recently published experimental data on equivalent systems. In particular, α conformation-probabilities for all of the dipeptides remain much lower than either β or P(II) propensities. This result is in agreement with most experimental data for dipeptides. However, it is also in contrast with most simulation studies performed so far with other force fields, where α conformations result even more probable than P(II) or β ones. We discuss the behavior of the OPLS-AA force field, which can be useful for the improvement of this model in reproducing the recent experimental observations on amino acid dipeptides.

Sidechain conformational dependence of hydrogen exchange in model peptides

Peptide hydrogens that are exposed to solvent in protein X-ray structures exhibit a billion-fold range in hydroxide-catalyzed exchange rates, and these rates have previously been shown to be predictable by continuum dielectric methods to within a factor of 7, based on single protein conformations. When using a protein coil library to model the Boltzmann-weighted conformational distribution for the various N-acetyl-[X-Ala]-N-methylamides and N-acetyl-[Ala-Y]-N-methylamides, the acidity of the central amide in the individual conformers of each peptide spans nearly a million-fold range. Nevertheless, population averaging of these conformer acidities predicts the standard sidechain-dependent hydrogen exchange correction factors for nonpolar model peptides to within a factor of 30% (10(0.11)) with a correlation coefficient r=0.91. Comparison with the analogous continuum dielectric calculations for the other N-acetyl-[X-Y]-N-methylamides indicates that deviations from the isolated residue hypothesis of classical polymer theory predict appreciable errors in the exchange rates for conformationally disordered peptides when the standard sidechain-dependent hydrogen exchange rate correction factors are assumed to be independently additive. Although electronic polarizability generally dominates the dielectric shielding for the approximately 10ps lifetime of peptide ionization, evidence is presented for modest contributions from rapid intrarotamer conformational reorganization of Asn and Gln sidechains.