Unsuspected pathway of the allosteric transition in hemoglobin

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Elucidating Dynamics of Adenylate Kinase from Enzyme Opening to Ligand Release

This study explores ligand-driven conformational changes in adenylate kinase (AK), which is known for its open-to-close conformational transitions upon ligand binding and release. By utilizing string free energy simulations, we determine the free energy profiles for both enzyme opening and ligand release and compare them with profiles from the apoenzyme. Results reveal a three-step ligand release process, which initiates with the opening of the adenosine triphosphate-binding subdomain (ATP lid), followed by ligand release and concomitant opening of the adenosine monophosphate-binding subdomain (AMP lid). The ligands then transition to nonspecific positions before complete dissociation. In these processes, the first step is energetically driven by ATP lid opening, whereas the second step is driven by ATP release. In contrast, the AMP lid opening and its ligand release make minor contributions to the total free energy for enzyme opening. Regarding the ligand binding mechanism, our results suggest that AMP lid closure occurs via an induced-fit mechanism triggered by AMP binding, whereas ATP lid closure follows conformational selection. This difference in the closure mechanisms provides an explanation with implications for the debate on ligand-driven conformational changes of AK. Additionally, we determine an X-ray structure of an AK variant that exhibits significant rearrangements in the stacking of catalytic arginines, explaining its reduced catalytic activity. In the context of apoenzyme opening, the sequence of events is different. Here, the AMP lid opens first while the ATP lid remains closed, and the free energy associated with ATP lid opening varies with orientation, aligning with the reported AK opening and closing rate heterogeneity. Finally, this study, in conjunction with our previous research, provides a comprehensive view of the intricate interplay between various structural elements, ligands, and catalytic residues that collectively contribute to the robust catalytic power of the enzyme.

Carbon monoxide separation: past, present and future

Large amounts of carbon monoxide are produced by industrial processes such as biomass gasification and steel manufacturing. The CO present in vent streams is often burnt, this produces a large amount of CO2, e.g., oxidation of CO from metallurgic flue gasses is solely responsible for 2.7% of manmade CO2 emissions. The separation of N2 from CO due to their very similar physical properties is very challenging, meaning that numerous energy-intensive steps are required for CO separation, making the CO separation from many process streams uneconomical in spite of CO being a valuable building block in the production of major chemicals through C1 chemistry and the production of linear hydrocarbons by the Fischer-Tropsch process. The development of suitable processes for the separation of carbon monoxide has both industrial and environmental significance. Especially since CO is a main product of electrocatalytic CO2 reduction, an emerging sustainable technology to enable carbon neutrality. This technology also requires an energy-efficient separation process. Therefore, there is a great need to develop energy efficient CO separation processes adequate for these different process streams. As such the urgency of separating carbon monoxide is gaining greater recognition, with research in the field becoming more and more crucial. This review details the principles on which CO separation is based and provides an overview of currently commercialised CO separation processes and their limitations. Adsorption is identified as a technology with the potential for CO separation with high selectivity and energy efficiency. We review the research efforts, mainly seen in the last decades, in developing new materials for CO separation via ad/bsorption and membrane technology. We have geared our review to both traditional CO sources and emerging CO sources, including CO production from CO2 conversion. To that end, a variety of emerging processes as potential CO2-to-CO technologies are discussed and, specifically, the need for CO capture after electrochemical CO2 reduction is highlighted, which is still underexposed in the available literature. Altogether, we aim to highlight the knowledge gaps that could guide future research to improve CO separation performance for industrial implementation.

Water dynamics around T0 vs R4 of hemoglobin from local hydrophobicity analysis

The local hydration around tetrameric hemoglobin (Hb) in its T₀ and R₄ conformational substates is analyzed based on molecular dynamics simulations. Analysis of the local hydrophobicity (LH) for all residues at the α₁β₂ and α₂β₁ interfaces, responsible for the quaternary T → R transition, which is encoded in the Monod-Wyman-Changeux model, as well as comparison with earlier computations of the solvent accessible surface area, makes clear that the two quantities measure different aspects of hydration. Local hydrophobicity quantifies the presence and structure of water molecules at the interface, whereas "buried surface" reports on the available space for solvent. For simulations with Hb frozen in its T₀ and R₄ states, the correlation coefficient between LH and buried surface is 0.36 and 0.44, respectively, but it increases considerably if the 95% confidence interval is used. The LH with Hb frozen and flexible changes little for most residues at the interfaces but is significantly altered for a few select ones: Thr41α, Tyr42α, Tyr140α, Trp37β, Glu101β (for T₀) and Thr38α, Tyr42α, Tyr140α (for R₄). The number of water molecules at the interface is found to increase by ∼25% for T₀ → R₄, which is consistent with earlier measurements. Since hydration is found to be essential to protein function, it is clear that hydration also plays an essential role in allostery.

Structural Basis of Sequential and Concerted Cooperativity

