Folding of peptide fragments comprising the complete sequence of proteins. Models for initiation of protein folding. I. Myohemerythrin.

In an attempt to delineate potential folding initiation sites for different protein structural motifs, we have synthesized series of peptides that span the entire length of the polypeptide chain of two proteins, and examined their conformational preferences in aqueous solution using proton nuclear magnetic resonance and circular dichroism spectroscopy. We describe here the behavior of peptides derived from a simple four-helix bundle protein, myohemerythrin. The peptides correspond to the sequences of the four long helices (the A, B, C and D helices), the N- and C-terminal loops and the connecting sequences between the helices. The peptides corresponding to the helices of the folded protein all exhibit preferences for helix-like conformations in solution. The conformational ensembles of the A- and D-helix peptides contain ordered helical forms, as shown by extensive series of medium-range nuclear Overhauser effect connectivities, while the B- and C-helix peptides exhibit conformational preferences for nascent helix. All four peptides adopt ordered helical conformations in mixtures of trifluoroethanol and water. The terminal and interconnecting loop peptides also appear to contain appreciable populations of conformers with backbone phi and psi angles in the alpha-region and include highly populated hydrophobic cluster and/or turn conformations in some cases. Trifluoroethanol is unable to drive these peptides towards helical conformations. Overall, the peptide fragments of myohemerythrin have a marked preference towards secondary structure formation in aqueous solution. In contrast, peptide fragments derived from the beta-sandwich protein plastocyanin are relatively devoid of secondary structure in aqueous solution (see accompanying paper). These results suggest that the two different protein structural motifs may require different propensities for formation of local elements of secondary structure to initiate folding, and that there is a prepartitioning of conformational space determined by the local amino acid sequence that is different for the helical and beta-sandwich structural motifs.     

Crystallographic analysis of synechocystis cyanoglobin reveals the structural changes accompanying ligand binding in a hexacoordinate hemoglobin

The crystal structures of cyanide and azide-bound forms of the truncated hemoglobin from Synechocystis are presented at 1.8 angstroms resolution. A comparison with the structure of the endogenously liganded protein reveals a conformational shift unprecedented in hemoglobins, and provides the first picture of a hexacoordinate hemoglobin in both the bis-histidyl and the exogenously coordinated states. The structural changes between the different conformations are confined to two regions of the protein; the B helix, and the E helix, including the EF loop. A molecular "hinge" controlling movement of the E helix is observed in the EF loop, which is composed of three principal structural elements: Arg64, the heme-d-propionate, and a three-residue extension of the F helix. Additional features of the structural transition between the two protein conformations are discussed as they relate to the complex ligand-binding behavior observed in hexacoordinate hemoglobins, and the potential physiological function of this class of proteins.     

Helix-forming tendencies of nonpolar amino acids predicted by Monte Carlo simulated annealing

Monte Carlo simulated annealing is applied to the study of the α-helix-forming tendencies of seven nonpolar amino acids, Ala, Leu, Met, Phe, Ile, Val, and Gly. Homooligomers of 10 amino acids are used and the helix tendency is calculated by folding α-helicies from completely random initial conformations. The results of the simulation imply that Met, Ala, and Leu are helix formers and that Val, Ile, and Gly are helix breakers, while Phe comes in between the two groups. The differences between helix formers and breakers turned out to be large in agreement with the recent experiments with short peptides. It is argued from the energy distributions of the obtained conformations that the helix tendency is small for the helix breakers because of steric hindrance of side chains. Homoglycine is shown to favor a random coil conformation. The β-strand tendencies of the same homooligomers are also considered, and they are shown to agree with the frequencies of amino acids in β-sheet from the protein data base. © 1994 Wiley-Liss, Inc.     

Hydration of the left spiral of the poly-L-proline type. Study by the Monte Carlo method]

The paper exhibits results of hydration shell Monte Carlo calculations in poly-L-proline II and extended helix conformation and in alpha-helical and beta-structural conformations for comparison. It was found that left-handed helix of poly-L-proline II type as well as epsilon-helix are characterized by very favorable hydration. Therefore this conformation has preference as compared to other standard conformations of the main polypeptide chain. This determined inevitability of cold denaturation of protein.     

Structural basis of ligand binding and release in insect pheromone-binding proteins: NMR structure of Antheraea polyphemus PBP1 at pH 4.5

