Importance of the single amino acid potential in water for secondary and tertiary structures of proteins

The structure of a protein molecule is considered to be primarily determined by the inter-amino-acid nonbonded interactions, such as hydrogen bonds. However, the conformational space of the polypeptide chain should be simultaneously restricted by the intrinsic conformational preferences of the individual amino acids. We present here precise single amino acid potential (SAAP) surfaces for glycine (For-Gly-NH(2)) and alanine (For-Ala-NH(2)) in water (epsilon = 78.39) and ether (epsilon = 4.335), which were calculated at the HF/6-31+G(d,p) level applying the self-consistent isodensity polarizable continuum model (SCIPCM) reaction field with geometry optimization in the corresponding solvents. The obtained Ramachandran potential surfaces in water showed distinct potential wells in the alpha- and beta-regions. The profiles were in almost perfect agreement with the Ramachandran plots of glycine and alanine residues in folded proteins, suggesting the Boltzmann distributions on the SAAP surfaces. Molecular simulations of polyalanines (For-Ala(n)-NH(2); n = 3-5) by using the SAAP force field equipped with the SCIPCM potentials revealed that the polyalanines readily form 3(10)-helical structures in water but not in vacuo. In ether (hydrophobic environments), the helical structures were relatively stable, but the most stable structure was assigned to a different one. These results indicated that the intrinsic conformational preferences of the individual amino acids (i.e., the SAAPs) in water are of significant importance not only for describing conformations of a polypeptide chain in the random coil state but also for understanding the folding to the secondary and tertiary structures.     

The SAAP force field. A simple approach to a new all-atom protein force field by using single amino acid potential (SAAP) functions in various solvents

A simple strategy to compose a new all-atom protein force field (named as the SAAP force field), which utilizes the single amino acid potential (SAAP) functions obtained in various solvents by ab initio molecular orbital calculation applying the isodensity polarizable continuum model (IPCM), is presented. We considered that the total energy function of a protein force field (E(TOTAL)) is divided into three components; a single amino acid potential term (E(SAAP)), an interamino acid nonbonded interaction term (E(INTER)), and a miscellaneous term (E(OTHERS)), which is ignored (or considered to be constant) at the current version of the force field. The E(INTER) term consists of electrostatic interactions (E(ES')) and van der Waals interactions (E(LJ')). Despite simplicity, the SAAP force field implicitly involves the correlation among individual terms of the Lifson's potential function within a single amino acid unit and can treat solvent effects unambiguously by choosing the SAAP function in an appropriate solvent and the dielectric constant (D) of medium. Application of the SAAP force field to the Monte Carlo simulation of For-Ala(2)-NH(2) in vacuo reasonably reproduced the results of the extensive conformational search by ab initio molecular orbital calculation. In addition, the preliminary Monte Carlo simulations for For-Gly(10)-NH(2) and For-Ala(10)-NH(2) showed reversible transitions from the extended to the pseudosecondary structures in water (D = 78.39) as well as in ether (D = 4.335). The result suggested that the new approach is efficient for fast modeling of protein structures in various environments. Decomposition analysis of the total energy function (E(TOTAL)) by using the SAAP force field suggested that conformational propensities of single amino acids (i.e., the E(SAAP) term) may play definitive roles on the topologies of protein secondary structures.     

Multicanonical Monte Carlo calculation of the free‐energy map of the base–amino acid interaction

