Asymmetric Cooperativity in a Symmetric Tetramer: Human Hemoglobin

Probing High-order Protein Complexes Using Native Mass Spectrometry and Hydrogen/Deuterium Exchange Mass Spectrometry: A Case Study Using Fresh and Commercial Hemoglobin Samples

Native mass spectrometry (MS) analysis of protein complexes is highly susceptible to matrix effect, and addressing this predicament using buffer exchange is a common approach. Nevertheless, optimization of the buffer exchange protocol is not trivial. With the use of hemoglobin (Hb) as the model entity, it was discovered that the native mass spectrum of protein assembly is highly dependent on the buffer-exchange protocol. Given the dependence of native MS on the purification protocol, this work attempts to use hydrogen/deuterium exchange mass spectrometry (HDX-MS) for comparative studies of hemoglobin complexes in untreated fresh and commercial samples. The information obtained from the HDX study was found to correlate well with the native mass spectrometry analysis of the properly buffer-exchanged Hb samples. Both native MS and HDX-MS showed that the fresh Hb sample has retained the expected tetrameric structure, whereas the commercial Hb has largely been denatured to the dimeric form. These findings prove the complementarity of native MS and HDX-MS in the analysis of high-order protein complexes and stress the necessity to validate the integrity of the high-order structures of the proteins prior to the use of the protein samples for other biomedical studies.

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.

Slow dynamics measured by phosphorescence lifetime reveals global conformational changes in human adult hemoglobin induced by allosteric effectors

The mechanism underlying allostery in hemoglobin (Hb) is still not completely understood. Various models describing the action of allosteric effectors on Hb function have been published in the literature. It has also been reported that some allosteric effectors-such as chloride ions, inositol hexaphosphate, 2,3-diphospho-glycerate and bezafibrate-considerably lower the oxygen affinity of Hb. In this context, an important question is the extent to which these changes influence the conformational dynamics of the protein. Earlier, we elaborated a challenging method based on phosphorescence quenching, which makes characterizing protein-internal dynamics possible in the ms time range. The experimental technique involves phosphorescence lifetime measurements in thermal equilibrium at varied temperatures from 10 K up to 273 K, based on the signal of Zn-protoporphyrin substituted for the heme in the β-subunits of Hb. The thermal activation of protein dynamics was observed by the enhancement of phosphorescence quenching attributed to O2 diffusion. It was shown that the thermal activation of protein matrix dynamics was clearly distinguishable from the dynamic activation of the aqueous solvent, and was therefore highly specific for the protein. In the present work, the same method was used to study the changes in the parameters of the dynamic activation of human HbA induced by binding allosteric effectors. We interpreted the phenomenon as phase transition between two states. The fitting of this model to lifetime data yielded the change of energy and entropy in the activation process and the quenching rate in the dynamically activated state. The fitted parameters were particularly sensitive to the presence of allosteric effectors and could be interpreted in line with results from earlier experimental studies. The results suggest that allosteric effectors are tightly coupled to the dynamics of the whole protein, and thus underline the importance of global dynamics in the regulation of Hb function.

Structure and function of H+/K+ pump mutants reveal Na+/K+ pump mechanisms

Ion-transport mechanisms evolve by changing ion-selectivity, such as switching from Na⁺ to H⁺ selectivity in secondary-active transporters or P-type-ATPases. Here we study primary-active transport via P-type ATPases using functional and structural analyses to demonstrate that four simultaneous residue substitutions transform the non-gastric H⁺/K⁺ pump, a strict H⁺-dependent electroneutral P-type ATPase, into a bona fide Na⁺-dependent electrogenic Na⁺/K⁺ pump. Conversion of a H⁺-dependent primary-active transporter into a Na⁺-dependent one provides a prototype for similar studies of ion-transport proteins. Moreover, we solve the structures of the wild-type non-gastric H⁺/K⁺ pump, a suitable drug target to treat cystic fibrosis, and of its Na⁺/K⁺ pump-mimicking mutant in two major conformations, providing insight on how Na⁺ binding drives a concerted mechanism leading to Na⁺/K⁺ pump phosphorylation.

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.

Gastric proton pump with two occluded K+ engineered with sodium pump-mimetic mutations

