Experimental resolution of early steps in protein folding: testing molecular dynamics simulations

Time-resolved Tyr fluorescence spectroscopy coupled with a laser-induced temperature-jump (T-jump) was employed to follow the folding relaxation dynamics of the B-domain of Staphylococcal protein A. The single Tyr is located in helix 1 (H1) and is a sensitive probe of the structure of this helix and the overall helical bundle structure. The results from this study were compared to those from a complementary infrared T-jump study on this protein [Vu, D. M.; Myers, J. K.; Oas, T. G.; Dyer, R. B. Biochemistry 2004, 43, 3582]. Both methods detect a microsecond process that follows the cooperative relaxation of the helical bundle core. However, a fast process (10-7 s) that follows the relaxation of the individual helices was observed only with the infrared probe. Thus, fast formation of H1 is not observed, but rather H1 forms in the microsecond phase, concomitantly with the docking to (and stabilization by) the other two helices to form the helical bundle structure. This observation validates the results of several previous molecular dynamics simulations that predict H1 formation only in the final assembly of the helix bundle.     

Electron density redistribution accounts for half the cooperativity of alpha helix formation

The energy of alpha helix formation is well known to be highly cooperative, but the origin and relative importance of the contributions to helical cooperativity have been unclear. Here we separate the energy of helix formation into short range and long range components by using two series of helical dimers of variable length. In one dimer series two monomeric helices interact by forming hydrogen bonds, while in the other they are coupled only through long range, primarily electrostatic interactions. Using Density Functional Theory, we find that approximately half of the cooperativity of helix formation is due to electrostatic interactions between residues, while the other half is due to nonadditive many-body effects brought about by redistribution of electron density with helix length.     

Peptoid oligomers with alpha-chiral, aromatic side chains: sequence requirements for the formation of stable peptoid helices.

The achiral backbone of oligo-N-substituted glycines or "peptoids" lacks hydrogen-bond donors, effectively preventing formation of the regular, intrachain hydrogen bonds that stabilize peptide alpha-helical structures. Yet, when peptoids are N-substituted with alpha-chiral, aromatic side chains, oligomers with as few as five residues form stable, chiral, polyproline-like helices in either organic or aqueous solution. The adoption of chiral secondary structure in peptoid oligomers is primarily driven by the steric influence of these bulky, chiral side chains. Interestingly, peptoid helices of this class exhibit intense circular dichroism (CD) spectra that closely resemble those of peptide alpha-helices. Here, we have taken advantage of this distinctive spectroscopic signature to investigate sequence-related factors that favor and disfavor stable formation of peptoid helices of this class, through a comparison of more than 30 different heterooligomers with mixed chiral and achiral side chains. For this family of peptoids, we observe that a composition of at least 50% alpha-chiral, aromatic residues is necessary for the formation of stable helical structure in hexameric sequences. Moreover, both CD and 1H-13C HSQC NMR studies reveal that these short peptoid helices are stabilized by the placement of an alpha-chiral, aromatic residue on the carboxy terminus. Additional stabilization can be provided by the presence of an "aromatic face" on the helix, which can be patterned by positioning aromatic residues with three-fold periodicity in the sequence. Extending heterooligomer chain length beyond 12-15 residues minimizes the impact of the placement, but not the percentage, of alpha-chiral aromatic side chains on overall helical stability. In light of these new data, we discuss implications for the design of helical, biomimetic peptoids based on this structural motif.     

Facile synthesis of pH‐responsive glycopolypeptides with adjustable sugar density

Well‐defined pH‐responsive glycopolypeptides were prepared by polymer‐analogous aqueous amide coupling of d ‐glucosamine to poly(α,l ‐glutamic acid) (PGA) using the coupling agent 4‐(4,6‐dimethoxy‐1,3,5‐triazin‐2‐yl)‐4‐methylmorpholinium chloride (DMT‐MM) without any organic solvents, additives, or buffers. Degrees of substitution (DS) up to 80% can be achieved, and the DS is adjustable by the molar ratio of DMT‐MM to PGA repeating units. Successful glycosylation of both low MW and high MW PGA was confirmed by ¹H NMR and FTIR spectroscopy as well as by an enhanced solubility at low pH. CD spectroscopy revealed that glycosylated PGAs with a DS up to 0.63 are able to undergo a pH‐responsive and reversible helix‐coil transition. However, for polymers with higher DS no transition occurs. A comparison with PGAs functionalized with monoethanolamine showed that the low helicity at high DS is not a steric effect due to the bulky sugar moieties, but a solvation effect. Preliminary turbidimetric tests with the lectin Concanavalin A indicate a biological activity of these glycosylated polypeptides. © 2013 Wiley Periodicals, Inc. J. Polym. Sci., Part A: Polym. Chem. 2013 , 51, 3925–3931     

