The construction of new proteins: V. A template-assembled synthetic protein (TASP) containing both a 4-helix bundle and beta-barrel-like structure

The construction of a template-assembled synthetic protein (TASP) designed to contain both a 4-helix bundle and a beta-barrel as two folding "domains" is described. For the de novo design of proteins, amphiphilic helices (alpha) and beta-sheets (beta) are covalently attached to a template peptide (T) carrying functional side chains suitably oriented to promote intramolecular folding of the secondary structure blocks into a characteristic packing arrangement, i.e., T8-(4 alpha)(4 beta). The design of this new macromolecule was assisted by computer modeling, which suggested a low-energy conformation with tight hydrophobic packing of the secondary structure subunits. Solid-phase synthesis of the "two-domain" TASP molecule was achieved using orthogonal protection techniques. The solution properties as well as circular dichroism (CD) and infrared spectroscopy (IR) data under various experimental conditions are consistent with the folded conformation suggested by modeling.     

Synthesis of antiparallel 4α-helix bundle TASP by chemoselective ligation

Summary The use of chemoselective ligation methods and orthogonal protection techniques allows access to Template-Assembled Synthetic Protein (TASP) molecules exhibiting a large variety of packing topologies. This is demonstrated for the synthesis of an antiparallel 4-helical bundle TASP by condensing amphiphilic peptide blocks, containing aldehyde functions at the C- or N-terminus, to a selectively addressable topological template via oxime bond formation. The resulting antiparallel 4-helix TASP is obtained in high yield and shows a template-induced helical conformation.     

Total chemical synthesis, characterization, and immunological properties of an MHC class I model using the TASP concept for protein de novo design.

The design, total chemical synthesis, and immunological properties of a four-alpha-helix bundle template-assembled synthetic protein (TASP) mimicking some of the structural features of the major histocompatibility complex (MHC) class I is described. In a first approach, the native sequence 58-74 of the alpha 1 heavy chain domain of HLA-A2 was modeled in order to increase helix stability and amphiphilicity of the 17-mer peptide, preserving the residues for potential T-cell receptor (TcR) binding properties. According to the TASP concept, these helical segments were covalently attached to a cyclic template molecule designed for the induction of a four-helix-bundle topology of the assembled peptide blocks. After extensive HPLC purification, stepwise solid-phase synthesis resulted in a TASP molecule of high chemical purity as demonstrated by analytical HPLC, mass spectrometry, and amino acid analysis. CD spectroscopic investigations are consistent with the onset of a partial alpha-helical conformation in aqueous buffer as well as in TFE. Antibodies raised directly against this four-alpha-helix bundle TASP molecule (without prior conjugation to a carrier molecule) were detected by ELISA. Flow cytometry studies showed that these antibodies recognize the native MHC class I molecule on the surface of HLA-A2-positive cells. The results indicate that the TASP approach represents a versatile tool for mimicking conformational epitopes.     

Characterization of the backbone geometry of protein native state structures

We present a simple method for analyzing the geometry of noncovalent residue-residue interactions stabilizing the protein structure, which takes into account the constraints on the local backbone geometry. We find that the principal geometrical constraints are amino acid aspecific and are associated with hydrogen bond formation in helices and sheets. In contrast, amino acid residues in nonhelical and nonextended conformations, which make noncovalent interactions stabilizing the protein tertiary structure, display greater flexibility. We apply the method to an analysis of the packing of helices in helical bundle proteins requiring an efficient packing of amino acid side-chains of the interacting helices.     

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.     

Electrostatic stabilization in four-helix bundle proteins.

