Analysis of peptide design in four-, five-, and six-helix bundle template assembled synthetic protein molecules

Four-, five-, and six-helix bundle template assembled synthetic proteins (TASPs) have been synthesized using disulfide bonds between cavitand templates and peptides, and characterized in terms of stability and structural specificity. The peptide sequence (CGGGEELLKKLEE LLKKG) used was originally designed for a four-helix bundle. The TASPs were analyzed using CD spectroscopy, chemical denaturation studies, NMR spectroscopy, sedimentation equilibria studies, and hydrophobic dye binding studies to determine the effect of a single peptide sequence when incorporated into bundles with different numbers of helices. If the design was indeed idealized for a four-helix bundle, then the five- and six-helix bundles should be less stable and manifest lower conformational specificity. The TASPs all demonstrated high stability and cooperative unfolding. However, the four-helix bundle was found to be significantly more stable and nativelike compared to the five- and six-helix bundles. This suggests that the peptide sequence is specific to the four-helix bundle, as designed. This result demonstrates the ability to design de novo proteins with specified structure, not just generic stability.     

Experimental and theoretical studies of the three-dimensional structure of human interleukin-4

The structure of human interleukin 4 (IL-4) was predicted utilizing a series of experimental and theoretical techniques. Circular Dichroism (CD) spectroscopy indicated that IL-4 belonged to the all alpha-helix class of protein structures. Secondary structure prediction, site-directed mutagenesis, and CD spectroscopy suggested a predominantly alpha-helical structure, consistent with a four-helix bundle structural motif. A human/mouse IL-4 chimera was constructed to qualitatively evaluate alternative secondary structure predictions. The four predicted helices were assembled into tertiary structures using established algorithms. The mapping of three disulfide bridges in IL-4 provided additional constraints on possible tertiary structures. Using accessible surface contact area as a criterion, the most suitable structures were right handed all antiparallel four-helix bundles with two overhand loop connections. Successful loop closure and incorporation of the three disulfide constraints were possible while maintaining the expected shape, solvent accessibility, and steric interactions between loops and helices. Lastly, energy minimization was used to regularize the chain.     

Relative stability of de novo four-helix bundle proteins: insights from coarse grained molecular simulations

We use a recently developed coarse-grained computational model to investigate the relative stability of two different sets of de novo designed four-helix bundle proteins. Our simulations suggest a possible explanation for the experimentally observed increase in stability of the four-helix bundles with increasing sequence length. In details, we show that both short subsequences composed only by polar residues and additional nonpolar residues inserted, via different point mutations in ad hoc positions, seem to play a significant role in stabilizing the four-helix bundle conformation in the longer sequences. Finally, we propose an additional mutation that rescues a short amino acid sequence that would otherwise adopt a compact misfolded state. Our work suggests that simple computational models can be used as a complementary tool in the design process of de novo proteins.     

Identification of a stability determinant on the edge of the Tet repressor four-helix bundle dimerization motif

Isofunctional tetracycline repressor (TetR) proteins isolated from different bacteria show a sequence identity between 38 and 88% of the residues. Their active state is a homodimer formed by a four-alpha-helix bundle as the main interaction motif. We utilize this sequence variation of isofunctional proteins to determine residues contributing to the stability of the four-helix bundle. The thermodynamic stabilities of two TetR proteins with 63% sequence identity were determined by urea-induced reversible denaturation followed by fluorescence and circular dichroism. Both methods yield identical results. The deltaG(o)U (H2O) values are 60 and 75 kJ x mol(-1). We have constructed TetR hybrid proteins derived from these wild types to identify the determinant leading to the 15 kJ x mol(-1) stability difference. Successive size reduction of the exchanged portion yielded two single residues affecting the overall protein stability. The P184Q exchange leads to a more stable protein, whereas the G181D exchange located at the solvent's exposed edge of the four-helix bundle is solely responsible for the reduced stability. Additional mutants based on crystal structures of TetR do not reveal any hint for steric interference of the Asp181 side chain with neighboring residues. Thus, this is an example for the role played by surface-exposed turn residues for the stability of four-helix bundles. We assume that the larger conformational flexibility of Gly and the reduction of the negative surface charge could favor formation of the turn on the edge of the four-helix bundle.     