Structural and spectroscopic studies of peptoid oligomers with alpha-chiral aliphatic side chains

Substantial progress has been made in the synthesis and characterization of various oligomeric molecules capable of autonomous folding to well-defined, repetitive secondary structures. It is now possible to investigate sequence-structure relationships and the driving forces for folding in these systems. Here, we present detailed analysis by X-ray crystallography, NMR, and circular dichroism (CD) of the helical structures formed by N-substituted glycine (or "peptoid") oligomers with alpha-chiral, aliphatic side chains. The X-ray crystal structure of a N-(1-cyclohexylethyl)glycine pentamer, the first reported for any peptoid, shows a helix with cis-amide bonds, approximately 3 residues per turn, and a pitch of approximately 6.7 A. The backbone dihedral angles of this pentamer are similar to those of a polyproline type I peptide helix, in agreement with prior modeling predictions. This crystal structure likely represents the major solution conformers, since the CD spectra of analogous peptoid hexamers, dodecamers, and pentadecamers, composed entirely of either (S)-N-(1-cyclohexylethyl)glycine or (S)-N-(sec-butyl)glycine monomers, also have features similar to those of the polyproline type I helix. Furthermore, this crystal structure is similar to a solution NMR structure previously described for a peptoid pentamer comprised of chiral, aromatic side chains, which suggests that peptoids containing either aromatic or aliphatic alpha-chiral side chains adopt fundamentally similar helical structures in solution, despite distinct CD spectra. The elucidation of detailed structural information for peptoid helices with alpha-chiral aliphatic side chains will facilitate the mimicry of biomolecules, such as transmembrane protein domains, in a distinctly stable form.     

The dynamics of the B-A transition of natural DNA double helices

The dynamics of the B-A transition of DNA double helices with different GC contents and various chain lengths has been characterized by an electric field pulse technique. The field-induced B-A reaction is separated from orientation effects using the magic angle technique. Amplitudes reflecting the B-A reaction are observed selectively in the limited range of ethanol contents, where CD spectra demonstrate the B-A transition. The maximum amplitude appears at 1-2% higher ethanol content than the center of the B-A transition observed by CD because electric field pulses induce a relatively large perturbation from the A- toward the B-form. The relaxation curves measured after pulse termination reflect a spectrum of up to three relaxation processes. For DNA's with approximately 50% GC, the main part of the amplitude ( approximately 75%) is associated with time constants of approximately 2 micros, and another major component appears with time constants of 50-100 micros. These relaxation effects have been observed for DNA samples with 859, 2629, 7160, and 48501 bp. The time constant associated with the main amplitude increases with decreasing GC content from approximately 2 micros at 50% GC to approximately 3 mus at 41% GC and approximately 10 micros at 0% GC at the center of the B-A transition. Model calculations on the kinetics of cooperative linear Ising lattices predict the appearance of a distinct maximum of the mean relaxation time at the center of the transition. The absence of such maximum in our experimental data indicates a low cooperativity of the B-A transition with a nucleation parameter of approximately 0.1. The rate of the B-A transition is lower by approximately 3 orders of magnitude than that predicted by molecular dynamics simulations.     

Formation Sites and Microscopic Conformation Study on the Konjac Glucomannan Molecular Helices

<正>In this work,the formation sites,helical parameters and hydrogen bond positions of Konjac glucomannan molecular helices were investigated using molecular dynamic simulation method.To our interest,the KGM chain is mainly composed by local left and right helix structures. The formation sites of KGM chain might locate at the chain-segments containing acetyl groups,and the left helix is the favorable conformation of KGM.Temperature-dependent molecule conformation study indicates that the right helix is dominant when the temperature is lower than 343 K,above which,however,the left helix is dominating(right helix disappears).In addition,intramolecular hydrogen bonds in the left helix can be found at the-OH groups on C(2),C(4)and C(6)of mannose residues;comparably,the intramolecular hydrogen bonds in the right helix can be mainly observed at the-OH groups on C(4)and C(6)of the mannose residues and C(3)of the glucose residues.In conclusion,molecular dynamic simulation is an efficient method for the microscopic conformation study of glucomannan molecular helices.     

The effects of individual amino acids on the fast folding dynamics of alpha-helical peptides

Nanosecond temperature jump experiments coupled to time-resolved infrared spectroscopy were carried out on a series of alanine-based peptides containing different guest amino acids to study the effects of residues with different helix propensities on the helix-coil dynamics.     

