Monte carlo simulations of protein-like heteropolymers.

Properties of a simple model of polypeptide chains were studied by the means of the Monte Carlo method. The chains were built on the (310) hybrid lattice. The residues interacted with long-range potential. There were two kinds of residues: hydrophobic and hydrophilic forming a typical helical pattern -HHPPHPP-. Short range potential was used to prefer helical conformations of the chain. It was found that at low temperatures the model chain formes dense and partially ordered structures (non-unique). The presence of the local potential led to an increase of helicity. The effect of the interplay between the two potentials was studied. After the collapse of the chain further annealing caused rearrangement of helical structures. Dynamic properties of the chain at low temperature depended strongly on the local chain ordering.     

Amino acid sequences of alpha-helical segments from S-carbosymethylkerateine-A. Complete sequence of a type-I segment.

The amino acid sequence of a type-I helical segment from the low-sulphur protein (S-carboxymethylkerateine-A) of wool was determined by combining automatic and manual-sequencing data. Whereas in the type-II helical segment most of the cationic groups occur in pairs, 11 of the 22 anionic residues in the sequence of the type-I segment were situated next to a second anionic residue. This suggests possible interactions between type-I and type-II helical segments in alpha-keratin. As observed with the sequence of a type-II helical segment a model constructed on 3.6 residues per turn of helix shows a line of hydrophobic residues along the helix, thereby supporting the physicochemical evidence that the molecule is predominantly helical and forms part of a coiled-coil structure. Examination of the sequence data by predictive methods indicates the possibilty of extensive sections of alpha-helix interspersed with discontinuities. The molecule contains a number of regions with peptide sequences identical with those found by other workers after enzymic digestion of fractions from oxidized wool.     

Role of recurrent hydrophobic residues in catalysis of helix formation by T cell-presented peptides in the presence of lipid vesicles

We tested the hypothesis that the recurrence of hydrophobic amino acids in a polypeptide at positions falling in an axial, hydrophobic strip if the sequence were coiled as an alpha helix, can lead to helical nucleation on a hydrophobic surface. The hydrophobic surface could anchor such residues, whereas the peptide sequence grows in a helical configuration that is stabilized by hydrogen bonds among carbonyl and amido NH groups along the peptidyl backbone of the helix, and by other intercycle interactions among amino acid side chains. Such bound, helical structures might protect peptides from proteases and/or facilitate transport to a MHC-containing compartment and thus be reflected in the selection of T cell-presented segments. Helical structure in a series of HPLC-purified peptides was estimated from circular dichroism measurements in: 1) 0.01 M phosphate buffer, pH 7.0, 2) that buffer with 45% trifluoroethanol (TFE), and 3) that buffer with di-O-hexadecyl phosphatidylcholine vesicles. By decreasing the dielectric constant of the buffer, TFE enhances intrapeptide interactions generally, whereas the lipid vesicles only provide a surface for hydrophobic interactions. The peptides varied in their strip-of-helix hydrophobicity indices (SOHHI; the mean Kyte-Doolittle hydrophobicities of residues in an axial strip of an alpha helix) and in proline content. Structural order for peptides with helical circular dichroism spectra was estimated as percentage helicity from circular dichroism theta 222 nm values and peptide concentration. A prototypic alpha helical peptide with three cycles plus two amino acids and an axial hydrophobic strip of four leucyl residues (SOHHI = 3.8) was disordered in phosphate buffer, 58% helical in that buffer with 48% TFE, and 36% helical in that buffer with vesicles. Percentage helicity in the presence of vesicles of the subset of peptides without proline followed their SOHHI values. Peptides with multiple prolyl residues had circular dichroism spectra with strong signals, but since they did not have altered spectra in the presence of vesicles relative to phosphate buffer alone, the hydrophobic surface of the vesicle did not appear to stabilize those structures.     

Evaluation of lipid‐binding properties of the N‐terminal helical segments in human apolipoprotein A‐I using fragment peptides