The gastric H⁺,K⁺-ATPase mediates electroneutral exchange of 1H⁺/1K⁺ per ATP hydrolysed across the membrane. Previous structural analysis of the K⁺-occluded E2-P transition state of H⁺,K⁺-ATPase showed a single bound K⁺ at cation-binding site II, in marked contrast to the two K⁺ ions occluded at sites I and II of the closely-related Na⁺,K⁺-ATPase which mediates electrogenic 3Na⁺/2K⁺ translocation across the membrane. The molecular basis of the different K⁺ stoichiometry between these K⁺-counter-transporting pumps is elusive. We show a series of crystal structures and a cryo-EM structure of H⁺,K⁺-ATPase mutants with changes in the vicinity of site I, based on the structure of the sodium pump. Our step-wise and tailored construction of the mutants finally gave a two-K⁺ bound H⁺,K⁺-ATPase, achieved by five mutations, including amino acids directly coordinating K⁺ (Lys791Ser, Glu820Asp), indirectly contributing to cation-binding site formation (Tyr340Asn, Glu936Val), and allosterically stabilizing K⁺-occluded conformation (Tyr799Trp). This quintuple mutant in the K⁺-occluded E2-P state unambiguously shows two separate densities at the cation-binding site in its 2.6 Å resolution cryo-EM structure. These results offer new insights into how two closely-related cation pumps specify the number of K⁺ accommodated at their cation-binding site.

The gastric H⁺,K⁺-ATPase is a proton pump that creates the acidic environment of the stomach lumen, maintaining high proton gradient across the gastric mucosa cell membrane. Here, structural analysis of rationally designed H⁺,K⁺-ATPase mutants provides insight into this and other P-type ATPases cation binding stoichiometry and mechanisms.

Interrogating the Quaternary Structure of Noncanonical Hemoglobin Complexes by Electrospray Mass Spectrometry and Collision-Induced Dissociation

Various activation methods are available for the fragmentation of gaseous protein complexes produced by electrospray ionization (ESI). Such experiments can potentially yield insights into quaternary structure. Collision-induced dissociation (CID) is the most widely used fragmentation technique. Unfortunately, CID of protein complexes is dominated by the ejection of highly charged monomers, a process that does not yield any structural insights. Using hemoglobin (Hb) as a model system, this work examines under what conditions CID generates structurally informative subcomplexes. Native ESI mainly produced tetrameric Hb ions. In addition, "noncanonical" hexameric and octameric complexes were observed. CID of all these species [(αβ)₂, (αβ)₃, and (αβ)₄] predominantly generated highly charged monomers. In addition, we observed hexamer → tetramer + dimer dissociation, implying that hexamers have a tetramer··dimer architecture. Similarly, the observation of octamer → two tetramer dissociation revealed that octamers have a tetramer··tetramer composition. Gas-phase candidate structures of Hb assemblies were produced by molecular dynamics (MD) simulations. Ion mobility spectrometry was used to identify the most likely candidates. Our data reveal that the capability of CID to produce structurally informative subcomplexes depends on the fate of protein-protein interfaces after transfer into the gas phase. Collapse of low affinity interfaces conjoins the corresponding subunits and favors CID via monomer ejection. Structurally informative subcomplexes are formed only if low affinity interfaces do not undergo a major collapse. However, even in these favorable cases CID is still dominated by monomer ejection, requiring careful analysis of the experimental data for the identification of structurally informative subcomplexes.

Structural Energy Landscapes and Plasticity of the Microstates of Apo Escherichia coli cAMP Receptor Protein

The theory for allostery has evolved to a modern energy landscape ensemble theory, the major feature of which is the existence of multiple microstates in equilibrium. The properties of microstates are not well defined due to their transient nature. Characterization of apo protein microstates is important because the specific complex of the ligand-bound microstate defines the biological function. The information needed to link biological function and structure is a quantitative correlation of the energy landscapes between the apo and holo protein states. We employed the Escherichia coli cAMP receptor protein (CRP) system to test the features embedded in the ensemble theory because multiple crystalline apo and holo structures are available. Small angle X-ray scattering data eliminated one of the three apo states but not the other two. We defined the underlying energy landscape differences among the apo microstates by employing the computation algorithm COREX/BEST. The same connectivity patterns among residues in apo CRP are retained upon binding of cAMP. The microstates of apo CRP differ from one another by minor structural perturbations, resulting in changes in the energy landscapes of the various domains of CRP. Using the differences in energy landscapes among these apo states, we computed the cAMP binding energetics that were compared with solution biophysical results. Only one of the three apo microstates yielded data consistent with the solution data. The relative magnitude of changes in energy landscapes embedded in various apo microstates apparently defines the ultimate outcome of the cooperativity of binding.

