Conformational Manifold of α-Aminoisobutyric Acid (Aib) Containing Alanine-Based Tripeptides in Aqueous Solution Explored by Vibrational Spectroscopy, Electronic Circular Dichroism Spectroscopy, and Molecular Dynamics Simulations

Strategic Single-Residue Substitution in the Antimicrobial Peptide Esc(1-21) Confers Activity against Staphylococcus aureus, Including Drug-Resistant and Biofilm Phenotype

Staphylococcus aureus, a bacterium resistant to multiple drugs, is a significant cause of illness and death worldwide. Antimicrobial peptides (AMPs) provide an excellent potential strategy to cope with this threat. Recently, we characterized a derivative of the frog-skin AMP esculentin-1a, Esc(1-21) (1) that is endowed with potent activity against Gram-negative bacteria but poor efficacy against Gram-positive strains. In this study, three analogues of peptide 1 were designed by replacing Gly8 with α-aminoisobutyric acid (Aib), Pro, and dPro (2-4, respectively). The single substitution Gly8 → Aib8 in peptide 2 makes it active against the planktonic form of Gram-positive bacterial strains, especially Staphylococcus aureus, including multidrug-resistant clinical isolates, with an improved biostability without resulting in cytotoxicity to mammalian cells. Moreover, peptide 2 showed a higher antibiofilm activity than peptide 1 against both reference and clinical isolates of S. aureus. Peptide 2 was also able to induce rapid bacterial killing, suggesting a membrane-perturbing mechanism of action. Structural analysis of the most active peptide 2 evidenced that the improved biological activity of peptide 2 is the consequence of a combination of higher biostability, higher α helical content, and ability to reduce membrane fluidity and to adopt a distorted helix, bent in correspondence of Aib8. Overall, this study has shown how a strategic single amino acid substitution is sufficient to enlarge the spectrum of activity of the original peptide 1, and improve its biological properties for therapeutic purposes, thus paving the way to optimize AMPs for the development of new broad-spectrum anti-infective agents.

Development of Orally Bioavailable Peptides Targeting an Intracellular Protein: From a Hit to a Clinical KRAS Inhibitor

Cyclic peptides as a therapeutic modality are attracting a lot of attention due to their potential for oral absorption and accessibility to intracellular tough targets. Here, starting with a drug-like hit discovered using an mRNA display library, we describe a chemical optimization that led to the orally available clinical compound known as LUNA18, an 11-mer cyclic peptide inhibitor for the intracellular tough target RAS. The key findings are as follows: (i) two peptide side chains were identified that each increase RAS affinity over 10-fold; (ii) physico-chemical properties (PCP) including Clog P can be adjusted by side-chain modification to increase membrane permeability; (iii) restriction of cyclic peptide conformation works effectively to adjust PCP and improve bio-activity; (iv) cellular efficacy was observed in peptides with a permeability of around 0.4 × 10^-6 cm/s or more in a Caco-2 permeability assay; and (v) while keeping the cyclic peptide's main-chain conformation, we found one example where the RAS protein structure was changed dramatically through induced-fit to our peptide side chain. This study demonstrates how the chemical optimization of bio-active peptides can be achieved without scaffold hopping, much like the processes for small molecule drug discovery that are guided by Lipinski's rule of five. Our approach provides a versatile new strategy for generating peptide drugs starting from drug-like hits.

Rational design and topochemical synthesis of polymorphs of a polymer

Packing a polymer in different ways can give polymorphs of the polymer having different properties. β-Turn forming peptides such as 2-aminoisobutyric acid (Aib)-rich peptides adopt several conformations by varying the dihedral angles. Aiming at this, a β-turn-forming peptide monomer would give different polymorphs and these polymorphs upon topochemical polymerization would yield polymorphs of the polymer, we designed an Aib-rich monomer N₃-(Aib)₃-NHCH₂-C[triple bond, length as m-dash]CH. This monomer crystallizes as two polymorphs and one hydrate. In all forms, the peptide adopts β-turn conformations and arranges in a head-to-tail manner with their azide and alkyne units proximally placed in a ready-to-react alignment. On heating, both the polymorphs undergo topochemical azide-alkyne cycloaddition polymerization. Polymorph I polymerized in a single-crystal-to-single-crystal (SCSC) fashion and the single-crystal X-ray diffraction analysis of the polymer revealed its screw-sense reversing helical structure. Polymorph II maintains its crystallinity during polymerization but gradually becomes amorphous upon storage. The hydrate III undergoes a dehydrative transition to polymorph II. Nanoindentation studies revealed that different polymorphs of the monomer and the corresponding polymers exhibited different mechanical properties, in accordance with their crystal packing. This work demonstrates the promising future of the marriage of polymorphism and topochemistry for obtaining polymorphs of polymers.

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.