Although the N‐terminal region in human apolipoprotein (apo) A‐I is thought to stabilize the lipid‐free structure of the protein, its role in lipid binding is unknown. Using synthetic fragment peptides, we examined the lipid‐binding properties of the first 43 residues (1–43) of apoA‐I in comparison with residues 44–65 and 220–241, which have strong lipid affinity in the molecule. Circular dichroism measurements demonstrated that peptides corresponding to each segment have potential propensity to form α‐helical structure in trifluoroethanol. Spectroscopic and thermodynamic measurements revealed that apoA‐I (1–43) peptide has the strong ability to bind to lipid vesicles and to form α‐helical structure comparable to apoA‐I (220–241) peptide. Substitution of Tyr‐18 located at the center of the most hydrophobic region in residues 1–43 with a helix‐breaking proline resulted in the impaired lipid binding, indicating that the α‐helical structure in this region is required to trigger the lipid binding. In contrast, apoA‐I (44–65) peptide exhibited a lower propensity to form α‐helical structure upon binding to lipid, and apoA‐I (44–65/S55P) peptide exhibited diminished, but not completely impaired, lipid binding, suggesting that the central region of residues 44–65 is not pivotally involved in the formation of the α‐helical structure and lipid binding. These results indicate that the most N‐terminal region of apoA‐I molecule, residues 1–43, contributes to the lipid interaction of apoA‐I through the hydrophobic helical residues. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

NMR structures of the C-terminal segment of surfactant protein B in detergent micelles and hexafluoro-2-propanol

Although the membrane-associated surfactant protein B (SP-B) is an essential component of lung surfactant, which is itself essential for life, the molecular basis for its activity is not understood. SP-B's biophysical functions can be partially mimicked by subfragments of the protein, including the C-terminus. We have used NMR to determine the structure of a C-terminal fragment of human SP-B that includes residues 63-78. Structure determination was performed both in the fluorinated alcohol hexafluoro-2-propanol (HFIP) and in sodium dodecyl sulfate (SDS) micelles. In both solvents, residues 68-78 take on an amphipathic helical structure, in agreement with predictions made by comparison to homologous saposin family proteins. In HFIP, the five N-terminal residues of the peptide are largely unstructured, while in SDS micelles, these residues take on a well-defined compact conformation. Differences in helical residue side chain positioning between the two solvents were also found, with better agreement between the structures for the hydrophobic face than the hydrophilic face. A paramagnetic probe was used to investigate the position of the peptide within the SDS micelles and indicated that the peptide is located at the water interface with the hydrophobic face of the helix oriented inward, the hydrophilic face of the helix oriented outward, and the N-terminal residues even farther from the micelle center than those on the hydrophilic face of the alpha-helix. Interactions of basic residues of SP-B with anionic lipid headgroups are known to have an impact on function, and these studies demonstrate structural ramifications of such interactions via the differences observed between the peptide structures determined in HFIP and SDS.     

Molecular determinants of the hanatoxin binding in voltage-gated K+-channel drk1

The carboxyl terminus of S3 segment (S3(C)) in voltage-gated potassium channels was proposed to bear the binding site for gating modifier toxins like Hanatoxin and a helical secondary structural arrangement was suggested. Due to the lack of complete structure in high resolution for such a channel molecule, no further direct experimental data to elucidate the mechanism for their binding conformations could thus far be derived. In order to examine the putative three-dimensional structure of S3(C) and to illustrate the residues required for Hanatoxin binding, molecular simulation and docking were performed, based on the solution structure of Hanatoxin and the structural information from lysine-scanning results for S3(C) fragment. From our results, it is indicated that both hydrophobic and electrostatic interactions are utilized to stabilize the toxin binding. Detailed docking residues and appropriate orientation for binding regarding hydrophobic/-philic environments are also described. Compared with the functional data proposed by previous studies, the helical structural arrangement for the C-terminus of S3 segment in voltage-gated potassium channels can therefore be further emphasized.     

NMR studies of the low-density lipoprotein receptor-binding peptide of apolipoprotein E bound to dodecylphosphocholine micelles.

Circular dichroism and NMR spectroscopy have been used to determine the structure of the low-density lipoprotein (LDL) receptor-binding peptide, comprising residues 130-152, of the human apolipoprotein E. This peptide has little persistent three-dimensional structure in solution, but when bound to micelles of dodecylphosphocholine (DPC) it adopts a predominantly alpha-helical structure. The three-dimensional structure of the DPC-bound peptide has been determined by using 1H-NMR spectroscopy: the structure derived from NOE-based distance constraints and restrained molecular dynamics is largely helical. The derived phi and psi angle order parameters show that the helical structure is well defined but with some flexibility that causes the structures not to be superimposable over the full peptide length. Deuterium exchange experiments suggest that many peptide amide groups are readily accessible to the solvent, but those associated with hydrophobic residues exchange more slowly, and this helix is thus likely to be positioned on the surface of the DPC micelles. In this conformation the peptide has one hydrophobic face and two that are rich in basic amino acid side chains. The solvent-exposed face of the peptide contains residues previously shown to be involved in binding to the LDL receptor.     

The structure of Escherichia coli heat-stable enterotoxin b by nuclear magnetic resonance and circular dichroism.

