Plant hemoglobins: a molecular fossil record for the evolution of oxygen transport

The evolution of oxygen transport hemoglobins occurred on at least two independent occasions. The earliest event led to myoglobin and red blood cell hemoglobin in animals. In plants, oxygen transport "leghemoglobins" evolved much more recently. In both events, pentacoordinate heme sites capable of inert oxygen transfer evolved from hexacoordinate hemoglobins that have unrelated functions. High sequence homology between hexacoordinate and pentacoordinate hemoglobins in plants has poised them for potential structural analysis leading to a molecular understanding of this important evolutionary event. However, the lack of a plant hexacoordinate hemoglobin structure in the exogenously ligand-bound form has prevented such comparison. Here we report the crystal structure of the cyanide-bound hexacoordinate hemoglobin from barley. This presents the first opportunity to examine conformational changes in plant hexacoordinate hemoglobins upon exogenous ligand binding, and reveals structural mechanisms for stabilizing the high-energy pentacoordinate heme conformation critical to the evolution of reversible oxygen binding hemoglobins.     

The role of dynamics in allosteric regulation

The biomolecular conformational changes often associated with allostery are, by definition, dynamic processes. Recent publications have disclosed the role of pre-existing equilibria of conformational substates in this process. In addition, the role of dynamics as an entropic carrier of free energy of allostery has been investigated. Recent work thus shows that dynamics is pivotal to allostery, and that it constitutes much more than just the move from the 'T'-state to the 'R'-state. Emerging computational studies have described the actual pathways of allosteric change.     

Contrasting roles of dynamics in protein allostery: NMR and structural studies of CheY and the third PDZ domain from PSD-95

Allosteric regulation is a ubiquitous phenomenon exploited in biological processes to control cells in a myriad of ways. It is also of emerging interest in the design of functional proteins and therapeutics. Even though allostery was proposed over 50 years ago and has been studied intensively from a structural perspective, many key details of allosteric mechanisms remain mysterious. Over the last decade significant attention has been paid to the “dynamic component” of allostery, as opposed to the analysis of rigid structures. Nuclear magnetic resonance spectroscopy and its ability to detect conformationally dynamic processes at atomic resolution have played an important role in expanding our understanding of allosteric mechanisms and opening up new questions. This article focuses on work that highlights how protein dynamics can factor into allosteric processes in distinct ways. Two cases are contrasted. The first considers the “traditionally allosteric” protein CheY, which undergoes a conformational change as a key element of its allostery. The second considers the more rarely observed “dynamic allostery” in a PDZ domain, in which allosteric behavior arises from changes in internal structural dynamics. Interestingly, the dynamic processes in these two contrasting examples occur on different timescales. In the case of the PDZ domain, subsequent experimental and computational work is reviewed to reveal a more complete picture of this interesting case of allostery.     

Penta- and hexa-coordinate ferric hemoglobins display distinct pH titration profiles measured by Soret peak shifts

Hemoglobins with diverse characteristics have been identified in all kingdoms of life. Their ubiquitous presence indicates that these proteins play important roles in physiology, though function for all hemoglobins are not yet established with certainty. Their physiological role may depend on their ability to bind ligands, which in turn is dictated by their heme chemistry. However, we have an incomplete understanding of the mechanism of ligand binding for these newly discovered hemoglobins and the measurement of their kinetic parameters depend on their coordination at the heme iron. To gain insights into their functional role, it is important to categorize the new hemoglobins into either penta- or hexa-coordinated varieties. We demonstrate that simple pH titration and absorbance measurements can determine the coordination state of heme iron atom in ferric hemoglobins, thus providing unambiguous information about the classification of new globins. This method is rapid, sensitive and requires low concentration of protein. Penta- and hexa-coordinate hemoglobins displayed distinct pH titration profiles as observed in a variety of hemoglobins. The pentacoordinate distal histidine mutant proteins of hexacoordinate hemoglobins and ligand-bound hexacoordinate forms of pentacoordinate hemoglobins reverse the pH titration profiles, thus validating the sensitivity of this spectroscopic technique.     

Allostery: An Overview of Its History,Concepts, Methods,and Applications

The concept of allostery has evolved in the past century. In this Editorial, we briefly overview the history of allostery, from the pre-allostery nomenclature era starting with the Bohr effect (1904) to the birth of allostery by Monod and Jacob (1961). We describe the evolution of the allostery concept, from a conformational change in a two-state model (1965, 1966) to dynamic allostery in the ensemble model (1999); from multi-subunit (1965) proteins to all proteins (2004). We highlight the current available methods to study allostery and their applications in studies of conformational mechanisms, disease, and allosteric drug discovery. We outline the challenges and future directions that we foresee. Altogether, this Editorial narrates the history of this fundamental concept in the life sciences, its significance, methodologies to detect and predict it, and its application in a broad range of living systems.     