The specific interactions between base pairs and amino acids were studied by the multicanonical Monte Carlo method. We sampled numerous interaction configurations and side‐chain conformations of the amino acid by the multicanonical algorithm, and calculated the free energies of the interactions between an amino acid at given C_α positions and a fixed base pair. The contour maps of free energy derived from this calculation represent the preferred C_α position of the amino acid around the base, and these maps of various combinations of bases and amino acids can be used to quantify the specificity of intrinsic base–amino acid interactions. Similarly, enthalpy and entropy maps will provide further details of the specific interactions. We have also calculated the free‐energy map of the orientations of the C_α C_β bond vector, which indicates the preferential orientation of the amino acid against the base. We compared the results obtained by the multicanonical method with those of the exhaustive sampling and canonical Monte Carlo methods. The free‐energy map of the base–amino acid interaction obtained by the multicanonical simulation method was nearly identical to the accurate result derived from the exhaustive sampling method. This indicates that a single multicanonical Monte Carlo simulation can produce an accurate free‐energy map. Multicanonical Monte Carlo sampling produced free‐energy maps that were more accurate than those produced by canonical Monte Carlo sampling. Thus, the multicanonical Monte Carlo method can serve as a powerful tool for estimating the free‐energy landscape of base–amino acid interactions and for elucidating the mechanism by which amino acids of proteins recognize particular DNA base pairs. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 954–962, 2000     

Electron capture dissociation product ion abundances at the X amino acid in RAAAA-X-AAAAK peptides correlate with amino acid polarity and radical stability

We present mechanistic studies aimed at improving the understanding of the product ion formation rules in electron capture dissociation (ECD) of peptides and proteins in Fourier transform ion cyclotron resonance mass spectrometry. In particular, we attempted to quantify the recently reported general correlation of ECD product ion abundance (PIA) with amino acid hydrophobicity. The results obtained on a series of model H-RAAAAXAAAAK-OH peptides confirm a direct correlation of ECD PIA with X amino acid hydrophobicity and polarity. The correlation factor (R) exceeds 0.9 for 12 amino acids (Ile, Val, His, Asn, Asp, Glu, Gln, Ser, Thr, Gly, Cys, and Ala). The deviation of ECD PIA for seven outliers (Pro is not taken into consideration) is explained by their specific radical stabilization properties (Phe, Trp, Tyr, Met, and Leu) and amino acid basicity (Lys, Arg). Phosphorylation of Ser, Thr, and Tyr decreases the efficiency of ECD around phosphorylated residues, as expected. The systematic arrangement of amino acids reported here indicates a possible route toward development of a predictive model for quantitative electron capture/transfer dissociation tandem mass spectrometry, with possible applications in proteomics.     

A genetic algorithm encoded with the structural information of amino acids and dipeptides for efficient conformational searches of oligopeptides

The genetic algorithm (GA) is an intelligent approach for finding minima in a highly dimensional parametric space. However, the success of GA searches for low energy conformations of biomolecules is rather limited so far. Herein an improved GA scheme is proposed for the conformational search of oligopeptides. A systematic analysis of the backbone dihedral angles of conformations of amino acids (AAs) and dipeptides is performed. The structural information is used to design a new encoding scheme to improve the efficiency of GA search. Local geometry optimizations based on the energy calculations by the density functional theory are employed to safeguard the quality and reliability of the GA structures. The GA scheme is applied to the conformational searches of Lys, Arg, Met‐Gly, Lys‐Gly, and Phe‐Gly‐Gly representative of AAs, dipeptides, and tripeptides with complicated side chains. Comparison with the best literature results shows that the new GA method is both highly efficient and reliable by providing the most complete set of the low energy conformations. Moreover, the computational cost of the GA method increases only moderately with the complexity of the molecule. The GA scheme is valuable for the study of the conformations and properties of oligopeptides. © 2016 Wiley Periodicals, Inc.     

Amino acids of ricin and its polypeptides

Ricin and its corresponding polypeptides (A & B chain) were purified from castor seed. The molecular weight of ricin subunits were 29,000 and 28,000 daltons. The amino acids in ricin determined were Asp45 The22 Ser40 Glu53 Cys4 Gly96 His5 Ile21 Leu33 Lys20 Met4 Phe13 Pro37 Tyr11 Ala45 Val23 Arg20 indicating that ricin contains approximately 516 amino acid residues. The amino acids of the two subunits of ricin A and B chains were Asp23 The12 Ser21 Glu29 Cys2 Gly48 His3 Ile12, Leu17 Lys10 Met2 Phe6 Pro17 Tyr7 Ala35 Val13 Arg13 while in B chain the amino acids were Asp22 The10 Ser19 Glu25 Cys2 Gly47 His1 Ile10, Leu15 Lys11 Met1 Phe7 Pro6 Tyr5 Ala32Val11 Arg10. The total helical content of ricin came around 53.6% which is a new observation.     