Structural characterization by computer experiments of the lipid-free LDL-receptor-binding domain of apolipoprotein E.

The structure and dynamics of the lipid-free LDL-receptor-binding domain of apolipoprotein E (apoE-RBD) has been investigated by Molecular Dynamics Simulations. ApoE-RBD in its monomeric lipid-free form is a singular four-helix bundle made up of four elongated amphipathic helices. Analysis of one 1.5 ns molecular dynamics trajectory of apoE-RBD performed in water indicates that the lipid-free domain adopts a structure that exhibits characteristics found in native proteins: it has very stable helices and presents a compact structure. Yet its interior exhibits a larger number of transient atomic-size cavities relative to that found in other proteins of similar size and its apolar side chains are more mobile. The latter features distinguish the elongated four-helix bundle as a slightly disordered structure, which shows a structural likeness with some de novo designed four-helix bundle proteins and shares with the latter a leucine-rich residue composition. We anticipate that these unique properties compared with other native helix bundles may be related to the postulated ability of apoE-RBD to undergo an opening of its bundle upon interaction with phospholipids. The distribution of empty cavities computed along the trajectory in the interface regions between the different pairs of helices reveals that the tertiary contacts in one of the interfaces are weaker suggesting that this particular interface could be more easily ruptured upon lipid association.     

Dynamic Monte Carlo simulations of globular protein folding/unfolding pathways. II. Alpha-helical motifs

Dynamic Monte Carlo simulations of the folding pathways of alpha-helical protein motifs have been undertaken in the context of a diamond lattice model of globular proteins. The first question addressed in the nature of the assembly process of an alpha-helical hairpin. While the hairpin could, in principle, be formed via the diffusion-collision-adhesion of isolated performed helices, this is not the dominant mechanism of assembly found in the simulations. Rather, the helices that form native hairpins are constructed on-site, with folding initiating at or near the turn in almost all cases. Next, the folding/unfolding pathways of four-helix bundles having tight bends and one and two long loops in the native state are explored. Once again, an on-site construction mechanism of folding obtains, with a hairpin forming first, followed by the formation of a three-helix bundle, and finally the fourth helix of the native bundle assembles. Unfolding is essentially the reverse of folding. A simplified analytic theory is developed that reproduces the equilibrium folding transitions obtained from the simulations remarkably well and, for the dominant folding pathway, correctly identifies the intermediates seen in the simulations. The analytic theory provides the free energy along the reaction co-ordinate and identifies the transition state for all three motifs as being quite close to the native state, with three of the four helices assembled, and approximately one turn of the fourth helix in place. The transition state is separated from the native conformation by a free-energy barrier of mainly energetic origin and from the denatured state by a barrier of mainly entropic origin. The general features of the folding pathway seen in all variants of the model four-helix bundles are similar to those observed in the folding of beta-barrel, Greek key proteins; this suggests that many of the qualitative aspects of folding are invariant to the particular native state topology and secondary structure.     

A comparison of three- and four-helix bundle TASP molecules.

We have designed, synthesized and characterized three- and four-helix bundle template-assembled synthetic proteins CTASPs). The TASPs were synthesized using disulphide bonds between the peptides and either the cyclotribenzylene (CTB) template, or the cavitand (BOWL) template, to form the three- and four helix bundles, respectively. The TASPs were constructed using peptides that were linked via their N-termini (peptide CGGGEELLKKXEELLKKG, where X = L, I, Nle or V), or via their C-termini (peptide GEELLKKLEELLKKGGGC). Each TASP was assayed for its structure, stability, 'native-like' characteristics and whether it was a monomer in solution. All TASPs were found to be highly helical, and highly resistant to chemical denaturation using guanidine hydrochloride (GnHCl). Analysis of the GnHCl-induced unfolding curves of the different TASPs demonstrated stability differences based on the number of helices in the bundle, the end of the helix that was attached to the template, and the identity of the core amino acid. The TASPs all had molten-globule structure, which is (generally) consistent with a degenerate sequence in the core. The four-helix bundle TASPs appeared to be monomers in solution, whereas there is some evidence that the three-helix bundle TASPs are weakly self-associating.     

