De novo design of transmembrane helix–helix interactions and measurement of stability in a biological membrane

Membrane proteins regulate a large number of cellular functions, and have great potential as tools for manipulation of biological systems. Developing these tools requires a robust and quantitative understanding of membrane protein folding and interactions within the bilayer. With this in mind, we have designed a series of proteins to probe the net thermodynamic contribution of well-known sequence motifs to transmembrane helix-helix association in a biological membrane. The proteins were designed from first principles (de novo) using current knowledge about membrane insertion and stabilizing interaction motifs. A simple poly-Leu “scaffold” was decorated with individual helix interaction motifs (G-XXX-G, polar residues, heptad repeat) to create transmembrane helix–helix interactions of increasing strength. The GALLEX assay, an in vivo assay for measurement of transmembrane helix self-association, was combined with computational methods to characterize the relative strength and mode of interaction for each sequence. In addition, the apparent free energy contribution (ΔΔG^app) of each motif to transmembrane helix self-association was measured in a biological membrane, results that are the first of their kind for these de novo designed sequences, and suggest that the free energy barrier to overcoming weak association is quite small (< 1.4 kcal mol^− 1) in a natural membrane. By quantifying and rationalizing the contribution of key motifs to transmembrane helix association, our work offers a route to direct the design of novel sequences for use in biotechnology or synthetic biology (e.g. molecular switches) and to predict the effects of sequence modification in known transmembrane domains (for control of cellular processes).     

Selection of inhibitory peptides for Aurora-A kinase from a phage-displayed library of helix–loop–helix peptides

Conformationally constrained peptide libraries have been made by grafting randomized amino acid sequences onto a rigid scaffold derived from natural proteins. Here, as a library scaffold, we propose a de novo designed helix–loop–helix motif. We constructed a peptide library of the loop region and screened against Aurora-A, which is a member of the Aurora family of serine/threonine protein kinases, to successfully isolate the inhibitory peptides. A semi-rational strategy, which combines phage-displayed libraries and de novo designed peptides, would provide a new way to generate selective peptide inhibitors for the protein kinase family.     

Conformational and dynamics simulation study of antimicrobial peptide hedistin—heterogeneity of its helix–turn–helix motif

Hedistin is an antimicrobial peptide isolated from the coelomocytes of Nereis diversicolor, possessing activity against a large spectrum of bacteria including the methicillin resistant Staphylococcus aureus and Vibrio alginolyticus. The three-dimensional structure of hedistin in both aqueous solution and deuterated dodecylphosphocholine (DPC) micelles was examined using circular dichroism (CD) and nuclear magnetic resonance (NMR) techniques. And, the early events of the antibacterial process of hedistin were simulated using palmitoyl-oleoyl-phophatidylcholine (POPC) lipid bilayers and molecular dynamics (MD) simulation methods. Hedistin lacks secondary structure in aqueous solution, however, in DPC micelles, it features with a heterogeneous helix–turn–helix moiety and exhibits obvious amphipathic nature. The turn region (residues Val9–Thr12) in the moiety is a four-residue hinge, lying in between the first N-terminal α-helix (residues Leu5–Lys8) and the second α-helix (residues Val13–Ala17) regions and causing an ～ 120° angle between the axes of the two helices. The segmental and nonlinear nature of hedistin structure is referred to as the heterogeneity of its helix–turn–helix motif which was found to be corresponding to a kind of discrete dynamics behavior, herein coined as its dynamical heterogeneity, at the early stage (0–50 ns) of the MD simulations. That is, the first helix segment, prior to (at 310 K) or following (at 363 K) the second helix, binds to the lipid head-group region and subsequently permeates into the hydrophobic lipid tail region, and the hinge is the last portion entering the lipid environment. This result implies that hedistin may adopt a “carpet” model action when disrupting bacterial membrane.     

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.     

The transmembrane domain of caveolin-1 exhibits a helix–break–helix structure

