Design and engineering of an O(2) transport protein

Structural analysis of heme proteins: implications for design and prediction

BACKGROUND

Heme is an essential molecule and plays vital roles in many biological processes. The structural determination of a large number of heme proteins has made it possible to study the detailed chemical and structural properties of heme binding environment. Knowledge of these characteristics can provide valuable guidelines in the design of novel heme proteins and help us predict unknown heme binding proteins.

RESULTS

In this paper, we constructed a non-redundant dataset of 125 heme-binding protein chains and found that these heme proteins encompass at least 31 different structural folds with all-α class as the dominating scaffold. Heme binding pockets are enriched in aromatic and non-polar amino acids with fewer charged residues. The differences between apo and holo forms of heme proteins in terms of the structure and the binding pockets have been investigated. In most cases the proteins undergo small conformational changes upon heme binding. We also examined the CP (cysteine-proline) heme regulatory motifs and demonstrated that the conserved dipeptide has structural implications in protein-heme interactions.

CONCLUSIONS

Our analysis revealed that heme binding pockets show special features and that most of the heme proteins undergo small conformational changes after heme binding, suggesting the apo structures can be used for structure-based heme protein prediction and as scaffolds for future heme protein design.

De novo design and molecular assembly of a transmembrane diporphyrin-binding protein complex

