Expression and purification of recombinant hemoglobin in Escherichia coli

Imperfect symmetry facilitated the evolution of specificity and high-order stoichiometry in vertebrate hemoglobin

Many proteins form paralogous multimers - molecular complexes in which evolutionarily related proteins are arranged into specific quaternary structures. Little is known about the mechanisms by which they acquired their stoichiometry (the number of total subunits in the complex) and heterospecificity (the preference of subunits for their paralogs rather than other copies of the same protein). Here we use ancestral protein reconstruction and biochemical experiments to study historical increases in stoichiometry and specificity during the evolution of vertebrate hemoglobin (Hb), a α2β2 heterotetramer that evolved from a homodimeric ancestor after a gene duplication. We show that the mechanisms for this evolutionary transition was simple. One hydrophobic substitution in subunit β after the gene duplication was sufficient to cause the ancestral dimer to homotetramerize with high affinity across a new interface. During this same interval, a single-residue deletion in subunit α at the older interface conferred specificity for the heterotetrameric form and the trans-orientation of subunits within it. These sudden transitions in stoichiometry and specificity were possible because the interfaces in Hb are isologous - involving the same surface patch on interacting subunits, rotated 180° relative to each other - but the symmetry is slightly imperfect. This architecture amplifies the impacts of individual mutations on stoichiometry and specificity, especially in higher-order complexes, and allows single substitutions to differentially affect heteromeric vs homomeric interactions. Many multimers are isologous, and symmetry in proteins is always imperfect; our findings therefore suggest that elaborate and specific molecular complexes may often evolve via simple genetic and physical mechanisms.

High-level expression of leghemoglobin in Kluyveromyces marxianus by remodeling the heme metabolism pathway

Soy leghemoglobin, when bound to heme, imparts a meat-like color and flavor and can serve as a substitute for animal-derived proteins. Enhancing cellular heme synthesis improves the recombinant expression of leghemoglobin in yeast. To achieve high-level expression of leghemoglobin A (LBA) in Kluyveromyces marxianus, a food-safe yeast, large-scale heme synthesis modules were transferred into K. marxianus using yeast artificial chromosomes (KmYACs). These modules contained up to 8 native and heterologous genes to promote the supply of heme precursors and downstream synthesis. Next, eight genes inhibiting heme or LBA synthesis were individually or combinatorially deleted, with the lsc1Δssn3Δ mutant yielding the best results. Subsequently, heme synthesis modules were combined with the lsc1Δssn3Δ mutant. In the resulting strains, the module genes were all actively expressed. Among these module genes, heterologous S. cerevisiae genes in the downstream heme synthesis pathway significantly enhanced the expression of their counterparts in K. marxianus, resulting in high heme content and LBA yield. After optimizing the medium recipe by adjusting the concentrations of glucose, glycine, and FeSO4·7H2O, a heme content of 66.32 mg/L and an intracellular LBA titer of 7.27 g/L were achieved in the engineered strain in a 5 L fermentor. This represents the highest intracellular expression of leghemoglobin in microorganisms to date. The leghemoglobin produced by K. marxianus can be utilized as a safe ingredient for plant-based protein products.

Evolution and molecular basis of a novel allosteric property of crocodilian hemoglobin

The extraordinary breath-hold diving capacity of crocodilians has been ascribed to a unique mode of allosterically regulating hemoglobin (Hb)-oxygenation in circulating red blood cells. We investigated the origin and mechanistic basis of this novel biochemical phenomenon by performing directed mutagenesis experiments on resurrected ancestral Hbs. Comparisons of Hb function between the common ancestor of archosaurs (the group that includes crocodilians and birds) and the last common ancestor of modern crocodilians revealed that regulation of Hb-O2 affinity via allosteric binding of bicarbonate ions represents a croc-specific innovation that evolved in combination with the loss of allosteric regulation by ATP binding. Mutagenesis experiments revealed that evolution of the novel allosteric function in crocodilians and the concomitant loss of ancestral function were not mechanistically coupled and were caused by different sets of substitutions. The gain of bicarbonate sensitivity in crocodilian Hb involved the direct effect of few amino acid substitutions at key sites in combination with indirect effects of numerous other substitutions at structurally disparate sites. Such indirect interaction effects suggest that evolution of the novel protein function was conditional on neutral mutations that produced no adaptive benefit when they first arose but that contributed to a permissive background for subsequent function-altering mutations at other sites. Due to the context dependence of causative substitutions, the unique allosteric properties of crocodilian Hb cannot be easily transplanted into divergent homologs of other species.

Introducing Negatively Charged Residues on the Surface of Fetal Hemoglobin Improves Yields in Escherichia coli

Fetal hemoglobin (HbF) has been developed into an important alternative protein for oxygen therapeutics. Such applications require extensive amounts of proteins, which only can be achieved via recombinant means. However, the expression of vertebrate hemoglobins in heterologous hosts is far from trivial. There are several issues that need to be dealt with. These include, among others, the solubility of the globin chains, equimolar expression of the globin chains, and access to high levels of free heme. In this study, we examined the impact of introducing negative charges on the surface of HbF. Three different HbF mutants were examined, carrying four additional negative charges on the α-subunit (rHbFα4), two additional negative charges on the γ-subunit (rHbFγ2) or a combination of these (rHbFα4/γ2). The increase in negative surface charge in these HbF mutants required the development of an alternate initial capture step in the downstream purification procedures. For the rHbFα4 mutant, we achieved a significantly enhanced yield of purified HbF with no apparent adverse effects on Hb functionality. However, the presence of non-functional Hb portions in the rHbFγ2 and rHbFα4/γ2 samples reduced the yields significantly for those mutants and indicated an imbalanced expression/association of globin chains. Furthermore, the autoxidation studies indicated that the rHbFγ2 and rHbFα4/γ2 mutants also were less oxidatively stable than rHbFα4 and wt rHbF. The study further verified the need for an improved flask culture protocol by optimizing cultivation parameters to enable yield-improving qualities of surface-located mutations.