Allostery is a property of biological macromolecules featuring cooperative ligand binding and regulation of ligand affinity by effectors. The definition was introduced by Monod and Jacob in 1963, and formally developed as the "concerted model" by Monod, Wyman, and Changeux in 1965. Since its inception, this model of cooperativity was seen as distinct from and not reducible to the "sequential model" originally formulated by Pauling in 1935, which was developed further by Koshland, Nemethy, and Filmer in 1966. However, it is difficult to decide which model is more appropriate from equilibrium or kinetics measurements alone. In this paper, we examine several cooperative proteins whose functional behavior, whether sequential or concerted, is established, and offer a combined approach based on functional and structural analysis. We find that isologous, mostly helical interfaces are common in cooperative proteins regardless of their mechanism. On the other hand, the relative contribution of tertiary and quaternary structural changes, as well as the asymmetry in the liganded state, may help distinguish between the two mechanisms.

Analysis of Key Genes Responsible for Low Urea Production in Saccharomyces cerevisiae JH301

There is a potential safety risk with ethyl carbamate (EC) in Hongqu Huangjiu production; 90% of the EC in rice wine is produced by the reaction of the urea with the alcohol of Saccharomyces cerevisiae. In our previous experiments, we screened and obtained a S. cerevisiae strain JH301 that offered low urea production. However, the key genes responsible for low urea production of strain JH301 remain unclear. Here, the whole genome sequencing of S. cerevisiae strain JH301 was accomplished via a next-generation high-throughput sequencing and long-read sequencing technology. There are six main pathways related to the urea metabolism of strain JH301 based on KEGG pathway mapping. Three species-specific genes are related to the urea metabolism pathways and were found in comparative genome analysis between strains JH301 and S288c during Hongqu Huangjiu production for the first time. Finally, the ARG80 gene was found to be likely a key gene responsible for low urea production of S. cerevisiae strain JH301, as determined by PCR and qRT-PCR check analyses from DNA and RNA levers. In conclusion, the results are useful for a scientific understanding of the mechanism of low urea production by Saccharomyces cerevisiae during Hongqu Huangjiu fermentation. It also is important to control the urea and EC contents in Hongqu Huangjiu production.

The functional role of the hemoglobin-water interface

The interface between hemoglobin (Hb) and its environment, in particular water, is of great physiological relevance. Here, results from in vitro, in vivo, and computational experiments (molecular dynamics simulations) are summarized and put into perspective. One of the main findings from the computations is that the stability of the deoxy, ligand-free T-state (T₀) can be stabilized relative to the deoxy R-state (R₀) only in sufficiently large simulation boxes for the hydrophobic effect to manifest itself. This effect directly influences protein stability and is operative also under physiological conditions. Furthermore, molecular simulations provide a dynamical interpretation of the Perutz model for Hb function. Results from experiments using higher protein concentrations and realistic cellular environments are also discussed. One of the next great challenges for computational studies, which as we show is likely to be taken up in the near future, is to provide a molecular-level understanding of the dynamics of proteins in such crowded environments.

Methemoglobin formation in mutant hemoglobin α chains: electron transfer parameters and rates

Hemoglobin-mediated transport of dioxygen (O₂) critically depends on the stability of the reduced (Fe²⁺) form of the heme cofactors. Some protein mutations stabilize the oxidized (Fe³⁺) state (methemoglobin, Hb M), causing methemoglobinemia, and can be lethal above 30%. The majority of the analyses of factors influencing Hb oxidation are retrospective and give insights only for inner-sphere mutations of heme (His58, His87). Herein, we report the first all-atom molecular dynamics simulations on both redox states and calculations of the Marcus electron transfer (ET) parameters for the α chain Hb oxidation and reduction rates for Hb M. The Hb wild-type (WT) and most of the studied α chain variants maintain globin structure except the Hb M Iwate (H87Y). The mutants forming Hb M tend to have lower redox potentials and thus stabilize the oxidized (Fe³⁺) state (in particular, the Hb Miyagi variant with K61E mutation). Solvent reorganization (λ_solv 73-96%) makes major contributions to reorganization free energy, whereas protein reorganization (λ_prot) accounts for 27-30% except for the Miyagi and J-Buda variants (λ_prot ∼4%). Analysis of heme-solvent H-bonding interactions among variants provide insights into the role of Lys61 residue in stabilizing the Fe²⁺ state. Semiclassical Marcus ET theory-based calculations predict experimental k_ET for the Cyt b5-Hb complex and provide insights into relative reduction rates for Hb M in Hb variants. Thus, our methodology provides a rationale for the effect of mutations on the structure, stability, and Hb oxidation reduction rates and has potential for identification of mutations that result in methemoglobinemia.

Effect of the abolition of intersubunit salt bridges on allosteric protein structural dynamics