Unique supramolecular assembly of a synthetic achiral α, γ-hybrid tripeptide

An achiral tripeptide, namely, Boc-γ-Abu-m-ABA-Aib-OMe (γ-Abu: γ−amino butyric acid; m-ABA: meta-aminobenzoic acid) was synthesized by solution phase procedure. The α, γ-hybrid peptide was designed in such a way that two dissimilar γ−amino acids, one flexible and another rigid, were positioned sidewise along with α-amino isobutyric acid (Aib) as C-terminal residue. The single crystal X-ray diffraction analysis revealed that two kinks were generated around centrally placed m-ABA. Interestingly, the peptide self-assembled via three intermolecular N–H···O and one intermolecular C–H···O hydrogen bonding interactions to supramlecular helical architecture.

Cα-Methyl-l-valine: A Preferential Choice over α-Aminoisobutyric Acid for Designing Right-Handed α-Helical Scaffolds

In synthetic peptides containing Gly and coded α-amino acids, one of the most common practices to enhance their helical extent is to incorporate a large number of l-Ala residues along with noncoded, strongly foldameric α-aminoisobutyric acid (Aib) units. Earlier studies have established that Aib-based peptides, with propensity for both the 3₁₀- and α-helices, have a tendency to form ordered three-dimensional structure that is much stronger than that exhibited by their l-Ala rich counterparts. However, the achiral nature of Aib induces an inherent, equal preference for the right- and left-handed helical conformations as found in Aib homopeptide stretches. This property poses challenges in the analysis of a model peptide helical conformation based on chirospectroscopic techniques like electronic circular dichroism (ECD), a very important tool for assigning secondary structures. To overcome such ambiguity, we have synthesized and investigated a thermally stable 14-mer peptide in which each of the Aib residues of our previously designed and reported analogue ABGY (where B stands for Aib) is replaced by C^α-methyl-l-valine (L-AMV). Analysis of the results described here from complementary ECD and ¹H nuclear magnetic resonance spectroscopic techniques in a variety of environments firmly establishes that the L-AMV-containing peptide exhibits a significantly stronger preference compared to that of its Aib parent in terms of conferring α-helical character. Furthermore, being a chiral α-amino acid, L-AMV shows an intrinsic, extremely strong bias for a quite specific (right-handed) screw sense. These findings emphasize the relevance of L-AMV as a more appropriate unit for the design of right-handed α-helical peptide models that may be utilized as conformationally constrained scaffolds.

Constrained Peptides with Fine-Tuned Flexibility Inhibit NF-Y Transcription Factor Assembly

Protein complex formation depends on the interplay between preorganization and flexibility of the binding epitopes involved. The design of epitope mimetics typically focuses on stabilizing a particular bioactive conformation, often without considering conformational dynamics, which limits the potential of peptidomimetics against challenging targets such as transcription factors. We developed a peptide-derived inhibitor of the NF-Y transcription factor by first constraining the conformation of an epitope through hydrocarbon stapling and then fine-tuning its flexibility. In the initial set of constrained peptides, a single non-interacting α-methyl group was observed to have a detrimental effect on complex stability. Biophysical characterization revealed how this methyl group affects the conformation of the peptide in its bound state. Adaption of the methylation pattern resulted in a peptide that inhibits transcription factor assembly and subsequent recruitment to the target DNA.

New 19-Residue Peptaibols from Trichoderma Clade Viride