Allostery and conformational free energy changes in human tryptophanyl‐tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of allostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side‐chains and then to compute the ligand‐induced population shift. Finally, we obtain the free‐energy landscape of the protein in equilibrium, characterizing the free‐energy minima accessed by the protein complexes. We have chosen human tryptophanyl‐tRNA synthetase (hTrpRS), a protein responsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Allostery in pharmacology: thermodynamics, evolution and design

This review focuses on basic models of allostery, the ambiguous application of the allosteric term in pharmacology illustrated by receptors, the role of thermodynamics in allosteric mechanisms, evolution and design of allostery. The initial step of ligand activation is closure of the agonist-binding cavity. Large entropy increases accompany the agonist-elicited conformational changes of pentameric ligand-gated ion channels due to cavity closure and rearrangement of transmembrane helices. The effects of point mutations on thermodynamic parameters of binding and function can reveal energetic coupling of neighbouring (and distant) amino acid residues in activation. High-order double-mutant cycle analysis and rate-equilibrium linear free-energy relationships can identify the trajectory and conformational spread of activation.Protein assembly and allostery can be deduced from colocalization and physicochemical principles. Molecular evolution has led from homooligomerization of protomers to heterotropic cooperativity and to allosteric regulation. Examples are discussed such as similar paths of protein (dis)assembly and evolution, irreversible evolution, statistical analysis of sequence homology revealing coevolution, different impacts of adaptation and evolution on hemoglobin, and the flagellar motor switch of bacteria. The driving force of dynamic allostery is associated with funnel-like free energy landscapes of protein binding and shifts in conformational fluctuations upon binding. Allostery can be designed based on our increasing knowledge of natural allosteric mechanisms and evolution. The allosteric principle has been applied for various bio/macro/molecular and signal transduction systems as well as in cognitive sciences.     

Structural and functional properties of class 1 plant hemoglobins

Nonsymbiotic class 1 plant hemoglobins are induced under hypoxia. Structurally they are protein dimers consisting of two identical subunits, each containing heme iron in a weak hexacoordinate state. The weak hexacoordination of heme-iron binding to the distal histidine results in an extremely high avidity to oxygen, with a dissociation constant in the nanomolar range. This low dissociation constant is due to rapid oxygen binding resulting in protein conformational changes that slow dissociation from the heme site. Class 1 hemoglobins are characterized by an increased rate of Fe3(+) reduction which is likely mediated by cysteine residue. This cysteine can form a reversible covalent bond between two monomers as shown by mass spectrometry analysis and, in addition to its structural role, prevents the molecule from autoxidation. The structural properties of class 1 hemoglobins allow them to serve as soluble electron transport proteins in the enzymatic system scavenging nitric oxide produced in low oxygen via reduction of nitrite. During oxygenation of nitric oxide to nitrate, oxidized ferric hemoglobin is formed (methemoglobin), which can be reduced by an associated reductase. The identified candidate for this reduction is monodehydroascorbate reductase. It is suggested that hemoglobin functions as a terminal electron acceptor during the hypoxic turnover of nitrogen, the process aided by its extremely high affinity for oxygen.     

Structural analysis in nonsymbiotic hemoglobins: What can we learn from inner cavities?

Plants contain three classes of hemoglobins which are not associated with nitrogen fixing bacteria, and have been accordingly termed nonsymbiotic hemoglobins. The function of nonsymbiotic hemoglobins is as yet mostly unknown. A NO dioxygenase activity has been proposed and demonstrated for some of them in vitro. In this context, a sound molecular mechanism that relates the structure with the biological activity is crucial to suggest a given physiological role. Insight into such a mechanism is now facilitated by recent progress made in both experimental and computational techniques. These studies have highlighted a number of key structural features implicated in the function of nonsymbiotic hemoglobins. The bis-histidyl hexacoordination of the heme in both its ferric and ferrous states provides a powerful and general tool to modulate reactivity, protein dynamics, and shape of the cavities. In addition, the specific arrangement of distal cavity residues provides effective protection against autoxidation. Inspection of the static crystal structures available for both liganded and unliganded states seems unsufficient to explain the function of these proteins. Function appears to be intimately linked with protein flexibility, which influences the dynamical behavior of inner cavities, capable of delivering apolar reactants to the reaction site, and removing charged reaction products. In this mini review, we demonstrate how the integration of information derived from experimental assays and computational studies is valuable and can shed light into the linkage between structural plasticity of nonsymbiotic hemoglobins and their biological role.     