A salt bridge, one of the representative structural factors established by non-covalent interactions, plays a crucial role in stabilizing the structure and regulating the protein function, but its role in dynamic processes has been elusive. Here, to scrutinize the structural and functional roles of the salt bridge in the process of performing the protein function, we investigated the effects of salt bridges on the allosteric structural transition of homodimeric hemoglobin (HbI) by applying time-resolved X-ray solution scattering (TRXSS) to the K30D mutant, in which the interfacial salt bridges of the wild type (WT) are abolished. The TRXSS data of K30D are consistent with the kinetic model that requires one monomer intermediate in addition to three structurally distinct dimer intermediates (I1, I2, and I3) observed in WT and other mutants. The kinetic and structural analyses show that K30D has an accelerated biphasic transition from I2 to I3 by more than nine times compared to WT and lacks significant structural changes in the transition from R-like I2 to T-like I3 observed in WT, unveiling that the loss of the salt bridges interrupts the R-T allosteric transition of HbI. Besides, the correlation between the bimolecular CO recombination rates in K30D, WT, and other mutants reveals that the bimolecular CO recombination is abnormally decelerated in K30D, indicating that the salt bridges also affect the cooperative ligand binding in HbI. These comparisons of the structural dynamics and kinetics of K30D and WT show that the interfacial salt bridges not only assist the physical connection of two subunits but also play a critical role in the global structural signal transduction of one subunit to the other subunit via a series of well-organized structural transitions.

Water Dynamics Around Proteins: T- and R-States of Hemoglobin and Melittin

The water dynamics, as characterized by the local hydrophobicity (LH), is investigated for tetrameric hemoglobin (Hb) and dimeric melittin. For the T₀ to R₀ transition in Hb, it is found that LH provides additional molecular-level insight into the Perutz mechanism, i.e., the breaking and formation of salt bridges at the α₁/β₂ and α₂/β₁ interface is accompanied by changes in LH. For Hb in cubic water boxes with 90 and 120 Å edge length it is observed that following a decrease in LH as a consequence of reduced water density or change of water orientation at the protein/water interface the α/β interfaces are destabilized; this is a hallmark of the Perutz stereochemical model for the T to R transition in Hb. The present work thus provides a dynamical view of the classical structural model relevant to the molecular foundations of Hb function. For dimeric melittin, earlier results by Cheng and Rossky [ Nature 1998, 392, 696-699] are confirmed and interpreted on the basis of LH from simulations in which the protein structure is frozen. For the flexible melittin dimer, the changes in the local hydration can be as much as 30% greater than for the rigid dimer, reflecting the fact that protein and water dynamics are coupled.

Functional (un)cooperativity in elevator transport proteins

The activity of enzymes is subject to regulation at multiple levels. Cooperativity, the interconnected behavior of active sites within a protein complex, directly affects protein activity. Cooperativity is a mode of regulation that requires neither extrinsic factors nor protein modifications. Instead, it allows enzymes themselves to modulate reaction rates. Cooperativity is an important regulatory mechanism in soluble proteins, but also examples of cooperative membrane proteins have been described. In this review, we summarize the current knowledge on interprotomer cooperativity in elevator-type proteins, a class of membrane transporters characterized by large rigid-body movements perpendicular to the membrane, and highlight well-studied examples and experimental approaches.

More than a Confinement: “Soft” and “Hard” Enzyme Entrapment Modulates Biological Catalyst Function

Catalysis makes chemical and biochemical reactions kinetically accessible. From a technological point of view, organic, inorganic, and biochemical catalysis is relevant for several applications, from industrial synthesis to biomedical, material, and food sciences. A heterogeneous catalyst, i.e., a catalyst confined in a different phase with respect to the reagents’ phase, requires either its physical confinement in an immobilization matrix or its physical adsorption on a surface. In this review, we will focus on the immobilization of biological catalysts, i.e., enzymes, by comparing hard and soft immobilization matrices and their effect on the modulation of the catalysts’ function. Indeed, unlike smaller molecules, the catalytic activity of protein catalysts depends on their structure, conformation, local environment, and dynamics, properties that can be strongly affected by the immobilization matrices, which, therefore, not only provide physical confinement, but also modulate catalysis.

Phthalides serve as potent modulators to boost fetal hemoglobin induction therapy for β-hemoglobinopathies

Fetal hemoglobin (HbF) induction therapy has become the most promising strategy for treating β-hemoglobinopathies, including sickle-cell diseases and β-thalassemia. However, subtle but critical structural difference exists between HbF and normal adult hemoglobin (HbA), which inevitably leads to reduced binding of the endogenous modulator 2,3-bisphosphoglycerate (2,3-BPG) to HbF and thus increased oxygen affinity and decreased oxygen transport efficiency of HbF. We combined the oxygen equilibrium experiments, resonance Raman (RR) spectroscopy, and molecular docking modeling, and we discuss 2 phthalides, z-butylidenephthalide and z-ligustilide, that can effectively lower the oxygen affinity of HbF. They adjust it to a level closer to that of HbA and make it a more satisfactory oxygen carrier for adults. From the oxygen equilibrium curve measurements, we show that the 2 phthalides are more effective than 2,3-BPG for modulating HbF. The RR spectra show that phthalides allosterically stabilize the oxygenated HbF in the low oxygen affinity conformation, and the molecular docking modeling reveals that the 2 chosen phthalides interact with HbF via the cleft around the γ1/γ2 interface with a binding strength ∼1.6 times stronger than that of 2,3-BPG. We discuss the implications of z-butylidenephthalide and z-ligustilide in boosting the efficacy of HbF induction therapy to mitigate the clinical severities of β-hemoglobinopathies.

Progress in Allosteric Database