Evolved increases in hemoglobin-oxygen affinity and the Bohr effect coincided with the aquatic specialization of penguins

In diving birds like penguins, physiologic considerations suggest that increased hemoglobin (Hb)-O₂ affinity may improve pulmonary O₂ extraction and enhance dive capacity. We integrated experimental tests on whole-blood and native Hbs of penguins with protein engineering experiments on reconstructed ancestral Hbs. The experiments involving ancestral protein resurrection enabled us to test for evolved changes in Hb function in the stem lineage of penguins after divergence from their closest nondiving relatives. We demonstrate that penguins evolved an increased Hb-O₂ affinity in conjunction with a greatly augmented Bohr effect (i.e., reduction in Hb-O₂ affinity at low pH) that should maximize pulmonary O₂ extraction without compromising O₂ delivery at systemic capillaries.

Dive capacities of air-breathing vertebrates are dictated by onboard O₂ stores, suggesting that physiologic specialization of diving birds such as penguins may have involved adaptive changes in convective O₂ transport. It has been hypothesized that increased hemoglobin (Hb)-O₂ affinity improves pulmonary O₂ extraction and enhances the capacity for breath-hold diving. To investigate evolved changes in Hb function associated with the aquatic specialization of penguins, we integrated comparative measurements of whole-blood and purified native Hb with protein engineering experiments based on site-directed mutagenesis. We reconstructed and resurrected ancestral Hb representing the common ancestor of penguins and the more ancient ancestor shared by penguins and their closest nondiving relatives (order Procellariiformes, which includes albatrosses, shearwaters, petrels, and storm petrels). These two ancestors bracket the phylogenetic interval in which penguin-specific changes in Hb function would have evolved. The experiments revealed that penguins evolved a derived increase in Hb-O₂ affinity and a greatly augmented Bohr effect (i.e., reduced Hb-O₂ affinity at low pH). Although an increased Hb-O₂ affinity reduces the gradient for O₂ diffusion from systemic capillaries to metabolizing cells, this can be compensated by a concomitant enhancement of the Bohr effect, thereby promoting O₂ unloading in acidified tissues. We suggest that the evolved increase in Hb-O₂ affinity in combination with the augmented Bohr effect maximizes both O₂ extraction from the lungs and O₂ unloading from the blood, allowing penguins to fully utilize their onboard O₂ stores and maximize underwater foraging time.

Effect of NH2-terminal acetylation on the oxygenation properties of vertebrate haemoglobin

In vertebrate haemoglobin (Hb), the NH2-terminal residues of the α- and β-chain subunits are thought to play an important role in the allosteric binding of protons (Bohr effect), CO2 (as carbamino derivatives), chloride ions, and organic phosphates. Accordingly, acetylation of the α- and/or β-chain NH2-termini may have significant effects on the oxygenation properties of Hb. Here we investigate the effect of NH2-terminal acetylation by using a newly developed expression plasmid system that enables us to compare recombinantly expressed Hbs that are structurally identical except for the presence or absence of NH2-terminal acetyl groups. Experiments with native and recombinant Hbs of representative vertebrates reveal that NH2-terminal acetylation does not impair the Bohr effect, nor does it significantly diminish responsiveness to allosteric cofactors, such as chloride ions or organic phosphates. These results suggest that observed variation in the oxygenation properties of vertebrate Hbs is principally explained by amino acid divergence in the constituent globin chains rather than post-translational modifications of the globin chain NH2-termini.

Recent Advances in the Microbial Synthesis of Hemoglobin

Hemoglobin is a cofactor-containing protein with heme that plays important roles in transporting and storing oxygen. Hemoglobins have been widely applied as acellular oxygen carriers, bioavailable iron-supplying agents, and food-grade coloring and flavoring agents. To meet increasing demands and overcome the drawbacks of chemical extraction, the biosynthesis of hemoglobin has become an attractive alternative. Several hemoglobins have recently been synthesized by various microorganisms through metabolic engineering and synthetic biology. In this review, we summarize the novel strategies that have been used to biosynthesize hemoglobin. These strategies can also serve as references for producing other heme-binding proteins.

Synthesis of Recombinant Human Hemoglobin With NH2 -Terminal Acetylation in Escherichia coli

