Design and synthesis of a globin fold

Comparison of Evolutionary Relationships between Branchiostoma floridae, Ciona intestinalis, and Homo sapiens Globins Provide Evidence of Gene Co-Option and Convergent Evolution

Globins have been studied as model proteins to elucidate the principles of protein evolution. This was achieved by understanding the relationship between amino acid sequence, three-dimensional structure, physicochemical properties, and physiological function. Previous molecular phylogenies of chordate globin genes revealed the monophyletic evolution of urochordate globins and suggested convergent evolution. However, to provide evidence of convergent evolution, it is necessary to determine the physicochemical and functional similarities between vertebrates and urochordate globins. In this study, we determined the expression patterns of Ciona globin genes using real-time RT-PCR. Two genes (Gb-1 and Gb-2) were predominantly expressed in the branchial sac, heart, and hemocytes and were induced under hypoxia. Combined with the sequence analysis, our findings suggest that Gb-1/-2 correspond to vertebrate hemoglobin-α/-β. However, we did not find a robust similarity between Gb-3, Gb-4, and vertebrate globins. These results suggested that, even though Ciona globins obtained their unique functions differently from vertebrate globins, the two of them shared some physicochemical features and physiological functions. Our findings offer a good example for understanding the molecular mechanisms underlying gene co-option and convergence, which could lead to evolutionary innovations.

Architectural digest: Thermodynamic stability and domain structure of a consensus monomeric globin

Artificial proteins representing the consensus of a set of homologous sequences have attracted attention for their increased thermodynamic stability and conserved activity. Here, we applied the consensus approach to a b-type heme-binding protein to inspect the contribution of a dissociable cofactor to enhanced stability and the chemical consequences of creating a generic heme environment. We targeted the group 1 truncated hemoglobin (TrHb1) subfamily of proteins for their small size (∼120 residues) and ease of characterization. The primary structure, derived from a curated set of ∼300 representative sequences, yielded a highly soluble consensus globin (cGlbN) enriched in acidic residues. Optical and NMR spectroscopies revealed high-affinity heme binding in the expected site and in two orientations. At neutral pH, proximal and distal iron coordination was achieved with a pair of histidine residues, as observed in some natural TrHb1s, and with labile ligation on the distal side. As opposed to studied TrHb1s, which undergo additional folding upon heme binding, cGlbN displayed the same extent of secondary structure whether the heme was associated with the protein or not. Denaturation required guanidine hydrochloride and showed that apo- and holoprotein unfolded in two transitions-the first (occurring with a midpoint of ∼2 M) was shifted to higher denaturant concentration in the holoprotein (∼3.7 M) and reflected stabilization due to heme binding, while the second transition (∼6.2 M) was common to both forms. Thus, the consensus sequence stabilized the protein but exposed the existence of two separately cooperative subdomains within the globin architecture, masked as one single domain in TrHb1s with typical stabilities. The results suggested ways in which specific chemical or thermodynamic features may be controlled in artificial heme proteins.

Structural basis for heme detoxification by an ATP-binding cassette-type efflux pump in gram-positive pathogenic bacteria

Bacterial pathogens acquire heme from the host hemoglobin as an iron nutrient for their virulence and proliferation in blood. Concurrently, they encounter cytotoxic-free heme that escapes the heme-acquisition process. To overcome this toxicity, many gram-positive bacteria employ an ATP-binding cassette heme-dedicated efflux pump, HrtBA in the cytoplasmic membranes. Although genetic analyses have suggested that HrtBA expels heme from the bacterial membranes, the molecular mechanism of heme efflux remains elusive due to the lack of protein studies. Here, we show the biochemical properties and crystal structures of Corynebacterium diphtheriae HrtBA, alone and in complex with heme or an ATP analog, and we reveal how HrtBA extracts heme from the membrane and releases it. HrtBA consists of two cytoplasmic HrtA ATPase subunits and two transmembrane HrtB permease subunits. A heme-binding site is formed in the HrtB dimer and is laterally accessible to heme in the outer leaflet of the membrane. The heme-binding site captures heme from the membrane using a glutamate residue of either subunit as an axial ligand and sequesters the heme within the rearranged transmembrane helix bundle. By ATP-driven HrtA dimerization, the heme-binding site is squeezed to extrude the bound heme. The mechanism sheds light on the detoxification of membrane-bound heme in this bacterium.

Protein-Embedded Metalloporphyrin Arrays Templated by Circularly Permuted Tobacco Mosaic Virus Coat Proteins