The de novo design of membrane proteins remains difficult despite recent advances in understanding the factors that drive membrane protein folding and association. We have designed a membrane protein PRIME (PoRphyrins In MEmbrane) that positions two non-natural iron diphenylporphyrins (Fe(III)DPP's) sufficiently close to provide a multicentered pathway for transmembrane electron transfer. Computational methods previously used for the design of multiporphyrin water-soluble helical proteins were extended to this membrane target. Four helices were arranged in a D(2)-symmetrical bundle to bind two Fe(II/III) diphenylporphyrins in a bis-His geometry further stabilized by second-shell hydrogen bonds. UV-vis absorbance, CD spectroscopy, analytical ultracentrifugation, redox potentiometry, and EPR demonstrate that PRIME binds the cofactor with high affinity and specificity in the expected geometry.

Engineering enzymes

Fundamental research into bioinorganic catalysis of the kind presented at this Faraday Discussion has the potential to turn inspiration drawn from impressive natural energy and chemical transformations into artificial catalyst constructions useful to mankind. Creating bio-inspired artificial constructions requires a level of understanding well beyond simple description of structures and mechanisms of natural enzymes. To be useful, such description must be augmented by a practical sense of structural and energetic engineering tolerances of the mechanism. Significant barriers to achieving an engineering understanding of enzyme mechanisms arise from natural protein complexity. In certain cases we can surmount these barriers to understanding, such as natural electron tunneling, coupling of electron tunneling to light capture and proton exchange as well as simpler bond breaking redox catalysis. Hope for similar solutions of more complex bioinorganic enzymes is indicated in several papers presented in this Discussion. Armed with an engineering understanding of mechanism, the current serious frustrations to successful creation of functional artificial proteins that are rooted in protein complexity can fall away. Here we discuss the genetic and biological roots of protein complexity and show how to dodge and minimize the effects of complexity. In the best-understood cases, artificial enzymes can be designed from scratch using the simplest of protein scaffolds.

Heme proteins--diversity in structural characteristics, function, and folding

The characteristics of heme prosthetic groups and their binding sites have been analyzed in detail in a data set of nonhomologous heme proteins. Variations in the shape, volume, and chemical composition of the binding site, in the mode of heme binding and in the number and nature of heme-protein interactions are found to result in significantly different heme environments in proteins with different functions in biology. Differences are also seen in the properties of the apo states of the proteins. The apo states of proteins that bind heme permanently in their functional form show some disorder, ranging from local unfolding in the heme binding pocket to complete unfolding to give a random coil. In contrast, proteins that bind heme transiently are fully folded in their apo and holo states, presumably allowing both apo and holo forms to remain biologically active resisting aggregation or proteolysis. The principles identified here provide a framework for the design of de novo proteins that will exhibit tight heme ligand binding and for the identification of the function of structural genomic target proteins with heme ligands.

Cytosolic disulfide bond formation in cells infected with large nucleocytoplasmic DNA viruses

Proteins that have evolved to contain stabilizing disulfide bonds generally fold in a membrane-delimited compartment in the cell [i.e., the endoplasmic reticulum (ER) or the mitochondrial intermembrane space (IMS)]. These compartments contain sulfhydryl oxidase enzymes that catalyze the pairing and oxidation of cysteine residues. In contrast, most proteins in a healthy cytosol are maintained in reduced form through surveillance by NADPH-dependent reductases and the lack of sulfhydryl oxidases. Nevertheless, one of the core functionalities that unify the broad and diverse set of nucleocytoplasmic large DNA viruses (NCLDVs) is the ability to catalyze disulfide formation in the cytosol. The substrates of this activity are proteins that contribute to the assembly, structure, and infectivity of the virions. If the last common ancestor of NCLDVs was present during eukaryogenesis as has been proposed, it is interesting to speculate that viral disulfide bond formation pathways may have predated oxidative protein folding in intracellular organelles.

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.

Metal-directed protein self-assembly

Proteins are nature's premier building blocks for constructing sophisticated nanoscale architectures that carry out complex tasks and chemical transformations. Some 70%-80% of all proteins are thought to be permanently oligomeric; that is, they are composed of multiple proteins that are held together in precise spatial organization through noncovalent interactions. Although it is of great fundamental interest to understand the physicochemical basis of protein self-assembly, the mastery of protein-protein interactions (PPIs) would also allow access to novel biomaterials with nature's favorite and most versatile building block. In this Account, we describe a new approach we have developed with this possibility in mind, metal-directed protein self-assembly (MDPSA), which utilizes the strength, directionality, and selectivity of metal-ligand interactions to control PPIs. At its core, MDPSA is inspired by supramolecular coordination chemistry, which exploits metal coordination for the self-assembly of small molecules into discrete, more-or-less predictable higher order structures. Proteins, however, are not exactly small molecules or simple metal ligands: they feature extensive, heterogeneous surfaces that can interact with each other and with metal ions in unpredictable ways. We begin by first describing the challenges of using entire proteins as molecular building blocks. We follow with an examination of our work on a model protein (cytochrome cb(562)), highlighting challenges toward establishing ground rules for MDPSA as well as progress in overcoming these challenges. Proteins are also nature's metal ligands of choice. In MDPSA, once metal ions guide proteins into forming large assemblies, they are by definition embedded within extensive interfaces formed between protein surfaces. These complex surfaces make an inorganic chemist's life somewhat difficult, yet they also provide a wide platform to modulate the metal coordination environment through distant, noncovalent interactions, exactly as natural metalloproteins and enzymes do. We describe our computational and experimental efforts toward restructuring the noncovalent interaction network formed between proteins surrounding the interfacial metal centers. This approach, of metal templating followed by the redesign of protein interfaces (metal-templated interface redesign, MeTIR), not only provides a route to engineer de novo PPIs and novel metal coordination environments but also suggests possible parallels with the evolution of metalloproteins.

Roles of glutamates and metal ions in a rationally designed nitric oxide reductase based on myoglobin

A structural and functional model of bacterial nitric oxide reductase (NOR) has been designed by introducing two glutamates (Glu) and three histidines (His) in sperm whale myoglobin. X-ray structural data indicate that the three His and one Glu (V68E) residues bind iron, mimicking the putative Fe(B) site in NOR, while the second Glu (I107E) interacts with a water molecule and forms a hydrogen bonding network in the designed protein. Unlike the first Glu (V68E), which lowered the heme reduction potential by approximately 110 mV, the second Glu has little effect on the heme potential, suggesting that the negatively charged Glu has a different role in redox tuning. More importantly, introducing the second Glu resulted in a approximately 100% increase in NOR activity, suggesting the importance of a hydrogen bonding network in facilitating proton delivery during NOR reactivity. In addition, EPR and X-ray structural studies indicate that the designed protein binds iron, copper, or zinc in the Fe(B) site, each with different effects on the structures and NOR activities, suggesting that both redox activity and an intermediate five-coordinate heme-NO species are important for high NOR activity. The designed protein offers an excellent model for NOR and demonstrates the power of using designed proteins as a simpler and more well-defined system to address important chemical and biological issues.

Building fibrous biomaterials from alpha-helical and collagen-like coiled-coil peptides

Over the decade and a half, interest has soared in the development of peptide-based biomaterials and their potential applications in biotechnology. This review outlines the advances during this time in the construction of biomaterials based on the alpha-helical coiled-coil and collagen-like peptides. These structures and the resulting designs are distinct from the more commonly used beta-structured peptides, which often lead to hydrogels comprising amyloid-like fibrils. The review covers basic design rules for these helical assemblies, and the various fibrous biomaterials that can be accomplished with them, which include rigid structures with straight, branched, or networked structures, decorated and functionalized systems, and most recently flexible fibers and entangled hydrogel networks. This plethora of alpha-helix-based biomaterials, together with more recent collagen-like assemblies, that are emerging from various laboratories complement those developed using beta-structured peptides, and open exciting new avenues for biomaterials research and potential new application areas.

Standardized evaluation of protein stability

We compare mean force potential values of a large series of PDB models of proteins and peptides and find that, either as monomers or polymers, proteins longer than 200-250 residues have equivalent MFP values that are averaged to -65+/-3 kcal/aa. This value is named the standard or stability value. The standard value is reached irrespective of sequences and 3D folds. Peptides are too short to follow the rule and frequently exist as populations of conformers; one exception is peptides in amyloid fibrils. Fibrils surpass the standard value in accordance with their uppermost stability. In parallel, we calculate median MFP values of amino acids in stably folded PDB models of proteins: median values vary from -25 for Gly to -115 kcal/aa for Trp. These median values are used to score primary sequences of proteins: all sequences converge to a mean value of -63.5+/-2.5 kcal/aa, i.e., only 1.5 kcal less than the folded model standard. Sequences from unfolded proteins have lower values. This supports the conclusion that sequences carry in an important message and more specifically that diversity of amino acids in sequences is mandatory for stability. We also use the median amino acid MFP to score residue stability in 3D folds. This demonstrates that 3D folds are compromises between fragments of high and fragments of low scores and that functional residues are often but not always in the extreme score values. The approach opens to possibilities of evaluating any 3D model and of detecting functional residues and should help in conducting mutation assays.

Metal templated design of protein interfaces

Metal coordination is a key structural and functional component of a large fraction of proteins. Given this dual role we considered the possibility that metal coordination may have played a templating role in the early evolution of protein folds and complexes. We describe here a rational design approach, Metal Templated Interface Redesign (MeTIR), that mimics the time course of a hypothetical evolutionary pathway for the formation of stable protein assemblies through an initial metal coordination event. Using a folded monomeric protein, cytochrome cb(562), as a building block we show that its non-self-associating surface can be made self-associating through a minimal number of mutations that enable Zn coordination. The protein interfaces in the resulting Zn-directed, D(2)-symmetrical tetramer are subsequently redesigned, yielding unique protein architectures that self-assemble in the presence or absence of metals. Aside from its evolutionary implications, MeTIR provides a route to engineer de novo protein interfaces and metal coordination environments that can be tuned through the extensive noncovalent bonding interactions in these interfaces.

Designing artificial enzymes by intuition and computation

The rational design of artificial enzymes, either by applying physico-chemical intuition of protein structure and function or with the aid of computational methods, is a promising area of research with the potential to tremendously impact medicine, industrial chemistry and energy production. Designed proteins also provide a powerful platform for dissecting enzyme mechanisms of natural systems. Artificial enzymes have come a long way from simple α-helical peptide catalysts to proteins that facilitate multistep chemical reactions designed by state-of-the-art computational methods. Looking forward, we examine strategies employed by natural enzymes that could be used to improve the speed and selectivity of artificial catalysts.

Rational design of a structural and functional nitric oxide reductase

Protein design provides an ultimate test of our knowledge about proteins and allows the creation of novel enzymes for biotechnological applications. While progress has been made in designing proteins that mimic native proteins structurally^1–3, it is more difficult to design functional proteins^4–8. In comparison to recent successes in designing non-metalloproteins^4,6,7,9,10, it is even more challenging to rationally design metalloproteins that reproduce both the structure and function of native metalloenzymes^5,8,11–20, since protein metal binding sites are much more varied than non-metal containing sites, in terms of different metal ion oxidation states, preferred geometry and metal ion ligand donor sets. Because of their variability, it has been difficult to predict metal binding site properties in silico, as many of the parameters for metal binding sites, such as force fields are ill-defined. Therefore, the successful design of a structural and functional metalloprotein will greatly advance the field of protein design and our understanding of enzymes. Here, we report a successful, rational design of a structural and functional model of a metalloprotein, nitric oxide reductase (NOR), by introducing three histidines and one glutamate, predicted as ligands in the active site of NOR, into the distal pocket of myoglobin. A crystal structure of the designed protein confirms that the minimized computer model contains a heme/non-heme Fe_B center that is remarkably similar to that in the crystal structure. This designed protein also exhibits NOR activity. This is the first designed protein that models both the structure and function of NOR, offering insight that the active site glutamate is required for both iron binding and activity. These results show that structural and functional metalloproteins can be rationally designed in silico.

Mimicking the evolution of a thermally stable monomeric four-helix bundle by fusion of four identical single-helix peptides

Internal symmetry is a common feature of the tertiary structures of proteins and protein domains. Probably, because the genes of homo-oligomeric proteins duplicated and fused, their evolutionary descendants are proteins with internal symmetry. To identify any advantages that cause monomeric proteins with internal symmetry to be selected evolutionarily, we characterized some of the physical properties of a recombinant protein with a sequence consisting of two tandemly fused copies of the Escherichia coli Lac repressor C-terminal alpha-helix. This polypeptide exists in solution mainly as dimer that likely maintains a four-helix bundle motif. Thermal unfolding experiments demonstrate that the protein is considerably more stable at elevated temperatures than is a homotetramer consisting of four non-covalently associated copies of a 21-residue polypeptide similar in sequence to that of the Lac repressor C-terminal alpha-helix. A tandem duplication of our helix-loop-helix polypeptide yields an even more thermally stable protein. Our results exemplify the concept that fusion of non-covalently assembled polypeptide chains leads to enhanced protein stability. Herein, we discuss how our work relates to the evolutionary selective-advantages realized when symmetrical homo-oligomers evolve into monomers. Moreover, our thermally stable single-chain four-helix bundle protein may provide a robust scaffold for development of new biomaterials.

Design of a novel heme protein with a non-heme globin scaffold

A binding site for iron protoporphyrin IX (heme) was designed and embedded in a photosynthetic non-heme protein, phycocyanin, which forms a globin-like backbone structure, called a globin fold, but lacks sequence similarity to the globin family containing myoglobins and hemoglobins. Based on the structural alignment of the phycocyanin and myoglobin molecules, the proximal and distal His residues were repositioned in the phycocyanin sequence for heme ligation. The heme-binding pocket was created around the His residues by several residue replacements in the phycocyanin core. The synthesized phycocyanin variant, designated as HPY, bound one heme per protein molecule and showed spectroscopic features characteristic of six-coordinated heme proteins. The heme-binding HPY exhibited redox activity with an electrochemical midpoint potential of -130 mV against the standard hydrogen electrode, which was approximately 200 mV lower than the potential of natural myoglobins but 50 mV higher than the typical values of designed heme proteins with four-helix bundle or globin scaffolds. HPY also displayed native-like folding properties, in contrast to these designed heme proteins. However, the bound ferrous heme of HPY was quickly reoxidized by air and did not stably bind O(2), unlike the natural globins. The present results demonstrated that the globin fold of a non-globin protein is suitable for binding heme but is not sufficient for the reversible O(2) binding and myoglobin functions. The comparison of HPY with the natural globins may yield new insights into the essential features for realizing the natural heme protein functions.

How does a simplified-sequence protein fold?

To investigate a putatively primordial protein we have simplified the sequence of a 56-residue alpha/beta fold (the immunoglobulin-binding domain of protein G) by replacing it with polyalanine, polythreonine, and diglycine segments at regions of the sequence that in the folded structure are alpha-helical, beta-strand, and turns, respectively. Remarkably, multiple folding and unfolding events are observed in a 15-micros molecular dynamics simulation at 330 K. The most stable state (populated at approximately 20%) of the simplified-sequence variant of protein G has the same alpha/beta topology as the wild-type but shows the characteristics of a molten globule, i.e., loose contacts among side chains and lack of a specific hydrophobic core. The unfolded state is heterogeneous and includes a variety of alpha/beta topologies but also fully alpha-helical and fully beta-sheet structures. Transitions within the denatured state are very fast, and the molten-globule state is reached in <1 micros by a framework mechanism of folding with multiple pathways. The native structure of the wild-type is more rigid than the molten-globule conformation of the simplified-sequence variant. The difference in structural stability and the very fast folding of the simplified protein suggest that evolution has enriched the primordial alphabet of amino acids mainly to optimize protein function by stabilization of a unique structure with specific tertiary interactions.

Design of functional metalloproteins

Metalloproteins catalyse some of the most complex and important processes in nature, such as photosynthesis and water oxidation. An ultimate test of our knowledge of how metalloproteins work is to design new metalloproteins. Doing so not only can reveal hidden structural features that may be missing from studies of native metalloproteins and their variants, but also can result in new metalloenzymes for biotechnological and pharmaceutical applications. Although it is much more challenging to design metalloproteins than non-metalloproteins, much progress has been made in this area, particularly in functional design, owing to recent advances in areas such as computational and structural biology.

Dimer interface migration in a viral sulfhydryl oxidase

Large double-stranded DNA viruses, including poxviruses and mimiviruses, encode enzymes to catalyze the formation of disulfide bonds in viral proteins produced in the cell cytosol, an atypical location for oxidative protein folding. These viral disulfide catalysts belong to a family of sulfhydryl oxidases that are dimers of a small five-helix fold containing a Cys-X-X-Cys motif juxtaposed to a flavin adenine dinucleotide cofactor. We report that the sulfhydryl oxidase pB119L from African swine fever virus (ASFV) uses for self-assembly surface different from that observed in homologs from mammals, plants, and fungi. Within a protein family, different packing interfaces for the same oligomerization state are extremely rare. We find that the alternate dimerization mode seen in ASFV pB119L is not characteristic of all viral sulfhydryl oxidases, as the flavin-binding domain from a mimivirus sulfhydryl oxidase assumes the same dimer structure as the known eukaryotic enzymes. ASFV pB119L demonstrates the potential of large double-stranded DNA viruses, which have faster mutation rates than their hosts and the tendency to incorporate host genes, to pioneer new protein folds and self-assembly modes.

Multiple pathways guide oxygen diffusion into flavoenzyme active sites

Dioxygen (O(2)) and other gas molecules have a fundamental role in a variety of enzymatic reactions. However, it is only poorly understood which O(2) uptake mechanism enzymes employ to promote efficient catalysis and how general this is. We investigated O(2) diffusion pathways into monooxygenase and oxidase flavoenzymes, using an integrated computational and experimental approach. Enhanced-statistics molecular dynamics simulations reveal spontaneous protein-guided O(2) diffusion from the bulk solvent to preorganized protein cavities. The predicted protein-guided diffusion paths and the importance of key cavity residues for oxygen diffusion were verified by combining site-directed mutagenesis, rapid kinetics experiments, and high-resolution X-ray structures. This study indicates that monooxygenase and oxidase flavoenzymes employ multiple funnel-shaped diffusion pathways to absorb O(2) from the solvent and direct it to the reacting C4a atom of the flavin cofactor. The difference in O(2) reactivity among dehydrogenases, monooxygenases, and oxidases ultimately resides in the fine modulation of the local environment embedding the reactive locus of the flavin.