Solution Structure of a Designed Four-α-Helix Bundle Maquette Scaffold

Catalysis and Electron Transfer in De Novo Designed Metalloproteins

One of the hallmark advances in our understanding of metalloprotein function is showcased in our ability to design new, non-native, catalytically active protein scaffolds. This review highlights progress and milestone achievements in the field of de novo metalloprotein design focused on reports from the past decade with special emphasis on de novo designs couched within common subfields of bioinorganic study: heme binding proteins, monometal- and dimetal-containing catalytic sites, and metal-containing electron transfer sites. Within each subfield, we highlight several of what we have identified as significant and important contributions to either our understanding of that subfield or de novo metalloprotein design as a discipline. These reports are placed in context both historically and scientifically. General suggestions for future directions that we feel will be important to advance our understanding or accelerate discovery are discussed.

De Novo Design, Solution Characterization, and Crystallographic Structure of an Abiological Mn-Porphyrin-Binding Protein Capable of Stabilizing a Mn(V) Species

De novo protein design offers the opportunity to test our understanding of how metalloproteins perform difficult transformations. Attaining high-resolution structural information is critical to understanding how such designs function. There have been many successes in the design of porphyrin-binding proteins; however, crystallographic characterization has been elusive, limiting what can be learned from such studies as well as the extension to new functions. Moreover, formation of highly oxidizing high-valent intermediates poses design challenges that have not been previously implemented: (1) purposeful design of substrate/oxidant access to the binding site and (2) limiting deleterious oxidation of the protein scaffold. Here we report the first crystallographically characterized porphyrin-binding protein that was programmed to not only bind a synthetic Mn-porphyrin but also maintain binding site access to form high-valent oxidation states. We explicitly designed a binding site with accessibility to dioxygen units in the open coordination site of the Mn center. In solution, the protein is capable of accessing a high-valent Mn(V)-oxo species which can transfer an O atom to a thioether substrate. The crystallographic structure is within 0.6 Å of the design and indeed contained an aquo ligand with a second water molecule stabilized by hydrogen bonding to a Gln side chain in the active site, offering a structural explanation for the observed reactivity.

De novo design of a hyperstable non-natural protein-ligand complex with sub-Å accuracy

Protein catalysis requires the atomic-level orchestration of side chains, substrates and cofactors, and yet the ability to design a small-molecule-binding protein entirely from first principles with a precisely predetermined structure has not been demonstrated. Here we report the design of a novel protein, PS1, that binds a highly electron-deficient non-natural porphyrin at temperatures up to 100 °C. The high-resolution structure of holo-PS1 is in sub-Å agreement with the design. The structure of apo-PS1 retains the remote core packing of the holoprotein, with a flexible binding region that is predisposed to ligand binding with the desired geometry. Our results illustrate the unification of core packing and binding-site definition as a central principle of ligand-binding protein design.

De Novo Construction of Redox Active Proteins

Relatively simple principles can be used to plan and construct de novo proteins that bind redox cofactors and participate in a range of electron-transfer reactions analogous to those seen in natural oxidoreductase proteins. These designed redox proteins are called maquettes. Hydrophobic/hydrophilic binary patterning of heptad repeats of amino acids linked together in a single-chain self-assemble into 4-alpha-helix bundles. These bundles form a robust and adaptable frame for uncovering the default properties of protein embedded cofactors independent of the complexities introduced by generations of natural selection and allow us to better understand what factors can be exploited by man or nature to manipulate the physical chemical properties of these cofactors. Anchoring of redox cofactors such as hemes, light active tetrapyrroles, FeS clusters, and flavins by His and Cys residues allow cofactors to be placed at positions in which electron-tunneling rates between cofactors within or between proteins can be predicted in advance. The modularity of heptad repeat designs facilitates the construction of electron-transfer chains and novel combinations of redox cofactors and new redox cofactor assisted functions. Developing de novo designs that can support cofactor incorporation upon expression in a cell is needed to support a synthetic biology advance that integrates with natural bioenergetic pathways.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.

Creating a bio-hybrid signal transduction pathway: opening a new channel of communication between cells and machines

Manipulation of signal transduction pathways presents a viable mechanism to interface cells with electronics. In this work, we present a two-step signal transduction pathway involving cellular and electronic transduction elements. In order to circumvent many of the conventional difficulties encountered when harnessing chemical signalling for the purpose of electronics communication, gaseous nitric oxide (NO) was selected as the signalling molecule. By genetic engineering of the nitric oxide synthase protein eNOS and insertion of light-oxygen-voltage (LOV) domains, we have created a photoactive version of the protein. The novel chimeric eNOS was found to be capable of producing NO in response to excitation by visible light. By coupling these mutant cells to a surface modified platinum electrode, it was possible to convert an optical signal into a chemical one, followed by subsequent conversion of the chemical signal into an electrical output.

Manipulating cofactor binding thermodynamics in an artificial oxygen transport protein

We report the mutational analysis of an artificial oxygen transport protein, HP7, which operates via a mechanism akin to that of human neuroglobin and cytoglobin. This protein destabilizes one of two heme-ligating histidine residues by coupling histidine side chain ligation with the burial of three charged glutamate residues on the same helix. Replacement of these glutamate residues with alanine, which is uncharged, increases the affinity of the distal histidine ligand by a factor of 13. Paradoxically, it also decreases heme binding affinity by a factor of 5 in the reduced state and 60 in the oxidized state. Application of a three-state binding model, in which an initial pentacoordinate binding event is followed by a protein conformational change to hexacoordinate, provides insight into the mechanism of this seemingly counterintuitive result: the initial pentacoordinate encounter complex is significantly destabilized by the loss of the glutamate side chains, and the increased affinity for the distal histidine only partially compensates for that. These results point to the importance of considering each oxidation and conformational state in the design of functional artificial proteins.

The crystal structure of the calcium-bound con-G[Q6A] peptide reveals a novel metal-dependent helical trimer