The NMR structure of the Antheraea polyphemus pheromone-binding protein 1 at pH 4.5, ApolPBP1A, was determined at 20 degrees C. The structure consists of six alpha-helices, which are arranged in a globular fold that encapsulates a central helix alpha7 formed by the C-terminal polypeptide segment 131-142. The 3D arrangement of these helices is anchored by the three disulfide bonds 19-54, 50-108 and 97-117, which were identified by NMR. Superposition of the ApolPBP1A structure with the structure of the homologous pheromone-binding protein of Bombyx mori at pH 4.5, BmorPBPA, yielded an rmsd of 1.7 A calculated for the backbone heavy-atoms N, Calpha and C' of residues 10-142. In contrast, the present ApolPBP1A structure is different from a recently proposed molecular model for a low-pH form of ApolPBP1 that does not contain the central helix alpha7. ApolPBP1 exhibits a pH-dependent transition between two different globular conformations in slow exchange on the NMR chemical shift timescale similar to BmorPBP, suggesting that the two proteins use the same mechanism of ligand binding and ejection. The extensive sequence homology observed for pheromone-binding proteins from moth species further implies that the previously proposed mechanism of ligand ejection involving the insertion of a C-terminal helix into the pheromone-binding site is a general feature of pheromone signaling in moths.     

Alpha-synuclein populates both elongated and broken helix states on small unilamellar vesicles

The misfolding of the protein α-synuclein (αS) has been implicated in the molecular chain of events leading to Parkinson disease. Physiologically, αS undergoes a transition from a random coil to helical conformation upon encountering synaptic vesicle membranes. On analogous small unilamellar vesicles (SUVs), the conformation of αS is dominated by a single elongated αS helix. However, alternative broken helix states have been postulated, mandating experimental clarification. Here, the upper limit for the free energy difference between elongated and broken helix conformations on SUVs resembling synaptic vesicles was determined to be 1.2 ± 0.4 kcal/mol, which amounts to a population ratio of 7.6:1 between both states (12% broken helices). In response to helix breaks at different positions, αS rearranged in an opportunistic manner, thereby minimizing helix abrogations to as little as one to two turns. Enthalpy and entropy measurements of gel state SUV-αS interactions indicated that broken helix states retain the ability to relieve membrane-packing stress. Thus, broken helix states are a distinct physiological feature of the vesicle-bound αS state, making it a "checkered" protein of multiple parallel conformations. A continuous interconversion between structural states may contribute to pathological αS misfolding.     

Folding propensities of peptide fragments of myoglobin.

Myoglobin has been studied extensively as a paradigm for protein folding. As part of an ongoing study of potential folding initiation sites in myoglobin, we have synthetized a series of peptides covering the entire sequence of sperm whale myoglobin. We report here on the conformation preferences of a series of peptides that cover the region from the A helix to the FG turn. Structural propensities were determined using circular dichroism and nuclear magnetic resonance spectroscopy in aqueous solution, trifluoroethanol, and methanol. Peptides corresponding to helical regions in the native protein, namely the B, C, D, and E helices, populate the alpha region of (phi, psi) space in water solution but show no measurable helix formation except in the presence of trifluoroethanol. The F-helix sequence has a much lower propensity to populate helical conformations even in TFE. Despite several attempts, we were not successful in synthesizing a peptide corresponding to the A-helix region that was soluble in water. A peptide termed the AB domain was constructed spanning the A- and B-helix sequences. The AB domain is not soluble in water, but shows extensive helix formation throughout the peptide when dissolved in methanol, with a break in the helix at a site close to the A-B helix junction in the intact folded myoglobin protein. With the exception of one local preference for a turn conformation stabilized by hydrophobic interactions, the peptides corresponding to turns in the folded protein do not measurably populate beta-turn conformations in water, and the addition of trifluoroethanol does not enhance the formation of either helical or turn structure. In contrast to the series of peptides described here, either studies of peptides from the GH region of myoglobin show a marked tendency to populate helical structures (H), nascent helical structures (G), or turn conformations (GH peptide) in water solution. This region, together with the A-helix and part of the B-helix, has been shown to participate in an early folding intermediate. The complete analysis of conformational properties of isolated myoglobin peptides supports the hypothesis that spontaneous secondary structure formation in local regions of the polypeptide may play an important role in the initiation of protein folding.     

Not too loose, not too tight--just right. Biphasic control of the Tsr HAMP domain