The development of new technologies for the efficient expression of recombinant hemoglobin (rHb) is of interest for experimental studies of protein biochemistry and the development of cell-free blood substitutes in transfusion medicine. Expression of rHb in Escherichia coli host cells has numerous advantages, but one disadvantage of using prokaryotic systems to express eukaryotic proteins is that they are incapable of performing post-translational modifications such as NH₂ -terminal acetylation. One possible solution is to coexpress additional enzymes that can perform the necessary modifications in the host cells. Here, we report a new method for synthesizing human rHb with proper NH₂ -terminal acetylation. Mass spectrometry experiments involving native and recombinant human Hb confirmed the efficacy of the new technique in producing correctly acetylated globin chains. Finally, functional experiments provided insights into the effects of NH₂ -terminal acetylation on O₂ binding properties. © 2020 Wiley Periodicals LLC. Basic Protocol 1: Gene synthesis and cloning the cassette to the expression plasmid Basic Protocol 2: Selection of E. coli expression strains for coexpression Basic Protocol 3: Large-scale recombinant hemoglobin expression and purification Support Protocol 1: Measuring O₂ equilibration curves Support Protocol 2: Mass spectrometry to confirm NH₂ -terminal acetylation.

Biochemical pedomorphosis and genetic assimilation in the hypoxia adaptation of Tibetan antelope

Developmental shifts in stage-specific gene expression can provide a ready mechanism of phenotypic change by altering the rate or timing of ontogenetic events. We found that the high-altitude Tibetan antelope (Panthelops hodgsonii) has evolved an adaptive increase in blood-O2 affinity by truncating the ancestral ontogeny of globin gene expression such that a high-affinity juvenile hemoglobin isoform (isoHb) completely supplants the lower-affinity isoHb that is expressed in the adult red blood cells of other bovids. This juvenilization of blood properties represents a canalization of an acclimatization response to hypoxia that has been well documented in adult goats and sheep. We also found the genomic mechanism underlying this regulatory isoHb switch, revealing how a reversible acclimatization response became genetically assimilated as an irreversible adaptation to chronic hypoxia.

Origin of complexity in haemoglobin evolution

Most proteins associate into multimeric complexes with specific architectures1,2, which often have functional properties such as cooperative ligand binding or allosteric regulation3. No detailed knowledge is available about how any multimer and its functions arose during evolution. Here we use ancestral protein reconstruction and biophysical assays to elucidate the origins of vertebrate haemoglobin, a heterotetramer of paralogous α- and β-subunits that mediates respiratory oxygen transport and exchange by cooperatively binding oxygen with moderate affinity. We show that modern haemoglobin evolved from an ancient monomer and characterize the historical 'missing link' through which the modern tetramer evolved-a noncooperative homodimer with high oxygen affinity that existed before the gene duplication that generated distinct α- and β-subunits. Reintroducing just two post-duplication historical substitutions into the ancestral protein is sufficient to cause strong tetramerization by creating favourable contacts with more ancient residues on the opposing subunit. These surface substitutions markedly reduce oxygen affinity and even confer cooperativity, because an ancient linkage between the oxygen binding site and the multimerization interface was already an intrinsic feature of the protein's structure. Our findings establish that evolution can produce new complex molecular structures and functions via simple genetic mechanisms that recruit existing biophysical features into higher-level architectures.

Adaptive Changes in Hemoglobin Function in High-Altitude Tibetan Canids Were Derived via Gene Conversion and Introgression

A key question in evolutionary biology concerns the relative importance of different sources of adaptive genetic variation, such as de novo mutations, standing variation, and introgressive hybridization. A corollary question concerns how allelic variants derived from these different sources may influence the molecular basis of phenotypic adaptation. Here, we use a protein-engineering approach to examine the phenotypic effect of putatively adaptive hemoglobin (Hb) mutations in the high-altitude Tibetan wolf that were selectively introgressed into the Tibetan mastiff, a high-altitude dog breed that is renowned for its hypoxia tolerance. Experiments revealed that the introgressed coding variants confer an increased Hb–O₂ affinity in conjunction with an enhanced Bohr effect. We also document that affinity-enhancing mutations in the β-globin gene of Tibetan wolf were originally derived via interparalog gene conversion from a tandemly linked β-globin pseudogene. Thus, affinity-enhancing mutations were introduced into the β-globin gene of Tibetan wolf via one form of intragenomic lateral transfer (ectopic gene conversion) and were subsequently introduced into the Tibetan mastiff genome via a second form of lateral transfer (introgression). Site-directed mutagenesis experiments revealed that the increased Hb–O₂ affinity requires a specific two-site combination of amino acid replacements, suggesting that the molecular underpinnings of Hb adaptation in Tibetan mastiff (involving mutations that arose in a nonexpressed gene and which originally fixed in Tibetan wolf) may be qualitatively distinct from functionally similar changes in protein function that could have evolved via sequential fixation of de novo mutations during the breed’s relatively short duration of residency at high altitude.

Emergence of a Chimeric Globin Pseudogene and Increased Hemoglobin Oxygen Affinity Underlie the Evolution of Aquatic Specializations in Sirenia

