Protein design, a minimalist approach

How do helix-helix interactions help determine the folds of membrane proteins? Perspectives from the study of homo-oligomeric helical bundles

The final, structure-determining step in the folding of membrane proteins involves the coalescence of preformed transmembrane helices to form the native tertiary structure. Here, we review recent studies on small peptide and protein systems that are providing quantitative data on the interactions that drive this process. Gel electrophoresis, analytical ultracentrifugation, and fluorescence resonance energy transfer (FRET) are useful methods for examining the assembly of homo-oligomeric transmembrane helical proteins. These methods have been used to study the assembly of the M2 proton channel from influenza A virus, glycophorin, phospholamban, and several designed membrane proteins-all of which have a single transmembrane helix that is sufficient for association into a transmembrane helical bundle. These systems are being studied to determine the relative thermodynamic contributions of van der Waals interactions, conformational entropy, and polar interactions in the stabilization of membrane proteins. Although the database of thermodynamic information is not yet large, a few generalities are beginning to emerge concerning the energetic differences between membrane and water-soluble proteins: the packing of apolar side chains in the interior of helical membrane proteins plays a smaller, but nevertheless significant, role in stabilizing their structure. Polar, hydrogen-bonded interactions occur less frequently, but, nevertheless, they often provide a strong driving force for folding helix-helix pairs in membrane proteins. These studies are laying the groundwork for the design of sequence motifs that dictate the association of membrane helices.

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.

Computational de novo design, and characterization of an A(2)B(2) diiron protein

Diiron proteins are found throughout nature and have a diverse range of functions; proteins in this class include methane monooxygenase, ribonucleotide reductase, Delta(9)-acyl carrier protein desaturase, rubrerythrin, hemerythrin, and the ferritins. Although each of these proteins has a very different overall fold, in every case the diiron active site is situated within a four-helix bundle. Additionally, nearly all of these proteins have a conserved Glu-Xxx-Xxx-His motif on two of the four helices with the Glu and His residues ligating the iron atoms. Intriguingly, subtle differences in the active site can result in a wide variety of functions. To probe the structural basis for this diversity, we designed an A(2)B(2) heterotetrameric four-helix bundle with an active site similar to those found in the naturally occurring diiron proteins. A novel computational approach was developed for the design, which considers the energy of not only the desired fold but also alternatively folded structures. Circular dichroism spectroscopy, analytical ultracentrifugation, and thermal unfolding studies indicate that the A and B peptides specifically associate to form an A(2)B(2) heterotetramer. Further, the protein binds Zn(II) and Co(II) in the expected manner and shows ferroxidase activity under single turnover conditions.

Towards the design and computational characterization of a membrane protein

The design of a transmembrane four-helix bundle is described. We start with an idealized four-helix bundle geometry, then use statistical information to build a plausible transmembrane bundle. Appropriate residues are chosen using database knowledge on the sequences of membrane helices and loops, then the packing of the bundle core is optimized, and favorable side chain rotamers from rotamer libraries are selected. Next, we use explicit physical knowledge from biomolecular simulation force fields and molecular dynamics simulations to test whether the designed structure is physically possible. These procedures test whether the designed protein will indeed be alpha-helical, well packed and stable over a time scale of several nanoseconds in a realistic lipid bilayer environment. We then test a modeling approach that does not include sophisticated database knowledge about proteins, but rather relies on applying our knowledge of the physics that governs protein motions. This independent validation of the design is based on simulated annealing and restrained molecular dynamics simulation in vacuo, comparable to procedures used to refine NMR and X-ray structures.

Increasing protein stability using a rational approach combining sequence homology and structural alignment: Stabilizing the WW domain

This study shows that a combination of sequence homology and structural information can be used to increase the stability of the WW domain by 2.5 kcal mol(-1) and increase the T(m) by 28 degrees C. Previous homology-based protein design efforts typically investigate positions with low sequence identity, whereas this study focuses on semi-conserved core residues and proximal residues, exploring their role(s) in mediating stabilizing interactions on the basis of structural considerations. The A20R and L30Y mutations allow increased hydrophobic interactions because of complimentary surfaces and an electrostatic interaction with a third residue adjacent to the ligand-binding hydrophobic cluster, increasing stability significantly beyond what additivity would predict for the single mutations. The D34T mutation situated in a pi-turn possibly disengages Asn31, allowing it to make up to three hydrogen bonds with the backbone in strand 1 and loop 2. The synergistic mutations A20R/L30Y in combination with the remotely located mutation D34T add together to create a hYap WW domain that is significantly more stable than any of the protein structures on which the design was based (Pin and FBP28 WW domains).

Design, synthesis, and characterization of a novel hemoprotein