Improved treatment of cyclic β‐amino acids and successful prediction of β‐peptide secondary structure using a modified force field: AMBER*C

We added parameters to the AMBER* force field to model cyclic β‐amino acid derivatives more accurately within the commonly used MacroModel program. In an effort to generate an improved treatment of cyclohexane and cyclopentane conformational preferences, carbon–carbon torsional parameters were modified and incorporated into a force field we call AMBER*C. Simulation of trans‐2‐aminocyclohexanecarboxylic acid (trans‐ACHC) and trans‐2‐aminocyclopentanecarboxylic acid (trans‐ACPC) derivatives using AMBER*C produces more realistic energy differences between (pseudo)diaxial and (pseudo)diequatorial conformations than does simulation using AMBER*. AMBER*C molecular dynamics simulations more accurately reproduce the experimental hydrogen‐bonding tendencies of simple diamide derivatives of trans‐ACHC and trans‐ACPC than do simulations using the AMBER* force field. More importantly, this modified force field allows accurate qualitative prediction of the helical secondary structures adopted by β‐amino acid homo‐oligomers. © 2000 John Wiley & Sons, Inc. J Comput Chem 21: 763–773, 2000     

The Molecular Dynamic Simulation of Zolpidem Interaction with Gamma Aminobutyric Acid Type A Receptor

Our goal was to generate the extracellular domain of gamma‐aminobutyric acid type A receptor (GABA_A receptor) by comparative modeling and to study the interaction of zolpidem with the GABA_A receptor. The modeling strategy was verified to provide reasonable 3‐dimensional coordinates. These coordinates helped to combine all the subunits well. The benzodiazepine (BZ) binding site was located in a binding pocket between the α₁ and γ₂ subunits of the GABA_A receptor. Zolpidem was selected to dock into the binding site. In our study, the residues of the binding pocket were suggested to be αHis129, αTyr187, αGly228, αThr234, αTyr237, γMet96, γPhe116, γSer130, γGly143, and γMet169. By the calculation of the docking module, the conformation of zolpidem docking in the BZ binding site was investigated. A hydrogen bond was found at γArg136 when zolpidem's conformation was in rank 2 of the docking score. The contracted binding pocket showed residues at αHis129, αTyr187, αGly228, αTyr237, γPhe116, and γMet169. Zolpidem docking in a contracted binding pocket might generate a hydrogen bond in α His 129.     

Determination of 13C isotopic enrichment of valine and threonine by GC-C-IRMS after formation of the N(O,S)-ethoxycarbonyl ethyl ester derivatives of the amino acids

We describe a new method of assessing, in a single run, ¹³C isotopic enrichment of both Val and Thr by gas chromatography–combustion–isotope-ratio mass spectrometry (GC–C–IRMS). This method characterised by a rapid one-step derivatisation procedure performed at room temperature to form the N(O,S)-ethoxycarbonyl ethyl ester derivatives, and a polar column for GC. The suitability of this method for Val and Thr in in-vivo samples (mucosal hydrolysate) was demonstrated by studying protein metabolism with two tracers (¹³C-valine or ¹³C-threonine). The intra-day and inter-day repeatability were both assessed either with standards or with in-vivo samples at natural abundance and at low ¹³C isotopic enrichment. For inter-day repeatability CVs were between 0.8 and 1.5% at natural abundance and lower than 5.5% at 0.112 and 0.190 atom% enrichment for Val and Thr, respectively. Overall isotopic precision was studied for eleven standard amino acid derivatives (those of Val, Ala, Leu, Iso, Gly, Pro, Asp, Thr, Ser, Met, and Phe) and was assessed at 0.32‰. The ¹³C isotopic measurement was then extended to the other amino acids (Ala, Val, Leu, Iso, Gly, Pro, Thr, and Phe) at natural abundance for in-vivo samples. The isotopic precision was better than 0.002 atom% per amino acid (for n = 4 rats). This analytical method was finally applied to an animal study to measure Thr utilization in protein synthesis.     