Caveolin is an integral membrane protein that is found in high abundance in caveolae. Both the N- and C- termini lie on the same side of the membrane, and the transmembrane domain has been postulated to form an unusual intra-membrane horseshoe configuration. To probe the structure of the transmembrane domain, we have prepared a construct of caveolin-1 that encompasses residues 96–136 (the entire intact transmembrane domain). Caveolin-1(96–136) was over-expressed and isotopically labeled in E. coli, purified to homogeneity, and incorporated into lyso-myristoylphosphatidylglycerol micelles. Circular dichroism and NMR spectroscopy reveal that the transmembrane domain of caveolin-1 is primarily α-helical (57–65%). Furthermore, chemical shift indexing reveals that the transmembrane domain has a helix–break–helix structure which could be critical for the formation of the intra-membrane horseshoe conformation predicted for caveolin-1. The break in the helix spans residues 108 to 110, and alanine scanning mutagenesis was carried out to probe the structural significance of these residues. Our results indicate that mutation of glycine 108 to alanine does not disrupt the structure, but mutation of isoleucine 109 and proline 110 to alanine dramatically alters the helix–break–helix structure. To explore the structural determinants further, additional mutagenesis was performed. Glycine 108 can be substituted with other small side chain amino acids (i.e. alanine), leucine 109 can be substituted with other β-branched amino acids (i.e. valine), and proline 110 cannot be substituted without disrupting the helix–break–helix structure.     

Kinetic aspects of the DNA helix—Cruciform transition

A kinetic model for the helix-cruciform transition is presented, mean lifetimes for the cruciform states are calculated and shown to be inconsistent with the notion of metastability.     

Contributions of cation-π interactions to the collagen triple helix stability

Cation–π interactions are found to be an important noncovalent force in proteins. Collagen is a right-handed triple helix composed of three left-handed PPII helices, in which (X–Y-Gly) repeats dominate in the sequence. Molecular modeling indicates that cation–π interactions could be formed between the X and Y positions in adjacent collagen strands. Here, we used a host–guest peptide system: (Pro-Hyp-Gly)₃-(Pro-Y-Gly-X-Hyp-Gly)-(Pro-Hyp-Gly)₃, where X is an aromatic residue and Y is a cationic residue, to study the cation–π interaction in the collagen triple helix. Circular dichroism (CD) measurements and T_m data analysis show that the cation–π interactions involving Arg have a larger contribution to the conformational stability than do those involving Lys, and Trp forms a weaker cation–π interaction with cationic residues than expected as a result of steric effects. The results also show that the formation of cation–π interactions between Arg and Phe depends on their relative positions in the strand. Moreover, the fluorinated and methylated Phe substitutions show that an electron-withdrawing or electron-donating substituent on the aromatic ring can modulate its π–electron density and the cation–π interaction in collagen. Our data demonstrate that the cation–π interaction could play an important role in stabilizing the collagen triple helix.     

Structural dynamics of a single-stranded RNA–helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3′ end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, ¹³C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A₆-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA–helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.     

Structure of D-DNA: 8-fold or 7-fold helix?

We have shown that both right- and left-handed uniform helical models (RU and LU models) could be built to give satisfactory agreement with the fibre diffraction data of poly[d(I-C)] in the D-form. Atomic coordinates of these two models are reported in the present work. Molecular transforms of these two models, as well as of the recently published Hoogsteen base-paired 7-fold helical structure of Drew and Dickerson, are given. In view of the work of Drew and Dickerson, attention is drawn to the presence of clear 004 and 008 reflections in the diffraction patterns of poly[d(I-C)] and poly[d(A-T)]. The available data strongly suggest an 8-fold helical structure for the D-form of DNA.     

Sequence preference in DNA binding: de novo designed helix–turn–helix metallopeptides recognize a family of DNA target sites

The DNA-binding behavior and target sequences of two designed metallopeptides have been investigated with an iterative electrophoresis mobility shift assay followed by PCR amplification, and by circular dichroism spectroscopy. Peptides P3W and P5b were designed based on the structural similarity of the helix–turn–helix motif of homeodomains and the EF-hand motifs of calmodulin, as previously described for P3W. Like P3W, P5b binds both Eu(III) (K _d=12.6±1.9 μM) and Ca(II) (K _d=70±8 μM) with reasonable affinity. Binding selection from a library of randomized 8-mer DNA oligonucleotide sequences identified one target family for CaP5b [5′-pur-T-pur-G-(G/C)-3′], and two target sites for CaP3W [5′-(A/T)-G-G-G-(T/C)-3′ and 5′-A-T-(G/T)-T-G-3′]. Circular dichroism studies indicate that unlike EuP3W, EuP5b is poorly folded in the absence of DNA. In the presence of DNA containing target-binding sites for both peptides, both EuP3W and EuP5b increase in helical content, in the latter case significantly. These results suggest that EuP5b binding to target DNA involves an induced-fit mechanism. These small chimeric metallopeptides have been found to bind selectively to DNA targets, analogous to natural protein–DNA interactions. This corroborates our earlier conclusions (J. Am. Chem. Soc. 125:6656, 2003) that sequence-preferential DNA cleavage by Ce(IV)P3W was due to sequence recognition. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Identification of helix capping and β-turn motifs from NMR chemical shifts