Trichoderma koningiopsis and T. gamsii belong to clade Viride of Trichoderma, the largest and most diverse group of this genus. They produce a wide range of bioactive secondary metabolites, including peptaibols with antibacterial, antifungal, and antiviral properties. The unusual amino acid residues of peptaibols, i.e., α-aminoisobutyric acid (Aib), isovaline (Iva), and the C-terminal 1,2-amino alcohol make them unique among peptides. In this study, the peptaibiomes of T. koningiopsis and T. gamsii were investigated by HPLC-ESI-MS. The examined strains appeared to produce 19-residue peptaibols, most of which are unknown from literature, but their amino acid sequences are similar to those of trikoningins, tricholongins, trichostrigocins, trichorzianins, and trichorzins. A new group of peptaibols detected in T. koningiopsis are described here under the name “Koningiopsin”. Trikoningin KA V, the closest peptaibol compound to the peptaibols produced by these two strains, was selected for structural investigation by short MD simulation, which revealed that many residues show high preference for left handed helix formation. The bioactivity of the peptaibol mixtures produced by T. koningiopsis and T. gamsii was tested on agar plates against bacteria, yeasts, and filamentous fungi. The results revealed characteristic differences in bioactivities towards the different groups of target microorganisms, which can be explained with the differences in their cell wall structures.

A spin-1 representation for dual-funnel energy landscapes

The interconversion between the left- and right-handed helical folds of a polypeptide defines a dual-funneled free energy landscape. In this context, the funnel minima are connected through a continuum of unfolded conformations, evocative of the classical helix-coil transition. Physical intuition and recent conjectures suggest that this landscape can be mapped by assigning a left- or right-handed helical state to each residue. We explore this possibility using all-atom replica exchange molecular dynamics and an Ising-like model, demonstrating that the energy landscape architecture is at odds with a two-state picture. A three-state model-left, right, and unstructured-can account for most key intermediates during chiral interconversion. Competing folds and excited conformational states still impose limitations on the scope of this approach. However, the improvement is stark: Moving from a two-state to a three-state model decreases the fit error from 1.6 kBT to 0.3 kBT along the left-to-right interconversion pathway.

Various modifications of the amphipathic dynorphin A pharmacophore for rat brain bradykinin receptors

As a unique endogenous opioid ligand, dynorphin A shows paradoxical neuroexcitatory effects at bradykinin receptors, and the effects are known to be amplified by the upregulation of dynorphin A under chronic pain and inflammatory conditions. In our earlier structure-activity relationship studies, the amphipathic dynorphin A fragment, [Des-Arg(7) ]-Dyn A-(4-11), was identified as a pharmacophore for the bradykinin receptors along with key structural features. Here, further modifications of the pharmacophore showed that the position of a Pro residue is also an important feature because of its role in making (or disrupting) a β-turn or 310 helix structure which is crucial for receptor recognition.

Peptide Backbone Composition and Protease Susceptibility: Impact of Modification Type, Position, and Tandem Substitution

The clinical utility of peptides is limited by their rapid degradation by endogenous proteases. Modification of the peptide backbone can generate functional analogues with enhanced proteolytic stability. Existing principles for the design of such oligomers have focused primarily on effective structural mimicry. A more robust strategy would incorporate a rational approach for engineering maximal proteolytic stability with minimal unnatural residue content. We report here the systematic comparison of the proteolytic resistance imparted by four backbone modifications commonly employed in the design of protease-stable analogues of peptides with complex folding patterns. The degree of protection was quantified as a function of modification type, position, and tandem substitution in the context of a long, unstructured host sequence and a canonical serine protease. These results promise to inform ongoing work to develop biostable mimics of increasingly complex peptides and proteins.

Conformational and thermodynamic properties of non-canonical α,α-dialkyl glycines in the peptaibol Alamethicin: molecular dynamics studies

In this work, we investigate the structure, dynamic and thermodynamic properties of noncanonical disubstituted amino acids (α,α-dialkyl glycines), also known as non-natural amino acids, in the peptaibol Alamethicin. The amino acids under study are Aib (α-amino isobutyric acid or α-methyl alanine), Deg (α,α-diethyl glycine), Dpg (α,α-dipropyl glycine), Dibg (α,α-di-isobutyl glycine), Dhg (α,α-dihexyl glycine), DΦg (α,α-diphenyl glycine), Dbzg (α,α-dibenzyl glycine), Ac6c (α,α-cyclohexyl glycine), and Dmg (α,α-dihydroxymethyl glycine). It is hypothesized that these amino acids are able to induce well-defined secondary structure in peptidomimetics. To test this hypothesis, new peptidomimetics of Alamethicin were constructed by replacing the native Aib positions of Alamethicin by one or more new α,α-dialkyl glycines. Dhg and Ac6c demonstrated the capacity to induce well-defined α-helical structures. Dhg and Ac6c also promote the thermodynamic stabilization of these peptides in a POPC model membrane and are better alternatives to the Aib in Alamethicin. These noncanonical amino acids also improved secondary structure properties, revealing preorganization in water and maintenance of α helical structure in POPC. We show that it is possible to optimize the helicity and thermodynamic properties of native Alamethicin, and we suggest that these amino acids could be incorporated in other peptides with similar structural effect.