Characterization of amino acids in silk sericin protein from Gonometa rufobrunnae by MEKC with phenyl isothiocyanate derivatization

An MEKC method was developed for the separation and characterization of phenyl-isothiocyanate (PITC)-labeled amino acids derived from Gonometa rufobrunnae silkworm after microdialysis sample cleanup. The influence of the buffer and SDS concentration on the resolution of the amino acids was investigated. A buffer system consisting of 25 mM phosphate, 10 mM borate buffer at pH 9.00, and 70 mM SDS showed the best results, with 13 PITC-amino acid derivatives being resolved out of 15 possible amino acids that were under study. Microdialysis sampling demonstrated its efficiency as a sample cleanup technique. Sericin protein from G. rufobrunnae was found to be characterized by at least 11 positively identified amino acids. These included His, Tyr, Ser, Ala, Phe, Lys, Gly, Arg, Cys, Glu, and Asp. Leu/Met and Val/Thr were coeluting pairs and hence could not be positively confirmed.     

Antihypertensive Effect of Angiotensin‐Converting Enzyme Inhibitory Peptide Obtained from Hen Ovotransferrin

Angiotensin I‐converting enzyme (ACE) inhibitory peptide was isolated from the hen ovotransferrin hydrolysate using chymotryptic hydrolysis by two steps of reverse‐phase high‐performance liquid chromatography. The amino sequence of this novel peptide was identified as Lys‐Val‐Arg‐Glu‐Gly‐Thr‐Thr‐Tyr that inhibited ACE activity in vitro in a concentration‐dependent manner with an effective concentration (IC₅₀) of 102.8 μM. Also, this inhibition was identified as noncompetitive using the Lineweaver‐Burk plot. Moreover, the antihypertensive action of this novel peptide was investigated by an intravenous injection into spontaneously hypertensive rats (SHR). A dose‐dependent reduction of systolic blood pressure by this peptide was observed after 40 min of treatment and it decreased the blood pressure markedly at the maximal dose (1 nmol/mL/kg). The maximal blood pressure lowering activity of this peptide was calculated as 163% of captopril (10 pmol/mL/kg) that was used as positive control. In conclusion, the obtained data suggests that Lys‐Val‐Arg‐Glu‐Gly‐Thr‐Thr‐Tyr has an ability to inhibit ACE activity and decrease the systolic blood pressure in hypertensive animals.     

Systematic study of the transfer of amino acids across the water/l,2-dichloroethane interface facilitated by dibenzo-18-crown-6

Failitated ion transfer reactions of 20 amino acids with dibenzo-18-crown-6 (DB18C6) at the water/1,2-dichloroethane (W/DCE) interfaces supported at the tips of micro- and nano-pipets were investigated systematically using cyclic voltammetry. It was found that there were only 10 amino acids, that is, Leu, Val, lle, Phe, Trp, Met, Ala, Gly, Cys, Gln (in brief), whose protonated forms as cations can give well-defined facilitated ion transfer voltammograms within the potential window, and the reaction pathway was proven to be consistent with the transfer by interfacial complexation/dissociation (TIC/TID) mechanisms. The association constants of DB 18C6 with different amino acids in the DCE (β°), and the kinetic parameters of reaction were evaluated based on the steady-state voltammetry of micro- or nano-pipets, respectively. The experimental results demonstrated that the selectivity of complexation of protonated amino acid by DB18C6 compared with that of alkali metal cations was low, which can be attributed to the vicinal effect arising from steric hindrance introduced by their side group and the steric bulk effect by lipophilic stabilization. Moreover, the association constants and the standard rate constants for different amino acids showed good correlations with their hydrophobicity (π), except Gly and Met, which inferred that the selectivity of such heterogeneous complex reaction for different amino acids with DB18C6, was not only affected by discrimination in binding these ions to the crown ether macro-cycle, but also significantly modified by the ion transfer Gibbs energy which was closely related to the structure of the transferred ions, protonated amino acids.     