An allosteric mechanism refers to the biological regulation process wherein macromolecules propagate the effect of ligand binding at one site to a spatially distant orthosteric locus, thus affecting activity. The theory has remained a trending topic in biology research for over 50 years, since the understanding of allostery is fundamental for gleaning numerous biological processes and developing new drug therapies. In the past two decades, the allosteric paradigm has evolved into more descriptive models, with ever-expanding amounts of experimental data pertaining to newly identified allosteric molecules. The AlloSteric Database (ASD, accessible at http://mdl.shsmu.edu.cn/ASD ), which is a comprehensive knowledge repository, has provided the public with integrated information encompassing allosteric proteins, modulators, sites, pathways, and networks to investigate allostery since 2009. In this chapter, we introduce the history and usage of the ASD and give attention to specific applications that have benefited from the ASD.

Gears-In-Motion: The Interplay of WW and PPIase Domains in Pin1

Pin1 belongs to the family of the peptidyl-prolyl cis-trans isomerase (PPIase), which is a class of enzymes that catalyze the cis/trans isomerization of the Proline residue. Pin1 is unique and only catalyzes the phosphorylated Serine/Threonine-Proline (S/T-P) motifs of a subset of proteins. Since the discovery of Pin1 as a key protein in cell cycle regulation, it has been implicated in numerous diseases, ranging from cancer to neurodegenerative diseases. The main features of Pin1 lies in its two main domains: the WW (two conserved tryptophan) domain and the PPIase domain. Despite extensive studies trying to understand the mechanisms of Pin1 functions, how these two domains contribute to the biological roles of Pin1 in cellular signaling requires more investigations. The WW domain of Pin1 is known to have a higher affinity to its substrate than that of the PPIase domain. Yet, the WW domain seems to prefer the trans configuration of phosphorylated S/T-P motif, while the PPIase catalyzes the cis to trans isomerasion. Such contradicting information has generated much confusion as to the actual mechanism of Pin1 function. In addition, dynamic allostery has been suggested to be important for Pin1 function. Henceforth, in this review, we will be looking at the progress made in understanding the function of Pin1, and how these understandings can aid us in overcoming the diseases implicated by Pin1 such as cancer during drug development.

Generating conformational transition paths with low potential-energy barriers for proteins

The knowledge of conformational transition paths in proteins can be useful for understanding protein mechanisms. Recently, we have introduced the As-Rigid-As-Possible (ARAP) interpolation method, for generating interpolation paths between two protein conformations. The method was shown to preserve well the rigidity of the initial conformation along the path. However, because the method is totally geometry-based, the generated paths may be inconsistent because the atom interactions are ignored. Therefore, in this article, we would like to introduce a new method to generate conformational transition paths with low potential-energy barriers for proteins. The method is composed of three processing stages. First, ARAP interpolation is used for generating an initial path. Then, the path conformations are enhanced by a clash remover. Finally, Nudged Elastic Band, a path-optimization method, is used to produce a low-energy path. Large energy reductions are found in the paths obtained from the method than in those obtained from the ARAP interpolation method alone. The results also show that ARAP interpolation is a good candidate for generating an initial path because it leads to lower potential-energy paths than two other common methods for path interpolation.

Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins

With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.

Tertiary dynamics of human adult hemoglobin fixed in R and T quaternary structures

Protein dynamics of human adult hemoglobin and its mutants restricted in R and T quaternary states following ligand photolysis were studied by time-resolved resonance Raman spectroscopy. In the time-resolved spectra, we observed spectral changes of in-plane stretching modes of heme and the iron-histidine stretching mode of the Fe-His bond for all the hemoglobin samples. The βD99N mutant, which adopts the R state in both the ligand-bound and the deoxy forms, showed similar temporal behaviors in time-resolved resonance Raman spectra as wild-type recombinant hemoglobin until 10 μs, consistent with the fact that the mutant undergoes only the tertiary structural changes in the R state. The βN102T mutant, which adopts the T state in both the ligand-bound and the deoxy forms, showed much slower tertiary structural changes, suggesting that the EF helical motion is decelerated by the change of the intersubunit interactions. The present data indicate that the allosteric kinetic response between the interhelical hydrogen bonds of the EF helices and the intersubunit hydrogen bonds is bidirectional. The implications of these results for understanding the allosteric pathway of Hb are discussed in detail.

Hinge-Shift Mechanism Modulates Allosteric Regulations in Human Pin1

Allostery, which is regulation from distant sites, plays a major role in biology. While traditional allostery is described in terms of conformational change upon ligand binding as an underlying principle, it is possible to have allosteric regulations without significant conformational change through modulating the conformational dynamics by altering the local effective elastic modulus of the protein upon ligand binding. Pin1 utilizes this dynamic allostery to regulate its function. It is a modular protein containing a WW domain and a larger peptidyl prolyl isomerase domain (PPIase) that isomerizes phosphoserine/threonine-proline (pS/TP) motifs. The WW domain serves as a docking module, whereas catalysis solely takes place within the PPIase domain. Here, we analyze the change in the dynamic flexibility profile of the PPIase domain upon ligand binding to the WW domain. Substrate binding to the WW domain induces the formation of a new rigid hinge site around the interface of the two domains and loosens the flexibility of a rigid site existing in the Apo form around the catalytic site. This hinge-shift mechanism enhances the dynamic coupling of the catalytic positions with the PPIase domain, where the rest of the domain can cooperatively respond to the local conformational changes around the catalytic site, leading to an increase in catalytic efficiency.

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.

Ab initio sampling of transition paths by conditioned Langevin dynamics

We propose a novel stochastic method to generate Brownian paths conditioned to start at an initial point and end at a given final point during a fixed time t_f under a given potential U(x). These paths are sampled with a probability given by the overdamped Langevin dynamics. We show that these paths can be exactly generated by a local stochastic partial differential equation. This equation cannot be solved in general but we present several approximations that are valid either in the low temperature regime or in the presence of barrier crossing. We show that this method warrants the generation of statistically independent transition paths. It is computationally very efficient. We illustrate the method first on two simple potentials, the two-dimensional Mueller potential and the Mexican hat potential, and then on the multi-dimensional problem of conformational transitions in proteins using the "Mixed Elastic Network Model" as a benchmark.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -“r” and “t”-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

Phthalide Derivatives from Angelica Sinensis Decrease Hemoglobin Oxygen Affinity: A New Allosteric-Modulating Mechanism and Potential Use as 2,3-BPG Functional Substitutes

Angelica sinensis (AS), one of the most versatile herbal medicines remains widely used due to its multi-faceted pharmacologic activities. Besides its traditional use as the blood-nourishing tonic, its anti-hypertensive, anti-cardiovascular, neuroprotective and anti-cancer effects have been reported. Albeit the significant therapeutic effects, how AS exerts such diverse efficacies from the molecular level remains elusive. Here we investigate the influences of AS and four representative phthalide derivatives from AS on the structure and function of hemoglobin (Hb). From the spectroscopy and oxygen equilibrium experiments, we show that AS and the chosen phthalides inhibited the oxygenated Hb from transforming into the high-affinity “relaxed” (R) state, decreasing Hb’s oxygen affinity. It reveals that phthalides cooperate with the endogenous Hb modulator, 2,3-bisphosphoglycerate (2,3-BPG) to synergetically regulate Hb allostery. From the docking modeling, phthalides appear to interact with Hb mainly through its α₁/α₂ interface, likely strengthening four (out of six) Hb “tense” (T) state stabilizing salt-bridges. A new allosteric-modulating mechanism is proposed to rationalize the capacity of phthalides to facilitate Hb oxygen transport, which may be inherently correlated with the therapeutic activities of AS. The potential of phthalides to serve as 2,3-BPG substitutes/supplements and their implications in the systemic biology and preventive medicine are discussed.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Conformational Transitions and Alternating-Access Mechanism in the Sarcoplasmic Reticulum Calcium Pump

Ion pumps are integral membrane proteins responsible for transporting ions against concentration gradients across biological membranes. Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA), a member of the P-type ATPases family, transports two calcium ions per hydrolyzed ATP molecule via an "alternating-access" mechanism. High-resolution crystallographic structures provide invaluable insight on the structural mechanism of the ion pumping process. However, to understand the molecular details of how ATP hydrolysis is coupled to calcium transport, it is necessary to gain knowledge about the conformational transition pathways connecting the crystallographically resolved conformations. Large-scale transitions in SERCA occur at time-scales beyond the current reach of unbiased molecular dynamics simulations. Here, we overcome this challenge by employing the string method, which represents a transition pathway as a chainofstates linking two conformational endpoints. Using a multiscale methodology, we have determined all-atom transition pathways for three main conformational transitions responsible for the alternating-access mechanism. The present pathways provide a clear chronology and ordering of the key events underlying the active transport of calcium ions by SERCA. Important conclusions are that the conformational transition that leads to occlusion with bound ATP and calcium is highly concerted and cooperative, the phosphorylation of Asp351 causes areorganization of the cytoplasmic domains that subsequently drives the opening of the luminal gate, and thereclosing of luminal gate induces a shift in the cytoplasmic domains that subsequently enables the dephosphorylation of Asp351-P. Formation of transient residue-residue contacts along the conformational transitions predicted by the computations provide an experimental route to test the general validity of the computational pathways.