Bohr equation and the lost allosteric Bohr effects in symmetry

Bohr, Hasselbalch and Krogh demonstrated a group of sigmoid curves under various carbon dioxide contents in 1904. Hill fitted these curves in 1910 with Hill equation without the physical meanings of Hill coefficient and dissociation constant. In 1965, Monod-Wyman-Changeux model (MWC) popularized the word “allostery” with 81 words of symmetry to define an orthosteric nature of cooperativity in a single and symmetric sigmoid curve. Paradoxically the MWC model didn’t quantify the homotropic Hill coefficient and confusingly described the symmetry of sigmoid shapes with three allosteric variables. A heterotropic Bohr equation, by clarifying the biophysical symmetry in allostery, suggests the solution of allosteric coefficients with only one Bohr variable. We reveal that the mathematical need of a fictional monomer by MWC model justify a symmetric logistic curve with a parabolic kernel of dissociation constant to model the 1904 sigmoid curves. The logistic-derived Bohr equation and its half-saturated P₅₀ equation successfully used the embedded P₅₀ values in the 1904 sigmoidal curves to quantify their hyperbolic conformational shifts and Hill coefficients (n) pending for a century. Both are the logarithmic functions of carbon dioxide. This truly quantitative Bohr equation digitizes the allosteric regulation of the orthosteric affinity by precisely cloning the original group of dissociation/association curves published in 1904. The Bohr equation honestly suggests that nature should have chosen the allosteric Bohr effects to modify hemoglobin to cope with the swift dynamic of gas exchange. The discovery of the Bohr function in Bohr equation challenges the feasibility of the orthosteric cooperativity of hemoglobin.

Combinatorial Genetics Reveals a Scaling Law for the Effects of Mutations on Splicing

Despite a wealth of molecular knowledge, quantitative laws for accurate prediction of biological phenomena remain rare. Alternative pre-mRNA splicing is an important regulated step in gene expression frequently perturbed in human disease. To understand the combined effects of mutations during evolution, we quantified the effects of all possible combinations of exonic mutations accumulated during the emergence of an alternatively spliced human exon. This revealed that mutation effects scale non-monotonically with the inclusion level of an exon, with each mutation having maximum effect at a predictable intermediate inclusion level. This scaling is observed genome-wide for cis and trans perturbations of splicing, including for natural and disease-associated variants. Mathematical modeling suggests that competition between alternative splice sites is sufficient to cause this non-linearity in the genotype-phenotype map. Combining the global scaling law with specific pairwise interactions between neighboring mutations allows accurate prediction of the effects of complex genotype changes involving >10 mutations.

Derivation and numerical profile analysis of a hierarchically formulated microscopic model of hemoglobin oxygen binding

To address complex thermodynamic systems with multiple interacting events, we have developed the concept of hierarchical thermodynamic interactions. In this study, this concept is extended to protein-ligand systems with similar but not identical protein subunits, and applied to the analysis of previously published NMR and UV-vis monitored hemoglobin oxygen binding data. Non-linear regression provided estimated errors for statistically significant parameters, but not for null (zero) valued parameters. A numerical/graphical profiling approach was therefore used to assess confidence intervals and correlations for both the statistically significant and nulled valued parameters in this model. Individual parameters were set to fixed values around their best-fit value, and the subset of statistically significant parameters re-minimized against hemoglobin oxygen binding data. Plots provide a graphical representation of parameter confidence intervals and correlations, and demonstrate how the two different data types - UV-vis and NMR - constrain the range of values for each parameter. This analysis further illustrates the value of hierarchically formulated models for the analysis of complex state systems, and illuminates the complexity of parameter space around the derived minimum microscopic model of hemoglobin oxygen binding.

Dissimilar flexibility of α and β subunits of human adult hemoglobin influences the protein dynamics and its alteration induced by allosteric effectors