Probing domain swapping for the neuronal SNARE complex with electron paramagnetic resonance

Highly conserved soluble N-ethylmaleimide sensitive factor attachment protein receptor (SNARE) proteins control membrane fusion at synapses. The target plasma membrane-associated SNARE proteins and the vesicle-associated SNARE protein assemble into a parallel four-helix bundle. Using a novel EPR approach, it is found that the SNARE four-helix bundles are interconnected via domain swapping that is achieved by substituting one of the two SNAP-25 helices with the identical helix from the second four-helical bundle. Domain swapping is likely to play a role in the multimerization of the SNARE complex that is required for successful membrane fusion. The new EPR application employed here should be useful to study other polymerizing proteins.     

Conformational flexibility of the acetylcholinesterase tetramer suggested by x-ray crystallography.

Acetylcholinesterase, a polymorphic enzyme, appears to form amphiphilic and nonamphiphilic tetramers from a single splice variant; this suggests discrete tetrameric arrangements where the amphipathic carboxyl-terminal sequences can be either buried or exposed. Two distinct, but related crystal structures of the soluble, trypsin-released tetramer of acetylcholinesterase from Electrophorus electricus were solved at 4.5 and 4.2 A resolution by molecular replacement. Resolution at these levels is sufficient to provide substantial information on the relative orientations of the subunits within the tetramer. The two structures, which show canonical homodimers of subunits assembled through four-helix bundles, reveal discrete geometries in the assembly of the dimers to form: (a) a loose, pseudo-square planar tetramer with antiparallel alignment of the two four-helix bundles and a large space in the center where the carboxyl-terminal sequences may be buried or (b) a compact, square nonplanar tetramer that may expose all four sequences on a single side. Comparison of these two structures points to significant conformational flexibility of the tetramer about the four-helix bundle axis and along the dimer-dimer interface. Hence, in solution, several conformational states of a flexible tetrameric arrangement of acetylcholinesterase catalytic subunits may exist to accommodate discrete carboxyl-terminal sequences of variable dimensions and amphipathicity.     

Computer simulations of the properties of the alpha2, alpha2C, and alpha2D de novo designed helical proteins

Reduced lattice models of the three de novo designed helical proteins alpha2, alpha2C, and alpha2D were studied. Low temperature stable folds were obtained for all three proteins. In all cases, the lowest energy folds were four-helix bundles. The folding pathway is qualitatively the same for all proteins studied. The energies of various topologies are similar, especially for the alpha2 polypeptide. The simulated crossover from molten globule to native-like behavior is very similar to that seen in experimental studies. Simulations on a reduced protein model reproduce most of the experimental properties of the alpha2, alpha2C, and alpha2D proteins. Stable four-helix bundle structures were obtained, with increasing native-like behavior on-going from alpha2 to alpha2D that mimics experiment.     

Design, synthesis and properties of synthetic chlorophyll proteins.

A chemoselective method is described for coupling chlorophyll derivatives with an aldehyde group to synthetic peptides or proteins modified with an aminoxyacetyl group at the epsilon-amino group of a lysine residue. Three template-assembled antiparallel four-helix bundles were synthesized for the ligation of one or two chlorophylls. This was achieved by coupling unprotected peptides to cysteine residues of a cyclic decapeptide by thioether formation. The amphiphilic helices were designed to form a hydrophobic pocket for the chlorophyll derivatives. Chlorophyll derivatives Zn-methyl-pheophorbide b and Zn-methyl-pyropheophorbide d were used. The aldehyde group of these chlorophyll derivatives was ligated to the modified lysine group to form an oxime bond. The peptide-chlorophyll conjugates were characterized by electrospray mass spectrometry, analytical HPLC, and UV/visible spectroscopy. Two four-helix bundle chlorophyll conjugates were further characterized by size-exclusion chromatography, circular dichroism, and resonance Raman spectroscopy.     