HAMP domains communicate between input and output signalling modules in a wide variety of bacterial sensor proteins. In the Tsr chemoreceptor, they convert a signal initiated by binding of serine to the periplasmic domain of the protein into regulation of receptor control of the CheA kinase, and ultimately of the direction of flagellar rotation. In this issue, Zhou et al. report an extensive mutational analysis of the Tsr HAMP domain that shows that it can assume a number of different signalling states, which presumably correspond to a variety of different conformations. The two conformational extremes of a tightly packed and a loosely packed HAMP four‐helix bundle support only low levels of CheA activity. Thus, Tsr HAMP does not function as a simple on‐off, two‐state device but rather as a dynamic structure with biphasic control. The normal physiological operating range of Tsr is proposed to be at intermediate degrees of packing of the HAMP four‐helix bundle, but HAMP domains in other proteins could occupy different portions of the conformational spectrum.     

Prediction,refinement, and persistency of transmembrane helix dimers in lipid bilayers using implicit and explicit solvent/lipid representations: Microsecond molecular dynamics simulations of ErbB1/B2 and EphA1

All‐atom simulations are carried out on ErbB1/B2 and EphA1 transmembrane helix dimers in lipid bilayers starting from their solution/DMPC bicelle NMR structures. Over the course of microsecond trajectories, the structures remain in close proximity to the initial configuration and satisfy the majority of experimental tertiary contact restraints. These results further validate CHARMM protein/lipid force fields and simulation protocols on Anton. Separately, dimer conformations are generated using replica exchange in conjunction with an implicit solvent and lipid representation. The implicit model requires further improvement, and this study investigates whether lengthy all‐atom molecular dynamics simulations can alleviate the shortcomings of the initial conditions. The simulations correct many of the deficiencies. For example, excessive helix twisting is eliminated over a period of hundreds of nanoseconds. The helix tilt, crossing angles, and dimer contacts approximate those of the NMR‐derived structure, although the detailed contact surface remains off‐set for one of two helices in both systems. Hence, even microsecond simulations are not long enough for extensive helix rotations. The alternate structures can be rationalized with reference to interaction motifs and may represent still sought after receptor states that are important in ErbB1/B2 and EphA1 signaling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Effect of Pseudorepeat Rearrangement on α-Synuclein Misfolding, Vesicle Binding, and Micelle Binding

The pathological and physiological hallmarks of the protein α-synuclein (aS) are its misfolding into cytotoxic aggregates and its binding to synaptic vesicles, respectively. Both events are mediated by seven 11-residue amphiphilic pseudorepeats and, most generally, involve a transition from intrinsically unstructured conformations to structured conformations. Based on aS interactions with aggregation-inhibiting small molecules, an aS variant termed shuffled α-synuclein (SaS), wherein the first six pseudorepeats had been rearranged, was introduced. Here, the effects of this rearrangement on misfolding, vesicle binding, and micelle binding are examined in reference to aS and β-synuclein to study the sequence characteristics underlying these processes. Fibrillization correlates with the distinct clustering of residues with high β-sheet propensities, while vesicle affinities depend on the mode of pseudorepeat interchange and loss. In the presence of micelles, the pseudorepeat region of SaS adopts an essentially continuous helix, whereas aS and β-synuclein encounter a distinct helix break, indicating that a more homogeneous distribution of surfactant affinities in SaS prevented the formation of an extensive helix break in the micelle-bound state. By demonstrating the importance of the distribution of β-sheet propensities and by revealing inhomogeneous aS surfactant affinities, the present study provides novel insights into two central themes of synuclein biology.     

Effect of hydrophobic surfactant proteins SP-B and SP-C on phospholipid monolayers. Protein structure studied using 2D IR and beta correlation analysis

We have applied two-dimensional infrared (2D IR) and betanu correlation spectroscopy to in-situ IR spectroscopy of pulmonary surfactant proteins SP-B and SP-C in lipid-protein monolayers at the air-water interface. For both SP-B and SP-C, a statistical windowed autocorrelation method identified two separate surface pressure regions that contained maximum amide I intensity changes: 4-25 mN/m and 25-40 mN/m. For SP-C, 2D IR and betanu correlation analyses of these regions indicated that SP-C adopts a variety of secondary structure conformations, including alpha-helix, beta-sheet, and an intermolecular aggregation of extended beta-sheet structure. The main alpha-helix band split into two peaks at high surface pressures, indicative of two different helix conformations. At low surface pressures, all conformations of the SP-C molecule reacted identically to increasing surface pressure and reoriented in phase with each other. Above 25 mN/m, however, the increasing surface pressure selectively affected the coexisting protein conformations, leading to an independent reorientation of the protein conformations. The asynchronous 2D IR spectrum of SP-B showed the presence of two alpha-helix components, consistent with two separate populations of alpha-helix in SP-B-a hydrophobic fraction associated with the lipid chains and a hydrophilic fraction parallel to the membrane surface. The distribution of correlation intensity between the two alpha-helix cross peaks indicated that the more hydrophobic helix fraction predominates at low surface pressures whereas the more hydrophilic helix fraction predominates at high surface pressures. The different SP-B secondary structures reacted identically to increasing surface pressure, leading to a reorientation of all SP-B subunits in phase with one another.     