As limits on O₂ availability during submergence impose severe constraints on aerobic respiration, the oxygen binding globin proteins of marine mammals are expected to have evolved under strong evolutionary pressures during their land-to-sea transition. Here, we address this question for the order Sirenia by retrieving, annotating, and performing detailed selection analyses on the globin repertoire of the extinct Steller’s sea cow (Hydrodamalis gigas), dugong (Dugong dugon), and Florida manatee (Trichechus manatus latirostris) in relation to their closest living terrestrial relatives (elephants and hyraxes). These analyses indicate most loci experienced elevated nucleotide substitution rates during their transition to a fully aquatic lifestyle. While most of these genes evolved under neutrality or strong purifying selection, the rate of nonsynonymous/synonymous replacements increased in two genes (Hbz-T1 and Hba-T1) that encode the α-type chains of hemoglobin (Hb) during each stage of life. Notably, the relaxed evolution of Hba-T1 is temporally coupled with the emergence of a chimeric pseudogene (Hba-T2/Hbq-ps) that contributed to the tandemly linked Hba-T1 of stem sirenians via interparalog gene conversion. Functional tests on recombinant Hb proteins from extant and ancestral sirenians further revealed that the molecular remodeling of Hba-T1 coincided with increased Hb–O₂ affinity in early sirenians. Available evidence suggests that this trait evolved to maximize O₂ extraction from finite lung stores and suppress tissue O₂ offloading, thereby facilitating the low metabolic intensities of extant sirenians. In contrast, the derived reduction in Hb–O₂ affinity in (sub)Arctic Steller’s sea cows is consistent with fueling increased thermogenesis by these once colossal marine herbivores.

The role of mutation bias in adaptive molecular evolution: insights from convergent changes in protein function

An underexplored question in evolutionary genetics concerns the extent to which mutational bias in the production of genetic variation influences outcomes and pathways of adaptive molecular evolution. In the genomes of at least some vertebrate taxa, an important form of mutation bias involves changes at CpG dinucleotides: if the DNA nucleotide cytosine (C) is immediately 5' to guanine (G) on the same coding strand, then-depending on methylation status-point mutations at both sites occur at an elevated rate relative to mutations at non-CpG sites. Here, we examine experimental data from case studies in which it has been possible to identify the causative substitutions that are responsible for adaptive changes in the functional properties of vertebrate haemoglobin (Hb). Specifically, we examine the molecular basis of convergent increases in Hb-O₂ affinity in high-altitude birds. Using a dataset of experimentally verified, affinity-enhancing mutations in the Hbs of highland avian taxa, we tested whether causative changes are enriched for mutations at CpG dinucleotides relative to the frequency of CpG mutations among all possible missense mutations. The tests revealed that a disproportionate number of causative amino acid replacements were attributable to CpG mutations, suggesting that mutation bias can influence outcomes of molecular adaptation. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Chromatographic separation of hemoglobin variants using robust molecularly imprinted polymers

Devising a robust, efficient and cost effective hemoglobin (Hb) purification strategy is one of the key challenges in the development of Hb-based blood substitutes. The aim of this study was to use molecularly imprinted polymers (MIPs) as a novel and efficient chromatographic resin to selectively recognize and purify different Hb variants. The results showed that the Hb-MIP material developed here could selectively recognize and purify various Hb directly from either crude E. coli extracts or human body fluids, such as blood plasma and cerebrospinal fluid (CSF), in one-step. The dynamic binding capacity at 10% breakthrough was around 7.4 mg mL^-1resin for adult Hb (HbA) and fetal Hb (HbF). This chromatographic material also allowed identification of changes related to amino acid substitutions on the Hb protein surface. For instance, when an additional lysine residue was introduced, the HbA αY42K mutant eluted later in an Hb-MIP column than wildtype HbA. Additional negative charges on the protein surface, such as aspartate, mitigated the interaction between the protein and imprinted polymers, and therefore an αA19D-αA12D HbF mutant eluted earlier, at -2.7 column volumes compared to wildtype HbF.

Molecular basis of hemoglobin adaptation in the high-flying bar-headed goose

During the adaptive evolution of a particular trait, some selectively fixed mutations may be directly causative and others may be purely compensatory. The relative contribution of these two classes of mutation to adaptive phenotypic evolution depends on the form and prevalence of mutational pleiotropy. To investigate the nature of adaptive substitutions and their pleiotropic effects, we used a protein engineering approach to characterize the molecular basis of hemoglobin (Hb) adaptation in the high-flying bar-headed goose (Anser indicus), a hypoxia-tolerant species renowned for its trans-Himalayan migratory flights. To test the effects of observed substitutions on evolutionarily relevant genetic backgrounds, we synthesized all possible genotypic intermediates in the line of descent connecting the wildtype bar-headed goose genotype with the most recent common ancestor of bar-headed goose and its lowland relatives. Site-directed mutagenesis experiments revealed one major-effect mutation that significantly increased Hb-O2 affinity on all possible genetic backgrounds. Two other mutations exhibited smaller average effect sizes and less additivity across backgrounds. One of the latter mutations produced a concomitant increase in the autoxidation rate, a deleterious side-effect that was fully compensated by a second-site mutation at a spatially proximal residue. The experiments revealed three key insights: (i) subtle, localized structural changes can produce large functional effects; (ii) relative effect sizes of function-altering mutations may depend on the sequential order in which they occur; and (iii) compensation of deleterious pleiotropic effects may play an important role in the adaptive evolution of protein function.

Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau

When different species experience similar selection pressures, the probability of evolving similar adaptive solutions may be influenced by legacies of evolutionary history, such as lineage-specific changes in genetic background. Here we test for adaptive convergence in hemoglobin (Hb) function among high-altitude passerine birds that are native to the Qinghai-Tibet Plateau, and we examine whether convergent increases in Hb-O₂ affinity have a similar molecular basis in different species. We documented that high-altitude parid and aegithalid species from the Qinghai-Tibet Plateau have evolved derived increases in Hb-O₂ affinity in comparison with their closest lowland relatives in East Asia. However, convergent increases in Hb-O₂ affinity and convergence in underlying functional mechanisms were seldom attributable to the same amino acid substitutions in different species. Using ancestral protein resurrection and site-directed mutagenesis, we experimentally confirmed two cases in which parallel substitutions contributed to convergent increases in Hb-O₂ affinity in codistributed high-altitude species. In one case involving the ground tit (Parus humilis) and gray-crested tit (Lophophanes dichrous), parallel amino acid replacements with affinity-enhancing effects were attributable to nonsynonymous substitutions at a CpG dinucleotide, suggesting a possible role for mutation bias in promoting recurrent changes at the same site. Overall, most altitude-related changes in Hb function were caused by divergent amino acid substitutions, and a select few were caused by parallel substitutions that produced similar phenotypic effects on the divergent genetic backgrounds of different species.