Bioenergetic processes in nature have relied on networks of cofactors for harvesting, storing, and transforming the energy from sunlight into chemical bonds. Models mimicking the structural arrangement and functional crosstalk of the cofactor arrays are important tools to understand the basic science of natural systems and to provide guidance for non-natural functional biomaterials. Here, we report an artificial multiheme system based on a circular permutant of the tobacco mosaic virus coat protein (cpTMV). The double disk assembly of cpTMV presents a gap region sandwiched by the two C₂-symmetrically related disks. Non-native bis-his coordination sites formed by the mutation of the residues in this gap region were computationally screened and experimentally tested. A cpTMV mutant Q101H was identified to create a circular assembly of 17 protein-embedded hemes. Biophysical characterization using X-ray crystallography, cyclic voltammetry, and electron paramagnetic resonance (EPR) suggested both structural and functional similarity to natural multiheme cytochrome c proteins. This protein framework offers many further engineering opportunities for tuning the redox properties of the cofactors and incorporating non-native components bearing varied porphyrin structures and metal centers. Emulating the electron transfer pathways in nature using a tunable artificial system can contribute to the development of photocatalytic materials and bioelectronics.

Functionalization of Bacterial Microcompartment Shell Proteins With Covalently Attached Heme

Heme is a versatile redox cofactor that has considerable potential for synthetic biology and bioelectronic applications. The capacity to functionalize non-heme-binding proteins with covalently bound heme moieties in vivo could expand the variety of bioelectronic materials, particularly if hemes could be attached at defined locations so as to facilitate position-sensitive processes like electron transfer. In this study, we utilized the cytochrome maturation system I to develop a simple approach that enables incorporation of hemes into the backbone of target proteins in vivo. We tested our methodology by targeting the self-assembling bacterial microcompartment shell proteins, and inserting functional hemes at multiple locations in the protein backbone. We found substitution of three amino acids on the target proteins promoted heme attachment with high occupancy. Spectroscopic measurements suggested these modified proteins covalently bind low-spin hemes, with relative low redox midpoint potentials (about -210 mV vs. SHE). Heme-modified shell proteins partially retained their self-assembly properties, including the capacity to hexamerize, and form inter-hexamer attachments. Heme-bound shell proteins demonstrated the capacity to integrate into higher-order shell assemblies, however, the structural features of these macromolecular complexes was sometimes altered. Altogether, we report a versatile strategy for generating electron-conductive cytochromes from structurally-defined proteins, and provide design considerations on how heme incorporation may interface with native assembly properties in engineered proteins.

Supramolecular polymer formation by a de novo hemoprotein with a synthetic diheme compound

Proteins are attractive materials for supramolecular chemistry due to their multifunctionality and self‐organization ability. In this work, we synthesized a diheme compound, in which two iron‐protoporphyrin IX molecules are associated via a linker chain, and introduced it into a de novo designed four‐helix bundle protein with two heme‐binding sites. The protein gradually bound the diheme compound by bis‐histidyl ligation and formed supramolecular polymers. Polymer formation was observed by atomic force microscopy (AFM), which revealed the highly branched, dendritic forms of the fibrous architecture. The present results may open a pathway toward nanowire construction with de novo heme‐proteins.

Design of a heme-binding peptide motif adopting a β-hairpin conformation

Heme-binding proteins constitute a large family of catalytic and transport proteins. Their widespread presence as globins and as essential oxygen and electron transporters, along with their diverse enzymatic functions, have made them targets for protein design. Most previously reported designs involved the use of α-helical scaffolds, and natural peptides also exhibit a strong preference for these scaffolds. However, the reason for this preference is not well-understood, in part because alternative protein designs, such as those with β-sheets or hairpins, are challenging to perform. Here, we report the computational design and experimental validation of a water-soluble heme-binding peptide, Pincer-1, composed of predominantly β-scaffold secondary structures. Such heme-binding proteins are rarely observed in nature, and by designing such a scaffold, we simultaneously increase the known fold space of heme-binding proteins and expand the limits of computational design methods. For a β-scaffold, two tryptophan zipper β-hairpins sandwiching a heme molecule were linked through an N-terminal cysteine disulfide bond. β-Hairpin orientations and residue selection were performed computationally. Heme binding was confirmed through absorbance experiments and surface plasmon resonance experiments (K_D = 730 ± 160 nm). CD and NMR experiments validated the β-hairpin topology of the designed peptide. Our results indicate that a helical scaffold is not essential for heme binding and reveal the first designed water-soluble, heme-binding β-hairpin peptide. This peptide could help expand the search for and design space to cytoplasmic heme-binding proteins.

