Deciphering the molecular code of hemoglobin allostery

Lanmodulin's EF 2-3 Domain: Insights from Infrared Spectroscopy and Simulations

Lanmodulins are small, ∼110-residue proteins with four EF-hand motifs that demonstrate a picomolar affinity for lanthanide ions, making them efficient in the recovery and separation of these technologically important metals. In this study, we examine the thermodynamic and structural complexities of lanthanide ion binding to a 41-residue domain, EF 2-3, that constitutes the two highest-affinity metal-binding sites in the lanmodulin protein from Methylorubrum extorquens. Using a combination of circular dichroism (CD) spectroscopy, isothermal titration calorimetry (ITC), two-dimensional infrared (2D IR) spectroscopy, and molecular dynamics (MD) simulations, we characterize the metal binding capabilities of EF 2-3. ITC demonstrates that binding occurs between peptide and lanthanides with conditional dissociation constants (Kd) in the range 20-30 μM, with no significant differences in the Kd values for La3+, Eu3+, and Tb3+ at pH 7.4. In addition, CD spectroscopy suggests that only one binding site of EF 2-3 undergoes a significant conformational change in the presence of lanthanides. 2D IR spectroscopy demonstrates the presence of both mono- and bidentate binding configurations in EF 2-3 with all three lanthanides. MD simulations, supported by Eu3+ luminescence measurements, explore these results, suggesting a competition between water-lanthanide and carboxylate-lanthanide interactions in the EF 2-3 domain. These results underscore the role of the core helical bundle of the protein architecture in influencing binding affinities and communication between the metal-binding sites in the full-length protein.

What is allosteric regulation? Exploring the exceptions that prove the rule!

"Allosteric" was first introduced to mean the other site (i.e., a site distinct from the active or orthosteric site), an adjective for "regulation" to imply a regulatory outcome resulting from ligand binding at another site. That original idea outlines a system with two ligand-binding events at two distinct locations on a macromolecule (originally a protein system), which defines a four-state energy cycle. An allosteric energy cycle provides a quantifiable allosteric coupling constant and focuses our attention on the unique properties of the four equilibrated protein complexes that constitute the energy cycle. Because many observed phenomena have been referenced as "allosteric regulation" in the literature, the goal of this work is to use literature examples to explore which systems are and are not consistent with the two-ligand thermodynamic energy cycle-based definition of allosteric regulation. We emphasize the need for consistent language so comparisons can be made among the ever-increasing number of allosteric systems. Building on the mutually exclusive natures of an energy cycle definition of allosteric regulation versus classic two-state models, we conclude our discussion by outlining how the often-proposed Rube-Goldberg-like mechanisms are likely inconsistent with an energy cycle definition of allosteric regulation.

Structural origin of cooperativity in human hemoglobin: a view from different roles of α and β subunits in the α2β2 tetramer

This mini-review, mainly based on our resonance Raman studies on the structural origin of cooperative O₂ binding in human adult hemoglobin (HbA), aims to answering why HbA is a tetramer consisting of two α and two β subunits. Here, we focus on the Fe-His bond, the sole coordination bond connecting heme to a globin. The Fe-His stretching frequencies reflect the O₂ affinity and also the magnitude of strain imposed through globin by inter-subunit interactions, which is the origin of cooperativity. Cooperativity was first explained by Monod, Wyman, and Changeux, referred to as the MWC theory, but later explained by the two tertiary states (TTS) theory. Here, we related the higher-order structures of globin observed mainly by vibrational spectroscopy to the MWC theory. It became clear from the recent spectroscopic studies, X-ray crystallographic analysis, and mutagenesis experiments that the Fe-His bonds exhibit different roles between the α and β subunits. The absence of the Fe-His bond in the α subunit in some mutant and artificial Hbs inhibits T to R quaternary structural change upon O₂ binding. However, its absence from the β subunit in mutant and artificial Hbs simply enhances the O₂ affinity of the α subunit. Accordingly, the inter-subunit interactions between α and β subunits are nonsymmetric but substantial for HbA to perform cooperative O₂ binding.

Impact of Conformational Substates and Energy Landscapes on Understanding Hemoglobin Kinetics and Function

