Molecular dynamics simulation of folding of a short helical peptide with many charged residues

A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.     

A structural study of model peptides derived from HIV-1 integrase central domain

The HIV-1 integrase (IN) catalyzes the integration of viral DNA in the human genome. In vitro the enzyme displays an equilibrium of monomers, dimers, tetramers and larger oligomers. However, its functional oligomeric form in vivo is not known. We report a study of the auto-associative properties of three peptides denoted K156, E156 and E159. These derive from the alpha4 helix of the IN catalytic core. The alpha4 helix is an amphipatic helix exposed at the surface of the protein and could be involved in the oligomerization process through its hydrophobic face. The peptides were obtained from the replacement of several amino acid residues by more helicogenic ones in the alpha4 helix peptide. K156 carries the basic residues Lys156 and Lys159, which have been shown important for the binding of IN to viral DNA. In E156 and E159 they are replaced with the acidic residue Glu. A fourth peptide K(E)156 obtained from the replacement of hydrophobic residues with Glu in K156 in order to abolish the auto-associative properties is used as a negative control. The capacity shown by peptides for alpha-helical formation is demonstrated by circular dichroism (CD) analysis performed in aqueous solution and in aqueous trifluoroethanol (TFE) mixtures. Both electrospray ionization mass spectrometry (ESI-MS) and glutaraldehyde chemical cross-linking show that peptides adopt different solvent-dependent equilibriums of monomers, dimers, trimers and tetramers. Oligomerization of peptides in aqueous solution is related to their ability to form helical structures. Addition of a small amount of TFE (<10%) stimulates helix stabilization and the interhelical hydrophobic contacts. Higher amounts of TFE alter the hydrophobic contacts and disrupt the oligomeric species. In addition to hydrophobic interactions, the patterns indicate that the biologically important Lys156 and Lys159 residues also participate in helix association. K(E)156 despite its ability to adopt a helical structure is unable to associate into oligomers, demonstrating the importance of hydrophobic contacts for oligomerization. Thus, the designed peptides provide us information on the functional properties of the alpha4 IN that seems to hold a dual role in DNA recognition and protein oligomerization.     

应用ABEEMσπ／MM模型探讨小α-螺旋在显性水中折叠／展开的可逆过程

应用ABEEMσπ／MM模型进行分子动力学模拟,研究了显性水溶液中小α-螺旋（短肽Ala5）折叠／展开的可逆过程．动力学分析显示,300K下α-螺旋可以保存2ns的时间,该结果支持Margulis等人的结论．每个结构与α-螺旋结构骨架重原子的均方根偏差的时间轨迹指出,“300K下螺旋成核现象在0．1ns内快速发生”的结论是不恰当的．通过对300、400和500K温度下的研究,首次定量地给出各温度下螺旋保存的时间分别为2ns、1～1．5ns和0．8ns,并且增加温度并不改变折叠／展开的方式,只是改变折叠／展开的速率．本文对“转化态集合”结构的分析表明,从螺旋到卷曲的转换,主要通过螺旋端的氢键断裂发生（92％）、尤其是C端的氢键断裂发生（50％）．氢键的破坏和形成在0．1ns的时间内完成．     

Host-guest chemistry in the gas phase: complex formation with 18-crown-6 enhances helicity of alanine-based peptides

The gas-phase helix propensities of alanine-based polypeptides are studied with different locations of a Lys residue and host-guest interactions with 18-Crown-6 (18C6). A series of model peptides Ac-Ala(9-n)-LysH(+)-Ala(n) (n = 0, 1, 3, 5, 7, and 9) is examined alone and with 18C6 using traveling wave ion mobility mass spectrometry combined with molecular dynamics (MD) simulations. The gas-phase helices are observed from the peptides whose Lys residue is located close to the C-terminus so that the Lys exerts its capping effect on the C-terminal carbonyl groups. The peptides, which interact with 18C6 in the gas phase, show enhanced helical propensities. The enhanced helicity of the peptide in the complex is attributed by isolation of the Lys butylammonium group from the helix backbone and the interaction of methylene groups of 18C6, which possess localized positive partial charges, with C-terminal carbonyl groups serving as a cap to stabilize the helix.     

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.     