Structural characterization of alpha-helices of implicitly solvated poly-alanine

The structural characteristics of alpha-helices in poly-alanine-based peptides have been investigated via molecular dynamics simulation with the goal of understanding the basic features of peptide simulations within the context of a model system, classical molecular dynamics with generalized Born (GB) solvation, and to shed insight into the formation and stabilization of alpha-helices in short peptides. The effects of peptide length, terminal charges, proline substitution, and temperature on the alpha-helical secondary structure have been studied. The simulations have shown that distinct secondary structure begins to develop in peptides with lengths approaching 10 residues while ambiguous structures occur in shorter peptides. The helical content of peptides with lengths > or =10 amino acids is observed to be nearly constant up to (Ala)(40). Interestingly, terminal charges and proline in the second position from the N-terminus alter the secondary structure locally with little effect on the overall alpha-helical content of the peptide. The free energy profile of helix formation was also investigated. A large increase in free energy accompanying the formation of helices with more than two consecutive hydrogen bonds in the (i, i + 4) pattern was observed while the free energy increases linearly with additional hydrogen bonds. Values for the change in enthalpy and entropy of helix nucleation and propagation are reported. Additionally the results obtained from the GB model are compared to explicit solvent simulations of two synthetic alanine-based peptides.     

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.     

Effects of As(III) binding on alpha-helical structure

As(III) displays a wide range of effects in cellular chemistry. Surprisingly, the structural consequences of arsenic binding to peptides and proteins are poorly understood. This study utilizes model alpha-helical peptides containing two cysteine (Cys) residues in various sequential arrangements and spatial locations to study the structural effects of arsenic binding. With i, and i + 1, i + 2, or i + 3 arrangements, CD spectroscopy shows that As(III) coordination causes helical destabilization when Cys residues are located at central or C-terminal regions of the helix. Interestingly, arsenic binding to i, i + 3 positions results in the elimination of helical structure and the formation of a relatively stable alternate fold. In contrast, helical stabilization is observed for peptides containing i, i + 4 Cys residues, with corresponding pseudo pairwise interaction energies (Delta G(pw) degrees) of -1.0 and -0.7 kcal/mol for C-terminal and central placements, respectively. Binding affinities and association rate constants show that As(III) binding is comparatively insensitive to the location of the Cys residues within these moderately stable helices. These data demonstrate that As(III) binding can be a significant modulator of helical secondary structure.     

Explanation of Ionic Sequences in Various Phenomena. VI. The Formation and Destruction of Helices in Bovine Plasma Albumin in Neutral and Acidic Solutions

Abstract

The structure of bovine plasma albumin (BPA) was examined by optical rotatory dispersion studies at both low (pH 1.5 and 2.0) and high (pH 9.0) pH values in various aqueous salt solutions. The resulting cationic sequences were compared to those observed by Pedersen for values of the sedimentation constant. At pH 9.0 the destruction of the “helix” produces an acidic sequence. The relative pH position of the “helix” transition, the fact that addition of salt increases the apparent helical content of BPA, and the observed acidic-type of sequence rule out the possibility of (1)ionic bonds between carboxylate and ?-amino groups, (2) hydrophobic bonds, or (3) hydrogen bonds between peptide linkages as major contributing forces in the formation of the helix. The stability of the “helix” in BPA between pH 3.0 and 9.0 must therefore be due to hydrogen bonds between carboxylate ions and hydroxyl groups such as those of serine, threonine, and tyrosine. Repulsive forces between the positively charged groups on BPA strengthen these bonds by preventing the expanded form of BPA from collapsing. At pH 2.0 two types of sequences were observed: The s⁰ _{20, w}, [α]_D and a₀ values gave an acidic-type cationic sequence. The b₀ (helix content), λ_c and [λ]₂₃₃ values gave essentially a nonpolar sequence. The nonpolar or hydrophobic salting-out sequences show that the formation of hydrophobic bonds at pH 2.0 hinders the formation of the helix or folded structure. The acidic sequences show that hydrogen bonds between carboxylic acid groups stabilize both the apparent helix or helices and the intermolecular aggregation of the BPA molecules. From a comparison of the S⁰ ₂₀,w values and the helical content of BPA at pH 9.0 it is also concluded that the formation of these apparent helices or folded structures expands or stiffens the BPA molecule.     

Peptoid oligomers with alpha-chiral, aromatic side chains: effects of chain length on secondary structure.