Hans Frauenfelder's discovery of conformational substates in studies of myoglobin carbon monoxide geminate rebinding kinetics at cryogenic temperatures (Austin RH, Beeson KW, Eisenstein L, Frauenfelder H, & Gunsalus IC (1975) Dynamics of Ligand Binding to Myoglobin. Biochemistry 14(24):5355-5373) followed by his introduction of energy landscape theory with Peter Wolynes (Frauenfelder H, Sligar SG, & Wolynes PG (1991) The Energy Landscapes and Motions of Proteins. Science 254(5038):1598-1603) marked the beginning of a new era in the physics and physical chemistry of proteins. Their work played a major role in demonstrating the power and importance of dynamics and of Kramers reaction rate theory for understanding protein function. The biggest impact of energy landscape theory has been in the protein folding field, which is well-known and has been documented in numerous articles and reviews, including a recent one of my own (Eaton WA (2021) Modern Kinetics and Mechanism of Protein Folding: a Retrospective. J. Phys. Chem. B. 125(14):3452-3467). Here I will describe the much less well-known impact of their modern view of proteins on both experimental and theoretical studies of hemoglobin kinetics and function. I will first describe how Frauenfelder's experiments motivated and influenced my own research on myoglobin, which were key ingredients to my work on understanding hemoglobin.

Photooxidase Mimicking with Adaptive Coordination Molecular Capsules

One important feature of enzyme catalysis is the induced-fit conformational change after binding substrates. Herein, we report a biomimetic water-soluble molecular capsule featuring adaptive structural change toward substrate binding, which offers an ideal platform for efficient photocatalysis. The molecular capsule was coordination-assembled from three anthracene-bridged bis-TPT [TPT = 2,4,6-tris(4-pyridyl)-1,3,5-triazine] ligands and six (bpy)Pd(NO3)2 (bpy = 2,2'-bipyridine). Once substrates bind to its hydrophobic cavity, this capsule would undergo quantitative capsule-to-bowl transformation. Visible-light absorption brought about by both the anthracene units and the charge-transfer absorption on the late-formed quintuple π-π stacked host-guest complex efficiently facilitates aerobic photooxidation for the sulfide guests by visible-light irradiation under mild conditions. Desired turnover numbers and product selectivity (sulfoxide over sulfone) have been achieved by the transformable nature of the catalyst and the hydrophilicity of the sulfoxide product. Such a photocatalytic process enabled by an adaptive coordination capsule and substrates as the allosteric effector paves the way for constructing artificial systems to mimic enzyme catalysis.

Impact of hemoglobin biophysical studies on molecular pathogenesis and drug therapy for sickle cell disease

Basic research on hemoglobin has been essential for understanding the origin and treatment of many hematological disorders due to abnormal hemoglobins. The most important of the hemoglobinopathies is sickle cell disease - Linus Pauling's "molecular disease" that gave birth to molecular medicine. In this review, I will describe the contributions of basic biophysical research on normal and sickle cell hemoglobin (HbS) to understanding the molecular pathogenesis of the disease and providing the conceptual basis for the various approaches to drug therapy that target HbS polymerization. Most prominent among these are the experimental results on the solubility of HbS as a function of oxygen saturation explained by the allosteric model of Monod, Wyman, and Changeux and the Gill-Wyman thermodynamic linkage relation between solubility and oxygen binding, the solubility of mixtures of HbS with normal or fetal hemoglobin explained by Minton's thermodynamic model, and the highly unusual kinetics of HbS polymerization explained by a novel double nucleation mechanism that also accounts for the aggregation kinetics of the Alzheimer's peptide. The HbS polymerization kinetics are of great importance to understanding the pathophysiology and clinical course, as well as guiding drug development for treating this common and severe disease. The article focuses primarily on experimental and theoretical results from my lab, so it is not a comprehensive review of the subject.

Inter-Active Site Communication Mediated by the Dimer Interface β-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig₀), the saturated holoenzyme (lig₂), and the typically elusive singly bound (lig₁) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface β-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique β-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig₀ to lig₁ and ultimately lig₂. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Small molecule cores demonstrate non-competitive inhibition of lactate dehydrogenase

Lactate dehydrogenase (LDH) has recently garnered attention as an attractive target for cancer therapies, owing to the enzyme's critical role in cellular metabolism. Current inhibition strategies, employing substrate or cofactor analogues, are insufficiently specific for use as pharmaceutical agents. The possibility of allosteric inhibition of LDH was postulated on the basis of theoretical docking studies of a small molecule inhibitor to LDH. The present study examined structural analogues of this proposed inhibitor to gauge its potency and attempt to elucidate the molecular mechanism of action. These analogues display encouraging in vitro inhibition of porcine heart LDH, including micromolar K_i values and a maximum inhibition of up to 50% in the steady state. Furthermore, Michaelis-Menten kinetics and fluorescence data both suggest the simple, acetaminophen derivatives are non-competitive in binding to the enzyme. Kinetic comparisons of a panel of increasingly decorated structural analogues imply that the binding is specific, and the small molecule core provides a privileged scaffold for further pharmaceutical development of a novel, allosteric drug.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Experiments on Hemoglobin in Single Crystals and Silica Gels Distinguish among Allosteric Models