Extending the New Generation of Structure Predictors to Account for Dynamics and Allostery

Recent progress in structure-prediction methods that rely on deep learning suggests that the atomic structure of almost any protein may soon be predictable directly from its amino acid sequence. This much-awaited revolution was driven by substantial improvements in the reliability of methods for inferring the spatial distances between amino acid pairs from an analysis of homologous sequences. Improved reliability has been accompanied, however, by a reduced ability to detect amino acid relationships that are not due to direct spatial contacts, such as those that arise from protein dynamics or allostery. Given the central importance of dynamics and allostery to protein activity, we argue that an important future advance would extend modeling beyond predicting a single static structure. Here, we briefly review some of the developments that have led to the remarkable recent achievement in structure prediction and speculate what methods and sources of information may be leveraged in the future to develop a modeling framework that addresses protein dynamics and allostery.     

Protein dynamics and allostery: an NMR view

Allostery, the process by which distant sites within a protein system are energetically coupled, is an efficient and ubiquitous mechanism for activity regulation. A purely mechanical view of allostery invoking only structural changes has developed over the decades as the classical view of the phenomenon. However, a fast growing list of examples illustrate the intimate link between internal motions over a wide range of time scales and function in protein-ligand interactions. Proteins respond to perturbations by redistributing their motions and they use fluctuating conformational states for binding and conformational entropy as a carrier of allosteric energy to modulate association with ligands. In several cases allosteric interactions proceed with minimal or no structural changes. We discuss emerging paradigms for the central role of protein dynamics in allostery.     

Highly reactive cysteine residues in rodent hemoglobins

Two hemoglobins with cysteine residues highly reactive toward electrophiles have been identified and characterized. Cys-125beta of guinea pig hemoglobin has a low pK(a) and forms conjugates with electrophiles more quickly than glutathione and several orders of magnitude more quickly than other protein thiols. This cysteine is capable of intercepting benzoquinone, a known carcinogenic metabolite, before other protein nucleophiles can be modified. Cys-13beta of mouse hemoglobin was observed to conjugate with electrophiles as quickly as glutathione. The structural basis of reactivity is different in the two hemoglobins and is analyzed in terms of hydrogen-bonding, solvent accessibility, and helix-dipole contributions. Complementing a previously characterized highly reactive cysteine in rat hemoglobin, identification of these cysteines suggests that the reactivity of these hemoglobins could represent a common function as a detoxification sink against carcinogens.     

Ligand binding and allostery can emerge simultaneously

A heterotropic allosteric effect involves an effector molecule that is distinct from the substrate or ligand of the protein. How heterotropic allostery originates is an unanswered question. We have previously created several heterotropic allosteric enzymes by recombining the genes for TEM1 beta-lactamase (BLA) and maltose binding protein (MBP) to create BLAs that are positively or negatively regulated by maltose. We show here that one of these engineered enzymes has approximately 10(6) M(-1) affinity for Zn(2+), a property that neither of the parental proteins possesses. Furthermore, Zn(2+) is a negative effector that noncompetitively switches off beta-lactam hydrolysis activity. Mutagenesis experiments indicate that the Zn(2+)-binding site does not involve a histidine or a cysteine, which is atypical of natural Zn(2+)-binding sites. These studies also implicate helices 1 and 12 of the BLA domain in allosteric signal propagation. These results support a model for the evolution of heterotropic allostery in which effector affinity and allosteric signaling emerge simultaneously.     

Molecular dynamics simulations of aptamer-binding reveal generalized allostery in thrombin

Thrombin is an attractive target for antithrombotic therapy due to its central role in thrombosis and hemostasis as well as its role in inducing tumor growth, metastasis, and tumor invasion. The thrombin-binding DNA aptamer (TBA), is under investigation for anticoagulant drugs. Although aptamer binding experiments have been revealed various effects on thrombin’s enzymatic activities, the detailed picture of the thrombin’s allostery from TBA binding is still unclear. To investigate thrombin’s response to the aptamer-binding at the molecular level, we compare the mechanical properties and free energy landscapes of the free and aptamer-bound thrombin using microsecond-scale all-atom GPU-based molecular dynamics simulations. Our calculations on residue fluctuations and coupling illustrate the allosteric effects of aptamer-binding at the atomic level, highlighting the exosite II, 60s, γ and the sodium loops, and the alpha helix region in the light chains involved in the allosteric changes. This level of details clarifies the mechanisms of previous experimentally demonstrated phenomena, and provides a prediction of the reduced autolysis rate after aptamer-binding. The shifts in thrombin’s ensemble of conformations and free energy surfaces after aptamer-binding demonstrate that the presence of bound-aptamer restricts the conformational freedom of thrombin suggesting that conformational selection, i.e. generalized allostery, is the dominant mechanism of thrombin-aptamer binding. The profound perturbation on thrombin’s mechanical and thermodynamic properties due to the aptamer-binding, which was revealed comprehensively as a generalized allostery in this work, may be exploited in further drug discovery and development.     