Fast folding dynamics of α-helical peptides – Effect of solvent additives and pH

The fast folding/unfolding dynamics of an alanine-based α-helical model peptide was investigated using nanosecond temperature jump experiments in D₂O. Temperature jumps were induced either by indirect heating, using visible laser pulses and a heat-transducing triphenylmethane dye, or by direct heating, using IR laser pulses. Although direct heating improves the quality, reliability and speed of the experiments, the heat-transducing dye was shown to not affect the helix-coil dynamics of alanine-based model peptides, thus validating the indirect heating method for this class of peptides. On the other hand, the presence of high concentrations of salt ions was found to alter the helix folding dynamics by screening of charged residues. We conclude that electrostatic interactions between residues have significant effects on the dynamics of helix folding. Similarly, the solvent pH was observed to affect the dynamics of helix folding, even in a range where no protonation/deprotonation is expected to occur.     

Backbone cleavages of [M - H](-) anions derived from caerin 1 peptides and some synthetic modifications. Molecular recognition initiating internal cyclisation of Glu23

The collision induced spectra of [M - H](-) anions from of caerin 1 peptides and some synthetic modifications show the usual alpha, beta and beta' backbone cleavages together with Ser (epsilon,gamma) and Glu (gamma) cleavages which break the peptide backbone in the vicinity of those residues. All of these cleavages require the peptide backbone to be flexible. There is also a backbone cleavage of a type not observed before. This cleavage involves nucleophilic attack of the carboxylate anion of the Glu23 side chain at the backbone CH of Ile 21. We propose that this cleavage requires the caerin peptide to be in an alpha helical conformation (the 3D structure that this peptide adopts in solution) in order that the interacting groups are held in close proximity.     

N,N'-linked oligoureas as foldamers: chain length requirements for helix formation in protic solvent investigated by circular dichroism, NMR spectroscopy, and molecular dynamics

N,N'-linked oligoureas with proteinogenic side chains are peptide backbone mimetics belonging to the gamma-peptide lineage. In pyridine, heptamer 4 adopts a stable helical fold reminiscent of the 2.6(14) helical structure proposed for gamma-peptide foldamers. In the present study, we have used a combination of CD and NMR spectroscopies to correlate far-UV chiroptical properties and conformational preferences of oligoureas as a function of chain length from tetramer to nonamer. Both the intensity of the CD spectra and NMR chemical shift differences between alphaCH2 diastereotopic protons experienced a marked increase for oligomers between four and seven residues. No major change in CD spectra occurred between seven and nine residues, thus suggesting that seven residues could be the minimum length required for stabilizing a dominant conformation. Unexpectedly, in-depth NMR conformational investigation of heptamer 4 in CD3OH revealed that the 2.5 helix probably coexists with partially (un)folded conformations and that Z-E urea isomerization occurs, to some degree, along the backbone. Removing unfavorable electrostatic interactions at the amino terminal end of 4 and adding one H-bond acceptor by acylation with alkyl isocyanate (4 --> 7) was found to reinforce the 2.5 helical population. The stability of the 2.5 helical fold in MeOH is further discussed in light of unrestrained molecular dynamics (MD) simulation. Taken together, these new data provide additional insight into the folding propensity of oligoureas in protic solvent and should be of practical value for the design of helical bioactive oligoureas.     

Molecular dynamics simulations of the helical antimicrobial peptide ovispirin-1 in a zwitterionic dodecylphosphocholine micelle: insights into host-cell toxicity

We have carried out a 40-ns all-atom molecular dynamics simulation of the helical antimicrobial peptide ovispirin-1 (OVIS) in a zwitterionic diphosphocholine (DPC) micelle. The DPC micelle serves as an economical and effective model for a cellular membrane owing to the presence of a choline headgroup, which resembles those of membrane phospholipids. OVIS, which was initially placed along a micelle diameter, diffuses out to the water-DPC interface, and the simulation stabilizes to an interface-bound steady state in 40 ns. The helical content of the peptide marginally increases in the process. The final conformation, orientation, and the structure of OVIS are in excellent agreement with the experimentally observed properties of the peptide in the presence of lipid bilayers composed of 75% zwitterionic lipids. The amphipathic peptide binds to the micelle with its hydrophobic face buried in the micellar core and the polar side chains protruding into the aqueous phase. There is overwhelming evidence that points to the significant and indispensable participation of hydrophobic residues in binding to the zwitterionic interface. The simulation starts with a conformation that is unbiased toward the final experimentally known binding state of the peptide. The ability of the model to reproduce experimental binding states despite this starting conformation is encouraging.     