Intrinsic energy landscapes of amino acid side-chains

Amino acid side-chain conformational properties influence the overall structural and dynamic properties of proteins and, therefore, their biological functions. In this study, quantum mechanical (QM) potential energy surfaces for the rotation of side-chain χ(1) and χ(2) torsions in dipeptides in the alphaR, beta, and alphaL backbone conformations were calculated. The QM energy surfaces provide a broad view of the intrinsic conformational properties of each amino acid side-chain. The extent to which intrinsic energetics dictates side-chain orientation was studied through comparisons of the QM energy surfaces with χ(1) and χ(2) free energy surfaces from probability distributions obtained from a survey of high resolution crystal structures. In general, the survey probability maxima are centered in minima of the QM surfaces as expected for sp(3) (or sp(2) for χ(2) of Asn, Phe, Trp, and Tyr) atom centers with strong variations between amino acids occurring in the energies of the minima indicating intrinsic differences in rotamer preferences. High correlations between the QM and survey data were found for hydrophobic side-chains except Met, suggesting minimal influence of the protein and solution environments on their conformational distributions. Conversely, low correlations for polar or charged side-chains indicate a dominant role of the environment in stabilizing conformations that are not intrinsically favored. Data also link the presence of off-rotamers in His and Trp to favorable interactions with the backbone. Results also suggest that the intrinsic energetics of the side-chains of Phe and Tyr may play important roles in protein folding and stability. Analyses on whether intrinsic side-chain energetics can influence backbone preference identified a strong correlation for residues in the alphaL backbone conformation. It is suggested that this correlation reflects the intrinsic instability of the alphaL backbone such that assumption of this backbone conformation is facilitated by intrinsically favorable side-chain conformations. Together our results offer a broad overview of the conformational properties of amino acid side-chains and the QM data may be used as target data for force field optimization.     

Efficient Access to Enantiopure γ4‐Amino Acids with Proteinogenic Side‐Chains and Structural Investigation of γ4‐Asn and γ4‐Ser in Hybrid Peptide Helices

Hybrid peptides composed of α‐ and β‐amino acids have recently emerged as new class of peptide foldamers. Comparatively, γ‐ and hybrid γ‐peptides composed of γ⁴‐amino acids are less studied than their β‐counterparts. However, recent investigations reveal that γ⁴‐amino acids have a higher propensity to fold into ordered helical structures. As amino acid side‐chain functional groups play a crucial role in the biological context, the objective of this study was to investigate efficient synthesis of γ⁴‐residues with functional proteinogenic side‐chains and their structural analysis in hybrid‐peptide sequences. Here, the efficient and enantiopure synthesis of various N‐ and C‐terminal free‐γ⁴‐residues, starting from the benzyl esters (COOBzl) of N‐Cbz‐protected (E)‐α,β‐unsaturated γ‐amino acids through multiple hydrogenolysis and double‐bond reduction in a single‐pot catalytic hydrogenation is reported. The crystal conformations of eight unprotected γ⁴‐amino acids (γ⁴‐Val, γ⁴‐Leu, γ⁴‐Ile, γ⁴‐Thr(OtBu), γ⁴‐Tyr, γ⁴‐Asp(OtBu), γ⁴‐Glu(OtBu), and γ‐Aib) reveals that these amino acids adopted a helix favoring gauche conformations along the central C^γ? C^β bond. To study the behavior of γ⁴‐residues with functional side chains in peptide sequences, two short hybrid γ‐peptides P1 (Ac‐Aib‐γ⁴‐Asn‐Aib‐γ⁴‐Leu‐Aib‐γ⁴‐Leu‐CONH₂) and P2 (Ac‐Aib‐γ⁴‐Ser‐Aib‐γ⁴‐Val‐Aib‐γ⁴‐Val‐CONH₂) were designed, synthesized on solid phase, and their 12‐helical conformation in single crystals were studied. Remarkably, the γ⁴‐Asn residue in P1 facilitates the tetrameric helical aggregations through interhelical H bonding between the side‐chain amide groups. Furthermore, the hydroxyl side‐chain of γ⁴‐Ser in P2 is involved in the interhelical H bonding with the backbone amide group. In addition, the analysis of 87 γ⁴‐residues in peptide single‐crystals reveal that the γ⁴‐residues in 12‐helices are more ordered as compared with the 10/12‐ and 12/14‐helices.     