Oligomeric N-substituted glycines or "peptoids" with alpha-chiral, aromatic side chains can adopt stable helices in organic or aqueous solution, despite their lack of backbone chirality and their inability to form intrachain hydrogen bonds. Helical ordering appears to be stabilized by avoidance of steric clash as well as by electrostatic repulsion between backbone carbonyls and pi clouds of aromatic rings in the side chains. Interestingly, these peptoid helices exhibit intense circular dichroism (CD) spectra that closely resemble those of peptide alpha-helices. Here, we have utilized CD to systematically study the effects of oligomer length, concentration, and temperature on the chiral secondary structure of organosoluble peptoid homooligomers ranging from 3 to 20 (R)-N-(1-phenylethyl)glycine (Nrpe) monomers in length. We find that a striking evolution in CD spectral features occurs for Nrpe oligomers between 4 and 12 residues in length, which we attribute to a chain length-dependent population of alternate structured conformers having cis versus trans amide bonds. No significant changes are observed in CD spectra of oligomers between 13 and 20 monomers in length, suggesting a minimal chain length of about 13 residues for the formation of stable poly(Nrpe) helices. Moreover, no dependence of circular dichroism on concentration is observed for an Nrpe hexamer, providing evidence that these helices remain monomeric in solution. In light of these new data, we discuss chain length-related factors that stabilize organosoluble peptoid helices of this class, which are important for the design of helical, biomimetic peptoids sharing this structural motif.     

Conformational manifold of alpha-aminoisobutyric acid (Aib) containing alanine-based tripeptides in aqueous solution explored by vibrational spectroscopy, electronic circular dichroism spectroscopy, and molecular dynamics simulations

Replacement of the alpha-proton of an alanine residue to generate alpha-aminoisobutyric acid (Aib) in alanine-based oligopeptides favors the formation of a 3(10) helix when the length of the oligopeptide is about four to six residues. This research was aimed at experimentally identifying the structural impact of an individual Aib residue in an alanine context of short peptides in water and Aib's influence on the conformation of nearest-neighbor residues. The amide I band profile of the IR, isotropic and anisotropic Raman, and vibrational circular dichroism (VCD) spectra of Ac-Ala-Ala-Aib-OMe, Ac-Ala-Aib-Ala-OMe, and Ac-Aib-Ala-Ala-OMe were measured and analyzed in terms of different structural models by utilizing an algorithm that exploits the excitonic coupling between amide I' modes. The conformational search was guided by the respective 1H NMR and electronic circular dichroism spectra of the respective peptides, which were also recorded. From these analyses, all peptides adopted multiple conformations. Aib predominantly sampled the right-handed and left-handed 3(10)-helix region and to a minor extent the bridge region between the polyproline (PPII) and the helical regions of the Ramachandran plot. Generally, alanine showed the anticipated PPII propensity, but its conformational equilibrium was shifted towards helical conformations in Ac-Aib-Ala-Ala-OMe, indicating that Aib can induce helical conformations of neighboring residues positioned towards the C-terminal direction of the peptide. An energy landscape exploration by molecular dynamics simulations corroborated the results of the spectroscopic studies. They also revealed the dynamics and pathways of potential conformational transitions of the corresponding Aib residues.     

Mapping the orientation of helices in micelle-bound peptides by paramagnetic relaxation waves

Many antimicrobial peptides form alpha-helices when bound to a membrane. In addition, around 80% of residues in membrane-bound proteins are found in alpha-helical regions. The orientation and location of such helical peptides and proteins in the membrane are key factors determining their function and activity. Here we present a new solution state NMR method for obtaining the orientation of helical peptides in a membrane-mimetic environment (micelle-bound) without any chemical perturbation of the peptide-micelle system. By monitoring proton longitudinal relaxation rates upon addition of the freely water-soluble and inert paramagnetic probe Gd(DTPA-BMA) to an alpha-helical peptide, a wavelike pattern with a periodicity of 3.6 residues per turn is observed. The tilt and azimuth (rotation) angle of the helix determine the shape of this paramagnetic relaxation wave and can be obtained by least-square fitting of measured relaxation enhancements. Results are presented for the 15-residue antimicrobial peptide CM15 which forms an amphipathic helix almost parallel to the surface of the micelle. Thus, a few fast experiments enable the identification of helical regions and determination of the helix orientation within the micelle without the need for covalent modification, isotopic labeling, or sophisticated equipment. This approach opens a path toward the topology determination of alpha-helical membrane-proteins without the need for a complete NOE-based structure determination.     

Polarographic study on poly alpha-L-glutamic acid-Cd2+,Co2+ complexes in the helix-coil pH region