The heat-stable enterotoxin b (STb) is secreted by enterotoxigenic Escherichia coli that cause secretory diarrhea in animals and humans. It is a 48-amino acid peptide containing two disulfide bridges, between residues 10 and 48 and 21 and 36, which are crucial for its biological activity. Here, we report the solution structure of STb determined by two- and three-dimensional NMR methods. Approximate interproton distances derived from NOE data were used to construct structures of STb using distance-geometry and simulated annealing procedures. The NMR-derived structure shows that STb is helical between residues 10 and 22 and residues 38 and 44. The helical structure in the region 10-22 is amphipathic and exposes several polar residues to the solvent, some of which have been shown to be important in determining the toxicity of STb. The hydrophobic residues on the opposite face of this helix make contacts with the hydrophobic residues of the C-terminal helix. The loop region between residues 21 and 36 has another cluster of hydrophobic residues and exposes Arg 29 and Asp 30, which have been shown to be important for intestinal secretory activity. CD studies show that reduction of disulfide bridges results in a dramatic loss of structure, which correlates with loss of function. Reduced STb adopts a predominantly random-coil conformation. Chromatographic measurements of concentrations of native, fully reduced, and single-disulfide species in equilibrium mixtures of STb in redox buffers indicate that the formation of the two disulfide bonds in STb is only moderately cooperative. Similar measurements in the presence of 8 M urea suggest that the native secondary structure significantly stabilizes the disulfide bonds.     

The sequence TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase displays structural ambivalence and interconverts between alpha-helical and beta-hairpin conformations mediated by collapsed conformational states.

The peptide TGAAKAVALVL from glyceraldehyde-3-phosphate dehydrogenase adopts a helical conformation in the crystal structure and is a site for two hydrated helical segments, which are thought to be helical folding intermediates. Overlapping sequences of four to five residues from the peptide, sample both helical and strand conformations in known protein structures, which are dissimilar to glyceraldehyde-3-phosphate dehydrogenase suggesting that the peptide may have a structural ambivalence. Molecular dynamics simulations of the peptide sequence performed for a total simulation time of 1.2 micros, starting from the various initial conformations using GROMOS96 force field under NVT conditions, show that the peptide samples a large number of conformational forms with transitions from alpha-helix to beta-hairpin and vice versa. The peptide, therefore, displays a structural ambivalence. The mechanism from alpha-helix to beta-hairpin transition and vice versa reveals that the compact bends and turns conformational forms mediate such conformational transitions. These compact structures including helices and hairpins have similar hydrophobic radius of gyration (Rgh) values suggesting that similar hydrophobic interactions govern these conformational forms. The distribution of conformational energies is Gaussian with helix sampling lowest energy followed by the hairpins and coil. The lowest potential energy of the full helix may enable the peptide to take up helical conformation in the crystal structure of the glyceraldehyde-3-phosphate dehydrogenase, even though the peptide has a preference for hairpin too. The relevance of folding and unfolding events observed in our simulations to hydrophobic collapse model of protein folding are discussed.     

Detection of secondary structure elements in proteins by hydrophobic cluster analysis.

Hydrophobic cluster analysis (HCA) is a protein sequence comparison method based on alpha-helical representations of the sequences where the size, shape and orientation of the clusters of hydrophobic residues are primarily compared. The effectiveness of HCA has been suggested to originate from its potential ability to focus on the residues forming the hydrophobic core of globular proteins. We have addressed the robustness of the bidimensional representation used for HCA in its ability to detect the regular secondary structure elements of proteins. Various parameters have been studied such as those governing cluster size and limits, the hydrophobic residues constituting the clusters as well as the potential shift of the cluster positions with respect to the position of the regular secondary structure elements. The following results have been found to support the alpha-helical bidimensional representation used in HCA: (i) there is a positive correlation (clearly above background noise) between the hydrophobic clusters and the regular secondary structure elements in proteins; (ii) the hydrophobic clusters are centred on the regular secondary structure elements; (iii) the pitch of the helical representation which gives the best correspondence is that of an alpha-helix. The correspondence between hydrophobic clusters and regular secondary structure elements suggests a way to implement variable gap penalties during the automatic alignment of protein sequences.     

Evaluation of the protein interfaces that form an intermolecular four-helix bundle as studied by computer simulation