The ability to form and control both secondary structure and oligomerization in short peptides has proven to be challenging owing to the structural instability of such peptides. The conantokin peptides are a family of γ-carboxyglutamic acid containing peptides produced in the venoms of predatory sea snails of the Conus family. They are examples of short peptides that form stable helical structures, especially in the presence of divalent cations. Both monomeric and dimeric conantokin peptides have been identified and represent a new mechanism of helix association, "the metallozipper motif" that is devoid of a hydrophobic interface between monomers. In the present study, a parallel/antiparallel three-helix bundle was identified and its crystal structure determined at high resolution. The three helices are almost perfectly parallel and represent a novel helix-helix association. The trimer interface is dominated by metal chelation between the three helices, and contains no interfacial hydrophobic interactions. It is now possible to produce stable monomeric, dimeric, or trimeric metallozippers depending on the peptide sequence and metal ion. Such structures have important applications in protein design.

Understanding properties of cofactors in proteins: redox potentials of synthetic cytochromes b

Haehnel et al. synthesized 399 different artificial cytochrome b (aCb) models. They consist of a template-assisted four-helix bundle with one embedded heme group. Their redox potentials were measured and cover the range from -148 to -89 mV. No crystal structures of these aCb are available. Therefore, we use the chemical composition and general structural principles to generate atomic coordinates of 31 of these aCb mutants, which are chosen to cover the whole interval of redox potentials. We start by modeling the coordinates of one aCb from scratch. Its structure remains stable after energy minimization and during molecular dynamics simulation over 2 ns. Based on this structure, coordinates of the other 30 aCb mutants are modeled. The calculated redox potentials for these 31 aCb agree within 10 mV with the experimental values in terms of root mean square deviation. Analysis of the dependence of heme redox potential on protein environment shows that the shifts in redox potentials relative to the model systems in water are due to the low-dielectric medium of the protein and the protonation states of the heme propionic acid groups, which are influenced by the surrounding amino acids. Alternatively, we perform a blind prediction of the same redox potentials using an empirical approach based on a linear scoring function and reach a similar accuracy. Both methods are useful to understand and predict heme redox potentials. Based on the modeled structure we can understand the detailed structural differences between aCb mutants that give rise to shifts in heme redox potential. On the other hand, one can explore the correlation between sequence variations and aCb redox potentials more directly and on much larger scale using the empirical prediction scheme, which--thanks to its simplicity--is much faster.

Artificial photoactive proteins

Solar power is the most abundant source of renewable energy. In this respect, the goal of making photoactive proteins is to utilize this energy to generate an electron flow. Photosystems have provided the blueprint for making such systems, since they are capable of converting the energy of light into an electron flow using a series of redox cofactors. Protein tunes the redox potential of the cofactors and arranges them such that their distance and orientation are optimal for the creation of a stable charge separation. The aim of this review is to present an overview of the literature with regard to some elegant functional structures that protein designers have created by introducing cofactors and photoactivity into synthetic proteins.

A mixed-α,β miniprotein stereochemically reprogrammed to high-binding affinity for acetylcholine

The protein-structure space is limited to L configuration in the asymmetric alpha-amino acid structures; the function space, on other hand, seems limitless because of the chemical diversity in the amino acid side chain structures. The chemical diversity in side chain structure may be multiplied beneficially with the stereochemical diversity in main chain structure; thus, de novo protein design may be explored for customizing molecular structures stereochemically and molecular functions chemically. Illustrating de novo design in the structure space of L and D alphabet, canonical all-beta folds of poly-L structure were reprogrammed as bracelet, boat, and canoe-shaped molecules-the "boat" as a receptor-like pocket and the "canoe" as a metal-ion receptor-simply by mutating specific L-amino acid residues to the corresponding D stereochemical structure. Demonstrating customization of molecular shape stereochemically and function chemically, a 15-residue mixed-alpha, beta miniprotein of canonical poly-L structure is now reprogrammed stereochemically as a cup-shaped receptor for acetylcholine via cation-pi interaction with a triad of aromatic side chains placed in mimicry of the acetylcholine-receptor sites both natural and artificial. Evidence from CD, fluorescence, NMR, DSC, ITC, MD, and molecular-docking studies is presented to show that a rationally designed 15-residue mixed-L, D peptide is a cooperatively ordered molecular fold in the stereochemically specified molecular morphology, submicromolar in affinity of acetylcholine and thus an acetylcholine receptor exceptionally small and simple. .

Thermodynamic investigation into the mechanisms of proton-coupled electron transfer events in heme protein maquettes

To study the engineering requirements for proton pumping in energy-converting enzymes such as cytochrome c oxidase, the thermodynamics and mechanisms of proton-coupled electron transfer in two designed heme proteins are elucidated. Both heme protein maquettes chosen, heme b-[H10A24]2 and heme b-[delta7-His]2, are four-alpha-helix bundles that display pH-dependent heme midpoint potential modulations, or redox-Bohr effects. Detailed equilibrium binding studies of ferric and ferrous heme b with these maquettes allow the individual contributions of heme-protein association, iron-histidine ligation, and heme-protein electrostatics to be elucidated. These data demonstrate that the larger, less well-structured [H10A24]2 binds heme b in both oxidation states tighter than the smaller and more well-structured [Delta7-His]2 due to a stronger porphyrin-protein hydrophobic interaction. The 66 mV (1.5 kcal/mol) difference in their heme reduction potentials observed at pH 8.0 is due mostly to stabilization of ferrous heme in [H10A24]2 relative to [delta7-His]2. The data indicate that porphyrin-protein hydrophobic interactions and heme iron coordination are responsible for the Kd value of 37 nM for the heme b-[delta7-His]2 scaffold, while the affinity of heme b for [H10A24]2 is 20-fold tighter due to a combination of porphyrin-protein hydrophobic interactions, iron coordination, and electrostatic effects. The data also illustrate that the contribution of bis-His coordination to ferrous heme protein affinity is limited, <3.0 kcal/mol. The 1H+/1e- redox-Bohr effect of heme b-[H10A24]2 is due to the greater absolute stabilization of the ferric heme (4.1 kcal/mol) compared to the ferrous heme (1.4 kcal/mol) binding upon glutamic acid deprotonation, i.e., an electrostatic response mechanism. The 2H+/1e- redox-Bohr effect observed for heme b-[delta7-His]2 is due to histidine protonation and histidine dissociation of ferrous heme b upon reduction, i.e., a ligand loss mechanism. These results indicate that the contribution of porphyrin-protein hydrophobic interactions to heme affinity is critical to maintaining the heme bound in both oxidation states and eliciting an electrostatic response from these designed heme protein scaffolds.