Design and fine-tuning redox potentials of metalloproteins involved in electron transfer in bioenergetics

Redox potentials are a major contributor in controlling the electron transfer (ET) rates and thus regulating the ET processes in the bioenergetics. To maximize the efficiency of the ET process, one needs to master the art of tuning the redox potential, especially in metalloproteins, as they represent major classes of ET proteins. In this review, we first describe the importance of tuning the redox potential of ET centers and its role in regulating the ET in bioenergetic processes including photosynthesis and respiration. The main focus of this review is to summarize recent work in designing the ET centers, namely cupredoxins, cytochromes, and iron-sulfur proteins, and examples in design of protein networks involved these ET centers. We then discuss the factors that affect redox potentials of these ET centers including metal ion, the ligands to metal center and interactions beyond the primary ligand, especially non-covalent secondary coordination sphere interactions. We provide examples of strategies to fine-tune the redox potential using both natural and unnatural amino acids and native and nonnative cofactors. Several case studies are used to illustrate recent successes in this area. Outlooks for future endeavors are also provided. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Intermediate Tyrosyl Radical and Amyloid Structure in Peroxide-Activated Cytoglobin

We characterized the peroxidase mechanism of recombinant rat brain cytoglobin (Cygb) challenged by hydrogen peroxide, tert-butylhydroperoxide and by cumene hydroperoxide. The peroxidase mechanism of Cygb is similar to that of myoglobin. Cygb challenged by hydrogen peroxide is converted to a Fe⁴⁺ oxoferryl π cation, which is converted to Fe⁴⁺ oxoferryl and tyrosyl radical detected by direct continuous wave-electron paramagnetic resonance and by 3,5-dibromo-4-nitrosobenzene sulfonate spin trapping. When organic peroxides are used as substrates at initial reaction times, and given an excess of peroxide present, the EPR signals of the corresponding peroxyl radicals precede those of the direct tyrosyl radical. This result is consistent with the use of peroxide as a reducing agent for the recycling of Cygb high-valence species. Furthermore, we found that the Cygb oxidation by peroxides leads to the formation of amyloid fibrils. This result suggests that Cygb possibly participates in the development of degenerative diseases; our findings also support the possible biological role of Cygb related to peroxidase activity.

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.

Design and engineering of an O(2) transport protein

The principles of natural protein engineering are obscured by overlapping functions and complexity accumulated through natural selection and evolution. Completely artificial proteins offer a clean slate on which to define and test these protein engineering principles, while recreating and extending natural functions. Here we introduce this method with the design of an oxygen transport protein, akin to human neuroglobin. Beginning with a simple and unnatural helix-forming sequence with just three different amino acids, we assembled a four-helix bundle, positioned histidines to bis-histidine ligate haems, and exploited helical rotation and glutamate burial on haem binding to introduce distal histidine strain and facilitate O(2) binding. For stable oxygen binding without haem oxidation, water is excluded by simple packing of the protein interior and loops that reduce helical-interface mobility. O(2) affinities and exchange timescales match natural globins with distal histidines, with the remarkable exception that O(2) binds tighter than CO.

Design of λ Cro Fold: Solution Structure of a Monomeric Variant of the De Novo Protein

One of the classical DNA-binding proteins, bacteriophage lambda Cro, forms a homodimer with a unique fold of alpha-helices and beta-sheets. We have computationally designed an artificial sequence of 60 amino acid residues to stabilize the backbone tertiary structure of the lambda Cro dimer by simulated annealing using knowledge-based structure-sequence compatibility functions. The designed amino acid sequence has 25% identity with that of natural lambda Cro and preserves Phe58, which is important for formation of the stably folded structure of lambda Cro. The designed dimer protein and its monomeric variant, which was redesigned by the insertion of a beta-hairpin sequence at the C-terminal region to prevent dimerization, were synthesized and biochemically characterized to be well folded. The designed protein was monomeric under a wide range of protein concentrations and its solution structure was determined by NMR spectroscopy. The solved structure is similar to that of a monomeric variant of natural lambda Cro with a root-mean-square deviation of the polypeptide backbones at 2.1A and has a well-packed protein core. Thus, our knowledge-based functions provide approximate but essential relationships between amino acid sequences and protein structures, and are useful for finding novel sequences that are foldable into a given target structure.

The response of internal dynamics to hydrophobic core mutations in the SH3 domain from the Fyn tyrosine kinase