The general question by what mechanism an "effector" molecule and the hemes of hemoglobin interact over widely separated intramolecular distances to change the oxygen affinity has been extensively investigated, and still has remained of central interest. In the present work we were interested in clarifying the general role of the protein matrix and its dynamics in the regulation of human adult hemoglobin (HbA). We used a spectroscopy approach that yields the compressibility (κ) of the protein matrix around the hemes of the subunits in HbA and studied how the binding of heterotropic allosteric effectors modify this parameter. κ is directly related to the variance of volume fluctuation, therefore it characterizes the molecular dynamics of the protein structure. For the experiments the heme groups either in the α or in the β subunits of HbA were replaced by fluorescent Zn-protoporphyrinIX, and series of fluorescence line narrowed spectra were measured at varied pressures. The evaluation of the spectra yielded the compressibility that showed significant dynamic asymmetry between the subunits: κ of the α subunit was 0.17±0.05/GPa, while for the β subunit it was much higher, 0.36±0.07/GPa. The heterotropic effectors, chloride ions, inositol hexaphosphate and bezafibrate did not cause significant changes in κ of the α subunits, while in the β subunits the effectors lead to a significant reduction down to 0.15±0.04/GPa. We relate our results to structural data, to results of recent functional studies and to those of molecular dynamics simulations, and find good agreements. The observed asymmetry in the flexibility suggests a distinct role of the subunits in the regulation of Hb that results in the observed changes of the oxygen binding capability.

Receptor Activation of HIV-1 Env Leads to Asymmetric Exposure of the gp41 Trimer

Structural rearrangements of HIV-1 glycoprotein Env promote viral entry through membrane fusion. Env is a symmetric homotrimer with each protomer composed of surface subunit gp120 and transmembrane subunit gp41. Cellular CD4- and chemokine receptor-binding to gp120 coordinate conformational changes in gp41, first to an extended prehairpin intermediate (PHI) and, ultimately, into a fusogenic trimer-of-hairpins (TOH). HIV-1 fusion inhibitors target gp41 in the PHI and block TOH formation. To characterize structural transformations into and through the PHI, we employed asymmetric Env trimers containing both high and low affinity binding sites for individual fusion inhibitors. Asymmetry was achieved using engineered Env heterotrimers composed of protomers deficient in either CD4- or chemokine receptor-binding. Linking receptor engagement to inhibitor affinity allowed us to assess conformational changes of individual Env protomers in the context of a functioning trimer. We found that the transition into the PHI could occur symmetrically or asymmetrically depending on the stoichiometry of CD4 binding. Sequential engagement of multiple CD4s promoted progressive exposure of individual fusion inhibitor binding sites in a CD4-dependent fashion. By contrast, engagement of only a single CD4 molecule led to a delayed, but symmetric, exposure of the gp41 trimer. This complex coupling between Env-CD4 interaction and gp41 exposure explained the multiphasic fusion-inhibitor titration observed for a mutant Env homotrimer with a naturally asymmetric gp41. Our results suggest that the spatial and temporal exposure of gp41 can proceed in a nonconcerted, asymmetric manner depending on the number of CD4s that engage the Env trimer. The findings have important implications for the mechanism of viral membrane fusion and the development of vaccine candidates designed to elicit neutralizing antibodies targeting gp41 in the PHI.

For HIV, cellular invasion requires merging viral and cellular membranes, an event achieved through the activity of the viral fusion protein Env. Env consists of three gp120 and three gp41 subunits symmetrically arranged on the viral surface. The gp120 subunits bind cellular receptors, which, in turn, coordinate gp41 conformational changes that promote membrane fusion. Understanding these structural rearrangements illuminates the mechanism of viral membrane fusion, and also spurs development of targeted inhibitors of viral entry and vaccine candidates that elicit antiviral immune responses. In this study, we employed a novel strategy to investigate individual subunits in the context of functioning Env complexes. The strategy links distinct gp120-receptor interactions to conformational changes that expose specific gp41 subunits. We found that, despite the initial symmetric arrangement of its subunits, Env conformational changes most often proceed quite asymmetrically, leading to exposure of only one-third of the gp41 trimer for much of the fusion event. This finding might explain why attempts to elicit potent anti-HIV antibodies to a fully exposed gp41 trimer have been largely unsuccessful. The study gives us a glimpse of the early structural transitions leading to Env-mediated membrane fusion and provides a framework for interrogating the fusion proteins of other membrane-encapsulated viruses.

An Origin of Cooperative Oxygen Binding of Human Adult Hemoglobin: Different Roles of the α and β Subunits in the α2β2 Tetramer