Three-dimensional structure and dynamics of a de novo designed, amphiphilic, metallo-porphyrin-binding protein maquette at soft interfaces by molecular dynamics simulations

The three-dimensional structure and dynamics of de novo designed, amphiphilic four-helix bundle peptides (or "maquettes"), capable of binding metallo-porphyrin cofactors at selected locations along the length of the core of the bundle, are investigated via molecular dynamics simulations. The rapid evolution of the initial design to stable three-dimensional structures in the absence (apo-form) and presence (holo-form) of bound cofactors is described for the maquettes at two different soft interfaces between polar and nonpolar media. This comparison of the apo- versus holo-forms allows the investigation of the effects of cofactor incorporation on the structure of the four-helix bundle. The simulation results are in qualitative agreement with available experimental data describing the structures at lower resolution and limited dimension.

Nativelike Structure in Designed Four α-Helix Bundles Driven by Buried Polar Interactions

We show that a single internal polar interaction per helix is sufficient to engender structural specificity in that helix in helical bundle proteins. Furthermore, we use histidine-binding cofactors of different shapes which bind directly into the core, demonstrating that this structural specificity is not the result of a prescribed complimentary, "knobs in holes" core packing. We show that we can switch structural specificity of individual helices on and off by ligating cofactors, singly and in pairs, which bind either one or two histidine ligands. To our knowledge, this is the first demonstration of such extensive manipulation of protein structure by ligand binding, an important result of general interest to those working with self-assembled molecular systems. Finally, as these proteins were designed without the use of computational modeling, we not only demonstrate that designing a uniquely structured cofactor binding protein is not as difficult as is generally believed, we have determined why this is so: hydrophobic core complementarity, which is very difficult to design, is not necessary. Instead, a much simpler design process entails the creation of core polar interactions which themselves can drive conformational specificity.

Intelligent design: the de novo engineering of proteins with specified functions

One of the principal successes of de novo protein design has been the creation of small, robust protein-cofactor complexes which can serve as simplified models, or maquettes, of more complicated multicofactor protein complexes commonly found in nature. Different maquettes, generated by us and others, recreate a variety of aspects, or functional elements, recognized as parts of natural enzyme function. The current challenge is to both expand the palette of functional elements and combine and/or integrate them in recreating familiar enzyme activities or generating novel catalysis in the simplest protein scaffolds.

Exploring amino-acid radical chemistry: protein engineering and de novo design

Amino-acid radical enzymes are often highly complex structures containing multiple protein subunits and cofactors. These properties have in many cases hampered the detailed characterization of their amino-acid redox cofactors. To address this problem, a range of approaches has recently been developed in which a common strategy is to reduce the complexity of the radical-containing system. This work will be reviewed and it includes the light-induced generation of aromatic radicals in small-molecule and peptide systems. Natural redox proteins, including the blue copper protein azurin and a bacterial photosynthetic reaction center, have been engineered to introduce amino-acid radical chemistry. The redesign strategies to achieve this remarkable change in the properties of these proteins will be described. An additional approach to gain insights into the properties of amino-acid radicals is to synthesize de novo designed model proteins in which the redox chemistry of these species can be studied. Here we describe the design, synthesis and characteristics of monomeric three-helix bundle and four-helix bundle proteins designed to study the redox chemistry of tryptophan and tyrosine. This work demonstrates that de novo protein design combined with structural, electrochemical and quantum chemical analyses can provide detailed information on how the protein matrix tunes the thermodynamic properties of tryptophan.

Ab initio prediction of the three-dimensional structure of a de novo designed protein: A double-blind case study

Ab initio structure prediction and de novo protein design are two problems at the forefront of research in the fields of structural biology and chemistry. The goal of ab initio structure prediction of proteins is to correctly characterize the 3D structure of a protein using only the amino acid sequence as input. De novo protein design involves the production of novel protein sequences that adopt a desired fold. In this work, the results of a double-blind study are presented in which a new ab initio method was successfully used to predict the 3D structure of a protein designed through an experimental approach using binary patterned combinatorial libraries of de novo sequences. The predicted structure, which was produced before the experimental structure was known and without consideration of the design goals, and the final NMR analysis both characterize this protein as a 4-helix bundle. The similarity of these structures is evidenced by both small RMSD values between the coordinates of the two structures and a detailed analysis of the helical packing.

Specular neutron reflectivity and the structure of artificial protein maquettes vectorially oriented at interfaces

Artificial peptides can be designed to possess a variety of functionalities. If these peptides can be ordered in an ensemble, the functionality can impart macroscopic material properties to the ensemble. Neutron reflectivity is shown to be an effective probe of the intramolecular structures of such peptides vectorially oriented at an interface, key to ensuring that the designed molecular structures translate into the desired material properties of the interface. A model-independent method is utilized to analyze the neutron reflectivity from an alkylated, di- alpha -helical peptide, containing perdeuterated leucine residues at one or two pre-selected positions, in mixed Langmuir monolayers with a phospholipid. The results presented here are more definitive than prior work employing x-ray reflectivity. They show explicitly that the di-helical peptide retains its designed alpha -helical secondary structure at the interface, when oriented perpendicular to the interface at high surface pressure, with the helices projecting into the aqueous subphase without penetrating the layer of phospholipid headgroups.

Aromatic interactions in peptides: Impact on structure and function