Exploring three different expression systems for recombinant expression of globins: Escherichia coli, Pichia pastoris and Spodoptera frugiperda

Globins are among the best investigated proteins in biological and medical sciences and represent a prime tool for the study of the evolution of genes and the structure-function relationship of proteins. Here, we explore the recombinant expression of globins in three different expression systems: Escherichia coli, Pichia pastoris and the baculovirus infected Spodoptera frugiperda. We expressed two different human globin types in these three expression systems: I) the well-characterized neuroglobin and II) the uncharacterized, circular permutated globin domain of the large chimeric globin androglobin. It is clear from the literature that E.coli is the most used expression system for expression and purification of recombinant globins. However, the major disadvantage of E. coli is the formation of insoluble aggregates. We experienced that, for more complex multi-domain globins, like the chimeric globin androglobin, it is recommended to switch to a higher eukaryotic expression system.

Acclimation to prolonged hypoxia alters hemoglobin isoform expression and increases hemoglobin oxygen affinity and aerobic performance in a marine fish

Hemoglobin (Hb) multiplicity is common in fish, yet despite its ubiquitous nature, the functional significance is unclear. Here we explore the hypothesis that Hb multiplicity plays a role in hypoxia tolerance using the red drum (Sciaenops ocellatus). Red drum is an economically and ecologically important species native to coastal regions and estuaries of the Gulf of Mexico – habitats that routinely experience pronounced hypoxic events. Using a transcriptomic approach, we demonstrate that red drum red blood cells express 7 and 5 Hbα and Hbβ isoforms, respectively. Phylogenetic analysis grouped these isoforms into distinct isoHb clades, and provided evidence of lineage specific expression of particular isoHbs. In normoxia, three isoHbs predominated (Hbα-3.1, -3.2, and Hbβ-3.1). A three-week hypoxia acclimation (48 mmHg) resulted in significant up-regulation of Hbα-2, Hbα-3.2, and Hbβ-3.1, effectively switching the predominantly expressed isoforms. Changes in subunit expression were correlated with a decrease in non-stripped hemolysate P₅₀. Similarly, hypoxia acclimation resulted in a 20% reduction in whole animal critical oxygen threshold (P_crit). Hypoxia acclimation was not associated with changes in gill morphology, hematocrit, or relative ventricular mass. Overall, these data provide support for the hypothesis that Hb isoform switching can provide a physiological benefit to counteract environmental stress in fishes.

Stability-Mediated Epistasis Restricts Accessible Mutational Pathways in the Functional Evolution of Avian Hemoglobin

If the fitness effects of amino acid mutations are conditional on genetic background, then mutations can have different effects depending on the sequential order in which they occur during evolutionary transitions in protein function. A key question concerns the fraction of possible mutational pathways connecting alternative functional states that involve transient reductions in fitness. Here we examine the functional effects of multiple amino acid substitutions that contributed to an evolutionary transition in the oxygenation properties of avian hemoglobin (Hb). The set of causative changes included mutations at intradimer interfaces of the Hb tetramer. Replacements at such sites may be especially likely to have epistatic effects on Hb function since residues at intersubunit interfaces are enmeshed in networks of salt bridges and hydrogen bonds between like and unlike subunits; mutational reconfigurations of these atomic contacts can affect allosteric transitions in quaternary structure and the propensity for tetramer-dimer dissociation. We used ancestral protein resurrection in conjunction with a combinatorial protein engineering approach to synthesize genotypes representing the complete set of mutational intermediates in all possible forward pathways that connect functionally distinct ancestral and descendent genotypes. The experiments revealed that 1/2 of all possible forward pathways included mutational intermediates with aberrant functional properties because particular combinations of mutations promoted tetramer-dimer dissociation. The subset of mutational pathways with unstable intermediates may be selectively inaccessible, representing evolutionary roads not taken. The experimental results also demonstrate how epistasis for particular functional properties of proteins may be mediated indirectly by mutational effects on quaternary structural stability.

Predictable convergence in hemoglobin function has unpredictable molecular underpinnings

To investigate the predictability of genetic adaptation, we examined the molecular basis of convergence in hemoglobin function in comparisons involving 56 avian taxa that have contrasting altitudinal range limits. Convergent increases in hemoglobin-oxygen affinity were pervasive among high-altitude taxa, but few such changes were attributable to parallel amino acid substitutions at key residues. Thus, predictable changes in biochemical phenotype do not have a predictable molecular basis. Experiments involving resurrected ancestral proteins revealed that historical substitutions have context-dependent effects, indicating that possible adaptive solutions are contingent on prior history. Mutations that produce an adaptive change in one species may represent precluded possibilities in other species because of differences in genetic background.