Helix‐sheet packing in proteins

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Although much is known about the packing geometry observed between α‐helices and between β‐sheets, there has been little progress on characterizing helix–sheet interactions. We present an analysis of the conformation of αβ₂ motifs in proteins, corresponding to all occurrences of helices in contact with two strands that are hydrogen bonded. The geometry of the αβ₂ motif is characterized by the azimuthal angle θ between the helix axis and an average vector representing the two strands, the elevation angle ψ between the helix axis and the plane containing the two strands, and the distance D between the helix and the strands. We observe that the helix tends to align to the two strands, with a preference for an antiparallel orientation if the two strands are parallel; this preference is diminished for other topologies of the β‐sheet. Side‐chain packing at the interface between the helix and the strands is mostly hydrophobic, with a preference for aliphatic amino acids in the strand and aromatic amino acids in the helix. From the knowledge of the geometry and amino acid propensities of αβ₂ motifs in proteins, we have derived different statistical potentials that are shown to be efficient in picking native‐like conformations among a set of non‐native conformations in well‐known decoy datasets. The information on the geometry of αβ₂ motifs as well as the related statistical potentials have applications in the field of protein structure prediction. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Thermodynamics and mechanism of alpha helix initiation in alanine and valine peptides

We used molecular dynamics simulations to study the folding/unfolding of one of turn of an alpha helix in Ac-(Ala)3-NHMe and Ac-(Val)3-NHMe. Using specialized sampling techniques, we computed free energy surfaces as functions of a conformational coordinate that corresponds to alpha helices at small values and to extended conformations at large values. Analysis of the peptide conformations populated during the simulations showed that alpha helices, reverse turns, and extended conformations correspond to minima on the free energy surfaces of both peptides. The free energy difference between alpha helix and extended conformations, determined from the equilibrium constants for helix unfolding, is approximately -1 kcal/mol for Ac-(Ala)3-NHMe and -5 kcal/mol for Ac-(Val)3-NHMe. The mechanism observed in our simulations, which includes reverse turns as important intermediates along the helix folding/unfolding pathway, is consistent with a mechanism proposed previously. Our results predict that both peptides (but especially the Ala peptide) have a much larger equilibrium constant for helix initiation than is predicted by the helix-coil transition theory with the host-guest parameters. We also predict a much greater difference in the equilibrium constants than the theory predicts. Insofar as helix initiation is concerned, our results suggest that the large difference between the helical propensities of Ala and Val cannot be explained by simple concepts such as side-chain rotamer restriction or unfavorable steric interactions. Rather, the origin of the difference appears to be quite complicated because it involves subtle differences in the solvation of the two peptides. The two peptides have similar turn-extended equilibria but very different helix-turn equilibria, and the difference in helical propensities reflects the fact that the helix-turn equilibrium strongly favors the turns in Ac-(Val)3-NHMe, while it favors the helices in Ac-(Ala)3-NHMe. We also computed thermodynamic decompositions of the free energy surfaces, and these revealed that the helix-turn equilibria are vastly different primarily because the changes in peptide-water interactions that accompany helix-to-turn conformational changes are qualitatively different for the two peptides.     

Interactions between relay helix and Src homology 1 (SH1) domain helix drive the converter domain rotation during the recovery stroke of myosin II

Myosin motor protein exists in two alternative conformations, prerecovery state M* and postrecovery state M**, on adenosine triphosphate binding. The details of the M*-to-M** transition, known as the recovery stroke to reflect its role as the functional opposite of the force-generating power stroke, remain elusive. The defining feature of the postrecovery state is a kink in the relay helix, a key part of the protein involved in force generation. In this article, we determine the interactions that are responsible for the appearance of the kink. We design a series of computational models that contain three other segments, relay loop, converter domain, and Src homology 1 (SH1) domain helix, with which relay helix interacts and determine their structure in accurate replica exchange molecular dynamics simulations in explicit solvent. By conducting an exhaustive combinatorial search among different models, we find that: (1) the converter domain must be attached to the relay helix during the transition, so it does not interfere with other parts of the protein and (2) the structure of the relay helix is controlled by SH1 helix. The kink is strongly coupled to the position of SH1 helix. It arises as a result of direct interactions between SH1 and the relay helix and leads to a rotation of the C-terminal part of the relay helix, which is subsequently transmitted to the converter domain.     