Aromatic interactions, including pi-pi, cation-pi, aryl-sulfur, and carbohydrate-pi interactions, have been shown to be prevalent in proteins through protein structure analysis, suggesting that they are important contributors to protein structure. However, the magnitude and significance of aromatic interactions is not defined by such studies. Investigation of aromatic interactions in the context of structured peptides has complemented studies of protein structure and has provided a wealth of information regarding the role of aromatic interactions in protein structure and function. Recent advances in this area are reviewed.

A model membrane protein for binding volatile anesthetics

Earlier work demonstrated that a water-soluble four-helix bundle protein designed with a cavity in its nonpolar core is capable of binding the volatile anesthetic halothane with near-physiological affinity (0.7 mM Kd). To create a more relevant, model membrane protein receptor for studying the physicochemical specificity of anesthetic binding, we have synthesized a new protein that builds on the anesthetic-binding, hydrophilic four-helix bundle and incorporates a hydrophobic domain capable of ion-channel activity, resulting in an amphiphilic four-helix bundle that forms stable monolayers at the air/water interface. The affinity of the cavity within the core of the bundle for volatile anesthetic binding is decreased by a factor of 4-3.1 mM Kd as compared to its water-soluble counterpart. Nevertheless, the absence of the cavity within the otherwise identical amphiphilic peptide significantly decreases its affinity for halothane similar to its water-soluble counterpart. Specular x-ray reflectivity shows that the amphiphilic protein orients vectorially in Langmuir monolayers at higher surface pressure with its long axis perpendicular to the interface, and that it possesses a length consistent with its design. This provides a successful starting template for probing the nature of the anesthetic-peptide interaction, as well as a potential model system in structure/function correlation for understanding the anesthetic binding mechanism.

A Stable Miniature Protein with Oxaloacetate Decarboxylase Activity

An 18-residue miniature enzyme, Apoxaldie-1, has been designed, based on the known structure of the neurotoxic peptide apamin. Three lysine residues were introduced on the solvent-exposed face of the apamin alpha-helix to serve as an active site for decarboxylation of oxaloacetate. The oxidised form of Apoxaldie-1, in which two disulfide bonds stabilise the alpha-helix, formed spontaneously. CD spectroscopy measurements revealed that, in its oxidised form, Apoxaldie-1 adopted a stably folded structure, which was lost upon reduction of the disulfide bonds. Despite its small size and the absence of a designed binding pocket, Apoxaldie-1 displayed saturation kinetics in its oxidised form and catalysed the decarboxylation of oxaloacetate at a rate that was almost four orders of magnitude faster than that observed with n-butylamine. This rivals the performance of the best synthetic oxaloacetate decarboxylases reported to date. Unlike those, however, Apoxaldie-1 displayed significant stability. It maintained its secondary structure at temperatures in excess of 75 degrees C, in the presence of high concentrations of guanidinium chloride and at pH values as low as 2.2. Apamin-based catalysts have potential for the generation of miniature peptides that display activity under nonphysiological conditions.

Design and Synthesis of de Novo Cytochromes c

Natural c-type cytochromes are characterized by the consensus Cys-X-X-Cys-His heme-binding motif (where X is any amino acid) by which the heme is covalently attached to protein by the addition of the sulfhydryl groups of two cysteine residues to the vinyl groups of the heme. In this work, the consensus sequence was used for the heme-binding site of a designed four-helix bundle, and the apoproteins with either a histidine residue or a methionine residue positioned at the sixth coordination site were synthesized and reacted with iron protoporphyrin IX (protoheme) under mild reducing conditions in vitro. These polypeptides bound one heme per helix-loop-helix monomer via a single thioether bond and formed four-helix bundle dimers in the holo forms as designed. They exhibited visible absorption spectra characteristic of c-type cytochromes, in which the absorption bands shifted to lower wavelengths in comparison with the b-type heme binding intermediates of the same proteins. Unexpectedly, the designed cytochromes c with bis-His-coordinated heme iron exhibited oxidation-reduction potentials similar to those of their b-type intermediates, which have no thioether bond. Furthermore, the cytochrome c with His and Met residues as the axial ligands exhibited redox potentials increased by only 15-30 mV in comparison with the cytochrome with the bis-His coordination. These results indicate that highly positive redox potentials of natural cytochromes c are not only due to the heme covalent structure, including the Met ligation, but also due to noncovalent and hydrophobic environments surrounding the heme. The covalent attachment of heme to the polypeptide in natural cytochromes c may contribute to their higher redox potentials by reducing the thermodynamic stability of the oxidized forms relatively against that of the reduced forms without the loss of heme.

Amphiphilic 4-helix bundles designed for biomolecular materials applications

Artificial peptides previously designed to possess alpha-helical bundle motifs have been either hydrophilic (i.e., soluble in polar media) or lipophilic (i.e., soluble in nonpolar media) in overall character. Realizations of these bioinspired bundles have succeeded in reproducing a variety of biomimetic functionality within the appropriate media. However, to translate their functionality into any biomolecular device applications at the macroscopic level, the bundles must be oriented in an ensemble, for example, at an interface. This goal has been realized in a new family of alpha-helical bundle peptides which are amphiphilic; namely, they assemble into 4-helix bundles with well-defined hydrophilic and hydrophobic domains. These peptides are capable of binding metalloporphyrin prosthetic groups at selected locations within these domains. We describe here the realization of one of the first members of this family, AP0, successfully designed for vectorial incorporation into soft interfaces between polar and nonpolar media.

The HP-1 maquette: from an apoprotein structure to a structured hemoprotein designed to promote redox-coupled proton exchange