Alteration of the α1β2/α2β1 subunit interface contributes to the increased hemoglobin-oxygen affinity of high-altitude deer mice

BACKGROUND

Deer mice (Peromyscus maniculatus) that are native to high altitudes in the Rocky Mountains have evolved hemoglobins with an increased oxygen-binding affinity relative to those of lowland conspecifics. To elucidate the molecular mechanisms responsible for the evolved increase in hemoglobin-oxygen affinity, the crystal structure of the highland hemoglobin variant was solved and compared with the previously reported structure for the lowland variant.

RESULTS

Highland hemoglobin yielded at least two crystal types, in which the longest axes were 507 and 230 Å. Using the smaller unit cell crystal, the structure was solved at 2.2 Å resolution. The asymmetric unit contained two tetrameric hemoglobin molecules.

CONCLUSIONS

The analyses revealed that αPro50 in the highland hemoglobin variant promoted a stable interaction between αHis45 and heme that was not seen in the αHis50 lowland variant. The αPro50 mutation also altered the nature of atomic contacts at the α1β2/α2β1 intersubunit interfaces. These results demonstrate how affinity-altering changes in intersubunit interactions can be produced by mutations at structurally remote sites.

Convergent Evolution of Hemoglobin Function in High-Altitude Andean Waterfowl Involves Limited Parallelism at the Molecular Sequence Level

A fundamental question in evolutionary genetics concerns the extent to which adaptive phenotypic convergence is attributable to convergent or parallel changes at the molecular sequence level. Here we report a comparative analysis of hemoglobin (Hb) function in eight phylogenetically replicated pairs of high- and low-altitude waterfowl taxa to test for convergence in the oxygenation properties of Hb, and to assess the extent to which convergence in biochemical phenotype is attributable to repeated amino acid replacements. Functional experiments on native Hb variants and protein engineering experiments based on site-directed mutagenesis revealed the phenotypic effects of specific amino acid replacements that were responsible for convergent increases in Hb-O₂ affinity in multiple high-altitude taxa. In six of the eight taxon pairs, high-altitude taxa evolved derived increases in Hb-O₂ affinity that were caused by a combination of unique replacements, parallel replacements (involving identical-by-state variants with independent mutational origins in different lineages), and collateral replacements (involving shared, identical-by-descent variants derived via introgressive hybridization). In genome scans of nucleotide differentiation involving high- and low-altitude populations of three separate species, function-altering amino acid polymorphisms in the globin genes emerged as highly significant outliers, providing independent evidence for adaptive divergence in Hb function. The experimental results demonstrate that convergent changes in protein function can occur through multiple historical paths, and can involve multiple possible mutations. Most cases of convergence in Hb function did not involve parallel substitutions and most parallel substitutions did not affect Hb-O₂ affinity, indicating that the repeatability of phenotypic evolution does not require parallelism at the molecular level.

The convergent evolution of similar traits in different species could be due to repeated changes at the genetic level or different changes that produce the same phenotypic effect. To investigate the extent to which convergence in phenotype is caused by repeated mutations, we investigated the molecular basis of convergent changes in the oxygenation properties of hemoglobin (Hb) in eight pairs of high- and low-altitude waterfowl taxa from the Andes. The results revealed that convergent increases in Hb-O₂ affinity in highland taxa involved a combination of unique and repeated amino acid replacements. However, convergent changes in Hb function generally did not involve parallel substitutions, indicating that repeatability in the evolution of protein function does not require repeatability at the sequence level.

Contribution of a mutational hot spot to hemoglobin adaptation in high-altitude Andean house wrens

A key question in evolutionary genetics is why certain mutations or certain types of mutation make disproportionate contributions to adaptive phenotypic evolution. In principle, the preferential fixation of particular mutations could stem directly from variation in the underlying rate of mutation to function-altering alleles. However, the influence of mutation bias on the genetic architecture of phenotypic evolution is difficult to evaluate because data on rates of mutation to function-altering alleles are seldom available. Here, we report the discovery that a single point mutation at a highly mutable site in the β(A)-globin gene has contributed to an evolutionary change in hemoglobin (Hb) function in high-altitude Andean house wrens (Troglodytes aedon). Results of experiments on native Hb variants and engineered, recombinant Hb mutants demonstrate that a nonsynonymous mutation at a CpG dinucleotide in the β(A)-globin gene is responsible for an evolved difference in Hb-O2 affinity between high- and low-altitude house wren populations. Moreover, patterns of genomic differentiation between high- and low-altitude populations suggest that altitudinal differentiation in allele frequencies at the causal amino acid polymorphism reflects a history of spatially varying selection. The experimental results highlight the influence of mutation rate on the genetic basis of phenotypic evolution by demonstrating that a large-effect allele at a highly mutable CpG site has promoted physiological differentiation in blood O2 transport capacity between house wren populations that are native to different elevations.

Oxygenation properties and isoform diversity of snake hemoglobins