Trapping quaternary structures of hemoglobin in single crystals or by encapsulation in silica gels has provided a demanding set of data to test statistical mechanical models of allostery. In this work, we compare the results of those experiments with predictions of the four major allosteric models for hemoglobin: the quaternary two-state model of Monod, Wyman, and Changeux; the tertiary two-state model of Henry et al., which is the simplest extension of the Monod-Wyman-Changeux model to include pre-equilibria of tertiary as well as quaternary conformations; the structure-based model of Szabo and Karplus; and the modification of the latter model by Lee and Karplus. We show that only the tertiary two-state model can provide a near quantitative explanation of the single-crystal and gel experimental results.

Identification of an L-phenylalanine binding site enhancing the cooperative responses of the calcium-sensing receptor to calcium

Functional positive cooperative activation of the extracellular calcium ([Ca(2+)]o)-sensing receptor (CaSR), a member of the family C G protein-coupled receptors, by [Ca(2+)]o or amino acids elicits intracellular Ca(2+) ([Ca(2+)]i) oscillations. Here, we report the central role of predicted Ca(2+)-binding site 1 within the hinge region of the extracellular domain (ECD) of CaSR and its interaction with other Ca(2+)-binding sites within the ECD in tuning functional positive homotropic cooperativity caused by changes in [Ca(2+)]o. Next, we identify an adjacent L-Phe-binding pocket that is responsible for positive heterotropic cooperativity between [Ca(2+)]o and L-Phe in eliciting CaSR-mediated [Ca(2+)]i oscillations. The heterocommunication between Ca(2+) and an amino acid globally enhances functional positive homotropic cooperative activation of CaSR in response to [Ca(2+)]o signaling by positively impacting multiple [Ca(2+)]o-binding sites within the ECD. Elucidation of the underlying mechanism provides important insights into the longstanding question of how the receptor transduces signals initiated by [Ca(2+)]o and amino acids into intracellular signaling events.

Probability fluxes and transition paths in a Markovian model describing complex subunit cooperativity in HCN2 channels

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are voltage-gated tetrameric cation channels that generate electrical rhythmicity in neurons and cardiomyocytes. Activation can be enhanced by the binding of adenosine-3',5'-cyclic monophosphate (cAMP) to an intracellular cyclic nucleotide binding domain. Based on previously determined rate constants for a complex Markovian model describing the gating of homotetrameric HCN2 channels, we analyzed probability fluxes within this model, including unidirectional probability fluxes and the probability flux along transition paths. The time-dependent probability fluxes quantify the contributions of all 13 transitions of the model to channel activation. The binding of the first, third and fourth ligand evoked robust channel opening whereas the binding of the second ligand obstructed channel opening similar to the empty channel. Analysis of the net probability fluxes in terms of the transition path theory revealed pronounced hysteresis for channel activation and deactivation. These results provide quantitative insight into the complex interaction of the four structurally equal subunits, leading to non-equality in their function.

A novel type of allosteric regulation: functional cooperativity in monomeric proteins

Cooperative functional properties and allosteric regulation in cytochromes P450 play an important role in xenobiotic metabolism and define one of the main mechanisms of drug-drug interactions. Recent experimental results suggest that ability to bind simultaneously two or more small organic molecules can be the essential feature of cytochrome P450 fold, and often results in rich and complex pattern of allosteric behavior. Manifestations of non-Michaelis kinetics include homotropic and heterotropic activation and inhibition effects depending on the stoichiometric ratios of substrate and effector, changes in the regio- and stereospecificity of catalytic transformations, and often give rise to the clinically important drug-drug interactions. In addition, functional response of P450 systems is modulated by the presence of specific and non-specific effector molecules, metal ions, membrane incorporation, formation of homo- and hetero-oligomers, and interactions with the protein redox partners. In this article we briefly overview the main factors contributing to the allosteric effects in cytochromes P450 with the main focus on the sources of cooperative behavior in xenobiotic metabolizing monomeric heme enzymes with their conformational flexibility and extremely broad substrate specificity. The novel mechanism of functional cooperativity in P450 enzymes does not require substantial binding cooperativity, rather it implies the presence of one or more binding sites with higher affinity than the single catalytically active site in the vicinity of the heme iron.

Kinetic trapping of a key hemoglobin intermediate