Here we describe a synthetic protein (6H7H) designed to bind four heme groups via bis-histidine axial ligation. The hemes are designed to bind perpendicular to another in an orientation that mimics the relative geometry of the two heme a groups in the active site of cytochrome c oxidase. Our newly developed protein-design program, called CORE, was implemented in the design of this novel hemoprotein. Heme titration studies resolved four distinct K(D) values (K(D1) = 80 nM, K(D2) = 18 nM, K(D3) > or = 3 mM, K(D4) < or = 570 nM, with K(D3) x K(D4) = 1700); positive cooperativity in binding between the first and second heme, as well as substantial positive cooperativity between the third and forth heme, was observed. Chemical and thermal denaturation studies reveal a stable protein with native-like properties. Visible circular dichroism spectroscopy of holo-6H7H indicates excitonic coupling between heme groups. Further electrochemical and spectroscopic characterization of the holo-protein support a structure that is consistent with the predefined target structure.

The effect of Zn(2+) on the secondary structure of a histidine-rich fusogenic peptide and its interaction with lipid membranes

Membrane fusion between uncharged lipid vesicles can be triggered by the peptide sequence 'B18' from the fertilization protein 'bindin', but it only proceeds efficiently in the presence of Zn(2+) ions. We studied (i) the interaction of Zn(2+) with the fusogenic peptide B18, (ii) the binding of B18 to 1-palmitoyl-2-oleoylglycero-3-phosphocholine (POPC), and (iii) the ternary system POPC/B18/Zn(2+). The complex formation of Zn(2+) with the central histidine-rich motif of B18 appears to shift the secondary structure away from a beta-sheet towards an alpha-helical conformation. Here we observe for the first time an essentially alpha-helical structure of the peptide when immersed in POPC bilayers which appears to represent its functional fusogenic state. Infrared linear dichroism suggests a peripheral, oblique insertion mode of B18, mediated by the hydrophobic patches along one side of the amphipathic peptide. Furthermore, the hydration level of the peptide is reduced, suggesting that the hydrophobic region of the bilayer is involved in the lipid/peptide interactions. The hydration capacity of the POPC/B18/Zn(2+) system is distinctly smaller than that of POPC/Zn(2+) without peptide. The accompanying decrease in the number of tightly bound water molecules per lipid can be interpreted as a reduction in the repulsive 'hydration' forces, which usually prevent the spontaneous fusion of lipid vesicles. Binding of the B18 peptide in the presence of Zn(2+) effectively renders the membrane surface more hydrophobic, thus allowing fusion to proceed.

Automatic protein design with all atom force-fields by exact and heuristic optimization

A fully automatic procedure for predicting the amino acid sequences compatible with a given target structure is described. It is based on the CHARMM package, and uses an all atom force-field and rotamer libraries to describe and evaluate side-chain types and conformations. Sequences are ranked by a quantity akin to the free energy of folding, which incorporates hydration effects. Exact (Branch and Bound) and heuristic optimisation procedures are used to identifying highly scoring sequences from an astronomical number of possibilities. These sequences include the minimum free energy sequence, as well as all amino acid sequences whose free energy lies within a specified window from the minimum. Several applications of our procedure are illustrated. Prediction of side-chain conformations for a set of ten proteins yields results comparable to those of established side-chain placement programs. Applications to sequence optimisation comprise the re-design of the protein cores of c-Crk SH3 domain, the B1 domain of protein G and Ubiquitin, and of surface residues of the SH3 domain. In all calculations, no restrictions are imposed on the amino acid composition and identical parameter settings are used for core and surface residues. The best scoring sequences for the protein cores are virtually identical to wild-type. They feature no more than one to three mutations in a total of 11-16 variable positions. Tests suggest that this is due to the balance between various contributions in the force-field rather than to overwhelming influence from packing constraints. The effectiveness of our force-field is further supported by the sequence predictions for surface residues of the SH3 domain. More mutations are predicted than in the core, seemingly in order to optimise the network of complementary interactions between polar and charged groups. This appears to be an important energetic requirement in absence of the partner molecules with which the SH3 domain interacts, which were not included in the calculations. Finally, a detailed comparison between the sequences generated by the heuristic and exact optimisation algorithms, commends a note of caution concerning the efficiency of heuristic procedures in exploring sequence space.

De novo design and structural characterization of proteins and metalloproteins

De novo protein design has recently emerged as an attractive approach for studying the structure and function of proteins. This approach critically tests our understanding of the principles of protein folding; only in de novo design must one truly confront the issue of how to specify a protein's fold and function. If we truly understand proteins, it should be possible to design receptors, enzymes, and ion channels from scratch. Further, as this understanding evolves and is further refined, it should be possible to design proteins and biomimetic polymers with properties unprecedented in nature.

SymROP: ROP protein with identical helices redesigned by all-atom contact analysis and molecular dynamics