We have used (15)N- and (2)H-NMR spin relaxation experiments to study the response of backbone and side-chain dynamics when a leucine or valine is substituted for a completely buried phenylalanine residue in the SH3 domain from the Fyn tyrosine kinase. Several residues show differences in the time scales and temperature dependences of internal motions when data for the three proteins are compared. Changes were also observed in the magnitude of dynamics, with the valine, and to a lesser extent leucine mutant, showing enhanced flexibility compared to the wild-type (WT) protein. The motions of many of the same amide and methyl groups are affected by both mutations, identifying a set of loci where dynamics are sensitive to interactions involving the targeted side chain. These results show that contacts within the hydrophobic core affect many aspects of internal mobility throughout the Fyn SH3 domain.

An in silico exploration of the neutral network in protein sequence space

Designating amino-acid sequences that fold into a common main-chain structure as "neutral sequences" for the structure, regardless of their function or stability, we investigated the distribution of neutral sequences in protein sequence space. For four distinct target structures (alpha, beta,alpha/beta and alpha+beta types) with the same chain length of 108, we generated the respective neutral sequences by using the inverse folding technique with a knowledge-based potential function. We assumed that neutral sequences for a protein structure have Z scores higher than or equal to fixed thresholds, where thresholds are defined as the Z score for the corresponding native sequence (case 1) or much greater Z score (case 2). An exploring walk simulation suggested that the neutral sequences mapped into the sequence space were connected with each other through straight neutral paths and formed an inherent neutral network over the sequence space. Through another exploring walk simulation, we investigated contiguous regions between or among the neutral networks for the distinct protein structures and obtained the following results. The closest approach distance between the two neutral networks ranged from 5 to 29 on the Hamming distance scale, showing a linear increase against the threshold values. The sequences located at the "interchange" regions between the two neutral networks have intermediate sequence-profile-scores for both corresponding structures. Introducing a "ball" in the sequence space that contains at least one neutral sequence for each of the four structures, we found that the minimal radius of the ball that is centered at an arbitrary position ranged from 35 to 50, while the minimal radius of the ball that is centered at a certain special position ranged from 20 to 30, in the Hamming distance scale. The relatively small Hamming distances (5-30) may support an evolution mechanism by transferring from a network for a structure to another network for a more beneficial structure via the interchange regions.

Knowledge-based potential functions in protein design

Predicting protein sequences that fold into specific native three-dimensional structures is a problem of great potential complexity. Although the complete solution is ultimately rooted in understanding the physical chemistry underlying the complex interactions between amino acid residues that determine protein stability, recent work shows that empirical information about these first principles is embedded in the statistics of protein sequence and structure databases. This review focuses on the use of 'knowledge-based' potentials derived from these databases in designing proteins. In addition, the data suggest how the study of these empirical potentials might impact our fundamental understanding of the energetic principles of protein structure.

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.

Knowledge-based potential defined for a rotamer library to design protein sequences

A knowledge-based potential for a rotamer library was developed to design protein sequences. Protein side-chain conformations are represented by 56 templates. Each of their fitness to a given structural site-environment is evaluated by a combined function of the three knowledge-based terms, i.e. two-body side-chain packing, one-body hydration and local conformation. The number of matches between the native sequence and the structural site-environment in the database and that of the virtually settled mismatches, counted in advance, were transformed into the energy scores. In the best-14 test (assessment for the reproduction ability of the native rotamer on its structural site within a quarter of 56 fitness rank positions), the structural stability analysis on mutants of human and T4 lysozymes and the inverse-folding search by a structure profile against the sequence database, this function performs better than the function deduced with the conventional normalization and our previously developed function. Targeting various structural motifs, de novo sequence design was conducted with the function. The sequences thus obtained exhibit reasonable molecular masses and hydrophobic/hydrophilic patterns similar to the native sequences of the target and act as if they were the homologs to the target proteins in BLASTP search. This significant improvement is discussed in terms of the reference state for normalization and the crucial role of short-range repulsion to prohibit residue bumps.

Designing protein beta-sheet surfaces by Z-score optimization

Studies of lattice models of proteins have suggested that the appropriate energy expression for protein design may include nonthermodynamic terms to accommodate negative design concerns. One method, developed in lattice model studies, maximizes a quantity known as the " Z-score," which compares the lowest energy sequence whose ground state structure is the target structure to an ensemble of random sequences. Here we show that, in certain circumstances, the technique can be applied to real proteins. The resulting energy expression is used to design the beta-sheet surfaces of two real proteins. We find experimentally that the designed proteins are stable and well folded, and in one case is even more thermostable than the wild type.