Role of Heme Pocket Water in Allosteric Regulation of Ligand Reactivity in Human Hemoglobin

Water molecules can enter the heme pockets of unliganded myoglobins and hemoglobins, hydrogen bond with the distal histidine, and introduce steric barriers to ligand binding. The spectrokinetics of photodissociated CO complexes of human hemoglobin and its isolated α and β chains were analyzed for the effect of heme hydration on ligand rebinding. A strong coupling was observed between heme hydration and quaternary state. This coupling may contribute significantly to the 20-60-fold difference between the R- and T-state bimolecular CO binding rate constants and thus to the modulation of ligand reactivity that is the hallmark of hemoglobin allostery. Heme hydration proceeded over the course of several kinetic phases in the tetramer, including the R to T quaternary transition. An initial 150 ns hydration phase increased the R-state distal pocket water occupancy, nw(R), to a level similar to that of the isolated α (∼60%) and β (∼10%) chains, resulting in a modest barrier to ligand binding. A subsequent phase, concurrent with the first step of the R → T transition, further increased the level of heme hydration, increasing the barrier. The final phase, concurrent with the final step of the allosteric transition, brought the water occupancy of the T-state tetramer, nw(T), even higher and close to full occupancy in both the α and β subunits (∼90%). This hydration level could present an even larger barrier to ligand binding and contribute significantly to the lower iron reactivity of the T state toward CO.