The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins

Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.

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.

(S)-N-[(4-{(S)-1-[2-(4-Meth-oxy-benz-amido)-2-methyl-propano-yl]pyrrolidine-2-carboxamido}-3,4,5,6-tetra-hydro-2H-pyran-4-yl)carbon-yl]proline dimethyl sulfoxide monosolvate (4-MeBz-Aib-Pro-Thp-Pro-OH)

The asymmetric unit of the title compound, C28H38N4O8·C2H6OS, contains one tetra-peptide and one disordered dimethyl sulfoxide (DMSO) mol-ecule. The central five-membered ring (Pro(2)) of the peptide mol-ecule has a disordered envelope conformation [occupancy ratio 0.879 (2):0.121 (2)] with the envelope flap atom, the central C atom of the three ring methylene groups, lying on alternate sides of the mean ring plane. The terminal five-membered ring (Pro(4)) also adopts an envelope conformation with the C atom of the methylene group closest to the carboxylic acid function as the envelope flap, and the six-membered tetra-hydro-pyrane ring shows a chair conformation. The tetra-peptide exists in a helical conformation, stabilized by an intra-molecular hydrogen bond between the amide N-H group of the heterocyclic α-amino acid Thp and the amide O atom of the 4-meth-oxy-benzoyl group. This inter-action has a graph set motif of S(10) and serves to maintain a fairly rigid β-turn structure. In the crystal, the terminal hy-droxy group forms a hydrogen bond with the amide O atom of Thp of a neighbouring mol-ecule, and the amide N-H group at the opposite end of the mol-ecule forms a hydrogen bond with the amide O atom of Thp of another neighbouring mol-ecule. The combination of both inter-molecular inter-actions links the mol-ecules into an extended three-dimensional framework.

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.

1H-Azepine-4-amino-4-carboxylic acid: a new α,α-disubstituted ornithine analogue capable of inducing helix conformations in short Ala-Aib pentapeptides

A very efficient synthesis of orthogonally protected 1H-azepine-4-amino-4-carboxylic acid, abbreviated as Azn, a conformationally restricted analogue of ornithine, was realized. It was obtained on a gram scale in good overall yield in five steps, three of which did not require isolation of the intermediates, starting from the readily available 1-amino-4-oxo-cyclohexane-4-carboxylic acid. Both enantiomers were used for the preparation of pentapeptide models containing Ala, Aib, and Azn. Conformational studies using both spectroscopic techniques (NMR, CD) and molecular dynamics on model 5-mer peptides showed that the (R)-Azn isomer possesses a marked helicogenic effect.

Structure analysis of unfolded peptides I: vibrational circular dichroism spectroscopy

Vibrational circular dichroism (VCD) spectroscopy is an invaluable spectroscopic techniques utilized to exploit the optical strength of vibrational transitions for structure analysis. In this chapter, we describe the protocol for measuring and self-consistently analyzing VCD and the corresponding FT-IR spectra of short peptides. This process involves the decomposition of the IR spectrum as well as simulations of the amide I band profiles in both spectra based on structural models of the peptides investigated. This type of spectral analysis should be complemented with similar investigations of Raman spectra, which are described in the subsequent chapter. The structural analysis of short, unfolded peptides described in this chapter can easily be extended for the analysis of longer unfolded peptides or even proteins. This is particularly important in view of the demonstrated biological relevance of intrinsically disordered peptides and proteins (IDPs).

N-(tert-butoxycarbonyl)-O-allyl-L-seryl-α-aminoisobutyryl-L-valine methyl ester: a protected tripeptide with an allylated serine residue