Synthetic heme-binding four-alpha-helix bundles show promise as working model systems, maquettes, for understanding heme cofactor-protein assembly and function in oxidoreductases. Despite successful inclusion of several key functional elements of natural proteins into a family of heme protein maquettes, the lack of 3D structures, due principally to conformational heterogeneity, has prevented them from achieving their full potential. We report here the design and synthesis of HP-1, a disulfide-bridged two-alpha-helix peptide that self-assembles to form an antiparallel twofold symmetric diheme four-alpha-helix bundle protein with a stable conformation on the NMR time-scale. The HP-1 design strategy began with the x-ray crystal structure of the apomaquette L31M, an apomaquette derived from the structurally heterogeneous tetraheme-binding H10H24 prototype. L31M was functionally redesigned to accommodate two hemes ligated to histidines and to retain the strong coupling of heme oxidation-reduction to glutamate acid-base transitions and proton exchange that was characterized in molten globule predecessors. Heme insertion was modeled with angular constraints statistically derived from natural proteins, and the pattern of hydrophobic and hydrophilic residues on each helix was then altered to account for this large structural reorganization. The transition to structured holomaquette involved the alteration of 6 of 31 residues in each of the four identical helices and, unlike our earlier efforts, required no design intermediates. Oxidation-reduction of both hemes displays an unusually low midpoint potential (-248 mV vs. normal hydrogen electrode at pH 9.0), which is strongly coupled to proton binding, as designed.

De novo designed cyclic-peptide heme complexes

The structural characterization of de novo designed metalloproteins together with determination of chemical reactivity can provide a detailed understanding of the relationship between protein structure and functional properties. Toward this goal, we have prepared a series of cyclic peptides that bind to water-soluble metalloporphyrins (FeIII and CoIII). Neutral and positively charged histidine-containing peptides bind with a high affinity, whereas anionic peptides bind only weakly to the negatively charged metalloporphyrin. Additionally, it was found that the peptide becomes helical only in the presence of the metalloporphyrin. CD experiments confirm that the metalloporphyrin binds specific cyclic peptides with high affinity and with isodichroic behavior. Thermal unfolding experiments show that the complex has "native-like" properties. Finally, NMR spectroscopy produced well dispersed spectra and experimental restraints that provide a high-resolution solution structure of the complexed peptide.

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.

Identification of amino acids involved in protein structural uniqueness: implication for de novo protein design

Structural uniqueness is characteristic of native proteins and is essential to express their biological functions. The major factors that bring about the uniqueness are specific interactions between hydrophobic residues and their unique packing in the protein core. To find the origin of the uniqueness in their amino acid sequences, we analyzed the distribution of the side chain rotational isomers (rotamers) of hydrophobic amino acids in protein tertiary structures and derived deltaS(contact), the conformational-entropy changes of side chains by residue-residue contacts in each secondary structure. The deltaS(contact) values indicate distinct tendencies of the residue pairs to restrict side chain conformation by inter-residue contacts. Of the hydrophobic residues in alpha-helices, aliphatic residues (Leu, Val, Ile) strongly restrict the side chain conformations of each other. In beta-sheets, Met is most strongly restricted by contact with Ile, whereas Leu, Val and Ile are less affected by other residues in contact than those in alpha-helices. In designed and native protein variants, deltaS(contact) was found to correlate with the folding-unfolding cooperativity. Thus, it can be used as a specificity parameter for designing artificial proteins with a unique structure.

Dissection of the pathway of molecular recognition by calmodulin

Amide hydrogen exchange has been used to examine the structural dynamics and energetics of the interaction of a peptide corresponding to the calmodulin-binding domain of smooth muscle myosin light chain kinase (smMLCKp) with calcium-saturated calmodulin. Heteronuclear NMR (15)N-(1)H correlation spectroscopy was used to quantify amide proton exchange rates of the uniformly (15)N-labeled domain bound to calmodulin. A key feature of a proposed model for molecular recognition by calmodulin [Ehrhardt et al. (1995) Biochemistry 34, 2731-2738] is tested by examination of the dependence of amide hydrogen exchange on applied hydrostatic pressure. Hydrogen exchange rates and corresponding protection factors (1/K(op)) for individual amide protons of the bound smMLCKp domain span 5 orders of magnitude at ambient pressure. Individual protection factors decrease significantly in a linear fashion with increasing hydrostatic pressure. A common pressure dependence is revealed by a constant large negative volume change across the residues comprising the core of the bound helical domain. The pattern of protection factors and their response to hydrostatic pressure is consistent with a structural reorganization that results in the concerted disruption of ion pairs between calmodulin and the bound domain. These observations reinforce a model for the molecular recognition pathway where formation of the initial encounter complex is followed by helix-coil transitions in the bound state and subsequent concerted formation of the extensive ion pair network defining the intermolecular contact surface between CaM and the target domain in the final, compact complex structure.

Oligomeric beta(beta)(alpha) miniprotein motifs: pivotal role of single hinge residue in determining the oligomeric state

The role of a single glycine hinge residue in the structure of BBAT1, a beta(beta)(alpha) peptide that forms a discrete homotrimeric structure in solution, was evaluated with 11 new peptide sequences which differ only in the identity of the residue at the hinge position. The integrity of the structure and oligomeric state of the peptides was evaluated by using a combination of analytical ultracentrifugation and circular dichroism spectroscopy. Initially, it was discovered that the glycine hinge adopts backbone dihedral angles favored in D-amino acids and that incorporation of D-alanine at the hinge position stabilizes the trimer species. Subsequently, the effect of the side chains of different D-amino acids at the hinge position was evaluated. While incorporation of polar amino acids led to a destabilization of the oligomeric form of the peptide, only peptides including D-Ser or D-Asp at the hinge position were able to achieve a discrete trimer species. Incorporation of hydrophobic amino acids D-Leu and D-Phe led to oligomerization beyond a trimer to a tetrameric form. The dramatic differences among the thermodynamic stabilities and oligomeric states of these peptides illustrates the pivotal role of the hinge residue in the oligomerization of the beta(beta)(alpha) peptides.

Hydrophobic modulation of heme properties in heme protein maquettes