The allosteric switching mechanism in bacteriophage MS2

We use all-atom simulations to elucidate the mechanisms underlying conformational switching and allostery within the coat protein of the bacteriophage MS2. Assembly of most icosahedral virus capsids requires that the capsid protein adopts different conformations at precise locations within the capsid. It has been shown that a 19 nucleotide stem loop (TR) from the MS2 genome acts as an allosteric effector, guiding conformational switching of the coat protein during capsid assembly. Since the principal conformational changes occur far from the TR binding site, it is important to understand the molecular mechanism underlying this allosteric communication. To this end, we use all-atom simulations with explicit water combined with a path sampling technique to sample the MS2 coat protein conformational transition, in the presence and absence of TR-binding. The calculations find that TR binding strongly alters the transition free energy profile, leading to a switch in the favored conformation. We discuss changes in molecular interactions responsible for this shift. We then identify networks of amino acids with correlated motions to reveal the mechanism by which effects of TR binding span the protein. We find that TR binding strongly affects residues located at the 5-fold and quasi-sixfold interfaces in the assembled capsid, suggesting a mechanism by which the TR binding could direct formation of the native capsid geometry. The analysis predicts amino acids whose substitution by mutagenesis could alter populations of the conformational substates or their transition rates.

The time scale of the quaternary structural changes in hemoglobin revealed using the transient grating technique

The quaternary structural transition between the R and T states of human hemoglobin was investigated using the transient grating technique. The results presented herein reveal that the quaternary structural change accompanied by the R-T transition occurs within a few microseconds.

From protein structure to function via single crystal optical spectroscopy

The more than 100,000 protein structures determined by X-ray crystallography provide a wealth of information for the characterization of biological processes at the molecular level. However, several crystallographic “artifacts,” including conformational selection, crystallization conditions and radiation damages, may affect the quality and the interpretation of the electron density maps, thus limiting the relevance of structure determinations. Moreover, for most of these structures, no functional data have been obtained in the crystalline state, thus posing serious questions on their validity in infereing protein mechanisms. In order to solve these issues, spectroscopic methods have been applied for the determination of equilibrium and kinetic properties of proteins in the crystalline state. These methods are UV-vis spectrophotometry, spectrofluorimetry, IR, EPR, Raman, and resonance Raman spectroscopy. Some of these approaches have been implemented with on-line instruments at X-ray synchrotron beamlines. Here, we provide an overview of investigations predominantly carried out in our laboratory by single crystal polarized absorption UV-vis microspectrophotometry, the most applied technique for the functional characterization of proteins in the crystalline state. Studies on hemoglobins, pyridoxal 5′-phosphate dependent enzymes and green fluorescent protein in the crystalline state have addressed key biological issues, leading to either straightforward structure-function correlations or limitations to structure-based mechanisms.

Two pathways mediate interdomain allosteric regulation in pin1

Allostery is an essential means for regulating biomolecular functions and provides unique opportunities for drug design, yet our ability to elucidate allosteric mechanisms remains limited. Here, based on extensive molecular dynamics simulations, we present an atomistic picture of the pathways mediating the allosteric regulation of the PPIase domain of Pin1 by its WW domain. Two pathways jointly propagate the action of substrate-WW binding to produce closure and rigidification of three PPIase catalytic-site loops. One pathway preexists in the apo protein, but remains dormant until substrate-WW binding completes the second. The reduction in conformational entropy and preorganization of the catalytic-site loops observed here may explain why substrate-WW binding enhances ligand affinity and catalytic activity of the PPIase domain and suggest a combination drug therapy for Pin1-related diseases. Whereas the traditional view of allostery has emphasized conformational transition, our study identifies a distinct role of conformational dynamics in eliciting allostery.

Differential control of heme reactivity in alpha and beta subunits of hemoglobin: a combined Raman spectroscopic and computational study

The use of hybrid hemoglobin (Hb), with mesoheme substituted for protoheme, allows separate monitoring of the α or β hemes along the allosteric pathway. Using resonance Raman (rR) spectroscopy in silica gel, which greatly slows protein motions, we have observed that the Fe-histidine stretching frequency, νFeHis, which is a monitor of heme reactivity, evolves between frequencies characteristic of the R and T states, for both α or β chains, prior to the quaternary R-T and T-R shifts. Computation of νFeHis, using QM/MM and the conformational search program PELE, produced remarkable agreement with experiment. Analysis of the PELE structures showed that the νFeHis shifts resulted from heme distortion and, in the α chain, Fe-His bond tilting. These results support the tertiary two-state model of ligand binding (Henry et al., Biophys. Chem. 2002, 98, 149). Experimentally, the νFeHis evolution is faster for β than for α chains, and pump-probe rR spectroscopy in solution reveals an inflection in the νFeHis time course at 3 μs for β but not for α hemes, an interval previously shown to be the first step in the R-T transition. In the α chain νFeHis dropped sharply at 20 μs, the final step in the R-T transition. The time courses are fully consistent with recent computational mapping of the R-T transition via conjugate peak refinement by Karplus and co-workers (Fischer et al., Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 5608). The effector molecule IHP was found to lower νFeHis selectively for α chains within the R state, and a binding site in the α1α2 cleft is suggested.