A heuristic approach to predicting the tertiary structure of bovine somatotropin

A combination of a heuristic approach and energy minimization was used to predict the three-dimensional structure of bovine somatotropin (bSt), also known as bovine growth hormone, a protein of 191 amino acids. The starting points for energy minimizations were generated from the following two types of inputs: (a) the amino acid sequence and (b) the heuristic inputs, which were derived according to physical, chemical, and biological principles by piecing together all useful information available. The predicted 3-D structure of the bSt molecule has all the features observed in four-helix bundle proteins. The four alpha-helices in bSt are intimately packed to form an assembly with an approximately square cross section. All the adjacent alpha-helices are antiparallel, with a somewhat tilted angle between each of the adjacent pairs so that the assembly of the four helices looks like a left-handed twisted bundle. There are two disulfide bonds in the bSt structure: one "hooking" the middle of a long loop with helix 4 so as to pull the long loop onto the surface of the helix bundle and the other "hooking" the C-terminal segment with the same helix so as to force the C-terminal segment to bend toward the helix bundle. As a consequence, a considerable part of the surface of the four-helix bundle is closely packed or intimately embraced by the loop segments. The predicted bSt structure has a hydrophobic core and a hydrophilic exterior surface. The energetic analysis of the predicted bSt structure indicates that the interaction between helices and loops plays a dominant role in stabilizing the four-helix bundle structure from the viewpoint of both electrostatic and nonbonded interactions. A technique called FOLD was meanwhile developed, by which one can fold a polypeptide chain into any shape as desired. This tool proved to be very useful during the heuristic model-building process.     

Stability of loops in the structure of lactose permease

Structural analysis of peptide fragments has provided useful information on the secondary structure of integral membrane proteins built from a helical bundle (up to seven transmembrane segments). Comparison of those results to recent X-ray crystallographic results showed agreement between the structures of the fragments and the structures of the intact proteins. Lactose permease of Escherichia coli (lac Y) offers an opportunity to test that hypothesis on a substantially larger integral membrane protein. Lac Y contains a bundle of 12 transmembrane segments connected by 11 loops. Eleven segments, each corresponding to one of the loops in this protein, were studied. Five of these segments form defined structures in solution as determined by multidimensional nuclear magnetic resonance. Four peptides form turns, and one peptide reveals the end of one of the transmembrane helices. These results suggest that some loops in helical bundles are stabilized by short-range interactions, particularly in smaller bundles, and such intrinsically stable loops may contribute to protein stability and influence the pathway of folding. Greater conformational flexibility may be found in large integral membrane proteins.     

New design of helix bundle peptide-polymer conjugates

We present a new design of peptide-polymer conjugates where a polymer chain is covalently linked to the side chain of a helix bundle-forming peptide. The effect of conjugated polymer chains on the peptide structure was examined using a de novo designed three-helix bundle and a photoactive four-helix bundle. Upon attachment of poly(ethylene glycol) to the exterior of the coiled-coil helix bundle, the peptide secondary structure was stabilized and the tertiary structure, that is, the coiled-coil helix bundle, was retained. When a heme-binding peptide as an example is used, the new peptide-polymer conjugate architecture also preserves the built-in functionalities within the interior of the helix bundle. It is expected that the conjugated polymer chains act to mediate the interactions between the helix bundle and its external environment. Thus, this new peptide-polymer conjugate design strategy may open new avenues to macroscopically assemble the helix bundles and may enable them to function in nonbiological environments.     

How cytochromes with different folds control heme redox potentials