Available data suggest that snake hemoglobins (Hbs) are characterized by a combination of unusual structural and functional properties relative to the Hbs of other amniote vertebrates, including oxygenation-linked tetramer-dimer dissociation. However, standardized comparative data are lacking for snake Hbs, and the Hb isoform composition of snake red blood cells has not been systematically characterized. Here we present the results of an integrated analysis of snake Hbs and the underlying α- and β-type globin genes to characterize 1) Hb isoform composition of definitive erythrocytes, and 2) the oxygenation properties of isolated isoforms as well as composite hemolysates. We used species from three families as subjects for experimental studies of Hb function: South American rattlesnake, Crotalus durissus (Viperidae); Indian python, Python molurus (Pythonidae); and yellow-bellied sea snake, Pelamis platura (Elapidae). We analyzed allosteric properties of snake Hbs in terms of the Monod-Wyman-Changeux model and Adair four-step thermodynamic model. Hbs from each of the three species exhibited high intrinsic O2 affinities, low cooperativities, small Bohr factors in the absence of phosphates, and high sensitivities to ATP. Oxygenation properties of the snake Hbs could be explained entirely by allosteric transitions in the quaternary structure of intact tetramers, suggesting that ligation-dependent dissociation of Hb tetramers into αβ-dimers is not a universal feature of snake Hbs. Surprisingly, the major Hb isoform of the South American rattlesnake is homologous to the minor HbD of other amniotes and, contrary to the pattern of Hb isoform differentiation in birds and turtles, exhibits a lower O2 affinity than the HbA isoform.

Intraspecific polymorphism, interspecific divergence, and the origins of function-altering mutations in deer mouse hemoglobin

Major challenges for illuminating the genetic basis of phenotypic evolution are to identify causative mutations, to quantify their functional effects, to trace their origins as new or preexisting variants, and to assess the manner in which segregating variation is transduced into species differences. Here, we report an experimental analysis of genetic variation in hemoglobin (Hb) function within and among species of Peromyscus mice that are native to different elevations. A multilocus survey of sequence variation in the duplicated HBA and HBB genes in Peromyscus maniculatus revealed that function-altering amino acid variants are widely shared among geographically disparate populations from different elevations, and numerous amino acid polymorphisms are also shared with closely related species. Variation in Hb-O₂ affinity within and among populations of P. maniculatus is attributable to numerous amino acid mutations that have individually small effects. One especially surprising feature of the Hb polymorphism in P. maniculatus is that an appreciable fraction of functional standing variation in the two transcriptionally active HBA paralogs is attributable to recurrent gene conversion from a tandemly linked HBA pseudogene. Moreover, transpecific polymorphism in the duplicated HBA genes is not solely attributable to incomplete lineage sorting or introgressive hybridization; instead, it is mainly attributable to recurrent interparalog gene conversion that has occurred independently in different species. Partly as a result of concerted evolution between tandemly duplicated globin genes, the same amino acid changes that contribute to variation in Hb function within P. maniculatus also contribute to divergence in Hb function among different species of Peromyscus. In the case of function-altering Hb mutations in Peromyscus, there is no qualitative or quantitative distinction between segregating variants within species and fixed differences between species.

Epistasis constrains mutational pathways of hemoglobin adaptation in high-altitude pikas

A fundamental question in evolutionary genetics concerns the roles of mutational pleiotropy and epistasis in shaping trajectories of protein evolution. This question can be addressed most directly by using site-directed mutagenesis to explore the mutational landscape of protein function in experimentally defined regions of sequence space. Here, we evaluate how pleiotropic trade-offs and epistatic interactions influence the accessibility of alternative mutational pathways during the adaptive evolution of hemoglobin (Hb) function in high-altitude pikas (Mammalia: Lagomorpha). By combining ancestral protein resurrection with a combinatorial protein-engineering approach, we examined the functional effects of sequential mutational steps in all possible pathways that produced an increased Hb–O₂ affinity. These experiments revealed that the effects of mutations on Hb–O₂ affinity are highly dependent on the temporal order in which they occur: Each of three β-chain substitutions produced a significant increase in Hb–O₂ affinity on the ancestral genetic background, but two of these substitutions produced opposite effects when they occurred as later steps in the pathway. The experiments revealed pervasive epistasis for Hb–O₂ affinity, but affinity-altering mutations produced no significant pleiotropic trade-offs. These results provide insights into the properties of adaptive substitutions in naturally evolved proteins and suggest that the accessibility of alternative mutational pathways may be more strongly constrained by sign epistasis for positively selected biochemical phenotypes than by antagonistic pleiotropy.

Integrating evolutionary and functional tests of adaptive hypotheses: a case study of altitudinal differentiation in hemoglobin function in an Andean Sparrow, Zonotrichia capensis

In air-breathing vertebrates, the physiologically optimal blood-O2 affinity is jointly determined by the prevailing partial pressure of atmospheric O2, the efficacy of pulmonary O2 transfer, and internal metabolic demands. Consequently, genetic variation in the oxygenation properties of hemoglobin (Hb) may be subject to spatially varying selection in species with broad elevational distributions. Here we report the results of a combined functional and evolutionary analysis of Hb polymorphism in the rufous-collared sparrow (Zonotrichia capensis), a species that is continuously distributed across a steep elevational gradient on the Pacific slope of the Peruvian Andes. We integrated a population genomic analysis that included all postnatally expressed Hb genes with functional studies of naturally occurring Hb variants, as well as recombinant Hb (rHb) mutants that were engineered through site-directed mutagenesis. We identified three clinally varying amino acid polymorphisms: Two in the α(A)-globin gene, which encodes the α-chain subunits of the major HbA isoform, and one in the α(D)-globin gene, which encodes the α-chain subunits of the minor HbD isoform. We then constructed and experimentally tested single- and double-mutant rHbs representing each of the alternative α(A)-globin genotypes that predominate at different elevations. Although the locus-specific patterns of altitudinal differentiation suggested a history of spatially varying selection acting on Hb polymorphism, the experimental tests demonstrated that the observed amino acid mutations have no discernible effect on respiratory properties of the HbA or HbD isoforms. These results highlight the importance of experimentally validating the hypothesized effects of genetic changes in protein function to avoid the pitfalls of adaptive storytelling.