Rational design of protein surface is important for creating higher order protein structures, but it is still challenging. In this study, we designed in silico the several binding interfaces on protein surfaces that allow a de novo protein–protein interaction to be formed. We used a computer simulation technique to find appropriate amino acid arrangements for the binding interface. The protein–protein interaction can be made by forming an intermolecular four-helix bundle structure, which is often found in naturally occurring protein subunit interfaces. As a model protein, we used a helical protein, YciF. Molecular dynamics simulation showed that a new protein–protein interaction is formed depending on the number of hydrophobic and charged amino acid residues present in the binding surfaces. However, too many hydrophobic amino acid residues present in the interface negatively affected on the binding. Finally, we found an appropriate arrangement of hydrophobic and charged amino acid residues that induces a protein–protein interaction through an intermolecular four-helix bundle formation.     

Trimerization specificity in HIV-1 gp41: analysis with a GCN4 leucine zipper model

The envelope glycoprotein of human immunodeficiency virus type 1 (HIV-1) consists of a complex of two noncovalently associated subunits, gp120 and gp41. Formation of gp120/gp41 oligomers is thought to be dependent on a 4-3 hydrophobic (heptad) repeat located in the amino-terminal region of the gp41 molecule. We have investigated the role of this heptad repeat in determining the oligomeric structure of gp41 by introducing its buried core residues into the first (a) and fourth (d) positions of the GCN4 leucine-zipper dimerization domain. The mutant peptides fold into trimeric, helical structures, as shown by circular dichroism and equilibrium sedimentation centrifugation. The 2.4 A resolution crystal structure of one such trimer reveals a parallel three-stranded, alpha-helical coiled coil. Thus, the buried core residues from the gp41 heptad repeat direct trimer formation. We suggest that the conserved amino-terminal heptad repeat within the gp41 ectodomain possesses trimerization specificity.     

Comparison of the conformation and electrostatic surface properties of magainin peptides bound to sodium dodecyl sulfate and dodecylphosphocholine micelles

The role played by noncovalent interactions in inducing a stable secondary structure onto the sodium dodecyl sulfate (SDS) and dodecylphosphocholine (DPC) micelle-bound conformations of (Ala(8,13,18))magainin 2 amide and the DPC micelle bound conformation of magainin 1 were determined. Two-dimensional NMR and molecular modeling investigations indicated that (Ala(8,13,18))magainin 2 amide bound to DPC micelles adopts a alpha-helical secondary structure involving residues 2-16. The four C-terminal residues converge to a lose beta-turn structure. (Ala(8,13,18))magainin 2 amide bound to SDS miscelles adopts a alpha-helical secondary structure involving residues 7-18. The C- and N-terminal residues exhibited a great deal of conformational flexibility. Magainin 1 bound to DPC micelles adopts a alpha-helical secondary structure involving residues 4-19. The C-terminal residues converge to a lose beta-turn structure. The results of this investigation indicate hydrophobic interactions are the major contributors to stabilizing the induced helical structure of the micelle-bound peptides. Electrostatic interactions between the polar head groups of the micelle and the cationic side chains of the peptides define the positions along the peptide backbone where the helical structures begin and end.     

Solution structure of the osteogenic 1-31 fragment of the human parathyroid hormone

The solution conformations of a selectively osteogenic 1-31 fragment of the human parathyroid hormone (hPTH), hPTH(1-31)NH(2), have been characterized by use of very high field NMR spectroscopy at 800 MHz. The combination of the CalphaH proton and (13)Calpha chemical shifts, (3)J(NH)(alpha) coupling constants, NH proton temperature coefficients, and backbone NOEs reveals that the hPTH(1-31)NH(2) peptide has well-formed helical structures localized in two distinct segments of the polypeptide backbone. There are also many characteristic NOEs defining specific side-chain/backbone and side-chain/side-chain contacts within both helical structures. The solution structure of hPTH(1-31)NH(2) contains a short N-terminal helical segment for residues 3-11, including the helix capping residues 3 and 11 and a long C-terminal helix for residues 16-30. The two helical structures are reinforced by well-defined capping motifs and side-chain packing interactions within and at both ends of these helices. On one face of the C-terminal helix, there are side-chain pairs of Glu22-Arg25, Glu22-Lys26, and Arg25-Gln29 that can form ion-pair and/or hydrogen bonding interactions. On the opposite face of this helix, there are characteristic hydrophobic interactions involving the aromatic side chain of Trp23 packing against the aliphatic side chains of Leu15, Leu24, Lys27, and Leu28. There is also a linear array of hydrophobic residues from Val2, to Leu7, to Leu11 and continuing on to residues His14 and Leu15 in the hinge region and to Trp23 in the C-terminal helix. Capping and hydrophobic interactions at the end of the N-terminal and at the beginning of the C-terminal helix appear to consolidate the helical structures into a V-shaped overall conformation for at least the folded population of the hPTH(1-31)NH(2) peptide. Stabilization of well-folded conformations in this linear 1-31 peptide fragment and possibly other analogues of human PTH may have a significant impact on the biological activities of the PTH peptides in general and specifically for the osteogenic/anabolic activities of bone-building PTH analogues.     