The electrochemical midpoint potentials (E(m)'s) of 13 cytochromes, in globin (c, c(2), c(551), c(553)), four-helix bundle (c', b(562)), alpha beta roll (b(5)), and beta sandwich (f) motifs, with E(m)'s spanning 450 mV were calculated with multiconformation continuum electrostatics (MCCE). MCCE calculates changes in oxidation free energy when a heme-axial ligand complex is moved from water into protein. Calculated and experimental E(m)'s are in good agreement for cytochromes with His-Met and bis-His ligated hemes, where microperoxidases provide reference E(m)'s. In all cytochromes, E(m)'s are raised by 130-260 mV relative to solvated hemes by the loss of reaction field (solvation) energy. However, there is no correlation between E(m) and heme surface exposure. Backbone amide dipoles in loops or helix termini near the axial ligands raise E(m)'s, but amides in helix bundles contribute little. Heme propionates lower E(m)'s. If the propionic acids are partially protonated in the reduced cytochrome, protons are released on heme oxidation, contributing to the pH dependence of the E(m). In all cytochromes studied except b(5)'s and low potential globins, buried side chains raise E(m)'s. MCCE samples ionizable group protonation states, heme redox states, and side chain rotamers simultaneously. Globins show the largest structural changes on heme oxidation and four-helix bundles the least. Given the calculated protein-induced E(m) shift and measured cytochrome E(m) the five-coordinate, His heme in c' is predicted to have a solution E(m) between that of isolated bis-His and His-Met hemes, while the reference E(m) for His-Ntr ligands in cytochrome f should be near that of His-Met hemes.     

Protein refolding versus aggregation: computer simulations on an intermediate-resolution protein model.

Computer simulations are performed on a system of eight model peptide chains to study how the competition between protein refolding and aggregation affects the optimal conditions for refolding of four-helix bundles. The discontinuous molecular dynamics algorithm is utilized along with an intermediate-resolution protein model that we developed for this work. Physically, the model is much more detailed than any model used to date for simulations of protein aggregation. Each model residue consists of a detailed, three-bead backbone and a simplified, single-bead side-chain. Excluded volume, hydrogen bond, and hydrophobic interactions are modeled with discontinuous (i.e. hard-sphere and square-well) potentials. Simulations efficiently sample conformational space, and complete folding trajectories from random initial configurations to two four-helix bundles are possible within two days on a single processor workstation. Folding of the bundles follows two main pathways, one through a trimeric intermediate and the other through an intermediate with two dimers. The proportion of trajectories that follow each route is significantly different for the eight-peptide system in this work than in a previously studied four-peptide system, which yields one four-helix bundle, suggesting, as our previous simulations have, that protein folding properties are strongly influenced by the presence of other proteins. Folding of the bundles is optimal within a fixed temperature range, with the high-temperature boundary a function of the complexity of the protein (or oligomer) to be folded and the low-temperature boundary a function of the complexity of the protein's environment. Above the optimal temperature range for folding, the model chains tend to unfold; below the optimal range, the model chains tend to aggregate. As has been seen previously, aggregates have substantial levels of native secondary structure, suggesting that aggregates are composed largely of partially folded intermediates, not denatured chains.     

The serine-rich domain from Crk-associated substrate (p130cas) is a four-helix bundle

p130(cas) (Crk-associated substrate) is a docking protein that is involved in assembly of focal adhesions and concomitant cellular signaling. It plays a role in physiological regulation of cell adhesion, migration, survival, and proliferation, as well as in oncogenic transformation. The molecule consists of multiple protein-protein interaction motifs, including a serine-rich region that is positioned between Crk and Src-binding sites. This study reports the first structure of a functional domain of Cas. The solution structure of the serine-rich region has been determined by NMR spectroscopy, demonstrating that this is a stable domain that folds as a four-helix bundle, a protein-interaction motif. The serine-rich region bears strong structural similarity to four-helix bundles found in other adhesion components like focal adhesion kinase, alpha-catenin, or vinculin. Potential sites for phosphorylation and interaction with the 14-3-3 family of cellular regulators are identified in the domain and characterized by site-directed mutagenesis and binding assays. Mapping the degree of amino acid conservation onto the molecular surface reveals a patch of invariant residues near the C terminus of the bundle, which may represent a previously unidentified site for protein interaction.     

Pierced Lasso Bundles Are a New Class of Knot-like Motifs