Human hemoglobin (Hb), which is an α2β2 tetramer and binds four O2 molecules, changes its O2-affinity from low to high as an increase of bound O2, that is characterized by 'cooperativity'. This property is indispensable for its function of O2 transfer from a lung to tissues and is accounted for in terms of T/R quaternary structure change, assuming the presence of a strain on the Fe-histidine (His) bond in the T state caused by the formation of hydrogen bonds at the subunit interfaces. However, the difference between the α and β subunits has been neglected. To investigate the different roles of the Fe-His(F8) bonds in the α and β subunits, we investigated cavity mutant Hbs in which the Fe-His(F8) in either α or β subunits was replaced by Fe-imidazole and F8-glycine. Thus, in cavity mutant Hbs, the movement of Fe upon O2-binding is detached from the movement of the F-helix, which is supposed to play a role of communication. Recombinant Hb (rHb)(αH87G), in which only the Fe-His in the α subunits is replaced by Fe-imidazole, showed a biphasic O2-binding with no cooperativity, indicating the coexistence of two independent hemes with different O2-affinities. In contrast, rHb(βH92G), in which only the Fe-His in the β subunits is replaced by Fe-imidazole, gave a simple high-affinity O2-binding curve with no cooperativity. Resonance Raman, 1H NMR, and near-UV circular dichroism measurements revealed that the quaternary structure change did not occur upon O2-binding to rHb(αH87G), but it did partially occur with O2-binding to rHb(βH92G). The quaternary structure of rHb(αH87G) appears to be frozen in T while its tertiary structure is changeable. Thus, the absence of the Fe-His bond in the α subunit inhibits the T to R quaternary structure change upon O2-binding, but its absence in the β subunit simply enhances the O2-affinity of α subunit.

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.

Structural basis for cooperative oxygen binding and bracelet-assisted assembly of Lumbricus terrestris hemoglobin

The iron-containing hemoglobins (Hbs) are essential proteins to serve as oxygen transporters in the blood. Among various kinds of Hbs, the earthworm Hbs are the champions in carrying oxygen due to not only their large size but also the unusually high cooperativity of ligand binding. However, the cooperative oxygen binding mechanisms are still mostly unknown. Here we report the cryo-electron microscopy structure of Lumbricus terrestris Hb in its native, oxygenated state at 9.1 Å resolution, showing remarkable differences from the carbon monoxide-binding X-ray structure. Our structural analysis first indicates that the cooperative ligand binding of L. terrestris Hb requires tertiary and quaternary transitions in the heme pocket and a global subunit movement facilitated by intra-ring and inter-ring contacts. Moreover, the additional sinusoidal bracelet provides the confirmation for the long-standing debate about the additional electron densities absent in the X-ray crystal structure.

Energy transfer in molecular devices

Protein machines often exhibit long-range interplay between different sites in order to achieve their biological tasks. We investigate and characterize the nonlinear energy localization and the basic mechanisms of energy transfer in protein devices. By studying two different model protein machines, with different biological functions, we show that genuinely nonlinear phenomena are responsible for energy transport between the different machine sites involved in the biological functions. The energy transfer turns out to be extremely efficient from an energetic point of view: by changing the energy initially provided to the model device, we identify a well defined range of energies where the time for the energy transport to occur is minimal and the amount of transferred energy is a maximum. Furthermore, by introducing an implicit solvent, we show that the energy is localized on the internal residues of the protein structure, thus minimizing the dissipation.

Structure and function of heme proteins in non-native states: a mini-review

Heme proteins perform various biological functions ranging from electron transfer, oxygen binding and transport, catalysis, to signaling. Although adopting proper native states is very important for these functions, progresses in representative heme proteins, including cytochrome c (cyt c), cytochrome b5 (cyt b5), myoglobin (Mb), neuroglobin (Ngb), cytochrome P450 (CYP) and heme-based sensor proteins such as CO sensor CooA, showed that various native functions, or new functions evolved, are also closely associated with non-native states. The structure and function relationship of heme proteins in non-native states is thus as important as that in native states for elucidating the precise roles of heme proteins in biological systems.

Tertiary and quaternary allostery in tetrameric hemoglobin from Scapharca inaequivalvis