Experience has shown that protein redesigns (using the backbone from a known protein structure) are far more likely to produce well-ordered, native-like structures than are true de novo designs. Therefore, to design a four-helix bundle made of identical short helices, we here proceed by an extensive redesign of the ROP protein. A fully symmetrical SymROP sequence derived from ROP was chosen by modeling ideal-geometry side chains, including hydrogens, while maintaining the "goodness-of-fit" of side-chain packing by calculating all-atom contact surfaces with the Reduce and Probe programs. To estimate the probable extent of backbone movement and side-chain mobility, restrained molecular dynamics simulations were compared for candidate sequences and controls, including substitution of Abu for all or half the core Ala residues. The resulting 17-residue designed sequence is 41% identical to the relevant regions in ROP. SymROP is intended for construction by the Template Assembled Synthetic Proteins approach, to control the bundle topology, to use short helices, and to allow blocked termini and unnatural amino acids. ROP protein has been a valuable system for studying helical protein structure because of its simplicity and regularity within a structure large enough to have a real hydrophobic core. The SymROP design carries that simplicity and regularity even further.

Synthesis and NMR solution structure of an alpha-helical hairpin stapled with two disulfide bridges

Helical coiled-coils and bundles are some of the most common structural motifs found in proteins. Design and synthesis of alpha-helical motifs may provide interesting scaffolds that can be useful as host structures to display functional sites, thus allowing the engineering of novel functional miniproteins. We have synthesized a 38-amino acid peptide, alpha2p8, encompassing the alpha-helical hairpin present in the structure of p8MTCP1, as an alpha-helical scaffold particularly promising for its stability and permissiveness of sequence mutations. The three-dimensional structure of this peptide has been solved using homonuclear two-dimensional NMR techniques at 600 MHz. After sequence specific assignment, a total of 285 distance and 29 dihedral restraints were collected. The solution structure of alpha2p8 is presented as a set of 30 DIANA structures, further refined by restrained molecular dynamics, using simulated annealing protocol with the AMBER force field. The RMSD values for the backbone and all heavy atoms are 0.65+/-0.25 and 1.51+/-0.21 A, respectively. Excised from its protein context, the alpha-hairpin keeps its native structure: an alpha-helical coiled-coil, similar to that found in superhelical structures, with two helices spanning residues 4-16 and 25-36, and linked by a short loop. This motif is stabilized by two interhelical disulfide bridges and several hydrophobic interactions at the helix interface, leaving most of its solvent-exposed surface available for mutation. This alpha-helical hairpin, easily amenable to synthetic chemistry and biological expression system, may represent a stable and versatile scaffold to display new functional sites and peptide libraries.

Synthesis and study of C(3)-symmetric hydropyran cyclooligolides with oriented aryl and alcohol appendages at 10 A spacing

Modular syntheses of C(3)-symmetric macrocycles with pendant aryl and hydroxymethyl groups are described. These functional groups, amenable to further elaboration, were installed early in each synthesis and carried through an iterative sequence of module coupling and macrolactonization. Association constants for macrolides 1a-c with alkali metal cation guests were determined, and sandwich-type complexes with Ba(2+) were confirmed for these macrocycles based on (1)H NMR studies, including Job plots. X-ray crystallographic data for macrolides 1a and 1c were obtained and are discussed in detail. These data provide support that the macrolides are structurally well-defined and preorganized for binding the potassium cation. Preparation of the tris(bromoacetylated) macrotriolide 43 exemplifies a functionalized platform suitable for elaboration with peptide or carbohydrate residues.

Diabetes-associated mutations in a beta-cell transcription factor destabilize an antiparallel "mini-zipper" in a dimerization interface

Maturity-onset diabetes of the young, a monogenic form of Type II diabetes mellitus, is most commonly caused by mutations in hepatic nuclear factor 1alpha (HNF-1alpha). Here, the dimerization motif of HNF-1alpha is shown to form an intermolecular four-helix bundle. One face contains an antiparallel coiled coil whereas the other contains splayed alpha-helices. The "mini-zipper" is complementary in structure and symmetry to the top surface of a transcriptional coactivator (dimerization cofactor of homeodomains). The bundle is destabilized by a subset of mutations associated with maturity-onset diabetes of the young. Impaired dimerization of a beta-cell transcription factor thus provides a molecular mechanism of metabolic deregulation in diabetes mellitus.

A new approach to the design of uniquely folded thermally stable proteins

A new computer program (CORE) is described that predicts core hydrophobic sequences of predetermined target protein structures. A novel scoring function is employed, which for the first time incorporates parameters directly correlated to free energies of unfolding (deltaGu), melting temperatures (Tm), and cooperativity. Metropolis-driven simulated annealing and low-temperature Monte Carlo sampling are used to optimize this score, generating sequences predicted to yield uniquely folded, stable proteins with cooperative unfolding transitions. The hydrophobic core residues of four natural proteins were predicted using CORE with the backbone structure and solvent exposed residues as input. In the two smaller proteins tested (Gbeta1, 11 core amino acids; 434 cro, 10 core amino acids), the native sequence was regenerated as well as the sequence of known thermally stable variants that exhibit cooperative denaturation transitions. Previously designed sequences of variants with lower thermal stability and weaker cooperativity were not predicted. In the two larger proteins tested (myoglobin, 32 core amino acids; methionine aminopeptidase, 63 core amino acids), sequences with corresponding side-chain conformations remarkably similar to that of native were predicted.