Using Markov state models to study self-assembly

Markov state models (MSMs) have been demonstrated to be a powerful method for computationally studying intramolecular processes such as protein folding and macromolecular conformational changes. In this article, we present a new approach to construct MSMs that is applicable to modeling a broad class of multi-molecular assembly reactions. Distinct structures formed during assembly are distinguished by their undirected graphs, which are defined by strong subunit interactions. Spatial inhomogeneities of free subunits are accounted for using a recently developed Gaussian-based signature. Simplifications to this state identification are also investigated. The feasibility of this approach is demonstrated on two different coarse-grained models for virus self-assembly. We find good agreement between the dynamics predicted by the MSMs and long, unbiased simulations, and that the MSMs can reduce overall simulation time by orders of magnitude.

A theoretical view of protein dynamics

Proteins are fascinating supramolecular structures, which are able to recognize ligands transforming binding information into chemical signals. They can transfer information across the cell, can catalyse complex chemical reactions, and are able to transform energy into work with much more efficiency than any human engine. The unique abilities of proteins are tightly coupled with their dynamic properties, which are coded in a complex way in the sequence and carefully refined by evolution. Despite its importance, our experimental knowledge of protein dynamics is still rather limited, and mostly derived from theoretical calculations. I will review here, in a systematic way, the current state-of-the-art theoretical approaches to the study of protein dynamics, emphasizing the most recent advances, examples of use and the expected lines of development in the near future.

A unified view of "how allostery works"

The question of how allostery works was posed almost 50 years ago. Since then it has been the focus of much effort. This is for two reasons: first, the intellectual curiosity of basic science and the desire to understand fundamental phenomena, and second, its vast practical importance. Allostery is at play in all processes in the living cell, and increasingly in drug discovery. Many models have been successfully formulated, and are able to describe allostery even in the absence of a detailed structural mechanism. However, conceptual schemes designed to qualitatively explain allosteric mechanisms usually lack a quantitative mathematical model, and are unable to link its thermodynamic and structural foundations. This hampers insight into oncogenic mutations in cancer progression and biased agonists' actions. Here, we describe how allostery works from three different standpoints: thermodynamics, free energy landscape of population shift, and structure; all with exactly the same allosteric descriptors. This results in a unified view which not only clarifies the elusive allosteric mechanism but also provides structural grasp of agonist-mediated signaling pathways, and guides allosteric drug discovery. Of note, the unified view reasons that allosteric coupling (or communication) does not determine the allosteric efficacy; however, a communication channel is what makes potential binding sites allosteric.

Reaction trajectory revealed by a joint analysis of protein data bank

Structural motions along a reaction pathway hold the secret about how a biological macromolecule functions. If each static structure were considered as a snapshot of the protein molecule in action, a large collection of structures would constitute a multidimensional conformational space of an enormous size. Here I present a joint analysis of hundreds of known structures of human hemoglobin in the Protein Data Bank. By applying singular value decomposition to distance matrices of these structures, I demonstrate that this large collection of structural snapshots, derived under a wide range of experimental conditions, arrange orderly along a reaction pathway. The structural motions along this extensive trajectory, including several helical transformations, arrive at a reverse engineered mechanism of the cooperative machinery (Ren, companion article), and shed light on pathological properties of the abnormal homotetrameric hemoglobins from α-thalassemia. This method of meta-analysis provides a general approach to structural dynamics based on static protein structures in this post genomics era.

Collective dynamics underlying allosteric transitions in hemoglobin

Hemoglobin is the prototypic allosteric protein. Still, its molecular allosteric mechanism is not fully understood. To elucidate the mechanism of cooperativity on an atomistic level, we developed a novel computational technique to analyse the coupling of tertiary and quaternary motions. From Molecular Dynamics simulations showing spontaneous quaternary transitions, we separated the transition trajectories into two orthogonal sets of motions: one consisting of intra-chain motions only (referred to as tertiary-only) and one consisting of global inter-chain motions only (referred to as quaternary-only). The two underlying subspaces are orthogonal by construction and their direct sum is the space of full motions. Using Functional Mode Analysis, we were able to identify a collective coordinate within the tertiary-only subspace that is correlated to the most dominant motion within the quaternary-only motions, hence providing direct insight into the allosteric coupling mechanism between tertiary and quaternary conformation changes. This coupling-motion is substantially different from tertiary structure changes between the crystallographic structures of the T- and R-state. We found that hemoglobin's allosteric mechanism of communication between subunits is equally based on hydrogen bonds and steric interactions. In addition, we were able to affect the T-to-R transition rates by choosing different histidine protonation states, thereby providing a possible atomistic explanation for the Bohr effect.

A gating mechanism of pentameric ligand-gated ion channels