The complete binding cascade of human hemoglobin consists of a series of partially ligated intermediates. The individual intermediate binding constants cannot be distinguished in O(2) binding curves, however, each constant can be determined from the O(2)-induced change in assembly constant for the α(2)β(2) tetramer from its constituent αβ dimers. The characterization of these O(2) binding constants has shown the Hb cascade to be asymmetric in nature, with binding dependent upon the specific distribution of O(2) among the four hemesites. A stopped-flow approach to measuring the dissociation constant of a key doubly ligated intermediate, that in which one dimer is oxygenated and the other is not, is described. The intermediate is transiently formed in the absence of O(2) and then allowed to dissociate in the presence of O(2). The free dimers thus released are trapped by the plasma protein haptoglobin, the rate limiting step being that of tetramer dissociation. The kinetic constant observed for the dissociation of this intermediate confirms the value for its equilibrium O(2) binding constant, previously determined under equilibrium conditions by subzero isoelectric focusing.

TOPOLOGICAL STUDY OF REACTIVE POTENTIAL-ENERGY SURFACES FROM MODEL DIABATIC STATES COUPLED TO SURROUNDING MEDIA

Reaction paths are studied with the help of diabatic potential-energy surfaces coupled in a generic external field. We show that all putative geometrical and topological features of two-dimensional (2D) potential-energy surfaces for an isomerization can be generated with a model consisting of strictly diabatic harmonic 2D wells coupled in a static (yet uniform) external field. For a large class of three-state models, we provide a phase diagram of possible topologies in the relevant parameter space for three-state models. By following the shift in diabatic electronic amplitudes (denoted by {|ck|2}) in coherent quantum states produced in the field, we assess whether the models within each phase can produce reaction paths visiting regions of configurational space that are structurally similar to the diabatic transition state.

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.

An investigation of the distal histidyl hydrogen bonds in oxyhemoglobin: effects of temperature, pH, and inositol hexaphosphate

On the basis of X-ray crystal structures and electron paramagnetic resonance (EPR) measurements, it has been inferred that the O(2) binding to hemoglobin is stabilized by the hydrogen bonds between the oxygen ligands and the distal histidines. Our previous study by multinuclear nuclear magnetic resonance (NMR) spectroscopy has provided the first direct evidence of such H-bonds in human normal adult oxyhemoglobin (HbO(2) A) in solution. Here, the NMR spectra of uniformly (15)N-labeled recombinant human Hb A (rHb A) and five mutant rHbs in the oxy form have been studied under various experimental conditions of pH and temperature and also in the presence of an organic phosphate, inositol hexaphosphate (IHP). We have found significant effects of pH and temperature on the strength of the H-bond markers, i.e., the cross-peaks for the side chains of the two distal histidyl residues, α58His and β63His, which form H-bonds with the O(2) ligands. At lower pH and/or higher temperature, the side chains of the distal histidines appear to be more mobile, and the exchange with water molecules in the distal heme pockets is faster. These changes in the stability of the H-bonds with pH and temperature are consistent with the changes in the O(2) affinity of Hb as a function of pH and temperature and are clearly illustrated by our NMR experiments. Our NMR results have also confirmed that this H-bond in the β-chain is weaker than that in the α-chain and is more sensitive to changes in pH and temperature. IHP has only a minor effect on these H-bond markers compared to the effects of pH and temperature. These H-bonds are sensitive to mutations in the distal heme pockets but not affected directly by the mutations in the quaternary interfaces, i.e., α(1)β(1) and/or α(1)β(2) subunit interface. These findings provide new insights regarding the roles of temperature, hydrogen ion, and organic phosphate in modulating the structure and function of hemoglobin in solution.

Cooperative properties of cytochromes P450

Cytochromes P450 form a large and important class of heme monooxygenases with a broad spectrum of substrates and corresponding functions, from steroid hormone biosynthesis to the metabolism of xenobiotics. Despite decades of study, the molecular mechanisms responsible for the complex non-Michaelis behavior observed with many members of this superfamily during metabolism, often termed 'cooperativity', remain to be fully elucidated. Although there is evidence that oligomerization may play an important role in defining the observed cooperativity, some monomeric cytochromes P450, particularly those involved in xenobiotic metabolism, also display this behavior due to their ability to simultaneously bind several substrate molecules. As a result, formation of distinct enzyme-substrate complexes with different stoichiometry and functional properties can give rise to homotropic and heterotropic cooperative behavior. This review aims to summarize the current understanding of cooperativity in cytochromes P450, with a focus on the nature of cooperative effects in monomeric enzymes.

The Hill coefficient: inadequate resolution of cooperativity in human hemoglobin