Charge substitutions generated by site-directed mutagenesis at the termini of adjacent anti-parallel alpha-helices in a four-helix bundle protein were used to determine a precise value for the contribution of indirect charge-charge interactions to overall protein stability, and to simulate the electrostatic effects of alpha-helix macrodipoles. Thermodynamic double mutant cycles were constructed to measure the interaction energy between such charges on adjacent anti-parallel helices in the four-helix bundle cytochrome b562 from Escherichia coli. Previously, theoretical calculations of helix macrodipole interactions using modeled four-helix bundle proteins have predicted values ranging over an order of magnitude from 0.2 to 2.5 kcal/mol. Our system represents the first experimental evidence for electrostatic interactions such as those between partial charges due to helix macrodipole charges. At the positions mutated, we have measured a favorable interaction energy of 0.6 kcal/mol between opposite charges simulating an anti-parallel helix pair. Pairs of negative or positive charges simulating a parallel orientation of helices produce an unfavorable interaction of similar magnitude. The interaction energies show a strong dependence upon ionic strength, consistent with an electrostatic effect. Indirect electrostatic contacts do appear to confer a limited stabilization upon the association of anti-parallel packing of helices, favoring this orientation by as much as 1 kcal/mol at 20 mM K phosphate.     

X-ray structure of a maquette scaffold

Maquettes are de novo designed mimicries of nature used to test the construction and engineering criteria of oxidoreductases. One type of scaffold used in maquette construction is a four-alpha-helical bundle. The sequence of the four-alpha-helix bundle maquettes follows a heptad repeat pattern typical of left-handed coiled-coils. Initial designs were molten globular due partly to the minimalist approach taken by the designers. Subsequent iterative redesign generated several structured scaffolds with similar heme binding properties. Variant [I(6)F(13)](2), a structured scaffold, was partially resolved with NMR spectroscopy and found to have a set of mobile inter-helical packing interfaces. Here, the X-ray structure of a similar peptide ([I(6)F(13)M(31)](2) i.e. ([CGGG EIWKL HEEFLKK FEELLKL HEERLKKM](2))(2) which we call L31M), has been solved using MAD phasing and refined to 2.8A resolution. The structure shows that the maquette scaffold is an anti-parallel four-helix bundle with "up-up-down-down" topology. No pre-formed heme-binding pocket exists in the protein scaffold. We report unexpected inter-helical crossing angles, residue positions and translations between the helices. The crossing angles between the parallel helices are -5 degrees rather than the expected +20 degrees for typical left-handed coiled-coils. Deviation of the scaffold from the design is likely due to the distribution and size of hydrophobic residues. The structure of L31M points out that four identical helices may interact differently in a bundle and heptad repeats with an alternating [HPPHHPP]/[HPPHHPH] (H: hydrophobic, P: polar) pattern are not a sufficient design criterion to generate left-hand coiled-coils.     

Template-assembled melittin: structural and functional characterization of a designed, synthetic channel-forming protein.

Template-assembled proteins (TASPs) comprising 4 peptide blocks, each of either the natural melittin sequence (melittin-TASP) or of a truncated melittin sequence (amino acids 6-26, melittin6-26-TASP), C-terminally linked to a (linear or cyclic) 10-amino acid template were synthesized and characterized, structurally by CD, by fluorescence spectroscopy, and by monolayer experiments, and functionally, by electrical conductance measurements on planar bilayers and release experiments on dye-loaded vesicles. Melittin-TASP and the truncated analogue preferentially adopt alpha-helical structures in methanol (56% and 52%, respectively) as in lipid membranes. Unlike in methanol, the melittin-TASP self-aggregates in water. On an air-water interface, the differently sized molecules can be self-assembled and compressed to a compact structure with a molecular area of around 600 A2, compatible with a 4-helix bundle preferentially oriented perpendicular to the interface. The proteins reveal a strong affinity for lipid membranes. A partition coefficient of 1.5 x 10(9) M-1 was evaluated from changes of the Trp fluorescence spectra of the TASP in water and in the lipid bilayer. In planar lipid bilayers, TASP molecules are able to form defined ion channels, exhibiting a small single-channel conductance of 7 pS (in 1 M NaCl). With increasing protein concentration in the lipid bilayer, additional, larger conductance states of up to 1 nS were observed. These states are likely to be formed by aggregated TASP structures as inferred from a strongly voltage-dependent channel activity on membranes of large area. In this respect, melittin-TASP reveals channel features of the native peptide, but with a considerably lower variation in the size of the channel states. Compared to the free peptide, template-assembled melittin has a much higher membrane activity: it is about 100 times more effective in channel formation and 20 times more effective in releasing dye molecules from lipid vesicles. This demonstrates that the lytic properties are not solely related to channel formation.     