The title compound [systematic name (6S,12S)-methyl 6-(allyloxymethyl)-12-isopropyl-2,2,9,9-tetramethyl-4,7,10-trioxo-3-oxa-5,8,11-triazatridecan-13-oate], C(21)H(37)N(3)O(7), containing the little studied O-allyl-L-serine residue [Ser(All)], crystallizes in the monoclinic space group C2 with one molecule in the asymmetric unit. The compound is an analogue of the Ser140-Val142 segment of the water channel aquaporin-4 (AQP4). It forms a distorted type-II β-turn with a P(II)-3(10L)-P(II) backbone conformation (P(II) is polyproline II). The overall backbone conformation is markedly different from that of the CO(Pro139)-Val142 stretch of rat AQP4, but is quite similar to the corresponding segment of human AQP4, despite significant differences at the level of the individual residues. The side chain of the Ser(All) residue adopts a gauche conformation relative to the backbone CO-C(α) and C(α)-N bonds. The H atoms of the two CH(2) groups in the Ser(All) side chain are almost eclipsed. The crystal packing of the title compound is divided into one-molecule-thick layers, each layer having a hydrophilic core and distinct hydrophobic interfaces on either side.

Sulfate ion interaction with 'anion recognition' short peptide motif at the N-terminus of an isolated helix: A conformational landscape

Anion-binding motifs in proteins are typically conserved in sequence and conformation. Crystal structural studies have shown that such motifs often occur in loop regions preceding a helix and interaction with the anions can induce their well defined conformational changes. In order to understand the properties of such motifs in isolation, we have synthesized an 18-residue chimeric polypeptide whose C-terminal part is a designed helix and its N-terminal consists of a C(alpha)NN anion binding structural motif containing residues Leu-Gly-Lys-Gln (residues 107-110 of protein DNA-glycosylase). We present evidence for the interaction of a sulfate (SO(4)(2-)) ion with the L-G-K-Q segment using complementary spectroscopic techniques. Moreover, upon interaction with SO(4)(2-) ion the N-terminal L-G-K-Q segment undergoes a non-helical to helical transition similar to what is observed in protein crystal structure. This work clearly demonstrates the "local" nature of anion binding and the accompanying conformational change that helps in understanding the influence of sequence/structural context of anion binding in proteins.

Conformational preferences of a short Aib/Ala-based water-soluble peptide as a function of temperature

The amino acid Aib predisposes a peptide to be helical with context-dependent preference for either 3(10)- or alpha- or a mixed helical conformation. Short peptides also show an inherent tendency to be unfolded. To characterize helical and unfolded states adopted by water-soluble Aib-containing peptides, the conformational preference of Ac-Ala-Aib-Ala-Lys-Ala-Aib-Lys-Ala-Lys-Ala-Aib-Tyr-NH(2) was determined by CD, NMR and MD simulations as a function of temperature. Temperature-dependent CD data indicated the contribution of two major components, each an admixture of helical and extended/polyproline II structures. Both right- and left-handed helical conformations were detected from deconvolution of CD data and (13)C NMR experiments. The presence of a helical backbone, more pronounced at the N-terminal, and a temperature-induced shift in alpha-helix/3(10)-helix equilibrium, more pronounced at the C-terminal, emerged from NMR data. Starting from polyproline II, the N-terminal of the peptide folded into a helical backbone in MD simulations within 5 ns at 60 degrees C. Longer simulations showed a mixed-helical backbone to be stable over the entire peptide at 5 degrees C while at 60 degrees C the mixed-helix was either stable at the N-terminus or occurred in short stretches through out the peptide, along with a significant population of polyproline II. Our results point towards conformational heterogeneity of water-soluble Aib-based peptide helices and the associated subtleties. The problem of analyzing CD and NMR data of both left- and right-handed helices are discussed, especially the validity of the ellipticity ratio [theta](222)/[theta](207), as a reporter of alpha-/3(10)- population ratio, in right- and left-handed helical mixtures.

Modeling of peptides containing D-amino acids: implications on cyclization