Folding kinetics of a naturally occurring helical peptide: implication of the folding speed limit of helical proteins

The folding mechanism and dynamics of a helical protein may strongly depend on how quickly its constituent alpha-helices can fold independently. Thus, our understanding of the protein folding problem may be greatly enhanced by a systematic survey of the folding rates of individual alpha-helical segments derived from their parent proteins. As a first step, we have studied the relaxation kinetics of the central helix (L9:41-74) of the ribosomal protein L9 from the bacterium Bacillus stearothermophilus , in response to a temperature-jump ( T-jump) using infrared spectroscopy. L9:41-74 has been shown to exhibit unusually high helicity in aqueous solution due to a series of side chain-side chain interactions, most of which are electrostatic in nature, while still remaining monomeric over a wide concentration range. Thus, this peptide represents an excellent model system not only for examining how the folding rate of naturally occurring helices differs from that of the widely studied alanine-based peptides, but also for estimating the folding speed limit of (small) helical proteins. Our results show that the T-jump induced relaxation rate of L9:41-74 is significantly slower than that of alanine-based peptides. For example, at 11 degrees C its relaxation time constant is about 2 micros, roughly seven times slower than that of SPE(5), an alanine-rich peptide of similar chain length. In addition, our results show that the folding rate of a truncated version of L9:41-74 is even slower. Taken together, these results suggest that individual alpha-helical segments in proteins may fold on a time scale that is significantly slower than the folding time of alanine-based peptides. Furthermore, we argue that the relaxation rate of L9:41-74 measured between 8 and 45 degrees C provides a realistic estimate of the ultimate folding rate of (small) helical proteins over this temperature range.     

Solvent effect on the folding dynamics and structure of E6-associated protein characterized from ab initio protein folding simulations

Solvent effect on protein conformation and folding mechanism of E6-associated protein (E6ap) peptide are investigated using a recently developed charge update scheme termed as adaptive hydrogen bond-specific charge (AHBC). On the basis of the close agreement between the calculated helix contents from AHBC simulations and experimental results, we observed based on the presented simulations that the two ends of the peptide may simultaneously take part in the formation of the helical structure at the early stage of folding and finally merge to form a helix with lowest backbone RMSD of about 0.9 A? in 40% 2,2,2-trifluoroethanol solution. However, in pure water, the folding may start at the center of the peptide sequence instead of at the two opposite ends. The analysis of the free energy landscape indicates that the solvent may determine the folding clusters of E6ap, which subsequently leads to the different final folded structure. The current study demonstrates new insight to the role of solvent in the determination of protein structure and folding dynamics.     

Exploring Helical Folding in Oligomers of Cyclopentane-Based ϵ-Amino Acids: A Computational Study

The conformational preferences of oligopeptides of an ϵ-amino acid (2-((1R,3S)-3-(aminomethyl)cyclopentyl)acetic acid, Amc₅a) with a cyclopentane substituent in the C^β−C^γ−C^δ sequence of the backbone were investigated using DFT methods in chloroform and water. The most preferred conformation of Amc₅a oligomers (dimer to hexamer) was the H₁₆ helical structure both in chloroform and water. Four residues were found to be sufficient to induce a substantial H₁₆ helix population in solution. The Amc₅a hexamer adopted a stable left-handed (M)-2.3₁₆ helical conformation with a rise of 4.8 Å per turn. The hexamer of Ampa (an analogue of Amc₅a with replacing cyclopentane by pyrrolidine) adopted the right-handed mixed (P)-2.9_18/16 helical conformation in chloroform and the (M)-2.4₁₆ helical conformation in water. Therefore, hexamers of ϵ-amino acid residues exhibited different preferences of helical structures depending on the substituents in peptide backbone and the solvent polarity as well as the chain length.     

Negative ion electrospray mass spectra of the maculatin peptides from the tree frogs Litoria genimaculata and Litoria eucnemis