Folding minimal sequences: the lower bound for sequence complexity of globular proteins

Alphabet size and informational entropy, two formal measures of sequence complexity, are herein applied to two prior studies on the folding of minimal proteins. These measures show a designed four-helix bundle to be unlike its natural counterparts but rather more like a coiled-coil dimer. Segments from a simplified sarc homology 3 domain and more than 2000000 segments from globular proteins both have lower bounds for alphabet size of 10 and for entropy near 2.9. These values are therefore suggested to be necessary and sufficient for folding into globular proteins having both rigid side chain packing and biological function.

The roles of side chain and backbone in protein structure probed with glycine- and sarcosine-rich synthetic leucine zipper peptides

The protein folding problem has long been a formidable challenge. Here we present a synthetic natural motif approach that exploits small preexisting structural models for the dissection of forces important in protein folding. An example for this approach is shown in the modification of a 31-residue leucine zipper peptide with the helix-breaking amino acid glycine and the hydrogen bond-breaking imino acid sarcosine. Circular dichroism and NMR experiments have shown that the glycine-modified leucine zipper peptide adopts a stable helical conformation similar to the native conformation while the sarcosine-modified leucine zipper peptide adopts a random coil conformation. These results provide valuable insight into the current controversy over the relative importance of long-range side chain-side chain interactions versus local backbone interactions in protein structure and suggest that the natural motif strategy may represent a useful model to study protein folding.

Molecular modeling and design of regioselectively addressable functionalized templates with rigidified three-dimensional structures

Extensive conformational analysis of a series of beta-alkyl substituted cyclopeptides-cyclo(Pro(1)-Xaa(2)-Nle(3)-Ala(4)-Nle(5)-Pro(6)-Xaa(7)- Nle(8)-Ala(9)-Nle(10)) and cyclo[Pro(1)-Xaa(2)-Nle(3)-(Cys(4)- Nle(5)-Pro(6)-Xaa(7)-Nle(8)-Cys(9))-Nle(10)] as well as their corresponding unsubstituted core structures cyclo(Pro(1)-Xaa(2)-Ala(3)-Ala(4)-Ala(5)-Pro(6)-Xaa(7)-Ala(8)-Ala(9)- Ala(10)) and cyclo(Pro(1)-Xaa(2)-Ala(3)-Cys(4)- Ala(5)-Pro(6)-Xaa(7)-Ala(8)-Cys(9)-Ala(10)) has been performed employing both the ECEPP/2 and the MAB force fields (Xaa = Gly, L-Ala, D-Ala, Aib, and D-Pro). Results show that (a) possible three-dimensional structures of the cyclo(Pro(1)-Gly(2)-Lys(3)-Ala(4)-Lys(5)-Pro(6)-Gly(7)-Lys(8)-Ala(9)- Lys(10)) molecule are not limited to a single extended "rectangular" conformation with all Lys side chains oriented at the same side of the molecule; (b) conformational equilibrium in monocyclic analogues obtained by replacements of conformationally flexible Gly residues for L-Ala, D-Ala, Aib, or D-Pro is not significantly shifted towards the target "rectangular" conformational type; and (c) introduction of disulfide bridges between positions 4 and 9 is a very powerful way to stabilize the target conformations in the resulting bicyclic molecules. These findings form the basis for further design of rigidified regioselectively addressable functionalized templates with many application areas ranging from biostructural to diagnostic purposes.

Miniaturized proteins: the backbone cyclic proteinomimetic approach

The field of proteinomimetics utilizes peptide-based molecules to mimic native protein functions. We describe a novel general method for mimicking proteins by small cyclic peptides for the purpose of drug design, and demonstrate its applicability on bovine pancreatic trypsin inhibitor (BPTI). These unique cyclic peptides, which both embody discontinuous residues of proteins in their bio-active conformation and ensure an induced fit, may overcome some of the pharmacological drawbacks attributed to proteins and peptides. This method, which we call the backbone cyclic (BC) proteinomimetic approach, combines backbone cyclization of peptides with a suitable selection method, cycloscan. Following this procedure, we have prepared a bicyclic nonapeptide, which mimics the binding region of BPTI. The X-ray crystal structure of the complex trypsin:mimetic, as well as kinetic studies, show that the BPTI mimetic binds to the specificity pocket of trypsin in a similar manner to BPTI. Inhibition measurements of various constructs revealed that backbone cyclization imposed the conformation crucial to binding.