Expression, stabilization and purification of membrane proteins via diverse protein synthesis systems and detergents involving cell-free associated with self-assembly peptide surfactants

G-protein coupled receptors (GPCRs) are involved in regulating most of physiological actions and metabolism in the bodies, which have become most frequently addressed therapeutic targets for various disorders and diseases. Purified GPCR-based drug discoveries have become routine that approaches to structural study, novel biophysical and biochemical function analyses. However, several bottlenecks that GPCR-directed drugs need to conquer the problems including overexpression, solubilization, and purification as well as stabilization. The breakthroughs are to obtain efficient protein yield and stabilize their functional conformation which are both urgently requiring of effective protein synthesis system methods and optimal surfactants. Cell-free protein synthesis system is superior to the high yields and post-translation modifications, and early signs of self-assembly peptide detergents also emerged to superiority in purification of membrane proteins. We herein focus several predominant protein synthesis systems and surfactants involving the novel peptide detergents, and uncover the advantages of cell-free protein synthesis system with self-assembling peptide detergents in purification of functional GPCRs. This review is useful to further study in membrane proteins as well as the new drug exploration.

Repeated elevational transitions in hemoglobin function during the evolution of Andean hummingbirds

Animals that sustain high levels of aerobic activity under hypoxic conditions (e.g., birds that fly at high altitude) face the physiological challenge of jointly optimizing blood-O2 affinity for O2 loading in the pulmonary circulation and O2 unloading in the systemic circulation. At high altitude, this challenge is especially acute for small endotherms like hummingbirds that have exceedingly high mass-specific metabolic rates. Here we report an experimental analysis of hemoglobin (Hb) function in South American hummingbirds that revealed a positive correlation between Hb-O2 affinity and native elevation. Protein engineering experiments and ancestral-state reconstructions revealed that this correlation is attributable to derived increases in Hb-O2 affinity in highland lineages, as well as derived reductions in Hb-O2 affinity in lowland lineages. Site-directed mutagenesis experiments demonstrated that repeated evolutionary transitions in biochemical phenotype are mainly attributable to repeated amino acid replacements at two epistatically interacting sites that alter the allosteric regulation of Hb-O2 affinity. These results demonstrate that repeated changes in biochemical phenotype involve parallelism at the molecular level, and that mutations with indirect, second-order effects on Hb allostery play key roles in biochemical adaptation.

Epistasis among adaptive mutations in deer mouse hemoglobin

Epistatic interactions between mutant sites in the same protein can exert a strong influence on pathways of molecular evolution. We performed protein engineering experiments that revealed pervasive epistasis among segregating amino acid variants that contribute to adaptive functional variation in deer mouse hemoglobin (Hb). Amino acid mutations increased or decreased Hb-O2 affinity depending on the allelic state of other sites. Structural analysis revealed that epistasis for Hb-O2 affinity and allosteric regulatory control is attributable to indirect interactions between structurally remote sites. The prevalence of sign epistasis for fitness-related biochemical phenotypes has important implications for the evolutionary dynamics of protein polymorphism in natural populations.

Deer mouse hemoglobin exhibits a lowered oxygen affinity owing to mobility of the E helix

The deer mouse, Peromyscus maniculatus, exhibits altitude-associated variation in hemoglobin oxygen affinity. To examine the structural basis of this functional variation, the structure of the hemoglobin was solved. Recombinant hemoglobin was expressed in Escherichia coli and was purified by ion-exchange chromatography. Recombinant hemoglobin was crystallized by the hanging-drop vapor-diffusion method using polyethylene glycol as a precipitant. The obtained orthorhombic crystal contained two subunits in the asymmetric unit. The refined structure was interpreted as the aquo-met form. Structural comparisons were performed among hemoglobins from deer mouse, house mouse and human. In contrast to human hemoglobin, deer mouse hemoglobin lacks the hydrogen bond between α1Trp14 in the A helix and α1Thr67 in the E helix owing to the Thr67Ala substitution. In addition, deer mouse hemoglobin has a unique hydrogen bond at the α1β1 interface between residues α1Cys34 and β1Ser128.

Crystallization of Chlorella deoxyuridine triphosphatase

Deoxyuridine triphosphatase (dUTPase) is a ubiquitous enzyme that has been widely studied owing to its function and evolutionary significance. The gene coding for the dUTPase from the Chlorella alga was codon-optimized and synthesized. The synthetic gene was expressed in Escherichia coli and recombinant core Chlorella dUTPase (chdUTPase) was purified. Crystallization of chdUTPase was performed by the repetitive hanging-drop vapor-diffusion method at 298 K with ammonium sulfate as the precipitant. In the presence of 2'-deoxyuridine-5'-[(α,β)-imido]triphosphate and magnesium, the enzyme produced die-shaped hexagonal R3 crystals with unit-cell parameters a = b = 66.9, c = 93.6 Å, γ = 120°. X-ray diffraction data for chdUTPase were collected to 1.6 Å resolution. The crystallization of chdUTPase with manganese resulted in very fragile clusters of needles.

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.