Collision-induced fragmentations of deprotonated maculatin 1 peptides provide significant sequencing information. When the peptide lacks those residues which can fragment through their alpha side chains (e.g. Thr, Ser, Glu and Gln in this study) the basic alpha and beta' backbone cleavages occur from the [Mbond;H](-) anion. When Thr, Ser, Glu and Gln are present, the ease of side-chain fragmentation of these residues is: Thr (loss of MeCHO) > Ser (CH(2)O) > Glu (H(2)O) > Gln (NH(3)). When one of more of these residues is (are) present, the alpha and beta' cleavages often occur from a fragment rather than the [Mbond;H](-) anion, e.g. for Thr, the [(Mbond;H)(-)bond;MeCHO](-) anion. These four residues also initiate gamma backbone cleavage reactions. The relative abundances of peaks resulting from gamma cleavage are Glu > Ser = Thr > Gln for maculatin 1 spectra. An unusual Gln19/Ile17 cyclisation/cleavage reaction occurs in maculatin spectra: the peptide [Mbond;H](-) anion must adopt a helical conformation in order for these two groups to interact. Analogous fragmentations have been reported previously in the negative ion spectra of the caerin 1 peptides.     

Molecular dynamics simulations of pro-apoptotic BH3 peptide helices in aqueous medium: relationship between helix stability and their binding affinities to the anti-apoptotic protein Bcl-X<Subscript>L</Subscript>

The B-cell lymphoma 2 (Bcl-2) family of proteins regulates the intrinsic pathway of apoptosis. Interactions between specific anti- and pro-apoptotic Bcl-2 proteins determine the fate of a cell. Anti-apoptotic Bcl-2 proteins have been shown to be over-expressed in certain cancers and they are attractive targets for developing anti-cancer drugs. Peptides from the BH3 region of pro-apoptotic proteins have been shown to interact with anti-apoptotic Bcl-2 proteins and induce biological activity similar to that observed in parent proteins. However, the specificity of BH3 peptides derived from different pro-apoptotic proteins differ for different anti-apoptotic Bcl-2 proteins. In this study, we have investigated the relationship between the stable helical nature of BH3 peptides and their affinities to Bcl-X_L, an anti-apoptotic Bcl-2 protein. We have carried out molecular dynamics simulations of six BH3 peptides derived from Bak, Bad and Bim pro-apoptotic proteins for a period of 50 ns each in aqueous medium. Due to the amphipathic nature of BH3 peptides, the hydrophobic residues on the hydrophobic face tend to cluster together in all BH3 peptides. While this process resulted in a complete loss of helical structure in 16-mer Bak and 16-mer Bad wild type peptides, stabilizing interactions in the hydrophilic face of the BH3 peptides and capping interactions helped to maintain partial helical character in 16-mer Bad mutant and 16-mer Bim peptides. The latter two 16-mer peptides exhibit higher affinity for Bcl-X_L. Similarly the longer BH3 peptides, 25-mer Bad and 33-mer Bim, also resulted in smaller and stable helical fragments and their helical conformation is stabilized by interactions between residues in the solvent-exposed hydrophilic half of the peptide. The stable nature of helical segment in a BH3 peptide can be directly correlated to its binding affinity and the helical region encompassed the highly conserved Leu residue. We propose that upon approaching the hydrophobic groove of anti-apoptotic proteins, a longer helix will be induced in high affinity BH3 peptides by extending the smaller stable helical segments around the conserved Leu residue in both N- and C-terminal regions. The results reported in this study will have implications in developing peptide-based inhibitors for anti-apoptotic Bcl-2 proteins.     

Structural Effects of L16Q,S20G,and L16Q‐S20G Mutations on hIAPP: A Comparative Molecular Dynamics Study