Extending the concept of template-assembled synthetic proteins

The creation of native-like macromolecules in copying nature's way represents a fascinating challenge in protein chemistry today. In the absence of a detailed knowledge of the complex folding pathway the ultimate goal in protein de novo design, the construction of artificial proteins with predetermined three-dimensional structure and tailor-made functions based on a defined, generally valid set of rules, appears to be still out of reach. With progress in synthesis strategies and biostructural characterization methods, topological templates have become a versatile tool for inducing and stabilizing secondary and tertiary structures, such as protein loops, beta-turns, alpha-helices, beta-sheets and a variety of folding motifs. In this article, we extend the concept of template-assembled synthetic proteins for the construction of protein-like topologies with multiply bridged, oligocyclic chain architectures termed locked-in tertiary folds that exhibit unique physicochemical and folding properties because of the highly confined conformational space. Furthermore, we show that some fundamental questions in protein assembly can be approached applying the template concept. Using covalent template trapping of self-associated peptide assemblies in aqueous solution the structural and physical forces guiding protein folding, supramolecular assembly and molecular recognition processes can be studied on a molecular level.

Surface binding of alamethicin stabilizes its helical structure: molecular dynamics simulations

Alamethicin is an amphipathic alpha-helical peptide that forms ion channels. An early event in channel formation is believed to be the binding of alamethicin to the surface of a lipid bilayer. Molecular dynamics simulations are used to compare the structural and dynamic properties of alamethicin in water and alamethicin bound to the surface of a phosphatidylcholine bilayer. The bilayer surface simulation corresponded to a loosely bound alamethicin molecule that interacted with lipid headgroups but did not penetrate the hydrophobic core of the bilayer. Both simulations started with the peptide molecule in an alpha-helical conformation and lasted 2 ns. In water, the helix started to unfold after approximately 300 ps and by the end of the simulation only the N-terminal region of the peptide remained alpha-helical and the molecule had collapsed into a more compact form. At the surface of the bilayer, loss of helicity was restricted to the C-terminal third of the molecule and the rod-shaped structure of the peptide was retained. In the surface simulation about 10% of the peptide/water H-bonds were replaced by peptide/lipid H-bonds. These simulations suggest that some degree of stabilization of an amphipathic alpha-helix occurs at a bilayer surface even without interactions between hydrophobic side chains and the acyl chain core of the bilayer.

Characterization of MB-1. A dimeric helical protein with a compact core

MB-1 is a de-novo protein designed to incorporate a large number of the nutritionally important amino acids methionine, lysine, leucine and threonine into a stable four-helix bundle protein. MB-1 has been expressed and purified from Escherichia coli, indicating it was resistant to intracellular proteases [Beauregard, M., Dupont, C., Teather, R.M. & Hefford, M.A. (1995) Bio/Technology 13, 974]. Here we report an analysis of the secondary, tertiary and quaternary structures in MB-1 using circular dichroism, fluorospectroscopy and size-exclusion chromatography. Our data indicate that the MB-1 structure is close to the target structure, an alpha-helical bundle, in many respects and is highly helical in solution. The single tyrosine incorporated into the designed protein as a spectrocopic probe of tertiary structure, is buried in a compact, folded core and becomes accessible on protein denaturation, as per design. Furthermore, MB-1 was found to be native-like in many respects: (a) protein denaturation induced by urea is cooperative and fully reversible; (b) its oligomeric state at moderate concentration is well defined; and (c) MB-1 has very low affinity for 8-anilino-1-naphthalenesulfonic acid (ANSA), leading to enhancement of ANSA fluorescence that resembles that of other native proteins. On the other hand, our analysis revealed two aspects that command further attention. The folding stability of MB-1 as assessed by urea and thermal denaturation is somewhat less than that found for natural globular proteins of similar size. Size-exclusion chromatography experiments and analysis of MB-1 denaturation indicate that MB-1 is dimeric, not monomeric as designed. In light of these results, the utility and the current limitations of our design approach are discussed.

Ultraviolet Raman examination of the environmental dependence of bombolitin I and bombolitin III secondary structure

Bombolitin I and III (BI and BIII) are small amphiphilic peptides isolated from bumblebee venom. Although they exist in predominately nonhelical conformations in dilute aqueous solutions, we demonstrate, using UV Raman spectroscopy, that they become predominately alpha-helical in solution at pH > 10, in high ionic strength solutions, and in the presence of trifluoroethanol (TFE) and dodecylphosphocholine (DPC) micelles. In this paper, we examine the effects of electrostatic and hydrophobic interactions that control folding of BI and BIII by systematically monitoring their secondary structures as a function of solution conditions. We determine the BI and BIII secondary structure contents by using the quantitative UV Raman methodology of Chi et al. (1998. Biochemistry. 37:2854-2864). Our findings suggest that the alpha-helix turn in BIII at neutral pH is stabilized by a salt bridge between residues Asp2 and Lys5. This initial alpha-helical turn results in different BI and BIII alpha-helical folding mechanisms observed in high pH and high salt concentrations: BIII folds from its single alpha-helix turn close to its N-terminal, whereas the BI alpha-helix probably nucleates within the C-terminal half. We also used quasielastic light scattering to demonstrate that the BI and BIII alpha-helix formation in 0.2 M Ca(ClO4)2 is accompanied by formation of trimers and hexamers, respectively.