Aliphatic poly(ester‐carbonate)s bearing amino groups and its RGD peptide grafting

This article deals with (1) synthesis of novel cyclic carbonate monomer (2‐oxo [1,3]dioxan‐5‐yl)carbamic acid benzyl ester (CAB) containing protected amino groups; (2) ring‐opening copolymerization of the cyclic monomer with L ‐lactide (LA) to provide novel degradable poly(ester‐carbonate)s with functional groups; (3) removal of the protective benzyloxycarbonyl (Cbz) groups by catalytic hydrogenation to afford the corresponding poly(ester‐co‐carbonate)s with free amino groups; (4) grafting of oligopeptide Gly‐Arg‐Gly‐Asp‐Ser‐Tyr (GRGDSY, abbreviated as RGD) onto the copolymer pendant amino groups in the presence of 1,1′‐carbonyldiimidazole (CDI). The structures of P(LA‐co‐CA/RGD) and its precursor were confirmed by ¹H NMR analysis. Cell experiments showed that P(LA‐co‐CA/RGD) had improved adhesion and proliferation behavior. Therefore, the novel RGD‐grafted block copolymer is promising for cell or tissue engineering applications. © Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 46: 7022–7032, 2008     

Solvation free energy of amino acids and side-chain analogues

The solvation free energies of amino acids and their side-chain analogues in water and cyclohexane are calculated by using Monte Carlo simulation. The molecular interactions are described by the OPLS-AA force field for the amino acids and the TIP4P model for water, and the free energies are determined by using the Bennett acceptance method. Results for the side-chain analogues in cyclohexane and in water are used to evaluate the performance of the force field for the van der Waals and the electrostatic interactions, respectively. Comparison of the calculated hydration free energies for the amino acid analogues and the full amino acids allows assessment of the additivity of the side chain contributions on the number of hydrating water molecules. The hydration free energies of neutral amino acids can be reasonably approximated by adding the contributions of their side chains to that of the hydration of glycine. However, significant nonadditivity in the free energy is found for the zwitterionic form of amino acids with polar side chains. In serine and threonine, intramolecular hydrogen bonds are formed between the polar side chains and backbone groups, leading to weaker solvation than for glycine. In contrast, such nonadditivity is not observed in tyrosine, in which the hydroxyl group is farther separated from, and therefore cannot form an intramolecular hydrogen bond with, the backbone. For histidine we find that a water molecule can form a bridge when the intramolecular hydrogen bond between the polar group and the backbone is broken.     

On the sodium and lithium ion affinities of some α-amino acids