Conservation of globin genes in the “living fossil” Latimeria chalumnae and reconstruction of the evolution of the vertebrate globin family

The (hemo-)globins are among the best-investigated proteins in biomedical sciences. These small heme-proteins play an important role in oxygen supply, but may also have other functions. In addition to well known hemoglobin and myoglobin, six other vertebrate globin types have been identified in recent years: neuroglobin, cytoglobin, globin E, globin X, globin Y, and androglobin. Analyses of the genome of the “living fossil” Latimeria chalumnae show that the coelacanth is the only known vertebrate that includes all eight globin types. Thus, Latimeria can also be considered as a “globin fossil”. Analyses of gene synteny and phylogenetic reconstructions allow us to trace the evolution and the functional changes of the vertebrate globin family. Neuroglobin and globin X diverged from the other globin types before the separation of Protostomia and Deuterostomia. The cytoglobins, which are unlikely to be involved in O₂ supply, form the earliest globin branch within the jawed vertebrates (Gnathostomata), but do not group with the agnathan hemoglobins, as it has been proposed before. There is strong evidence from phylogenetic reconstructions and gene synteny that the eye-specific globin E and muscle-specific myoglobin constitute a common clade, suggesting a similar role in intracellular O₂ supply. Latimeria possesses two α- and two β-hemoglobin chains, of which one α-chain emerged prior to the divergence of Actinopterygii and Sarcopterygii, but has been retained only in the coelacanth. Notably, the embryonic hemoglobin α-chains of Gnathostomata derive from a common ancestor, while the embryonic β-chains – with the exception of a more complex pattern in the coelacanth and amphibians – display a clade-specific evolution. Globin Y is associated with the hemoglobin gene cluster, but its phylogenetic position is not resolved. Our data show an early divergence of distinct globin types in the vertebrate evolution before the emergence of tetrapods. The subsequent loss of globins in certain taxa may be associated with changes in the oxygen-dependent metabolism. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.     

Modulation of allostery of pyruvate kinase by shifting of an ensemble of microstates

Since the introduction of the concepts of allostery about four decades ago, much advancement has been made in elucidating the structure-function correlation in allostery. However, there are still a number of issues that remain unresolved. In this review we used mammalian pyruvate kinase (PK) as a model system to understand the role of protein dynamics in modulating cooperativity. PK has a triosephosphate isomerase (TIM)(α/β)₈ barrel structural motif. PK is an ideal system to address basic questions regarding regulatory mechanisms about this common (α/β)₈ structural motif. The simplest model accounting for all of the solution thermodynamic and kinetic data on ligand-enzyme interactions involves two conformational states, inactive E^T and active E^R. These conformational states are represented by domain movements. Further studies provide the first evidence for a differential effect of ligand binding on the dynamics of the structural elements, not major secondary structural changes. These data are consistent with our model that allosteric regulation of PK is the consequence of perturbation of the distribution of an ensemble of states in which the inactive E^T and active E^R represent the two extreme end states. Sequence differences and ligands can modulate the distribution of states leading to alterations of functions. The future work includes: defining the network of functionally connected residues; elucidating the chemical principles governing the sequence differences which affect functions; and probing the nature of mutations on the stability of the secondary structural elements, which in turn modulate allostery.     

Unique features of the hemoglobin system of the Antarctic notothenioid fish Gobionotothen gibberifrons.