We have investigated the properties of the two hemes bound to histidine in the H10 positions of the uniquely structured apo form of the heme binding four-helix bundle protein maquette [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2) for the amino acids at positions 6 (I), 13 (F) and 24 (H), respectively. The primary structure of each alpha-helix, alpha-SH, in [I(6)F(13)H(24)](2) is Ac-CGGGEI(6)WKL.H(10)EEF(13)LKK.FEELLKL.H(24)EERLKK.L-CONH(2). In our nomenclature, [I(6)F(13)H(24)] represents the disulfide-bridged di-alpha-helical homodimer of this sequence, i.e., (alpha-SS-alpha), and [I(6)F(13)H(24)](2) represents the dimeric four helix bundle composed of two di-alpha-helical subunits, i.e., (alpha-SS-alpha)(2). We replaced the histidines at positions H24 in [I(6)F(13)H(24)](2) with hydrophobic amino acids incompetent for heme ligation. These maquette variants, [I(6)F(13)I(24)](2), [I(6)F(13)A(24)](2), and [I(6)F(13)F(24)](2), are distinguished from the tetraheme binding parent peptide, [I(6)F(13)H(24)](2), by a reduction in the heme:four-helix bundle stoichiometry from 4:1 to 2:1. Iterative redesign has identified phenylalanine as the optimal amino acid replacement for H24 in the context of apo state conformational specificity. Furthermore, the novel second generation diheme [I(6)F(13)F(24)](2) maquette was related to the first generation diheme [H10A24](2) prototype, [L(6)L(13)A(24)](2) in the present nomenclature, via a sequential path in sequence space to evaluate the effects of conservative hydrophobic amino acid changes on heme properties. Each of the disulfide-linked dipeptides studied was highly helical (>77% as determined from circular dichroism spectroscopy), self-associates in solution to form a dimer (as determined by size exclusion chromatography), is thermodynamically stable (-DeltaG(H)2(O) >18 kcal/mol), and possesses conformational specificity that NMR data indicate can vary from multistructured to single structured. Each peptide binds one heme with a dissociation constant, K(d1) value, tighter than 65 nM forming a series of monoheme maquettes. Addition of a second equivalent of heme results in heme binding with a K(d2) in the range of 35-800 nM forming the diheme maquette state. Single conservative amino acid changes between peptide sequences are responsible for up to 10-fold changes in K(d) values. The equilibrium reduction midpoint potential (E(m7.5)) determined in the monoheme state ranges from -156 to -210 mV vs SHE and in the diheme state ranges from -144 to -288 mV. An observed heme-heme electrostatic interaction (>70 mV) in the diheme state indicates a syn global topology of the di-alpha-helical monomers. The heme affinity and electrochemistry of the three H24 variants studied identify the tight binding sites (K(d1) and K(d2) values <200 nM) having the lower reduction midpoint potentials (E(m7.5) values of -155 and -260 mV) with the H10 bound hemes in the parent tetraheme state of [H10H24-L6I,L13F](2), here called [I(6)F(13)H(24)](2). The results of this study illustrate that conservative hydrophobic amino acid changes near the heme binding site can modulate the E(m) by up to +/-50 mV and the K(d) by an order of magnitude. Furthermore, the effects of multiple single amino acid changes on E(m) and K(d) do not appear to be additive.

Emerging principles of de novo catalyst design

Considerable progress has been made in the understanding of how to exploit hydrophobic and charge-charge interactions in forming binding sites for peptides and small molecules in folded polypeptide catalysts. This knowledge has enabled the introduction of feedback and control functions into catalytic cycles and the construction of folded polypeptide catalysts that follow saturation kinetics. Major advances have also been made in the design of metalloproteins and metallopeptides, especially with regards to understanding redox potential control.

Artificial cytochrome b: computer modeling and evaluation of redox potentials

We generated atomic coordinates of an artificial protein that was recently synthesized to model the central part of the native cytochrome b (Cb) subunit consisting of a four-helix bundle with two hemes. Since no X-ray structure is available, the structural elements of the artificial Cb were assembled from scratch using all known chemical and structural information available and avoiding strain as much as possible. Molecular dynamics (MD) simulations applied to this model protein exhibited root-mean-square deviations as small as those obtained from MD simulations starting with the crystal structure of the native Cb subunit. This demonstrates that the modeled structure of the artificial Cb is relatively rigid and strain-free. The model structure of the artificial Cb was used to determine the redox potentials of the two hemes by calculating the electrostatic energies from the solution of the linearized Poisson-Boltzmann equation (LPBE). The calculated redox potentials agree within 20 meV with the experimentally measured values. The dependence of the redox potentials of the hemes on the protein environment was analyzed. Accordingly, the total shift in the redox potentials is mainly due to the low dielectric medium of the protein, the protein backbone charges, and the salt bridges formed between the arginines and the propionic acid groups of the hemes. The difference in the shift of the redox potentials is due to the interactions with the hydrophilic side chains and the salt bridges formed with the propionic acids of the hemes. For comparison and to test the computational procedure, the redox potentials of the two hemes in the native Cb from the cytochrome bc(1) (Cbc(1)) complex were also calculated. Also in this case the computed redox potentials agree well with experiments.

Main chain and side chain dynamics of a heme protein: 15N and 2H NMR relaxation studies of R. capsulatus ferrocytochrome c2

A detailed characterization of the main chain and side chain dynamics in R. capsulatus ferrocytochrome c(2) derived from (2)H NMR relaxation of methyl group resonances is presented. (15)N relaxation measurements confirm earlier results indicating that R. capsulatus ferrocytochrome c(2) exhibits minor rotational anisotropy in solution. The current study is focused on the use of deuterium relaxation in side chain methyl groups, which has been shown to provide a detailed and accurate measure of internal dynamics. Results obtained indicate that the side chains of ferrocytochrome c(2) exhibit a wide range of motional amplitudes, but are more rigid than generally found in the interior of nonprosthetic group bearing globular proteins. This unusual rigidity is ascribed to the interactions of the protein with the large heme prosthetic group. This observation has significant implications for the potential of the heme-protein interface to modulate the redox properties of the protein and also points to the need for great precision in the design and engineering of heme proteins.

Proof of principle in a de novo designed protein maquette: an allosterically regulated, charge-activated conformational switch in a tetra-alpha-helix bundle