Cyclic peptides are therapeutically attractive due to their high bioavailability, potential selectivity, and scaffold novelty. Furthermore, the presence of D-residues induces conformational preferences not followed by peptides consisting of naturally abundant L-residues. Therefore, comprehending how amino acids induce turns in peptides, subsequently facilitating cyclization, is significant in peptide design. Here, we performed 20-ns explicit-solvent molecular dynamics simulations for three diastereomeric peptides with stereochemistries: LLLLL, LLLDL, and LDLDL. Experimentally LLLLL and LDLDL readily cyclize, whereas LLLDL cyclizes in low yield. Simulations at 310 K produced conformations with inter-terminal hydrogen bonds that correlated qualitatively with the experimental cyclization trend. Energies obtained for representative structures from quantum chemical (B3LYP/PCM/cc-pVTZ//HF/6-31G*) calculations predicted pseudo-cyclic and extended conformations as the most stable for LLLLL and LLLDL, respectively, in agreement with the experimental data. In contrast, the most stable conformer predicted for peptide LDLDL was not a pseudo-cyclic structure. Moreover, D-residues preferred the experimentally less populated alpha(L) rotamers even when simulations were performed at a higher temperature and with strategically selected starting conformations. Energies calculated with molecular mechanics were consistent only with peptide LLLLL. Thus, the conformational preferences obtained for the all L: -amino acid peptide were in agreement with the experimental observations. Moreover, refinement of the force field is expected to provide far-reaching conformational sampling of peptides containing D-residues to further develop force field-based conformational-searching methods.

Optimized molecular dynamics force fields applied to the helix-coil transition of polypeptides

Obtaining the correct balance of secondary structure propensities is a central priority in protein force-field development. Given that current force fields differ significantly in their alpha-helical propensities, a correction to match experimental results would be highly desirable. We have determined simple backbone energy corrections for two force fields to reproduce the fraction of helix measured in short peptides at 300 K. As validation, we show that the optimized force fields produce results in excellent agreement with nuclear magnetic resonance experiments for folded proteins and short peptides not used in the optimization. However, despite the agreement at ambient conditions, the dependence of the helix content on temperature is too weak, a problem shared with other force fields. A fit of the Lifson-Roig helix-coil theory shows that both the enthalpy and entropy of helix formation are too small: the helix extension parameter w agrees well with experiment, but its entropic and enthalpic components are both only about half the respective experimental estimates. Our structural and thermodynamic analyses point toward the physical origins of these shortcomings in current force fields, and suggest ways to address them in future force-field development.

Helix forming tendency of valine substituted poly-alanine: a molecular dynamics investigation

In this study, classical molecular dynamics simulations have been carried out on the valine (guest) substituted poly alanine (host) using the host-guest peptide approach to understand the role of valine in the formation and stabilization of helix. Valine has been substituted in the host peptide starting from N terminal to C terminal. Various structural parameters have been obtained from the molecular dynamics simulation to understand the tolerance of helical motif to valine. Depending on the position of valine in the host peptide, it stabilizes (or destabilizes) the formation of the helical structure. The substitution of valine in the poly alanine at some positions has no effect on the helix formation (deformation). It is interesting to observe the coexistence of 3 10 and alpha-helix in the peptides due to the dynamical nature of the hydrogen bonding interaction and sterical interactions.

Interaction of a tripeptide with cesium perfluorooctanoate micelles

The interaction of alanyl-phenylalanyl-alanine (Ala-Phe-Ala) with the micelles formed by cesium perfluorooctanoate (CsPFO) in water was studied in the isotropic phase by means of 1H NMR and by molecular dynamics (MD) simulations. Information on the location of the peptide was experimentally obtained from selective variations in Ala-Phe-Ala chemical shifts and from differential line broadening in the presence of the paramagnetic ion Mn2+. The peptide-micelle association constant was estimated analyzing the chemical shift variations of the most sensitive Ala-Phe-Ala resonances with the peptide concentration. MD simulations of Ala-Phe-Ala in the micellar environment confirmed the experimental observations, identifying the hydrogen bonding interactions of the different peptide moieties with the micelle, yielding a binding constant close to the experimental one. NOESY experiments suggest that the peptide in the micellar environment does not adopt a preferred conformation but is mainly unstructured. Details on the conformational behavior of the peptide in the micellar solution observed through MD were consistent with a different conformational equilibrium in the proximity of the micelle. Information on Ala-Phe-Ala dynamics was obtained from 1H T1 data and compared to MD simulation results on the overall tumbling motion.