The Hill coefficient nH is a dimensionless parameter that has long been used as a measure of the extent of cooperativity. Originally derived from the oxygen-binding curve of human hemoglobin (Hb) by A. V. Hill in 1910, and reinvented by J. Wyman several decades later, nH is indexed to the stoichiometry of ligation and is indirectly related to the overall cooperative free energy for binding all four oxygen ligands. However, the overall cooperative free energy of Hb ligation can be measured directly by experimental methods. The microscopic cooperative free energies that relate to energetic coupling between specific subunit pairs can also be experimentally determined, while the Hill coefficient is, by its nature, a macroscopic parameter that cannot detect differences among specific subunit-subunit couplings. Its continued use in studies of the mechanism of cooperativity in Hb is therefore of increasingly limited value.

Allostery: an illustrated definition for the 'second secret of life'

Although allosteric regulation is the 'second secret of life', the molecular mechanisms that give rise to allostery currently elude understanding. In my opinion, experimental progress is hampered by a commonly used but misleading definition of allostery as protein structural changes that are elicited by the binding of a single ligand. Allostery is more strictly defined in functional terms as a comparison of how one ligand binds in the absence, versus the presence, of a second ligand. Therefore, as each of the two binding events involves two protein complexes, a study of allostery must consider four complexes and not just two. Such a comparison can distinguish allosteric from non-allosteric protein changes, the importance of which is frequently overlooked. When a study of all four complexes is not feasible, an alternative, albeit limited, strategy can identify subsets of allosteric-specific changes.

Overview of the quantitation of protein interactions

The biological function of many proteins involves reversible interactions with other proteins, nucleic acids, or other non-protein ligands. Such interactions play many different roles in a wide range of cellular processes. A few examples are: (1) storing or transporting key metabolites (e.g., O(2) storage by myoglobin); (2) forming and maintaining the quaternary structure of multi-subunit enzymes; (3) specific binding and recognition events (antigen-antibody, hormone-receptor, transcription factor-promoter); and (4) self-assembly of large structures (microtubules, chromatin). Thus, the quantitative characterization of such interactions represents an important part of understanding the function of such proteins and their role in these cellular events. This unit sets the tone for the rest of the chapter, and gives important information necessary to understand some of the topics that will be covered in future supplements, such as sedimentation equilibrium (analytical and micro-preparative), surface plasmon resonance (SPR), size-exclusion chromatography (SEC) with on-line light scattering, and chemical cross-linking.

Protein dynamics explain the allosteric behaviors of hemoglobin

Bohr, Hasselbalch, and Krogh discovered homotropic and heterotropic allosteric behaviors of hemoglobin (Hb) in 1903/1904. A chronological description since then of selected principal models of the allosteric mechanism of Hb, such as the Adair scheme, the MWC two-state concerted model, the KNF induced-fit sequential model, the Perutz stereochemical model, the tertiary two-state model, and the global allostery model (an expanded MWC models), is concisely presented, followed by analysis and discussion of their limitations and deficiencies. The determination of X-ray crystallographic structures of deoxy- and ligated-Hb and the structure-based stereochemical model by Perutz are an epoch-making event in this history. However, his assignment of low-affinity deoxy- and high-affinity oxy-quaternary structures of Hb to the T- and R-states, respectively, though apparently reasonable, and as well as his hypothesis that the T-/R-quaternary structural transition regulates the oxygen-affinity, have created confusions and side-tracked studies of Hb on the structure-function relationship. The differences in static molecular structures of Hb between T(deoxy)- and R(oxy)-quaternary states reported in detail by Perutz and others are ligation-linked structural changes, but not related to the control/regulation of the oxygen-affinity. The oxygen-affinity (K(T) and K(R)) of Hb has been shown to be regulated by the heterotropic effector-linked tertiary structural changes without involving the T/R-quaternary changes. However, a recent high-resolution crystallographic analysis of Hb with different oxygen-affinities shows that static molecular structures of Hb determined by crystallography can neither identify the nature of the T(low-affinity) functional state nor decipher the mechanism by which Hb stores free energy in the T(low-affinity) functional state. Molecular dynamics simulations show that fluctuations of helices of oxy-Hb are increased upon de-oxygenation and/or binding 2,3-biphosphoglycerate. These are known to lower the oxygen-affinity of Hb. It is proposed that the coordination mode of the heme Fe with proximal and distal His is modulated by these helical fluctuations, resulting in the modulation of the oxygen-affinity of Hb. Therefore, it is proposed that the oxygen-affinity of Hb is regulated by pentanary (the 5th-order time-dependent or dynamic) tertiary structural changes rather than the T-/R-quaternary structural transitions in Hb. Homotropic and heterotropic allosteric effects of Hb are oxygen- and effector-linked, conformational entropy-driven entropy-enthalpy compensation phenomena and not much to do with static structural changes. The dynamic allostery model, which integrates these observations, provides the structural basis for the global allostery model (an expanded MWC model).

Hb(alphaalpha,betabeta): a novel fusion construct for a dimeric, four-domain hemoglobin