New understanding of the engineering and allosteric regulation of natural protein conformational switches (such as those that couple chemical and ionic signals, mechanical force, and electro/chemical free energy for biochemical activation, catalysis, and motion) can be derived from simple de novo designed synthetic protein models (maquettes). We demonstrate proof of principle of both reversible switch action and allosteric regulation in a tetra-alpha-helical bundle protein composed of two identical di-helical subunits containing heme coordinated at a specific position close to the disulfide loop region. Individual bundles assume one of two switch states related by large-scale mechanical changes: a syn-topology (helices of the different subunits parallel) or anti-topology (helices antiparallel). Both the spectral properties of a coproporphyrin probe appended to the loop region and the distance-dependent redox interaction between the hemes identify the topologies. Beginning from a syn-topology, introduction of ferric heme in each subunit (either binding or redox change) shifts the topological balance by 25-50-fold (1.9-2.3 kcal/mol) to an anti-dominance. Charge repulsion between the two internal cationic ferric hemes drives the syn- to anti-switch, as demonstrated in two ways. When fixed in the syn-topology, the second ferric heme binding is 25-80-fold (1.9-2.6 kcal/mol) weaker than the first, and adjacent heme redox potentials are split by 80 mV (1.85 kcal/mol), values that energetically match the shift in topological balance. Allosteric and cooperative regulation of the switch by ionic strength exploits the shielded charge interactions between the two hemes and the exposed, cooperative interactions between the coproporphyrin carboxylates.

Oligomerization of uniquely folded mini-protein motifs: development of a homotrimeric betabetaalpha peptide

The discovery of a discretely folded homotrimeric betabetaalpha motif (BBAT1) was recently reported (J. Am. Chem. Soc. 2001, 123, 1002-1003). Herein the design, synthesis, and analysis of a small library of peptides which led to the isolation of BBAT1 is described. betabetaalpha peptides based on the monomeric sequence of BBA5 (Folding Des. 1998, 120, 95-103) were synthesized to include the anthranilic acid/nitrotyrosine fluorescence quenching pair to rapidly detect interpeptide association. In the first generation of peptides synthesized, truncations in the loop region connecting the beta-hairpin to the alpha-helix revealed that a two-residue deletion in the loop promoted an interpeptide association as detected by fluorescence quenching. An additional library of 22 loop-truncated betabetaalpha peptides was subsequently synthesized to include a variety of sequence mutations in an effort to enhance the observed peptide-peptide binding. From the fluorescence quenching screen, peptide B2 was found to possess the strongest fluorescence-quenching response, indicative of a strong peptide-peptide association. Due the poor solubility of peptide B2, the S-methylated cysteine at position 9 in the loop was substituted with a glycine to generate peptide BBAT1 which possessed greatly improved water solubility and formed discrete trimers. The successful design of this oligomeric betabetaalpha structure will likely aid the design of more complex alpha-beta superstructures and further our understanding of the factors controlling protein-protein interactions at alpha-beta protein interfaces.

Designing metal-peptide models for protein structure and function

Recent progress in the rational design of metal sites within peptide model systems shows increasing control in the placement of metals within helical bundles and inclusion of sophisticated elements such as second-sphere ligand interactions. A crystallographically characterized two-metal peptide model for diiron proteins represents a major achievement in de novo design methodologies. Increasingly complex and robust models for electron transfer through and between helices, and electrode-supported electron-transfer peptides, have been constructed. Design elements for peptide-supported ferredoxins and mononuclear Fe(II) and Zn(II) sites have been refined.

Heme redox potential control in de novo designed four-alpha-helix bundle proteins

The effects of various mechanisms of metalloporphyrin reduction potential modulation were investigated experimentally using a robust, well-characterized heme protein maquette, synthetic protein scaffold H10A24 [¿CH(3)()CONH-CGGGELWKL.HEELLKK.FEELLKL.AEERLKK. L-CONH(2)()¿(2)](2). Removal of the iron porphyrin macrocycle from the high dielectric aqueous environment and sequestration within the hydrophobic core of the H10A24 maquette raises the equilibrium reduction midpoint potential by 36-138 mV depending on the hydrophobicity of the metalloporphyrin structure. By incorporating various natural and synthetic metalloporphyrins into a single protein scaffold, we demonstrate a 300-mV range in reduction potential modulation due to the electron-donating/withdrawing character of the peripheral macrocycle substituents. Solution pH is used to modulate the metalloporphyrin reduction potential by 160 mV, regardless of the macrocycle architecture, by controlling the protonation state of the glutamate involved in partial charge compensation of the ferric heme. Attempts to control the reduction potential by inserting charged amino acids into the hydrophobic core at close proximity to the metalloporphyrin lead to varied success, with H10A24-L13E lowering the E(m8.5) by 40 mV, H10A24-E11Q raising it by 50 mV, and H10A24-L13R remaining surprisingly unaltered. Modifying the charge of the adjacent metalloporphyrin, +1 for iron(III) protoporphyrin IX or neutral for zinc(II) protoporphyrin IX resulted in a loss of 70 mV [Fe(III)PPIX](+) - [Fe(III)PPIX](+) interaction observed in maquettes. Using these factors in combination, we illustrate a 435-mV variation of the metalloporphyrin reduction midpoint potential in a simple heme maquette relative to the about 800-mV range observed for natural cytochromes. Comparison between the reduction potentials of the heme maquettes and other de novo designed heme proteins reveals global trends in the E(m) values of synthetic cytochromes.

De novo design of helical bundles as models for understanding protein folding and function

De novo protein design has proven to be a powerful tool for understanding protein folding, structure, and function. In this Account, we highlight aspects of our research on the design of dimeric, four-helix bundles. Dimeric, four-helix bundles are found throughout nature, and the history of their design in our laboratory illustrates our hierarchic approach to protein design. This approach has been successfully applied to create a completely native-like protein. Structural and mutational analysis allowed us to explore the determinants of native protein structure. These determinants were then applied to the design of a dinuclear metal-binding protein that can now serve as a model for this important class of proteins.

Self-assembly of heme A and heme B in a designed four-helix bundle: implications for a cytochrome c oxidase maquette