Orientation of the pore-forming peptide GALA in POPC vesicles determined by a BODIPY-avidin/biotin binding assay

We determined the orientation of a biotinylated version of the pore-forming peptide GALA (WEAALAEALAEALAEHLAEALAEALEALAA) at pH 5.0 in large unilamellar phosphatidylcholine vesicles, using the enhancement of BODIPY-avidin fluorescence subsequent to its irreversible binding to a biotin moiety. GALA and its variants were biotinylated at the N- or C-terminus. BODIPY-avidin was either added externally or was pre-encapsulated in vesicles to assess the fraction of liposome-bound biotinylated GALA that exposed its labeled terminus to the external or internal side of the bilayer, respectively. Under conditions where most of the membrane-bound peptides were involved in transmembrane aggregates and formed aqueous pores (at a lipid/bound peptide molar ratio of 2500/1), the head-to-tail (N- to C-terminus) orientation of the membrane-inserted peptides was such that 3/4 of the peptides exposed their N-terminus on the inside of the vesicle and their C-terminus on the outside. Under conditions resulting in reduced pore formation (at higher lipid/peptide molar ratios), we observed an increase in the fraction of GALA termini exposed to the outside of the vesicle. These results are consistent with a model (Parente et al., Biochemistry, 29:8720, 1990) that requires a critical number of peptides (M) in an aggregate to form a transbilayer structure. When the peptides form an aggregate of size i, with i < M = 4 to 6, the orientation of the peptides is mostly parallel to the membrane surface, such that both termini of the biotinylated peptide are exposed to external BODIPY-avidin. This BODIPY-avidin/biotin binding assay should be useful to determine the orientation of other membrane-interacting molecules.

Peptide-bilayer interactions: simulations of dermaseptin B, an antimicrobial peptide

Dermaseptins, a family of antimicrobial peptides, are believed to act by forming amphipathic alpha-helices which associate with the cell membrane, leading to its permeabilisation and disruption. A simple mean field method is described for simulation of the interactions of peptides with lipid bilayers which includes an approximate representation of the electrostatic effects of the head-group region of the bilayer. Starting from an atomistic model of a PC phospholipid bilayer we calculate an average electrostatic potential along the bilayer normal. By combining the interaction of the peptide with this electrostatic potential and with the hydrophobic core of the membrane we arrive at a more complete description of peptide-bilayer energetics than would be obtained using sidechain hydrophobicities alone. Using this interaction potential in MD simulations of the frog skin peptide dermaseptin B reveals that the lipid bilayer stabilises the alpha-helical conformation of the peptide. This is in agreement with FTIR data. A surface associated orientation thus appears to be the most stable arrangement of the peptide, at least at zero ionic strength and without taking account of possible peptide-peptide interactions.

An autonomous folding unit mediates the assembly of two-stranded coiled coils

Subunit oligomerization of many proteins is mediated by coiled-coil domains. Although the basic features contributing to the thermodynamic stability of coiled coils are well understood, the mechanistic details of their assembly have not yet been dissected. Here we report a 13-residue sequence pattern that occurs with limited sequence variations in many two-stranded coiled coils and that is absolutely required for the assembly of the Dictyostelium discoideum actin-bundling protein cortexillin I and the yeast transcriptional activator GCN4. The functional relationship between coiled-coil "trigger" sequences was manifested by replacing the intrinsic trigger motif of GCN4 with the related sequence from cortexillin I. We demonstrate that these trigger sequences represent autonomous helical folding units that, in contrast to arbitrarily chosen heptad repeats, can mediate coiled-coil formation. Aside from being of general interest for protein folding, trigger motifs should be of particular importance in the protein de novo design.

Fully integrated biocatalytic electrodes based on bioaffinity interactions