PackHelix: a tool for helix-sheet packing during protein structure prediction

The three‐dimensional structure of a protein is organized around the packing of its secondary structure elements. Predicting the topology and constructing the geometry of structural motifs involving α‐helices and/or β‐strands are therefore key steps for accurate prediction of protein structure. While many efforts have focused on how to pack helices and on how to sample exhaustively the topologies and geometries of multiple strands forming a β‐sheet in a protein, there has been little progress on generating native‐like packings of helices on sheets. We describe a method that can generate the packing of multiple helices on a given β‐sheet for αβα sandwich type protein folds. This method mines the results of a statistical analysis of the conformations of αβ₂ motifs in protein structures to provide input values for the geometric attributes of the packing of a helix on a sheet. It then proceeds with a geometric builder that generates multiple arrangements of the helices on the sheet of interest by sampling through these values and performing consistency checks that guarantee proper loop geometry between the helices and the strands, minimal number of collisions between the helices, and proper formation of a hydrophobic core. The method is implemented as a module of ProteinShop. Our results show that it produces structures that are within 4–6 Å RMSD of the native one, regardless of the number of helices that need to be packed, though this number may increase if the protein has several helices between two consecutive strands in the sequence that pack on the sheet formed by these two strands. Proteins 2011; Published 2011 Wiley‐Liss, Inc.     

Mutational analyses of HAMP helices suggest a dynamic bundle model of input–output signalling in chemoreceptors

To test the gearbox model of HAMP signalling in the Escherichia coli serine receptor, Tsr, we generated a series of amino acid replacements at each residue of the AS1 and AS2 helices. The residues most critical for Tsr function defined hydrophobic packing faces consistent with a four-helix bundle. Suppression patterns of helix lesions conformed to the predicted packing layers in the bundle. Although the properties and patterns of most AS1 and AS2 lesions were consistent with both proposed gearbox structures, some mutational features specifically indicate the functional importance of an x-da bundle over an alternative a-d bundle. These genetic data suggest that HAMP signalling could simply involve changes in the stability of its x-da bundle. We propose that Tsr HAMP controls output signals by modulating destabilizing phase clashes between the AS2 helices and the adjoining kinase control helices. Our model further proposes that chemoeffectors regulate HAMP bundle stability through a control cable connection between the transmembrane segments and AS1 helices. Attractant stimuli, which cause inward piston displacements in chemoreceptors, should reduce cable tension, thereby stabilizing the HAMP bundle. This study shows how transmembrane signalling and HAMP input–output control could occur without the helix rotations central to the gearbox model.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Energetics and stability of transmembrane helix packing: a replica-exchange simulation with a knowledge-based membrane potential

The energetics and stability of the packing of transmembrane helices were investigated by Monte Carlo simulations with the replica-exchange method. The helices were modeled with a united atom representation, and the CHARMM19 force field was employed. Based on known experimental structures of membrane proteins, an implicit knowledge-based potential was developed to describe the helix-membrane interactions at the residue level, whose validity was tested through prediction of the orientations when single helices were inserted into a membrane. Two systems were studied in this article, namely the glycophorin A dimer, and helices A and B of Bacteriorhodopsin. For the glycophorin A dimer, the most stable structure (0.5 A away from the experimental structure) is mainly stabilized by the favorable helix-helix interactions, and has the most population regardless of the helix-membrane interaction. However, for helices A and B of Bacteriorhodopsin, it was found that the packing determined by helix-helix interactions is nonspecific, and a native-like structure (0.2 A from the experimental one) can be identified from several structural analogs as the most stable one only after applying the membrane potential. Our results suggest that the contribution from the helix-membrane interaction could be critical in the correct packing of transmembrane helices in the membrane.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Different secondary structure elements as scaffolds for protein folding transition states of two homologous four-helix bundles