The conformation change picture of human islet amyloid polypeptide (hIAPP) is outlined using molecular dynamics simulation, and the structural influences of L16Q, S20G, and L16Q‐S20G mutations on the conformation of hIAPP are analyzed. Particularly, the conformational changes of the amyloidogenic‐related regions of residues 15–17 and 20–29 are emphasized. Our studies find that, for WT hIAPP, residues 15–17 always maintain a stable α‐helix structure, residues 20–25 structurally fluctuate between turn and 5‐helix, and residues 26–29 mainly adopt coil and bend structures. The hydrogen bonds between the polar groups of hIAPP, long‐rang van der Waals forces between the residues, and hydrophobic interactions between the residues of hIAPP are important driving forces to maintain the secondary structure of hIAPP. The replacement of leucine16 by glutamine stabilizes the helix structure of residues 15–17 and 20–23 of hIAPP monomer, and the structure of residues 24–29 fluctuates between helix and turn. The relatively stable helix structures of residues 15–17 and 20–29 are supposed to be beneficial for L16Q hIAPP to resist the aggregation as observed in the experiment. The substitution of serine20 by glycine drives residues 15–17 and 20–29 of hIAPP to transform from helix structure to β‐strands or coil structures with higher extension and flexibility, which may promote the aggregation of hIAPP as the experiments reported. These results are significant to understand the aggregation mechanism of hIAPP monomer into the dimer, trimer, oligomers and fibrils associated with the type 2 diabetes at the atomic level.     

Mechanisms of amphipathic helical peptide denaturation by guanidinium chloride and urea: a molecular dynamics simulation study

Urea and GdmCl are widely used to denature proteins at high concentrations. Here, we used MD simulations to study the denaturation mechanisms of helical peptide in different concentrations of GdmCl and urea. It was found that the helical structure of the peptide in water simulation is disappeared after 5 ns while the helicity of the peptide is disappeared after 70 ns in 2 M urea and 25 ns in 1 M GdmCl. Surprisingly, this result shows that the helical structure in low concentration of denaturants is remained more with respect to that solvated in water. The present work strongly suggests that urea interact more preferentially to non-polar and aromatic side chains in 2 M urea; therefore, hydrophobic residues are in more favorable environment in 2 M urea. Our results also reveal that the hydrogen bonds between urea and the backbone is the dominant mechanism by which the peptide is destabilized in high concentration of urea. In 1 M and 2 M GdmCl, GdmCl molecules tend to engage in transient stacking interactions with aromatics and hydrophobic planar side chains that lead to displacement of water from the hydration surface, providing more favorable environment for them. This shows that accumulation of GdmCl around hydrophobic surfaces in 1 M and 2 M GdmCl solutions prevents proper solvation of the peptide at the beginning. In high GdmCl concentrations, water solvate the peptide better than 1 M and 2 M GdmCl. Therefore, our results strongly suggest that hydrogen bonds between water and the peptide are important factors in the destabilization of peptide in GdmCl solutions.     

Expected and unexpected results from combined beta-hairpin design elements

A model beta-hairpin dodecapeptide [EFGWVpGKWTIK] was designed by including a favorable D-ProGly Type II' beta-turn sequence and a Trp-zip interaction, while also incorporating a beta-strand unfavorable glycine residue in the N-terminal strand. This peptide is highly folded and monomeric in aqueous solution as determined by combined analysis with circular dichroism and 1H NMR spectroscopy. A peptide representing the folded conformation of the model beta-hairpin [cyclic(EFGWVpGKWTIKpG)] and a linear peptide representing the unfolded conformation [EFGWVPGKWTIK] yield unexpected relative deviations between the CD and 1H NMR spectroscopic results that are attributed to variations in the packing interactions of the aromatic side chains. Mutational analysis of the model beta-hairpin indicates that the Trp-zip interaction favors folding and stability relative to an alternate hydrophobic cluster between Trp and Tyr residues [EFGYVpGKWTIK]. The significance of select diagonal interactions in the model beta-hairpin was tested by rearranging the cross-strand hydrophobic interactions to provide a folded peptide [EWFGIpGKTYWK] displaying evidence of an unusual backbone conformation at the hydrophobic cluster. This unusual conformation does not appear to be a result of the glycine residue in the beta-strand, as replacement with a serine results in a peptide [EWFSIpGKTYWK] with a similar and seemingly characteristic CD spectrum. However, an alternate arrangement of hydrophobic residues with a Trp-zip interaction in a similar position to the parent beta-hairpin [EGFWVpGKWITK] results in a folded beta-hairpin conformation. The differences between side chain packing of these peptides precludes meaningful thermodynamic analysis and illustrates the caution necessary when interpreting beta-hairpin folding thermodynamics that are driven, at least in part, by aromatic cross strand interactions.     