Hemoglobin-based blood substitutes are one of the options available to derive a resuscitating fluid taking into account clinical and physiological demands. In this paper we investigated a novel protein, Hb(alphaalpha,betabeta) obtained as a combination of two homodimers alpha(2) and beta(2) both derived from a fusion gene containing two alfa chains or two beta chains, each respectively coupled via a specific linker. The construct here described is thus a novel heterodimeric hemoglobin carrying four heme groups. The protein cannot dissociate into dimers, as demonstrated by its absence of reactivity versus haptoglobin, and is expected to have a relatively long circulating half-life. The modification does not increase the autoxidation rate, but increases the oxygen affinity, due to a destabilization of the T quaternary state. Characterization of the biochemical properties of this protein in comparison with HbA is reported.

Ligand effects on the protein ensemble: unifying the descriptions of ligand binding, local conformational fluctuations, and protein stability

Detailed description of the structural and physical basis of allostery, cooperativity, and other manifestations of long-range communication between binding sites in proteins remains elusive. Here we describe an ensemble-based structural-thermodynamic model capable of treating explicitly the coupling between ligand binding reactions, local fluctuations in structure, and global conformational transitions. The H(+) binding reactions of staphylococcal nuclease and the effects of pH on its stability were used to illustrate the properties of proteins that can be described quantitatively with this model. Each microstate in the native ensemble was modeled to have dual structural character; some regions were treated as folded and retained the same atomic geometry as in the crystallographic structure while other regions were treated thermodynamically as if they were unfolded. Two sets of pK(a) values were used to describe the affinity of each H(+) binding site. One set, calculated with a standard continuum electrostatics method, describes H(+) binding to sites in folded parts of the protein. A second set of pK(a) values, obtained from model compounds in water, was used to describe H(+) binding to sites in unfolded regions. An empirical free energy function, parameterized to reproduce folding thermodynamics measured by differential scanning calorimetry, was used to calculate the probability of each microstate. The effects of pH on the distribution of microstates were determined by the H(+) binding properties of each microstate. The validity of the calculations was established by comparison with a number of different experimental observables.

Coarse-grained free energy functions for studying protein conformational changes: a double-well network model

In this work, a double-well network model (DWNM) is presented for generating a coarse-grained free energy function that can be used to study the transition between reference conformational states of a protein molecule. Compared to earlier work that uses a single, multidimensional double-well potential to connect two conformational states, the DWNM uses a set of interconnected double-well potentials for this purpose. The DWNM free energy function has multiple intermediate states and saddle points, and is hence a "rough" free energy landscape. In this implementation of the DWNM, the free energy function is reduced to an elastic-network model representation near the two reference states. The effects of free energy function roughness on the reaction pathways of protein conformational change is demonstrated by applying the DWNM to the conformational changes of two protein systems: the coil-to-helix transition of the DB-loop in G-actin and the open-to-closed transition of adenylate kinase. In both systems, the rough free energy function of the DWNM leads to the identification of distinct minimum free energy paths connecting two conformational states. These results indicate that while the elastic-network model captures the low-frequency vibrational motions of a protein, the roughness in the free energy function introduced by the DWNM can be used to characterize the transition mechanism between protein conformations.

Apparent specific volume of human hemoglobin: Effect of ligand state and contribution of heme

The apparent specific volumes of human deoxy-, oxy-, met-, and CN-met hemoglobin (Hb) were measured with a vibrating tube densitometer. The values were calculated from the difference in density between protein solutions and solvents with which they were in dialysis equilibrium. The results obtained were very similar to the value for horse HbCO often used for sedimentation studies of Hbs. The apparent specific volumes of oxy- and CN-metHb are approximately 0.0020 cm(3)/g higher than those of deoxy- and metHb. This small reproducible difference could be due either to a small conformational difference between the liganded and unliganded species or to different interactions with components of the solvent. On the basis of these results, a simple method for the determination of the contribution of the heme to the apparent specific volume is proposed. The contribution can be estimated from the difference between the measured volume of each molecular species and that calculated from the amino acid composition.

Structural Stability of Covalently Linked GroES Heptamer: Advantages in the Formation of Oligomeric Structure