The clam Scapharca inaequivalvis possesses two cooperative oxygen binding hemoglobins in its red cells: a homodimeric HbI and a heterotetrameric A2B2 HbII. Each AB dimeric half of HbII is assembled in a manner very similar to that of the well-studied HbI. This study presents crystal structures of HbII along with oxygen binding data both in the crystalline state and in wet nanoporous silica gels. Despite very similar ligand-linked structural transitions observed in HbI and HbII crystals, HbII in the crystal or encapsulated in silica gels apparently exhibits minimal cooperativity in oxygen binding, in contrast with the full cooperativity exhibited by HbI crystals. However, oxygen binding curves in the crystal indicate the presence of a significant functional inequivalence of A and B chains. When this inequivalence is taken into account, both crystal and R state gel functional data are consistent with the conservation of a tertiary contribution to cooperative oxygen binding, quantitatively similar to that measured for HbI, and are in keeping with the structural information. Furthermore, our results indicate that to fully express cooperative ligand binding, HbII requires quaternary transitions hampered by crystal lattice and gel encapsulation, revealing greater complexity in cooperative function than the direct communication across a dimeric interface observed in HbI.

The mechanism of allosteric coupling in choline kinase α1 revealed by the action of a rationally designed inhibitor

Applying a CHOK hold: Combined experimental and computational studies of the binding mode of a rationally designed inhibitor of the dimeric choline kinase α1 (CHOKα1) explain the molecular mechanism of negative cooperativity (see scheme) and how the monomers are connected. The results give insight into how the symmetry of the dimer can be partially conserved despite a lack of conservation in the static crystal structures.

Dynamic dissociating homo-oligomers and the control of protein function

Homo-oligomeric protein assemblies are known to participate in dynamic association/disassociation equilibria under native conditions, thus creating an equilibrium of assembly states. Such quaternary structure equilibria may be influenced in a physiologically significant manner either by covalent modification or by the non-covalent binding of ligands. This review follows the evolution of ideas about homo-oligomeric equilibria through the 20th and into the 21st centuries and the relationship of these equilibria to allosteric regulation by the non-covalent binding of ligands. A dynamic quaternary structure equilibria is described where the dissociated state can have alternate conformations that cannot reassociate to the original multimer; the alternate conformations dictate assembly to functionally distinct alternate multimers of finite stoichiometry. The functional distinction between different assemblies provides a mechanism for allostery. The requirement for dissociation distinguishes this morpheein model of allosteric regulation from the classical MWC concerted and KNF sequential models. These models are described alongside earlier dissociating allosteric models. The identification of proteins that exist as an equilibrium of diverse native quaternary structure assemblies has the potential to define new targets for allosteric modulation with significant consequences for further understanding and/or controlling protein structure and function. Thus, a rationale for identifying proteins that may use the morpheein model of allostery is presented and a selection of proteins for which published data suggests this mechanism may be operative are listed.

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.

Introduction: twenty five years of the Gibbs Conference on Biothermodynamics

In 2011, the Gibbs Conference on Biothermodynamics will celebrate its 25th anniversary. Since the inaugural meeting in 1987, it has brought together laboratories that lived, breathed and argued about the molecular logic of macromolecular machines. The participants have a deep commitment to understanding the nature of physico-chemical forces that govern regulation of biological systems, and share a passion for applying linkage theory. The collective goal is to understand how ligand binding, subunit assembly and conformational change drive what we observe as physiological processes such as regulated transport, enzyme cascades, gene regulation, membrane permeability, viral infection, intracellular trafficking and folding of macromolecules. In this special issue, articles by former organizers of the Gibbs Conference showcase the current breadth and depth of the field of biothermodynamics, and how thoroughly it is integrated with the study of macromolecular structures, computational modeling and physiological studies of human health and disease.

Exploring and exploiting allostery: Models, evolution, and drug targeting

The concept of allostery was elaborated almost 50years ago by Monod and coworkers to provide a framework for interpreting experimental studies on the regulation of protein function. In essence, binding of a ligand at an allosteric site affects the function at a distant site exploiting protein flexibility and reshaping protein energy landscape. Both monomeric and oligomeric proteins can be allosteric. In the past decades, the behavior of allosteric systems has been analyzed in many investigations while general theoretical models and variations thereof have been steadily proposed to interpret the experimental data. Allostery has been established as a fundamental mechanism of regulation in all organisms, governing a variety of processes that range from metabolic control to receptor function and from ligand transport to cell motility. A number of studies have shed light on how evolutionary pressures have favored and molded the development of allosteric features in specific macromolecular systems. The widespread occurrence of allostery has been recently exploited for the development and design of allosteric drugs that bind to either physiological or non-physiological allosteric sites leading to gain of function or loss of function. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

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.