Molecular dynamics simulations predict a tilted orientation for the helical region of dynorphin A(1-17) in dimyristoylphosphatidylcholine bilayers

The structural properties of the endogenous opioid peptide dynorphin A(1-17) (DynA), a potential analgesic, were studied with molecular dynamics simulations in dimyristoylphosphatidylcholine bilayers. Starting with the known NMR structure of the peptide in dodecylphosphocholine micelles, the N-terminal helical segment of DynA (encompassing residues 1-10) was initially inserted in the bilayer in a perpendicular orientation with respect to the membrane plane. Parallel simulations were carried out from two starting structures, systems A and B, that differ by 4 A in the vertical positioning of the peptide helix. The complex consisted of approximately 26,400 atoms (dynorphin + 86 lipids + approximately 5300 waters). After >2 ns of simulation, which included >1 ns of equilibration, the orientation of the helical segment of DynA had undergone a transition from parallel to tilted with respect to the bilayer normal in both the A and B systems. When the helix axis achieved a approximately 50 degrees angle with the bilayer normal, it remained stable for the next 1 ns of simulation. The two simulations with different starting points converged to the same final structure, with the helix inserted in the bilayer throughout the simulations. Analysis shows that the tilted orientation adopted by the N-terminal helix is due to specific interactions of residues in the DynA sequence with phospholipid headgroups, water, and the hydrocarbon chains. Key elements are the "snorkel model"-type interactions of arginine side chains, the stabilization of the N-terminal hydrophobic sequence in the lipid environment, and the specific interactions of the first residue, Tyr. Water penetration within the bilayer is facilitated by the immersed DynA, but it is not uniform around the surface of the helix. Many water molecules surround the arginine side chains, while water penetration near the helical surface formed by hydrophobic residues is negligible. A mechanism of receptor interaction is proposed for DynA, involving the tilted orientation observed from these simulations of the peptide in the lipid bilayer.     

Solution structure of prosurvival Mcl-1 and characterization of its binding by proapoptotic BH3-only ligands

The B cell lymphoma-2 (Bcl-2) homologs myeloid cell leukemia-1 (Mcl-1) and A1 are prosurvival factors that selectively bind a subset of proapoptotic Bcl homology (BH) 3-only proteins. To investigate the molecular basis of the selectivity, we determined the solution structure of the C-terminal Bcl-2-like domain of Mcl-1. This domain shares features expected of a prosurvival Bcl-2 protein, having a helical fold centered on a core hydrophobic helix and a surface-exposed hydrophobic groove for binding its cognate partners. A number of residues in the binding groove differentiate Mcl-1 from its homologs, and in contrast to other Bcl-2 homologs, Mcl-1 has a binding groove in a conformation intermediate between the open structures characterized by peptide complexes and the closed state observed in unliganded structures. Mutagenesis of potential binding site residues was used to probe the contributions of groove residues to the binding properties of Mcl-1. Although mutations in Mcl-1 had little impact on binding, a single mutation in the BH3-only ligand Bad enabled it to bind both Mcl-1 and A1 while retaining its binding to Bcl-2, Bcl-xL, and Bcl-w. Elucidating the selective action of certain BH3-only ligands is required for delineating their mode of action and will aid the search for effective BH3-mimetic drugs.     

Structural and biophysical properties of the integrin-associated cytoskeletal protein talin

Talin is a large cytoskeletal protein (2541 amino acid residues) which plays a key role in integrin-mediated events that are crucial for cell adhesion, migration, proliferation and survival. This review summarises recent work on the structure of talin and on some of the structurally better defined interactions with other proteins. The N-terminal talin head (approx. 50 kDa) consists of an atypical FERM domain linked to a long flexible rod (approx. 220 kDa) made up of a series of amphipathic helical bundle domains. The F3 FERM subdomain in the head binds the cytoplasmic tail of integrins, but this interaction can be inhibited by an interaction of F3 with a helical bundle in the talin rod, the so-called “autoinhibited form” of the molecule. The talin rod contains a second integrin-binding site, at least two actin-binding sites and a large number of binding sites for vinculin, which is important in reinforcing the initial integrin–actin link mediated by talin. The vinculin binding sites are defined by hydrophobic residues buried within helical bundles, and these must unfold to allow vinculin binding. Recent experiments suggest that this unfolding may be mediated by mechanical force exerted on the talin molecule by actomyosin contraction.