Interactions between poly alpha- L-glutamic acid (PGA) and metal ions Cd(2+), Co(2+) were studied by direct current polarography. The diffusion currents of these ions decreased sharply in the presence of PGA in the pH region from 5.0 through neutral. A corresponding increase in the helix content of the PGA-metal ion complex was revealed by CD measurements on the same solutions. Helix contents determined by polarography were in good agreement with those by CD in the neutral pH region. On the contrary, the decrease of current in lower acidic pH regions was independent of helix formation and suggested that metal ions coordinate to sporadically-dissociated carboxylate groups to cause aggregation of the intra and/or inter polymer chains. The diffusion current of the ions, therefore, is a parameter sensitive to the conformational changes of PGA from acidic through neutral pH region.     

Infrared study of the stability and folding kinetics of a 15-residue beta-hairpin

The thermal stability and folding kinetics of a 15-residue beta-hairpin (SESYINPDGTWTVTE) have been studied by using infrared (IR) spectroscopy coupled with laser-induced temperature-jump (T-jump) technique for rapid folding-unfolding initiation. An alternative method based on analyzing IR difference spectra was also introduced to obtain thermodynamic properties of beta-sheets, which complements the commonly used circular dichroism (CD) and fluorescence techniques. Equilibrium IR measurements indicate that the thermal unfolding of this beta-hairpin is fairly broad. However, it can be described by a two-state transition with a thermal melting temperature of approximately 29 degrees C. Time-resolved IR measurements following a T-jump, probed at 1634 cm(-1), indicate that the folding of this beta-hairpin follows first-order kinetics and is amazingly fast. At 300 K, the folding time is approximately 0.8 micros, which is only 2-3 times slower than that of alpha-helix formation. Additionally, the energetic barrier for folding is small (approximately 2 kcal mol(-1)). These results, in conjunction with results from other studies, support a view that the details of native contacts play a dominant role in the kinetics of beta-hairpin folding.     

Length-Dependent Collective Vibrational Dynamics in Alpha-Helices

Functions of protein molecules in nature are closely associated with their well-defined three-dimensional structures and dynamics in body fluid. So far, many efforts have been made to reveal the relation of protein structure, dynamics, and function, but the structural origin of protein dynamics, especially for secondary structures as building blocks of conformation transition, is still ambiguous. Here we theoretically uncover the collective vibrations of elastic poly-alanine α-helices and find vibration patterns that are distinctively different over residue numbers ranging from 20 to 80. Contrary to the decreasing vibration magnitude from ends to the middle region for short helices, the vibration magnitude for long helices takes the minimum at approximately 1/5 of helix length from ends but reaches a peak at the center. Further analysis indicates that major vibrational modes of helical structures strongly depend on their lengths, where the twist mode dominates in the vibrations of short helices while the bend mode dominates the long ones analogous to an elastic Euler beam. The helix-coil transition pathway is also affected by the alternation of the first-order mode in helices with different lengths. The dynamic properties of the helical polypeptides are promising to be harnessed for de novo design of protein-based materials and artificial biomolecules in clinical treatments.     

Two-dimensional infrared spectra of the 13C=18O isotopomers of alanine residues in an alpha-helix

The parameters needed to describe the two-dimensional infrared (2D IR) spectra of the isotopically labeled alpha-helix are presented. The 2D IR spectra in the amide-I' spectral region of a series of singly 13C=18O-labeled 25-residue alpha-helices were measured by three-pulse heterodyned spectral interferometry. The dependence of the spectra on the population time was measured. Individual isotopomer levels (residues 11-14) were clearly identified in 2D IR, downshifted by approximately 61 cm(-1) from the main helical band. By analyzing the line shapes of the 13C=18O diagonal peaks that appeared at approximately 1571.3 +/- 0.8 cm(-1) for all four labeled samples, we observed wider structural distributions for residues 14 and 11 than those for 12 and 13. A small fast component in the correlation function was used to estimate the dynamics of these distributions. In all cases, the v = 1 --> 2 transition showed a more Lorentzian-like line shape and also decayed faster than the v = 0 --> 1 transition, indicating that the population relaxation time of the v = 2 state was significantly faster than the v = 1 state. The amide transitions with naturally abundant 13C=16O appeared at approximately 1594 cm(-1), forming very weak and blurred cross-peaks with 13C=18O isotopomer modes. The effects of spectral interferences on the coherence time dependence of the detection frequency spectrum were also investigated. The methods of first moments and Wigner analysis were developed to circumvent the interference effects on the weak isotopomer transitions. The structural origin of the distributions for individual isotopomers was proposed to be an effect of nearby lysine residues on the intrahelical hydrogen-bond network.