Integrated bioelectrocatalytically active electrodes are assembled by the deposition of enzymes onto respective electrically contacted affinity matrices and further cross-linking of the enzyme monolayers. A catalyst-NAD(+)-dyad for the binding of the NAD(+)-dependent enzymes and cytochrome-like molecules for the binding of the heme-protein-dependent enzymes are used to construct integrated electrically contacted biocatalytic systems. NAD(+)-dependent lactate dehydrogenase (LDH) is assembled onto a pyrroloquinoline quinone-NAD+ monolayer. The redox-active monolayer is organized via covalent attachment of pyrroloquinoline quinone (PQQ) to a cystamine monolayer associated with a Au-electrode, followed by covalent linkage of N6-(2-aminoethyl)-NAD+ to the monolayer. The interface modified with the PQQ-NAD(+)-dyad provides temporary affinity binding for LDH and allows cross-linking of the enzyme monolayer. The cross-linked LDH is bioelectrocatalytically active towards oxidation of lactate. The bioelectrocatalyzed process involves the PQQ-mediated oxidation of the immobilized NADH. Integrated, electrically contacted bioelectrodes are produced by the affinity binding and further cross-linking of nitrate reductase (NR) (cytochrome-dependent, E.C. 1.9.6.1 from E. coli) or CoII-protoporphyrin IX reconstituted myoglobin (CoII-Mb) atop the microperoxidase-11 (MP-11) monolayer associated with a Au-electrode. The MP-11 monolayer provides an affinity interface for the temporary binding of the enzymes, that allows the cross-linkage of the enzyme molecules. The MP-11 assembly acts as electron transfer mediator for the reduction of the secondary enzyme layer. The integrated bioelectrodes consisting of NR and CoII-Mb show catalytic activities for NO3- reduction and acetylene-dicarboxylic acid hydrogenation, respectively. Two FeIII-protoporphyrin IX units are reconstituted into a four alpha-helix bundle de novo protein assembled as a monolayer on a Au-electrode. Vectorial electron transfer proceeds in the synthetic heme-protein monolayer. Cross-linking of an affinity complex generated between the FeIII-protoporphyrin IX reconstituted de novo protein monolayer and NR yields an integrated, electrically contacted enzyme electrode that stimulates the bioelectrocatalyzed reduction of nitrate.

Crystallization and stabilization of MB-1, a de novo designed protein for optimized feeding technology

Milk Bundle-1 is a de novo protein that was designed for application in agriculture. It has a high content of selected essential amino acids, and is intended to adopt an alpha-helical bundle fold. Crystallization experiments with MB-1 have been carried out on the ground and in reduced gravity on board Columbia orbiter during mission STS-80. Rather small crystals were obtained (< 0.05 mm) in both environments. Among other factors, the lack of stability of purified MB-1 has been detrimental to crystal growth. We report here on our progress with regard to optimizing crystal growth conditions, protein purification and protein stability. The first MB-1 mutant we present (MB-1-His) contains a poly-histidine tail, allowing the use of metal affinity chromatography for purification. MB-1-His has been found to keep its original mass for a month at room temperature, a spectacular improvement over MB-1. The other mutant (MB-1-Cys) was engineered to carry a cysteine residue on a solvent exposed face. The exposed cysteine binds readily to p-HMB, and allows for dimerization of MB-1-Cys. The dimer was found to be twice as stable as MB-1 during proteolytic degradation studies.

Design of a 20-amino acid, three-stranded beta-sheet protein

A 20-residue protein (named Betanova) forming a monomeric, three-stranded, antiparallel beta sheet was designed using a structural backbone template and an iterative hierarchical approach. Structural and physicochemical characterization show that the beta-sheet conformation is stabilized by specific tertiary interactions and that the protein exhibits a cooperative two-state folding-unfolding transition, which is a hallmark of natural proteins. The Betanova molecule constitutes a tractable model system to aid in the understanding of beta-sheet formation, including beta-sheet aggregation and amyloid fibril formation.

Addition of a poly-(6X) His tag to Milk Bundle-1 and purification using immobilized metal-affinity chromatography

Milk Bundle 1 (MB-1) is a de novo designed protein enriched in M, T, K, and L. Its future application is as a high-quality dietary protein source for ruminants. The protein is currently expressed in Escherichia coli and is being characterized to solve its folded conformation. MB-1 has marginal stability at room temperature, which has hindered our attempts at characterization. To increase the stability of the protein at room temperature, the purification procedure was examined and changed to hopefully increase its effectiveness. We describe here the production and purification of a new MB-1 with six His residues at the C-terminal end. This allows the new mutant (MB-1-His) to bind metal ions and to be purified with immobilized metal-affinity chromatography (IMAC). MB-1-His obtained using IMAC was purer on SDS-PAGE than both MB-1 or MB-1-His isolated using the current protocol. The IMAC protocol is more economical and more efficient; preliminary results show that the protein purified by this method is also quite stable at room temperature.

From coiled coils to small globular proteins: design of a native-like three-helix bundle

A monomolecular native-like three-helix bundle has been designed in an iterative process, beginning with a peptide that noncooperatively assembled into an antiparallel three-helix bundle. Three versions of the protein were designed in which specific interactions were incrementally added. The hydrodynamic and spectroscopic properties of the proteins were examined by size exclusion chromatography, sedimentation equilibrium, fluorescence spectroscopy, and NMR. The thermodynamics of folding were evaluated by monitoring the thermal and guanidine-induced unfolding transitions using far UV circular dichroism spectroscopy. The attainment of a unique, native-like state was achieved through the introduction of: (1) helix capping interactions; (2) electrostatic interactions between partially exposed charged residues; (3) a diverse collection of apolar side chains within the hydrophobic core.