Comparison of the folding processes for homologue proteins can provide valuable information about details in the interactions leading to the formation of the folding transition state. Here the folding kinetics of 18 variants of yACBP and 3 variants of bACBP have been studied by Phi-value analysis. In combination with Phi-values from previous work, detailed insight into the transition states for folding of both yACBP and bACBP has been obtained. Of the 16 sequence positions that have been studied in both yACBP and bACBP, 5 (V12, I/L27, Y73, V77, and L80) have high Phi-values and appear to be important for the transition state formation in both homologues. Y31, A34, and A69 have high Phi-values only in yACBP, while F5, A9, and I74 have high Phi-values only in bACBP. Thus, additional interactions between helices A2 and A4 appear to be important for the transition state of yACBP, whereas additional interactions between helices A1 and A4 appear to be important for the transition state of bACBP. To examine whether these differences could be assigned to different packing of the residues in the native state, a solution structure of yACBP was determined by NMR. Small changes in the packing of the hydrophobic side-chains, which strengthen the interactions between helices A2 and A4, are observed in yACBP relative to bACBP. It is suggested that different structure elements serve as scaffolds for the folding of the 2 ACBP homologues.     

Packing of transmembrane helices in bacteriorhodopsin folding: structure and thermodynamics

We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.     

Role of stabilization centers in 4 helix bundle proteins

Stabilization centers (SCs) were shown to play an important role in preventing decay of three-dimensional protein structures. These residue clusters, stabilized by cooperative long-range interactions, were proposed to serve as anchoring points for arranging secondary structure elements. In all-alpha proteins, SC elements appear less frequently than in all-beta, alpha/beta, and alpha+beta proteins suggesting that tertiary structure formation of all-alpha proteins is governed by different principles than in other protein classes. Here we analyzed the relation between the formation of stabilization centers and the inter-axial angles (Omega) of alpha-helices in 4 helix bundle proteins. In the distance range, where dipoles have dominant effect on the helix pair arrangement, those helix pairs, where residues from both helices participate in SC elements, appear as parallel more frequently than those helices where no SC elements are present. For SC containing helix pairs, the energetic difference between the parallel and anti-parallel states decreases considerably from 1.1 kcal/mol to 0.4 kcal/mol. Although the observed effect is weak for more distant helices, a competition between the SC element formation and the optimal dipole-dipole interaction of alpha-helices is proposed as a mechanism for tertiary structure formation in 4 helix bundle proteins. The SC-forming potential of different arrangements as well as the pitfalls of the SC definition are also discussed.     

Crystal structure of a designed, thermostable, heterotrimeric coiled coil

Electrostatic interactions are often critical for determining the specificity of protein-protein complexes. To study the role of electrostatic interactions for assembly of helical bundles, we previously designed a thermostable, heterotrimeric coiled coil, ABC, in which charged residues were employed to drive preferential association of three distinct, 34-residue helices. To investigate the basis for heterotrimer specificity, we have used multiwavelength anomalous diffraction (MAD) analysis to determine the 1.8 A resolution crystal structure of ABC. The structure shows that ABC forms a heterotrimeric coiled coil with the intended arrangement of parallel chains. Over half of the ion pairs engineered to restrict helix associations were apparent in the experimental electron density map. As seen in other trimeric coiled coils, ABC displays acute knobs-into-holes packing and a buried anion coordinated by core polar amino acids. These interactions validate the design strategy and illustrate how packing and polar contacts determine structural uniqueness.