In order to understand how inter-subunit association stabilizes oligomeric proteins, a single polypeptide chain variant of heptameric co-chaperonin GroES (tandem GroES) was constructed from Escherichia coli heptameric GroES by linking consecutively the C-terminal of one subunit to the N-terminal of the adjacent subunit with a small linker peptide. The tandem GroES (ESC7) showed properties similar to wild-type GroES in structural aspects and co-chaperonin activity. In unfolding and refolding equilibrium experiments using guanidine hydrochloride (Gdn-HCl) as a denaturant at a low protein concentration (50 microg ml(-1)), ESC7 showed a two-state transition with a greater resistance toward Gdn-HCl denaturation (Cm=1.95 M) compared to wild-type GroES (Cm=1.1 M). ESC7 was found to be about 10 kcal mol(-1) more stable than the wild-type GroES heptamer at 50 microg ml(-1). Kinetic unfolding and refolding experiments of ESC7 revealed that the increased stability was mainly attributed to a slower unfolding rate. Also a transient intermediate was detected in the refolding reaction. Interestingly, at the physiological GroES concentration (>1 mg ml(-1)), the free energy of unfolding for GroES heptamer exceeded that for ESC7. These results showed that at low protein concentrations (<1 mg ml(-1)), the covalent linking of subunits contributes to the stability but also complicates the refolding kinetics. At physiological concentrations of GroES, however, the oligomeric state is energetically preferred and the advantages of covalent linkage are lost. This finding highlights a possible advantage in transitioning from multi-domain proteins to oligomeric proteins with small subunits in order to improve structural and kinetic stabilities.

Electronic reorganization: Origin of sigma trans promotion effect

Binding of two ligands trans to each other by some transition metal complexes may be cooperative [Khoroshun et al., Mol Phys 2002, 100, 523]. Several interesting consequent effects include (i) inverse relationship between bond strength and binding affinity; (ii) smaller coordination barriers to formation of weaker bonds; (iii) enhancement of Lewis acidity with increased number of ligands. We describe a simple model, sigma trans promotion effect (TPE), which considers electronic reorganization between two Lewis structures, and predicts the above-mentioned effects. The applied result of present study is the unified perspective on several facts of heme chemistry. Particularly, we reiterate an important but often overlooked notion, developed previously within the spin pairing model [Drago and Corden, Acc Chem Res 1980, 13, 353], that, in hemoproteins, the proximal histidine and the distal ligand such as O2 or CO cooperate in promoting electronic reorganization. As a result, depopulation of dz2 orbital upon ligand binding contributes to the phenomenon of hemoglobin cooperativity. The presented density functional (B3LYP) calculations on realistic models, the processes of carbon monoxide binding by Fe(II) porphyrins and dinitrogen binding by triamido/triamidoamine Mo(III) complexes, particularly the evaluation of the coordination barriers due to spin-state change by location of the minima on seams of crossing, support the TPE model predictions. From a broader theoretical perspective, the present study would hopefully stimulate the development of much needed frameworks and tools for facile comparisons of wave functions and their properties between different geometries, species, and electronic states. Advancement of practical wave function comparisons may yield fresh qualitative perspectives on chemical reactivity, and promote better understanding of related concepts such as electronic reorganization.

Acid destabilization of the solution conformation of Bcl-xL does not drive its pH-dependent insertion into membranes

Regulation of programmed cell death by Bcl-xL is dependent on both its solution and integral membrane conformations. A conformational change from solution to membrane is also important in this regulation. This conformational change shows a pH-dependence similar to the translocation domain of diphtheria toxin, where an acid-induced molten globule conformation in the absence of lipid vesicles mediates the change from solution to membrane conformations. By contrast, Bcl-xL deltaTM in the absence of lipid vesicles exhibits no gross conformational changes upon acidification as observed by near- and far-UV circular dichroism spectropolarimetry. Additionally, no significant local conformational changes upon acidification were observed by heteronuclear NMR spectroscopy of Bcl-xL deltaTM. Under conditions that favor the solution conformation (pH 7.4), the free energy of folding for Bcl-xL deltaTM (deltaG(o)) was determined to be 15.8 kcal x mol(-1). Surprisingly, under conditions that favor a membrane conformation (pH 4.9), deltaG(o) was 14.6 kcal x mol(-1). These results differ from those obtained with many other membrane-insertable proteins where acid-induced destabilization is important. Therefore, other contributions must be necessary to destabilize the solution conformation Bcl-xL and favor the membrane conformation at pH 4.9. Such contributions might include the presence of a negatively charged membrane or an electrostatic potential across the membrane. Thus, for proteins that adopt both solution and membrane conformations, an obligatory molten globule intermediate may not be necessary. The absence of a molten globule intermediate might have evolved to protect Bcl-xL from intracellular proteases as it undergoes this conformational change essential for its activity.

NON-MICHAELIS-MENTEN KINETICS IN CYTOCHROME P450-CATALYZED REACTIONS