A differential scanning calorimetric study of the thermal unfolding of apo- and holo-cytochrome b562

Cytochrome b562 is a four-helix-bundle protein containing a non-covalently bound b-type heme prosthetic group. In the absence of heme, cytochrome b562 remains highly structured under native conditions. Here we report thermodynamic data for the thermal denaturation of the holo- and apoproteins as determined by differential scanning calorimetry. Thermal denaturation of holocytochrome b562 is a highly reversible process, and unexpectedly does not involve dissociation of the heme prosthetic group. Thermal denaturation of the corresponding apoprotein, with the heme group chemically removed, remains a cooperative, reversible process. Apocytochrome b562 is substantially destabilized relative to the holoprotein: the t1/2 is more than ten degrees lower, and enthalpy and heat capacity changes are about one-half of the holoprotein values. However, the energetic parameters of apocytochrome b562 denaturation are within the range of observed values for small proteins.

The construction of metal centers in proteins by rational design

Metalloprotein properties result from the interplay between coordination requirements of the metal center, protein stability, and modulation of the metal center by the surrounding protein matrix. Simple metal centers, which exercise control over the protein by affecting stability or enzyme activity, have been created by rational design. Complex centers, which require control by the protein matrix, have also been constructed.

Simplified proteins: minimalist solutions to the 'protein folding problem'

Recent research has suggested that stable, native proteins may be encoded by simple sequences of fewer than the full set of 20 proteogenic amino acids. Studies of the ability of simple amino acid sequences to encode stable, topologically complex, native conformations and to fold to these conformations in a biologically relevant time frame have provided insights into the sequence determinants of protein structure and folding kinetics. They may also have important implications for protein design and for theories of the origins of protein synthesis itself.

Effects of lanthanide binding on the stability of de novo designed alpha-helical coiled-coils

Effects of La3+ ion binding on the stability of de novo designed two-stranded alpha-helical coiled-coils were studied. The coiled-coils were composed of two 35-residue polypeptide chains based on the "native" heptad sequence Q(g)V(a)G(b)A(c)L(d)Q(e)K(f) and each contained a Cys residue at position 2a to allow formation of an interchain disulfide bridge. The effect of LaCl3 on the stability of five analogs containing two or three Glu substitutions per chain at heptad positions e and g was observed by urea denaturation at 20 degrees C. The analog E2(15,20), in which Glu residues are involved in interhelical i to i' + 5 repulsions, was stabilized relative to the control native peptide by addition of 50 mM LaCl3 to the buffer, whereas two analogs, in which Glu residues do not interact, were destabilized. These results suggest that LaCl3 may preferentially stabilize the folded state of E2(15,20) by the "bridging" of La3+ ions between two pairs of Glu residues usually involved in interhelical repulsions. Two analogs designed to contain two La3+ binding sites composed of three Glu residues each show greater stabilization by LaCl3 than E2(15,20) in the disulfide-bridged form. The apparent stabilization of E2(15,20) by La3+ binding was not observed with either Ca2+ or Mg2+, indicating that the effect is specific for trivalent versus divalent cations.

Solution structure of human p8MTCP1, a cysteine-rich protein encoded by the MTCP1 oncogene, reveals a new alpha-helical assembly motif

MTCP1 (for Mature-T-Cell Proliferation) is the first gene unequivocally identified in the group of uncommon leukemias with a mature phenotype. The three-dimensional solution structure of the human p8(MTCP1) protein encoded by the MTCP1 oncogene was determined by homonuclear proton two-dimensional NMR methods at 600 MHz. After sequence specific assignments, a total of 931 distance restraints and 57 dihedral restraints were collected. The location of the three previously unassigned disulfide bridges was determined from preliminary DIANA structures, using a statistical analysis of intercystinyl distances. The solution structure of p8(MTCP1) is presented as a set of 30 DIANA structures, further refined by restrained molecular dynamics using a simulated annealing protocol with the AMBER force field. The r.m.s.d. values with respect to the mean structure for the backbone and all heavy atoms for a family of 30 structures are 0.73(+/-0.28) and 1.17(+/-0.23) A, when the structured core of the protein (residues 5 to 63) is considered. The solution structure of p8(MTCP1) reveals an original scaffold consisting of three alpha helices, associated with a new cysteine motif. Two of the helices are covalently paired by two disulfide bridges, forming an alpha-hairpin which resembles an antiparallel coiled-coil. The third helix is oriented roughly parallel to the plane defined by the alpha-antiparallel motif and its axis forms an angle of approximately 60 degrees with respect to the main axis of this motif.