The decomposition of 59 different cluster ions (generated by fast atom bombardment) consisting of two different amino acids and a sodium ion was analysed. The only fragment ions of significant abundance could be assigned to sodium ion-bound amino acids. Assuming that the most abundant ion in the fragment ion spectrum corresponds to the amino acid with the highest sodium ion affinity (SIA), the 20 common α-amino acids could be ordered with increasing sodium ion affinity as follows: Gly, Ala, Cys, Val, (Leu, Ile), Ser, Met, Thr, (Phe, Pro), Asp, Tyr, (Glu, Lys), Trp, Asn, Gln, His, Arg. Quantitative determinations were carried out by comparison of the lithium ion affinity (LIA) of Ala with that of dimethylformamide (DMF) in a fragment ion scan of the ion-bound dimer Ala—Li⁺—DMF. LIA(Ala) was calculated from LIA(Ala) = LIA(DMF) – (1/C)ln[I(AlaLi⁺)/I(DMF—Li⁺)], where the constant C was estimated from measurements of proton-bound amine–amino acid clusters. From fragment ion analysis of nine other Li⁺-bound α-amino acid dimers, the following lithium ion affinities were obtained: Gly 51.0, Ala 52.6, Sar 53.5, α-aminobutyric acid 53.7, glycine methyl ester 54.7 and Val 54.8. SIA(Ala) was estimated to be 75% of the lithium ion affinity and from fragment ion analysis of ten Na⁺-bound α-amino acid dimers the following sodium ion affinities were obtained: Gly 37.9, Ala 39.4, α-aminobutyric acid 40.3, Val 41.0, glycine methylster 41.0 and Sar 41.2.     

Aza-amino acid scanning of secondary structure suited for solid-phase peptide synthesis with fmoc chemistry and aza-amino acids with heteroatomic side chains

Aza-peptides, peptide analogues in which the alpha-carbon of one or more of the amino acid residues is replaced with a nitrogen atom, exhibit a propensity for adopting beta-turn conformations. A general Fmoc-protection protocol for the stepwise solid-phase synthesis of aza-peptides has now been developed based on the activation of N'-alkyl fluoren-9-ylmethyl carbazates with phosgene for coupling the aza-amino acid residues. This method has proven effective for introducing aza-amino acid residues with aliphatic (Ala, Leu, Val, and Gly) and aromatic (Phe, Tyr, and Trp) side chains. Acid promoted loss of aromatic side chains was noted with aza-Trp and aza-Tyr residues during peptide cleavage and suppressed by temperature control in the case of the latter. In addition, aza-peptides with heteroatomic side chain residues (Lys, Orn, Arg, and Asp) were conveniently synthesized using this protocol. Partial aza-amino acid scans were performed on three biologically active peptides: the potent tetrapeptide melanocortin receptor agonist, Ac-His-d-Phe-Arg-Trp-NH2; the growth hormone secretagogue hexapeptide, GHRP-6, His-d-Trp-Ala-Trp-d-Phe-Lys-NH2; and the human calcitonin gene-related peptide (hCGRP) antagonist, FVPTDVGPFAF-NH2. This practical procedure for aza-amino acid scanning using Fmoc-based solid-phase synthesis should find general utility for probing the existence and importance of beta-turn conformations in bioactive peptides.     

Validation of the GROMOS 54A7 Force Field Regarding Mixed α/β‐Peptide Molecules

A molecular‐dynamics (MD) simulation study of two heptapeptides containing α‐ and β‐amino acid residues is presented. According to NMR experiments, the two peptides differ in dominant fold when solvated in MeOH: peptide 3 adopts predominantly β‐hairpin‐like conformations, while peptide 8 adopts a 14/15‐helical fold. The MD simulations largely reproduce the experimental data. Application of NOE atom? atom distance restraining improves the agreement with experimental data, but reduces the conformational sampling. Peptide 3 shows a variety of conformations, while still agreeing with the NOE and ³J‐coupling data, whereas the conformational ensemble of peptide 8 is dominated by one helical conformation. The results confirm the suitability of the GROMOS 54A7 force field for simulation or structure refinement of mixed α/β‐peptides in MeOH.