The cytochrome P450 monooxygenases (CYPs) are the dominant enzyme system responsible for xenobiotic detoxification and drug metabolism. Several CYP isoforms exhibit non-Michaelis-Menten, or “atypical,” steady state kinetic patterns. The allosteric kinetics confound prediction of drug metabolism and drug-drug interactions, and they challenge the theoretical paradigms of allosterism. Both homotropic and heterotropic ligand effects are now widely documented. It is becoming apparent that multiple ligands can simultaneously bind within the active sites of individual CYPs, and the kinetic parameters change with ligand occupancy. In fact, the functional effect of any specific ligand as an activator or inhibitor can be substrate dependent. Divergent approaches, including kinetic modeling and X-ray crystallography, are providing new information about how multiple ligand binding yields complex CYP kinetics.

Local conformational fluctuations can modulate the coupling between proton binding and global structural transitions in proteins

Local conformational fluctuations in proteins can affect the coupling between ligand binding and global structural transitions. This finding was established by monitoring quantitatively how the population distribution in the ensemble of microstates of staphylococcal nuclease was affected by proton binding. Analysis of acid unfolding and proton-binding data with an ensemble-based model suggests that local fluctuations: (i) can be effective modulators of ligand-binding affinities, (ii) are important determinants of the cooperativity of ligand-driven global structural transitions, and (iii) are well represented thermodynamically as local unfolding processes. These studies illustrate how an ensemble-based description of proteins can be used to describe quantitatively the interdependence of local conformational fluctuations, ligand-binding processes, and global structural transitions. This level of understanding of the relationship between conformation, energy, and dynamics is required for a detailed mechanistic understanding of allostery, cooperativity, and other complex functional and regulatory properties of macromolecules.

Ligand-linked stability of mutants of the C-domain of calmodulin

There is a necessary energetic linkage between ligand binding and stability in biological molecules. The critical glutamate in Site 4 was mutated to create two mutants of the C-domain of calmodulin yielding E140D and E140Q. These proteins were stably folded in the absence of calcium, but had dramatically impaired binding of calcium. We determined the stability of the mutant proteins in the absence and presence of calcium using urea-induced unfolding monitored by circular dichroism (CD) spectroscopy. These calcium-dependent unfolding curves were fit to models that allowed for linkage of stability to binding of a single calcium ion to the native and unfolded states. Simultaneous analysis of the unfolding profiles for each mutant yielded estimates for calcium-binding constants that were consistent with results from direct titrations monitored by fluorescence. Binding to the unfolded state was not an important energetic contributor to the ligand-linked stability of these mutants.

New insights into allosteric mechanisms from trapping unstable protein conformations in silica gels

To understand why the classical two-state allosteric model of Monod, Wyman, and Changeux explains cooperative oxygen binding by hemoglobin but does not explain changes in oxygen affinity by allosteric inhibitors, we have investigated the kinetic properties of unstable conformations transiently trapped by encapsulation in silica gels. Conformational trapping reveals that after nanosecond photodissociation of carbon monoxide a large fraction of the subunits of the T quaternary structure has kinetic properties almost identical to those of subunits of the R quaternary structure. Addition of allosteric inhibitors reduces both the fraction of R-like subunits and the oxygen affinity of the T quaternary structure. These kinetic and equilibrium results are readily explained by a recently proposed generalization of the Monod-Wyman-Changeux model in which a pre-equilibrium between two functionally different tertiary, rather than quaternary, conformations plays the central role.

The hemoglobin cyanomet ligation analogue and carbon monoxide induce similar allosteric mechanisms

Current thermodynamic models of protein cooperativity predicting sigmoidal ligand equilibrium curves differ in the assumptions regarding the structural/functional properties of the intermediate ligation states. Quantitative information on the intermediates cannot be extracted from the equilibrium curves, but must be obtained from direct studies of the intermediates. Since the intermediates are intrinsically unstable species, ligation analogues with reduced mobility are indispensable tools for cooperativity studies provided that the tertiary/quaternary changes triggered by the ligation analogue are similar to those observed using the physiological ligands. We demonstrate that the valency exchange reactions occurring in mixtures of deoxy and cyanomethemoglobin yield non-random distributions of deoxy/cyanomet intermediates that resemble those observed in the equilibrium with carbon monoxide. Previous and new data using the analogue, in agreement with the studies of the CO intermediates, indicate that the mechanism of hemoglobin cooperativity is neither purely concerted nor sequential nor combinatorial, but contains some elements of each of these models.

Allosteric proteins: lessons to be learned from the hemoglobin intermediates

Allosteric proteins, such as hemoglobin, are assemblies of functional units, which undergo quaternary structural transitions in response to concentration changes of a specific ligand. Functional properties of hemoglobin ligation intermediates indicate that the tertiary structural changes induced by the ligand do not promote an equilibrium of quaternary structures.