Heme A, a prosthetic group of cytochrome c oxidase [EC 1.9.3.1], has been introduced into two de novo designed four helix bundle proteins, [H10A24](2) and [H10H24](2), known to bind 2-4 equiv of heme B, respectively [Robertson, D. E., Farid, R. S., Moser, C. C., Mulholland, S. E., Pidikiti, R., Lear, J. D., Wand, A., J., DeGrado, W. F., and Dutton, P. L. (1994) Nature 368, 425-432]. [H10A24](2), [Ac-CGGGELWKL x HEELLKK x FEELLKL x AEERLKK x L-CONH(2)](2)(2), binds two heme A molecules per four-helix unit via bis-histidine ligation at the 10,10' positions with measured K(d) values of <0.1 and 5 nM, values much lower than those measured for heme B (K(d) values of 50 and 800 nM). The heme A-protein complex, [heme A-H10A24](2), exhibits well-defined absorption spectra in both the ferric and ferrous states, and an electron paramagnetic resonance spectrum characteristic of a low spin heme in the ferric form. A single midpoint redox potential (E(m8)) was determined for [heme A-H10A24](2) at -45 mV (vs SHE), which is significantly higher than that of the protein bound heme B (-130 and -200 mV). The observation of a single midpoint redox potential for [heme A-H10A24](2) and a pair of midpoints for [heme B-H10A24](2) indicates that the di-alpha-helical monomers are oriented in an anti topology (disulfides on opposite sides of bundle) in the former (lacking heme-heme electrostatic interaction) and syn in the latter. A mixture of global topologies was indicated by the potentiometric titration of the related [heme A-H10H24](2) which possess two distinct reduction potentials of +41 (31%) and -65 mV (69%). Self-assembly of the mixed cofactor heme A-heme B-[H10A24](2) was accomplished by addition of a single equivalent of each heme A and heme B to [H10A24](2). The single midpoint redox potential of heme B, E(m8) = -200 mV, together with the split midpoint redox potential of heme A in heme A-heme B-[H10A24](2), E(m8) = +28 mV (33%) and -65 mV (67%), indicated the existence of both syn and anti topologies of the two di-alpha-helical monomers in this four helix bundle. Synthesis of the mixed cofactor [heme A-heme B-H10H24](2) was accomplished by addition of a 2 equiv of each heme A and heme B to [H10H24](2) and potentiometry indicated the pair of hemes B resided in the 10,10' sites and heme A occupied the 24,24' sites. The results indicate that heme peripheral structure controls the orientation of the di-alpha-helical monomers in the four-helix bundle which are interchangeable between syn and anti topologies. In the reduced form, [heme A-H10A24](2), reacts quantitatively to form [carbonmonoxy-heme A-H10A24](2) as evidenced by optical spectroscopy. The synthetic [heme A-H10A24](2) can be enzymatically reduced by NAD(P)H with natural reductases under anaerobic conditions, and reversibly oxidized by dioxygen to the ferric form.

Retrostructural analysis of metalloproteins: application to the design of a minimal model for diiron proteins

De novo protein design provides an attractive approach for the construction of models to probe the features required for function of complex metalloproteins. The metal-binding sites of many metalloproteins lie between multiple elements of secondary structure, inviting a retrostructural approach to constructing minimal models of their active sites. The backbone geometries comprising the metal-binding sites of zinc fingers, diiron proteins, and rubredoxins may be described to within approximately 1 A rms deviation by using a simple geometric model with only six adjustable parameters. These geometric models provide excellent starting points for the design of metalloproteins, as illustrated in the construction of Due Ferro 1 (DF1), a minimal model for the Glu-Xxx-Xxx-His class of dinuclear metalloproteins. This protein was synthesized and structurally characterized as the di-Zn(II) complex by x-ray crystallography, by using data that extend to 2.5 A. This four-helix bundle protein is comprised of two noncovalently associated helix-loop-helix motifs. The dinuclear center is formed by two bridging Glu and two chelating Glu side chains, as well as two monodentate His ligands. The primary ligands are mostly buried in the protein interior, and their geometries are stabilized by a network of hydrogen bonds to second-shell ligands. In particular, a Tyr residue forms a hydrogen bond to a chelating Glu ligand, similar to a motif found in the diiron-containing R2 subunit of Escherichia coli ribonucleotide reductase and the ferritins. DF1 also binds cobalt and iron ions and should provide an attractive model for a variety of diiron proteins that use oxygen for processes including iron storage, radical formation, and hydrocarbon oxidation.

A designed four-alpha-helix bundle that binds the volatile general anesthetic halothane with high affinity

The structural features of volatile anesthetic binding sites on proteins are being examined with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work has suggested that introducing a cavity into the hydrophobic core improves anesthetic binding affinity. The more polarizable methionine side chain was substituted for a leucine, in an attempt to enhance the dispersion forces between the ligand and the protein. The resulting bundle variant has an improved affinity (K(d) = 0.20 +/- 0.01 mM) for halothane binding, compared with the leucine-containing bundle (K(d) = 0.69 +/- 0.06 mM). Photoaffinity labeling with (14)C-halothane reveals preferential labeling of the W15 residue in both peptides, supporting the view that fluorescence quenching by bound anesthetic reports both the binding energetics and the location of the ligand in the hydrophobic core. The rates of amide hydrogen exchange were similar for the two bundles, suggesting that differences in binding affinity were not due to changes in protein stability. Binding of halothane to both four-alpha-helix bundle proteins stabilized the native folded conformations. Molecular dynamics simulations of the bundles illustrate the existence of the hydrophobic core, containing both W15 residues. These results suggest that in addition to packing defects, enhanced dispersion forces may be important in providing higher affinity anesthetic binding sites. Alternatively, the effect of the methionine substitution on halothane binding energetics may reflect either improved access to the binding site or allosteric optimization of the dimensions of the binding pocket. Finally, preferential stabilization of folded protein conformations may represent a fundamental mechanism of inhaled anesthetic action.