Pentameric ligand-gated ion channels (pLGICs) play a central role in intercellular communication in the nervous system and are involved in fundamental processes such as attention, learning, and memory. They are oligomeric protein assemblies that convert a chemical signal into an ion flux through the postsynaptic membrane, but the molecular mechanism of gating ions has remained elusive. Here, we present atomistic molecular dynamics simulations of the prokaryotic channels from Gloeobacter violaceus (GLIC) and Erwinia chrysanthemi (ELIC), whose crystal structures are thought to represent the active and the resting states of pLGICs, respectively, and of the eukaryotic glutamate-gated chloride channel from Caenorhabditis elegans (GluCl), whose open-channel structure was determined complexed with the positive allosteric modulator ivermectin. Structural observables extracted from the trajectories of GLIC and ELIC are used as progress variables to analyze the time evolution of GluCl, which was simulated in the absence of ivermectin starting from the structure with bound ivermectin. The trajectory of GluCl with ivermectin removed shows a sequence of structural events that couple agonist unbinding from the extracellular domain to ion-pore closing in the transmembrane domain. Based on these results, we propose a structural mechanism for the allosteric communication leading to deactivation/activation of the GluCl channel. This model of gating emphasizes the coupling between the quaternary twisting and the opening/closing of the ion pore and is likely to apply to other members of the pLGIC family.

Mapping polymerization and allostery of hemoglobin S using point mutations

Hemoglobin is a complex system that undergoes conformational changes in response to oxygen, allosteric effectors, mutations, and environmental changes. Here, we study allostery and polymerization of hemoglobin and its variants by application of two previously described methods: (i) AllosMod for simulating allostery dynamics given two allosterically related input structures and (ii) a machine-learning method for dynamics- and structure-based prediction of the mutation impact on allostery (Weinkam et al. J. Mol. Biol. 2013, 425, 647-661), now applicable to systems with multiple coupled binding sites, such as hemoglobin. First, we predict the relative stabilities of substates and microstates of hemoglobin, which are determined primarily by entropy within our model. Next, we predict the impact of 866 annotated mutations on hemoglobin's oxygen binding equilibrium. We then discuss a subset of 30 mutations that occur in the presence of the sickle cell mutation and whose effects on polymerization have been measured. Seven of these HbS mutations occur in three predicted druggable binding pockets that might be exploited to directly inhibit polymerization; one of these binding pockets is not apparent in the crystal structure, but only in structures generated by AllosMod. For the 30 mutations, we predict that mutation-induced conformational changes within a single tetramer tend not to significantly impact polymerization; instead, these mutations more likely impact polymerization by directly perturbing a polymerization interface. Finally, our analysis of allostery allows us to hypothesize why hemoglobin evolved to have multiple subunits and a persistent low frequency sickle cell mutation.

The concept of allosteric interaction and its consequences for the chemistry of the brain

Throughout this Reflections article, I have tried to follow up on the genesis in the 1960s and subsequent evolution of the concept of allosteric interaction and to examine its consequences within the past decades, essentially in the field of the neuroscience. The main conclusion is that allosteric mechanisms built on similar structural principles operate in bacterial regulatory enzymes, gene repressors (and the related nuclear receptors), rhodopsin, G-protein-coupled receptors, neurotransmitter receptors, ion channels, and so on from prokaryotes up to the human brain yet with important features of their own. Thus, future research on these basic cybernetic sensors is expected to develop in two major directions: at the elementary level, toward the atomic structure and molecular dynamics of the conformational changes involved in signal recognition and transduction, but also at a higher level of organization, the contribution of allosteric mechanisms to the modulation of brain functions.

Comparative analysis of oxy-hemoglobin and aquomet-hemoglobin by hydrogen/deuterium exchange mass spectrometry

The function of hemoglobin (Hb) as oxygen transporter is mediated by reversible O2 binding to Fe(2+) heme in each of the α and β subunits. X-ray crystallography revealed different subunit arrangements in oxy-Hb and deoxy-Hb. The deoxy state is stabilized by additional contacts, causing a rigidification that results in strong protection against hydrogen/deuterium exchange (HDX). Aquomet-Hb is a dysfunctional degradation product with four water-bound Fe(3+) centers. Heme release from aquomet-Hb is relatively facile, triggering oxidative damage of membrane lipids. Aquomet-Hb crystallizes in virtually the same conformation as oxy-Hb. Hence, it is commonly implied that the solution-phase properties of aquomet-Hb should resemble those of the oxy state. This work compares the structural dynamics of oxy-Hb and aquomet-Hb by HDX mass spectrometry (MS). It is found that the aquomet state exhibits a solution-phase structure that is significantly more dynamic, as manifested by elevated HDX levels. These enhanced dynamics affect the aquomet α and β subunits in a different fashion. The latter undergoes global destabilization, whereas the former shows elevated HDX levels only in the heme binding region. It is proposed that these enhanced dynamics play a role in facilitating heme release from aquomet-Hb. Our findings should be of particular interest to the MS community because oxy-Hb and aquomet-Hb serve as widely used test analytes for probing the relationship between biomolecular structure in solution and in the gas phase. We are not aware of any prior comparative HDX/MS experiments on oxy-Hb and aquomet-Hb.