We present an empirical method for identification of distinct structural motifs in proteins on the basis of experimentally determined backbone and ¹³C^β chemical shifts. Elements identified include the N-terminal and C-terminal helix capping motifs and five types of β-turns: I, II, I′, II′ and VIII. Using a database of proteins of known structure, the NMR chemical shifts, together with the PDB-extracted amino acid preference of the helix capping and β-turn motifs are used as input data for training an artificial neural network algorithm, which outputs the statistical probability of finding each motif at any given position in the protein. The trained neural networks, contained in the MICS (motif identification from chemical shifts) program, also provide a confidence level for each of their predictions, and values ranging from ca 0.7–0.9 for the Matthews correlation coefficient of its predictions far exceed those attainable by sequence analysis. MICS is anticipated to be useful both in the conventional NMR structure determination process and for enhancing on-going efforts to determine protein structures solely on the basis of chemical shift information, where it can aid in identifying protein database fragments suitable for use in building such structures.     

Two subtypes of HLA-B51 differing by substitution at position 171 of the α2 helix

Newly defined antigens of the B5, B35 cross-reacting group have been found in Japanese and North American Indians. Nucleotide sequencing of the alleles encoding the Japanese B5.35 antigen and the variant B5 antigen from the Piman Indians show them to be identical. This new allele, B ^* 5102, differs from B ^* 5101 by a single nucleotide substitution that changes residue 171 from histidine to tyrosine. Residue 171, which is part of the ₂ helix, is believed to contribute directly to peptide interaction in the A pocket of the binding groove and is either histidine or tyrosine in all HLA-A, B, C heavy chains. Tyrosine 171 is shared by B ^* 5102, B ^* 3501, B ^* 3502, and B ^* 5301 and must be responsible for the serological cross-reactivities of these molecules not shared with B ^* 5101. Stimulation of lymphocytes from a B ^* 5101 positive donor with B ^* 5102 positive cells failed to generate cytotoxic T cells with specificity for the difference between these molecules. However, one out of five clones of cytotoxic T cells raised against B ^* 5101 failed to lyse targets expressing B ^* 5102. Substitution of histidine for tyrosine at residue 171 affected recognition of HLA-B35-restricted human minor histocompatibility antigen-specific T cell clones.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession number M68964.     

Transmembrane helix–helix interactions are modulated by the sequence context and by lipid bilayer properties

Folding of polytopic transmembrane proteins involves interactions of individual transmembrane helices, and multiple TM helix–helix interactions need to be controlled and aligned to result in the final TM protein structure. While defined interaction motifs, such as the GxxxG motif, might be critically involved in transmembrane helix–helix interactions, the sequence context as well as lipid bilayer properties significantly modulate the strength of a sequence specific transmembrane helix–helix interaction. Structures of 11 transmembrane helix dimers have been described today, and the influence of the sequence context as well as of the detergent and lipid environment on a sequence specific dimerization is discussed in light of the available structural information. This article is part of a Special Issue entitled: Protein Folding in Membranes.     

Hexafluoroisopropanol-induced helix–sheet transition of stem bromelain: Correlation to function

Stem bromelain is a proteolytic phytoprotein with a variety of therapeutic effects. Understanding its structural properties could provide insight into the mechanisms underlying its clinical utility. Stem bromelain was evaluated for its conformational and folding properties at the pH conditions it encounters when administered orally. It exists as a partially folded intermediate at pH 2.0. The conformational changes to this intermediate state were evaluated using fluorinated alcohols known to induce changes similar to those seen in vivo. Studies using circular dichroism, fluorescence emission spectroscopy, binding of the hydrophobic dye 1-anilino-8-naphthalene sulfonic acid and mass spectrometry indicate that treatment with 10–30% hexafluoroisopropanol induces the partially folded intermediate to adopt much of the native protein's secondary structure, but only a rudimentary tertiary structure, characteristic of the molten globule state. Addition of slightly higher concentrations of hexafluoroisopropanol caused transformation from an α-helix to a β-sheet and induced formation of a compact nonnative structure. This nonnative form was more inhibitory of cell survival than either the native or the partially folded intermediate forms, as measured by enhanced suppression of proliferative cues (e.g., extracellular-signal-regulated kinase) and initiation of apoptotic events. The nonnative form also showed better antitumorigenic properties, as evaluated using an induced two-stage mouse skin papilloma model. In contrast, the nonnative state showed only a fraction of the proteolytic activity of the native form. This study demonstrates that hexafluoroisopropanol can induce a conformational change in stem bromelain to a form with potentially useful therapeutic properties different from those of the native protein.