Molecular carpentry: piecing together helices and hairpins in designed peptides

The design of a peptide that contains two distinct elements of secondary structure, helix and beta-hairpin, is described. Two designed 17-residue peptides: Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-Gly-Gly-Leu-Phe-Val-D-Pro-Gly-Leu-Phe-Val-OMe (I) and Boc-Leu-Aib-Val-Ala-Leu-Aib-Val-Gly-Gly-Leu-Val-Val-D-Pro-Gly-Leu-Val-Val-OMe (II) have been conformationally characterized by NMR spectroscopy. Peptides I and II contain a seven-residue helical module at the N terminus and a eight-residue beta-hairpin module at the C terminus, which are connected by a conformationally flexible Gly-Gly segment. The choice of the secondary-structure modules is based upon prior crystallographic and spectroscopic analysis of the individual modules. Analysis of 500 MHz 1H NMR data, recorded as solutions in methanol, suggests that the observed pattern of chemical shifts, 3JHN CalphaH values, temperature coefficients of the NH chemical shifts, and backbone inter-residue nuclear Overhauser effects favor helical structures for residues 1-7 and beta-hairpin structures for residues 10-17. The spectroscopic data are compatible with termination of the helical segment by formation of a Schellman motif; this restricts Gly(8) to a left-handed alpha-helical conformation. Gly(9) is the only residue with multiple conformational possibilities in phi,psi space. Possible orientations of the two secondary-structure modules are considered. This study validates the use of stereochemically rigid peptide modules as prefabricated elements in the construction of synthetic protein mimics.     

Monitoring folding transitions of synthetic,branched-chain polypeptides by capillary zone electrophoresis

The coil/helix transition of a synthetic, branched-chain polymeric polypeptide (poly (Lys(Glu(1)-DL-Ala(3))EAK), 50-Lys residues long in the backbone, as a function of increasing molarities of methanol in solution, is here studied by both, circular dichroism (CD) and capillary zone electrophoresis. CD spectra showed that, at 75% v/v methanol, the transition from random coil to fully helical structure was obtained, in a pH 1.1 HCI solution in the presence of 20 mM NaCI. CZE studies, run in parallel, exhibited the classical unfolding to folding sigmoidal transition, with mid-point at 60% v/v methanol concentration, plateauing at ca. 80% v/v organic solvent. Surprisingly, though, such unfolding to folding transition was accompanied by an expansion, rather than a contraction, of the resulting ordered polypeptide. As the charge of the polypeptide (a pure polycation at a pH of 2.1 in CZE) was kept rigorously constant, a plot of the radius of the polymer along the sigmoidal transition clearly showed that the radius of gyration of the helical, structured polypeptide was in fact larger than that of the random coil. Such results were confirmed by molecular dynamics simulations, which indicated that the dimensions of such polypeptide, in alpha-helix configuration, were 8.5 nm (in length) and 3.2 nm (in diameter), whereas those of the corresponding random coil were 7.2 nm (in length) and 5.1 nm (length of shorter axis). It would thus appear that the randomized structure assumes the shape of a more compact object, roughly resembling a "rugby ball".     

Effect of thioxopeptide bonds on alpha-helix structure and stability

Thioxoamide (thioamide) bonds are nearly isosteric substitutions for amides but have altered hydrogen-bonding and photophysical properties. They are thus well-suited backbone modifications for physicochemical studies on peptides and proteins. The effect of thioxoamides on protein structure and stability has not been subject to detailed experimental investigations up to date. We used alanine-based model peptides to test the influence of single thioxoamide bonds on alpha-helix structure and stability. The results from circular dichroism measurements show that thioxoamides are strongly helix-destabilizing. The effect of an oxo-to-thioxoamide backbone substitution is of similar magnitude as an alanine-to-glycine substitution resulting in a helix destabilization of about 7 kJ/mol. NMR characterization of a helical peptide with a thioxopeptide bond near the N-terminus indicates that the thioxopeptide moiety is tolerated in helical structures. The thioxoamide group is engaged in an i, i+4 hydrogen bond, arguing against the formation of a 3(10)-helical structure as suggested for the N-termini of alpha-helices in general and for thioxopeptides in particular.