100ps molecular dynamic simulation of d(TATCACC)2.

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.     

Catalytic control in the EGF receptor and its connection to general kinase regulatory mechanisms

In contrast to the active conformations of protein kinases, which are essentially the same for all kinases, inactive kinase conformations are structurally diverse. Some inactive conformations are, however, observed repeatedly in different kinases, perhaps reflecting an important role in catalysis. In this review, we analyze one of these recurring conformations, first identified in CDK and Src kinases, which turned out to be central to understanding of how kinase domain of the EGF receptor is activated. This mechanism, which involves the stabilization of the active conformation of an α helix, has features in common with mechanisms operative in several other kinases.     

Towards understanding a molecular switch mechanism: thermodynamic and crystallographic studies of the signal transduction protein CheY

The signal transduction protein CheY displays an alpha/beta-parallel polypeptide folding, including a highly unstable helix alpha4 and a strongly charged active site. Helix alpha4 has been shown to adopt various positions and conformations in different crystal structures, suggesting that it is a mobile segment. Furthermore, the instability of this helix is believed to have functional significance because it is involved in protein-protein contacts with the transmitter protein kinase CheA, the target protein FliM and the phosphatase CheZ. The active site of CheY comprises a cluster of three aspartic acid residues and a lysine residue, all of which participate in the binding of the Mg(2+) needed for the protein activation. Two steps were followed to study the activation mechanism of CheY upon phosphorylation: first, we independently substituted the three aspartic acid residues in the active site with alanine; second, several mutations were designed in helix alpha 4, both to increase its level of stability and to improve its packing against the protein core. The structural and thermodynamic analysis of these mutant proteins provides further evidence of the connection between the active-site area and helix alpha 4, and helps to understand how small movements at the active site are transmitted and amplified to the protein surface.     

Molecular dynamics simulations of helix denaturation.

An understanding of the structural transitions that an alpha-helix undergoes will help to elucidate such motions in proteins and their role in protein folding. We present the results of molecular dynamics simulations to investigate these transitions in a short polyalanine peptide (13 residues) both in vacuo and in the presence of solvent. The denaturation of this peptide was monitored as a function of temperature (ranging from 5 to 200 degrees C). In vacuo, the helical state predominated at all temperatures, whereas in solution the helix melted with increasing temperature. The peptide was predominantly helical at low temperature in solution, while at intermediate temperatures the peptide spent the bulk of the time fluctuating between different conformations with intermediate amounts of helix, e.g. not completely helical nor entirely non-helical. Many of these conformations consisted of short helical segments with intervening non-helical residues. At high temperature the peptide unfolded and adopted various collapsed unstructured states. The intrahelical hydrogen bonds that break at high temperature were not fully compensated by hydrogen bonds with water molecules in the partially unfolded forms of the peptide. Increases in temperature disrupted both the helical structure and the peptide-water interactions. Water played a major but indirect role in facilitating unfolding, as opposed to specifically competing for the intrapeptide hydrogen bonds. The implications of our results to protein folding are discussed.     

Molecular dynamics simulations of an engineered T4 lysozyme exclude helix to sheet transition,and provide insights into long distance,intra‐protein switchable motion

An engineered variant of T4 lysozyme serves as a model for studying induced remote conformational changes in a full protein context. The design involves a duplicated surface helix, flanked by two loops, that switches between two different conformations spanning about 20 Å. Molecular dynamics simulations of the engineered protein, up to 1 μs, rule out α‐helix to β‐sheet transitions within the duplicated helix as suggested by others. These simulations highlight how the use of different force fields can lead to radical differences in the structure of the protein. In addition, Markov state modeling and transition path theory were employed to map a 6.6 μs simulation for possible early intermediate states and to provide insights into the onset of the switching motion. The putative intermediates involve the folding of one helical turn in the C‐terminal loop through energy driven, sequential rearrangement of nearby salt bridges around the key residue Arg63. These results provide a first step towards understanding the energetics and dynamics of a rather complicated intra‐protein motion.