A four-helix bundle is a well-characterized motif often used as a target for designed pharmaceutical therapeutics and nutritional supplements. Recently, we discovered a new structural complexity within this motif created by a disulphide bridge in the long-chain helical bundle cytokine leptin. When oxidized, leptin contains a disulphide bridge creating a covalent-loop through which part of the polypeptide chain is threaded (as seen in knotted proteins). We explored whether other proteins contain a similar intriguing knot-like structure as in leptin and discovered 11 structurally homologous proteins in the PDB. We call this new helical family class the Pierced Lasso Bundle (PLB) and the knot-like threaded structural motif a Pierced Lasso (PL). In the current study, we use structure-based simulation to investigate the threading/folding mechanisms for all the PLBs along with three unthreaded homologs as the covalent loop (or lasso) in leptin is important in folding dynamics and activity. We find that the presence of a small covalent loop leads to a mechanism where structural elements slipknot to thread through the covalent loop. Larger loops use a piercing mechanism where the free terminal plugs through the covalent loop. Remarkably, the position of the loop as well as its size influences the native state dynamics, which can impact receptor binding and biological activity. This previously unrecognized complexity of knot-like proteins within the helical bundle family comprises a completely new class within the knot family, and the hidden complexity we unraveled in the PLBs is expected to be found in other protein structures outside the four-helix bundles. The insights gained here provide critical new elements for future investigation of this emerging class of proteins, where function and the energetic landscape can be controlled by hidden topology, and should be take into account in ab initio predictions of newly identified protein targets.     

Central Region of Talin Has a Unique Fold That Binds Vinculin and Actin

Talin is an adaptor protein that couples integrins to F-actin. Structural studies show that the N-terminal talin head contains an atypical FERM domain, whereas the N- and C-terminal parts of the talin rod include a series of α-helical bundles. However, determining the structure of the central part of the rod has proved problematic. Residues 1359–1659 are homologous to the MESDc1 gene product, and we therefore expressed this region of talin in Escherichia coli. The crystal structure shows a unique fold comprised of a 5- and 4-helix bundle. The 5-helix bundle is composed of nonsequential helices due to insertion of the 4-helix bundle into the loop at the C terminus of helix α3. The linker connecting the bundles forms a two-stranded anti-parallel β-sheet likely limiting the relative movement of the two bundles. Because the 5-helix bundle contains the N and C termini of this module, we propose that it is linked by short loops to adjacent bundles, whereas the 4-helix bundle protrudes from the rod. This suggests the 4-helix bundle has a unique role, and its pI (7.8) is higher than other rod domains. Both helical bundles contain vinculin-binding sites but that in the isolated 5-helix bundle is cryptic, whereas that in the isolated 4-helix bundle is constitutively active. In contrast, both bundles are required for actin binding. Finally, we show that the MESDc1 protein, which is predicted to have a similar fold, is a novel actin-binding protein.     

Energetic approach to the folding of four alpha-helices connected sequentially

The packing of four alpha-helices, which each consist of 12 Ala residues and are sequentially connected to each other by a segment of 10 Ala residues, has been investigated by means of energy minimizations. For the lowest energy structure thus obtained, the following features have been found: (i) the four alpha-helices are intimately packed to form an assembly with an approximately square section; (ii) the distances of closest approach between two adjacent interhelix axes are 7.7 +/- 0.2 A and those between two diagonal interhelix axes are 11.2 +/- 0.2 A; (iii) the adjacent interhelix angles are -163 +/- 2 degrees; and (iv) the diagonal interhelix angles are 24 +/- 4 degrees. These results indicate that the polypeptide chain, driven by energetics (nonbonded and electrostatic interactions), is folded into a typical left-handed twisted four-helix bundle with an approximately 4-fold symmetric array, as observed in most four alpha-helix proteins. Furthermore, it has been found that the interaction between the loops formed by the connecting segments and the other part of molecule plays a significant role in stabilizing such a bundle structure. The technology developed here and the relevant knowledge obtained through this study are very useful for the study of modeling four-helix bundle proteins.