The hemolysate of the Antarctic teleost Gobionotothen gibberifrons (family Nototheniidae) contains two hemoglobins (Hb 1 and Hb 2). The concentration of Hb 2 (15-20% of the total hemoglobin content) is higher than that found in most cold-adapted Notothenioidei. Unlike the other Antarctic species so far examined having two hemoglobins, Hb 1 and Hb 2 do not have globin chains in common. Therefore this hemoglobin system is made of four globins (two alpha- and two beta-chains). The complete amino-acid sequence of the two hemoglobins (Hb 1, alpha2(1)beta2(1); Hb 2, alpha2(2)beta2(2)) has been established. The two hemoglobins have different functional properties. Hb 2 has lower oxygen affinity than Hb 1, and higher sensitivity to the modulatory effect of organophosphates. They also differ thermodynamically, as shown by the effects on the oxygen-binding properties brought about by temperature variations. The oxygen-transport system of G. gibberifrons, with two functionally distinct hemoglobins, suggests that the two components may have distinct physiological roles, in relation with life style and the environmental conditions which the fish may have to face. The unique features of the oxygen-transport system of this species are reflected in the phylogeny of the hemoglobin amino-acid sequences, which are intermediate between those of other fish of the family Nototheniidae and of species of the more advanced family Bathydraconidae.     

Interfacial and distal-heme pocket mutations exhibit additive effects on the structure and function of hemoglobin

Protein engineering strategies seek to develop a hemoglobin-based oxygen carrier with optimized functional properties, including (i) an appropriate O 2 affinity, (ii) high cooperativity, (iii) limited NO reactivity, and (iv) a diminished rate of auto-oxidation. The mutations alphaL29F, alphaL29W, alphaV96W and betaN108K individually impart some of these traits and in combinations produce hemoglobin molecules with interesting ligand-binding and allosteric properties. Studies of the ligand-binding properties and solution structures of single and multiple mutants have been performed. The aromatic side chains placed in the distal-heme pocket environment affect the intrinsic ligand-binding properties of the mutated subunit itself, beyond what can be explained by allostery, and these changes are accompanied by local structural perturbations. In contrast, hemoglobins with mutations in the alpha 1beta 1 and alpha 1beta 2 interfaces display functional properties of both "R"- and "T"-state tetramers because the equilibrium between quaternary states is altered. These mutations are accompanied by global structural perturbations, suggesting an indirect, allostery-driven cause for their effects. Combinations of the distal-heme pocket and interfacial mutations exhibit additive effects in both structural and functional properties, contribute to our understanding of allostery, and advance protein-engineering methods for manipulating the O 2 binding properties of the hemoglobin molecule.     

Globin gene family evolution and functional diversification in annelids

Globins are the most common type of oxygen-binding protein in annelids. In this paper, we show that circulating intracellular globin (Alvinella pompejana and Glycera dibranchiata), noncirculating intracellular globin (Arenicola marina myoglobin) and extracellular globin from various annelids share a similar gene structure, with two conserved introns at canonical positions B12.2 and G7.0. Despite sequence divergence between intracellular and extracellular globins, these data strongly suggest that these three globin types are derived from a common ancestral globin-like gene and evolved by duplication events leading to diversification of globin types and derived functions. A phylogenetic analysis shows a distinct evolutionary history of annelid extracellular hemoglobins with respect to intracellular annelid hemoglobins and mollusc and arthropod extracellular hemoglobins. In addition, dehaloperoxidase (DHP) from the annelid, Amphitrite ornata, surprisingly exhibits close phylogenetic relationships to some annelid intracellular globins. We have characterized the gene structure of A. ornata DHP to confirm assumptions about its homology with globins. It appears that it has the same intron position as in globin genes, suggesting a common ancestry with globins. In A. ornata, DHP may be a derived globin with an unusual enzymatic function.     

Monosaccharides bound to hemoglobins in normal and diabetic individuals. Evidence for glucose, mannose and galactose as sugars released by methanolysis of the different hemoglobin components

Direct evidence is given for the presence of glucose, mannose and galactose as the products of hydrolysis of hemoglobins A1a1, A1a2, A1b, A1c and A0. The presence of galactose cannot be explained by the earlier hypothesis of Amadori rearrangement and suggests the existence of further complex rearrangements. Monosaccharide content of the different hemoglobin components varies from 0.2-2.0 mol/mol of alpha beta dimer with an increase of 1.5-2.0-times in diabetic components. This increase is not accompanied by net charge differences, suggesting that additionally bound sugars are not responsible for the pI modification of these hemoglobins. The pattern of glucose, mannose and galactose ratio in normal individuals divides these hemoglobins into two classes, hemoglobins A1b, A1c and A0 (ratio 0.60:0.25:0.15) on one hand and hemoglobins A1a1 and A1a2 (ratio 0.40:0.40: 0.20) on the other. These findings suggest that diverse mechanisms for sugar binding might exist between these two classes of glycosylated hemoglobins. This difference disappears in diabetic components suggesting that the non-NH2-terminal sites are glycosylated in all components by a common mechanism. Increase in glucose at the expense of mannose and galactose, as observed in diabetics, could be an indicator of recent glycosylation.