Speed of intersubunit communication in proteins

Allosteric modulation of α1β3γ2 GABAA receptors by farnesol through the neurosteroid sites

In mammalian cells, all-trans farnesol, a 15-carbon isoprenol, is a product of the mevalonate pathway. It is the natural substrate of alcohol dehydrogenase and a substrate for CYP2E1, two enzymes implicated in ethanol metabolism. Studies have shown that farnesol is present in the human brain and inhibits voltage-gated Ca2+ channels at much lower concentrations than ethanol. Here we show that farnesol modulates the activity of γ-aminobutyric acid type A receptors (GABAARs), some of which also mediate the sedative activity of ethanol. Electrophysiology experiments performed in HEK cells expressing human α1β3γ2 or α6β3γ2 GABAARs revealed that farnesol increased chloride currents through positive allosteric modulation of these receptors and showed dependence on both the alcoholic functional group of farnesol and the length of the alkyl chain for activity. In silico studies using long-timescale unbiased all-atom molecular dynamics (MD) simulations of the human α1β3γ2 GABAA receptors revealed that farnesol modulates the channel by directly binding to the transmembrane neurosteroid-binding site, after partitioning into the surrounding membrane and reaching the receptor by lateral diffusion. Channel activation by farnesol was further characterized by several structural and dynamic variables, such as global twisting of the receptor's extracellular domain, tilting of the transmembrane M2 helices, radius, cross-sectional area, hydration status, and electrostatic potential of the channel pore. Our results expand the pharmacological activities of farnesol to yet another class of ion channels implicated in neurotransmission, thus providing a novel path for understanding and treatment of diseases involving GABAA receptor dysfunction.

A retrospective on statistical mechanical models for hemoglobin allostery

Understanding allosteric interactions in proteins has become one of the major research areas in protein science. The original aim of the famous theoretical model of Monod, Wyman, and Changeux (MWC) was to explain the regulation of enzymatic activity in biochemical pathways. However, its first successful quantitative application was to explain cooperative oxygen binding by hemoglobin, often called the "hydrogen molecule of biology." The combination of its original application and the enormous amount of research on hemoglobin has made it the paradigm for studies of allostery, especially for multi-subunit proteins, and for the development of statistical mechanical models to describe how structure determines function. This article is a historical account of the development of statistical mechanical models for hemoglobin to explain both the cooperative binding of oxygen (called homotropic effects by MWC) and how oxygen binding is affected by ligands that bind distant from the heme oxygen binding site (called heterotropic allosteric effects by MWC). This account makes clear the many remaining challenges for describing the relationship of structure to function for hemoglobin in terms of a satisfactory statistical mechanical model.

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

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

Role of hemoglobin structural-functional relationships in oxygen transport

The molecular mechanism of O₂ binding to hemoglobin (Hb) has been critically reviewed on the basis of the information built up in the last decades. It allows to describe in detail from the kinetic and thermodynamic viewpoint the process of O₂ uptake in the lungs and release to the tissues, casting some light on the physiological and pathological aspects of this process. The relevance of structural-functional relationships for O₂ binding is particularly outlined in the case of poorly vascularized tissues, such as retina, briefly discussing of strategies employed for optimization of oxygen supply to this type of tissues.

How fast are the motions of tertiary-structure elements in proteins?

Protein motions occur on multiple time and distance scales. Large-scale motions of protein tertiary-structure elements, i.e., domains, are particularly intriguing as they are essential for the catalytic activity of many enzymes and for the functional cycles of protein machines and motors. Theoretical estimates suggest that domain motions should be very fast, occurring on the nanosecond or microsecond time scales. Indeed, free-energy barriers for domain motions are likely to involve salt bridges, which can break in microseconds. Experimental methods that can directly probe domain motions on fast time scales have appeared only in recent years. This Perspective discusses briefly some of these techniques, including nuclear magnetic resonance and single-molecule fluorescence spectroscopies. We introduce a few recent studies that demonstrate ultrafast domain motions and discuss their potential roles. Particularly surprising is the observation of tertiary-structure element dynamics that are much faster than the functional cycles in some protein machines. These swift motions can be rationalized on a case-by-case basis. For example, fast domain closure in multi-substrate enzymes may be utilized to optimize relative substrate orientation. Whether a large mismatch in time scales of conformational dynamics vs functional cycles is a general design principle in proteins remains to be determined.

The Greenland shark Somniosus microcephalus-Hemoglobins and ligand-binding properties

A large amount of data is currently available on the adaptive mechanisms of polar bony fish hemoglobins, but structural information on those of cartilaginous species is scarce. This study presents the first characterisation of the hemoglobin system of one of the longest-living vertebrate species (392 ± 120 years), the Arctic shark Somniosus microcephalus. Three major hemoglobins are found in its red blood cells and are made of two copies of the same α globin combined with two copies of three very similar β subunits. The three hemoglobins show very similar oxygenation and carbonylation properties, which are unaffected by urea, a very important compound in marine elasmobranch physiology. They display identical electronic absorption and resonance Raman spectra, indicating that their heme-pocket structures are identical or highly similar. The quaternary transition equilibrium between the relaxed (R) and the tense (T) states is more dependent on physiological allosteric effectors than in human hemoglobin, as also demonstrated in polar teleost hemoglobins. Similar to other cartilaginous fishes, we found no evidence for functional differentiation among the three isoforms. The very similar ligand-binding properties suggest that regulatory control of O₂ transport may be at the cellular level and that it may involve changes in the cellular concentrations of allosteric effectors and/or variations of other systemic factors. The hemoglobins of this polar shark have evolved adaptive decreases in O₂ affinity in comparison to temperate sharks.

Direct observation of conformational population shifts in crystalline human hemoglobin

Although X-ray crystallography is the most commonly used technique for studying the molecular structure of proteins, it is not generally able to monitor the dynamic changes or global domain motions that often underlie allostery. These motions often prevent crystal growth or reduce crystal order. We have recently discovered a crystal form of human hemoglobin that contains three protein molecules allowed to express a full range of quaternary structures, whereas maintaining strong X-ray diffraction. Here we use this crystal form to investigate the effects of two allosteric effectors, phosphate and bezafibrate, by tracking the structures and functions of the three hemoglobin molecules following the addition of each effector. The X-ray analysis shows that the addition of either phosphate or bezafibrate not only induces conformational changes in a direction from a relaxed-state to a tense-state, but also within relaxed-state populations. The microspectrophotometric O₂ equilibrium measurements on the crystals demonstrate that the binding of each effector energetically stabilizes the lowest affinity conformer more strongly than the intermediate affinity one, thereby reducing the O₂ affinity of tense-state populations, and that the addition of bezafibrate causes an ∼5-fold decrease in the O₂ affinity of relaxed-state populations. These results show that the allosteric pathway of hemoglobin involves shifts of populations rather than a unidirectional conversion of one quaternary structure to another, and that minor conformers of hemoglobin may have a disproportionate effect on the overall O₂ affinity.

A mechanism for acetylcholine receptor gating based on structure, coupling, phi, and flip

Gupta et al. use single-channel electrophysiology to investigate the gating mechanism of acetylcholine receptor ion channels. They propose that channel opening starts at the M2–M3 linker and ligand-binding sites and proceeds through four brief intermediate conformations before ending with the collapse of a gate bubble.

Nicotinic acetylcholine receptors are allosteric proteins that generate membrane currents by isomerizing (“gating”) between resting and active conformations under the influence of neurotransmitters. Here, to explore the mechanisms that link the transmitter-binding sites (TBSs) with the distant gate, we use mutant cycle analyses to measure coupling between residue pairs, phi value analyses to sequence domain rearrangements, and current simulations to reproduce a microsecond shut component (“flip”) apparent in single-channel recordings. Significant interactions between amino acids separated by >15 Å are rare; an exception is between the αM2–M3 linkers and the TBSs that are ∼30 Å apart. Linker residues also make significant, local interactions within and between subunits. Phi value analyses indicate that without agonists, the linker is the first region in the protein to reach the gating transition state. Together, the phi pattern and flip component suggest that a complete, resting↔active allosteric transition involves passage through four brief intermediate states, with brief shut events arising from sojourns in all or a subset. We derive energy landscapes for gating with and without agonists, and propose a structure-based model in which resting→active starts with spontaneous rearrangements of the M2–M3 linkers and TBSs. These conformational changes stabilize a twisted extracellular domain to promote transmembrane helix tilting, gate dilation, and the formation of a “bubble” that collapses to initiate ion conduction. The energy landscapes suggest that twisting is the most energetically unfavorable step in the resting→active conformational change and that the rate-limiting step in the reverse process is bubble formation.

Motion of proximal histidine and structural allosteric transition in soluble guanylate cyclase

We investigated the changes of heme coordination in purified soluble guanylate cyclase (sGC) by time-resolved spectroscopy in a time range encompassing 11 orders of magnitude (from 1 ps to 0.2 s). After dissociation, NO either recombines geminately to the 4-coordinate (4c) heme (τG1 = 7.5 ps; 97 ± 1% of the population) or exits the heme pocket (3 ± 1%). The proximal His rebinds to the 4c heme with a 70-ps time constant. Then, NO is distributed in two approximately equal populations (1.5%). One geminately rebinds to the 5c heme (τG2 = 6.5 ns), whereas the other diffuses out to the solution, from where it rebinds bimolecularly (τ = 50 μs with [NO] = 200 μM) forming a 6c heme with a diffusion-limited rate constant of 2 × 10(8) M(-1)⋅s(-1). In both cases, the rebinding of NO induces the cleavage of the Fe-His bond that can be observed as an individual reaction step. Saliently, the time constant of bond cleavage differs depending on whether NO binds geminately or from solution (τ5C1 = 0.66 μs and τ5C2 = 10 ms, respectively). Because the same event occurs with rates separated by four orders of magnitude, this measurement implies that sGC is in different structural states in both cases, having different strain exerted on the Fe-His bond. We show here that this structural allosteric transition takes place in the range 1-50 μs. In this context, the detection of NO binding to the proximal side of sGC heme is discussed.

Experimental basis for a new allosteric model for multisubunit proteins

Monod, Wyman, and Changeux (MWC) explained allostery in multisubunit proteins with a widely applied theoretical model in which binding of small molecules, so-called allosteric effectors, affects reactivity by altering the equilibrium between more reactive (R) and less reactive (T) quaternary structures. In their model, each quaternary structure has a single reactivity. Here, we use silica gels to trap protein conformations and a new kind of laser photolysis experiment to show that hemoglobin, the paradigm of allostery, exhibits two ligand binding phases with the same fast and slow rates in both R and T quaternary structures. Allosteric effectors change the fraction of each phase but not the rates. These surprising results are readily explained by the simplest possible extension of the MWC model to include a preequilibrium between two tertiary conformations that have the same functional properties within each quaternary structure. They also have important implications for the long-standing question of a structural explanation for the difference in hemoglobin oxygen affinity of the two quaternary structures.

Fast protein folding kinetics

Fast-folding proteins have been a major focus of computational and experimental study because they are accessible to both techniques: they are small and fast enough to be reasonably simulated with current computational power, but have dynamics slow enough to be observed with specially developed experimental techniques. This coupled study of fast-folding proteins has provided insight into the mechanisms, which allow some proteins to find their native conformation well <1 ms and has uncovered examples of theoretically predicted phenomena such as downhill folding. The study of fast folders also informs our understanding of even 'slow' folding processes: fast folders are small; relatively simple protein domains and the principles that govern their folding also govern the folding of more complex systems. This review summarizes the major theoretical and experimental techniques used to study fast-folding proteins and provides an overview of the major findings of fast-folding research. Finally, we examine the themes that have emerged from studying fast folders and briefly summarize their application to protein folding in general, as well as some work that is left to do.

Ligand uptake modulation by internal water molecules and hydrophobic cavities in hemoglobins

Internal water molecules play an active role in ligand uptake regulation, since displacement of retained water molecules from protein surfaces or cavities by incoming ligands can promote favorable or disfavorable effects over the global binding process. Detection of these water molecules by X-ray crystallography is difficult given their positional disorder and low occupancy. In this work, we employ a combination of molecular dynamics simulations and ligand rebinding over a broad time range to shed light into the role of water molecules in ligand migration and binding. Computational studies on the unliganded structure of the thermostable truncated hemoglobin from Thermobifida fusca (Tf-trHbO) show that a water molecule is in the vicinity of the iron heme, stabilized by WG8 with the assistance of YCD1, exerting a steric hindrance for binding of an exogenous ligand. Mutation of WG8 to F results in a significantly lower stabilization of this water molecule and in subtle dynamical structural changes that favor ligand binding, as observed experimentally. Water is absent from the fully hydrophobic distal cavity of the triple mutant YB10F-YCD1F-WG8F (3F), due to the lack of residues capable of stabilizing it nearby the heme. In agreement with these effects on the barriers for ligand rebinding, over 97% of the photodissociated ligands are rebound within a few nanoseconds in the 3F mutant case. Our results demonstrate the specific involvement of water molecules in shaping the energetic barriers for ligand migration and binding.

Probing anisotropic structure changes in proteins with picosecond time-resolved small-angle X-ray scattering

We have exploited the principle of photoselection and the method of time-resolved small-angle X-ray scattering (SAXS) to investigate protein size and shape changes following photoactivation of photoactive yellow protein (PYP) in solution with ∼150 ps time resolution. This study partially overcomes the orientational average intrinsic to solution scattering methods and provides structural information at a higher level of detail. Photoactivation of the p-coumaric acid (pCA) chromophore in PYP produces a highly contorted, short-lived, red-shifted intermediate (pR0), and triggers prompt, protein compaction of approximately 0.3% along the direction defined by the electronic transition dipole moment of the chromophore. Contraction along this dimension is accompanied by expansion along the orthogonal directions, with the net protein volume change being approximately -0.25%. More than half the strain arising from formation of pR0 is relieved by the pR0 to pR1 structure transition (1.8 ± 0.2 ns), with the persistent strain presumably contributing to the driving force needed to generate the spectroscopically blue-shifted pB signaling state. The results reported here are consistent with the near-atomic resolution structural dynamics reported in a recent time-resolved Laue crystallography study of PYP crystals and suggest that the early time structural dynamics in the crystalline state carry over to proteins in solution.

Ligand migration through hemeprotein cavities: insights from laser flash photolysis and molecular dynamics simulations

The presence of cavities and tunnels in the interior of proteins, in conjunction with the structural plasticity arising from the coupling to the thermal fluctuations of the protein scaffold, has profound consequences on the pathways followed by ligands moving through the protein matrix. In this perspective we discuss how quantitative analysis of experimental rebinding kinetics from laser flash photolysis, trapping of unstable conformational states by embedding proteins within the nanopores of silica gels, and molecular simulations can synergistically converge to gain insight into the migration mechanism of ligands. We show how the evaluation of the free energy landscape for ligand diffusion based on the outcome of computational techniques can assist the definition of sound reaction schemes, leading to a comprehensive understanding of the broad range of chemical events and time scales that encompass the transport of small ligands in hemeproteins.

ATP regulation of the ligand-binding properties in temperate and cold-adapted haemoglobins. X-ray structure and ligand-binding kinetics in the sub-Antarctic fish Eleginops maclovinus

The major haemoglobin of the sub-Antarctic fish Eleginops maclovinus was structurally and functionally characterised with the aim to compare molecular environmental adaptations in the O(2)-transport system of sub-Antarctic fishes of the suborder Notothenioidei with those of their high-latitude relatives. Ligand-binding kinetics of the major haemoglobin of E. maclovinus indicated strong stabilisation of the liganded quaternary T state, enhanced in the presence of the physiological allosteric effector ATP, compared to that of high-Antarctic Trematomus bernacchii. The activation enthalpy for O(2) dissociation was dramatically lower than that in T. bernacchii haemoglobin, suggesting remarkable differences in temperature sensitivity and structural changes associated with O(2) release and exit from the protein. The haemoglobin functional properties, together with the X-ray structure of the CO form at 1.49 Å resolution, the first of a temperate notothenioid, strongly support the hypothesis that in E. maclovinus, whose life-style varies according to changes in habitat, the mechanisms that regulate O(2) affinity and the ATP-induced Root effect differ from those of high-Antarctic Notothenioids.

The Monod-Wyman-Changeux allosteric model accounts for the quaternary transition dynamics in wild type and a recombinant mutant human hemoglobin

The acknowledged success of the Monod-Wyman-Changeux (MWC) allosteric model stems from its efficacy in accounting for the functional behavior of many complex proteins starting with hemoglobin (the paradigmatic case) and extending to channels and receptors. The kinetic aspects of the allosteric model, however, have been often neglected, with the exception of hemoglobin and a few other proteins where conformational relaxations can be triggered by a short and intense laser pulse, and monitored by time-resolved optical spectroscopy. Only recently the application of time-resolved wide-angle X-ray scattering (TR-WAXS), a direct structurally sensitive technique, unveiled the time scale of hemoglobin quaternary structural transition. In order to test the generality of the MWC kinetic model, we carried out a TR-WAXS investigation in parallel on adult human hemoglobin and on a recombinant protein (HbYQ) carrying two mutations at the active site [Leu(B10)Tyr and His(E7)Gln]. HbYQ seemed an ideal test because, although exhibiting allosteric properties, its kinetic and structural properties are different from adult human hemoglobin. The structural dynamics of HbYQ unveiled by TR-WAXS can be quantitatively accounted for by the MWC kinetic model. Interestingly, the main structural change associated with the R-T allosteric transition (i.e., the relative rotation and translation of the dimers) is approximately 10-fold slower in HbYQ, and the drop in the allosteric transition rate with ligand saturation is steeper. Our results extend the general validity of the MWC kinetic model and reveal peculiar thermodynamic properties of HbYQ. A possible structural interpretation of the characteristic kinetic behavior of HbYQ is also discussed.

Following ligand migration pathways from picoseconds to milliseconds in type II truncated hemoglobin from Thermobifida fusca

CO recombination kinetics has been investigated in the type II truncated hemoglobin from Thermobifida fusca (Tf-trHb) over more than 10 time decades (from 1 ps to ∼100 ms) by combining femtosecond transient absorption, nanosecond laser flash photolysis and optoacoustic spectroscopy. Photolysis is followed by a rapid geminate recombination with a time constant of ∼2 ns representing almost 60% of the overall reaction. An additional, small amplitude geminate recombination was identified at ∼100 ns. Finally, CO pressure dependent measurements brought out the presence of two transient species in the second order rebinding phase, with time constants ranging from ∼3 to ∼100 ms. The available experimental evidence suggests that the two transients are due to the presence of two conformations which do not interconvert within the time frame of the experiment. Computational studies revealed that the plasticity of protein structure is able to define a branched pathway connecting the ligand binding site and the solvent. This allowed to build a kinetic model capable of describing the complete time course of the CO rebinding kinetics to Tf-trHb.

Heme reactivity is uncoupled from quaternary structure in gel-encapsulated hemoglobin: a resonance Raman spectroscopic study

Encapsulation of hemoglobin (Hb) in silica gel preserves structure and function but greatly slows protein motion, thereby providing access to intermediates along the allosteric pathway that are inaccessible in solution. Resonance Raman (RR) spectroscopy with visible and ultraviolet laser excitation provides probes of heme reactivity and of key tertiary and quaternary contacts. These probes were monitored in gels after deoxygenation of oxyHb and after CO binding to deoxyHb, which initiate conformational change in the R-T and T-R directions, respectively. The spectra establish that quaternary structure change in the gel takes a week or more but that the evolution of heme reactivity, as monitored by the Fe-histidine stretching vibration, ν(FeHis), is completed within two days, and is therefore uncoupled from the quaternary structure. Within each quaternary structure, the evolving ν(FeHis) frequencies span the full range of values between those previously associated with the high- and low-affinity end states, R and T. This result supports the tertiary two-state (TTS) model, in which the Hb subunits can adopt high- and low-affinity tertiary structures, r and t, within each quaternary state. The spectra also reveal different tertiary pathways, involving the breaking and reformation of E and F interhelical contacts in the R-T direction but not the T-R direction. In the latter, tertiary motions are restricted by the T quaternary contacts.

Unsuspected pathway of the allosteric transition in hemoglobin

Large conformational transitions play an essential role in the function of many proteins, but experiments do not provide the atomic details of the path followed in going from one end structure to the other. For the hemoglobin tetramer, the transition path between the unliganded (T) and tetraoxygenated (R) structures is not known, which limits our understanding of the cooperative mechanism in this classic allosteric system, where both tertiary and quaternary changes are involved. The conjugate peak refinement algorithm is used to compute an unbiased minimum energy path at atomic detail between the two end states. Although the results confirm some of the proposals of Perutz [Perutz MF (1970) Stereochemistry of cooperative effects in haemoglobin. Nature 228:726-734], the subunit motions do not follow the textbook description of a simple rotation of one αβ-dimer relative to the other. Instead, the path consists of two sequential quaternary rotations, each involving different subdomains and axes. The quaternary transitions are preceded and followed by phases of tertiary structural changes. The results explain the recent photodissociation measurements, which suggest that the quaternary transition has a fast (2 μs) as well as a slow (20 μs) component and provide a testable model for single molecule FRET experiments.

Unveiling the timescale of the R-T transition in human hemoglobin

Time-resolved wide-angle X-ray scattering, a recently developed technique allowing to probe global structural changes of proteins in solution, was used to investigate the kinetics of R-T quaternary transition in human hemoglobin and to systematically compare it to that obtained with time-resolved optical spectroscopy under nearly identical experimental conditions. Our data reveal that the main structural rearrangement associated with the R-T transition takes place approximately 2 mus after the photolysis of hemoglobin at room temperature and neutral pH. This finding suggests that the 20-mus step observed with time-resolved optical spectroscopy corresponds to a small and localized structural change.

Myoglobin strikes back

Over the last half century, myoglobin (Mb) has been an excellent model system to test a number of concepts, theories, and new experimental methods that proved valuable to investigate protein structure, function, evolution, and dynamics. Mb's function, most often considered just an oxygen repository, has considerably diversified over the last 15 years, especially because it was shown to have a role in the biochemistry of quenching and synthesizing nitric oxide in the red muscle, thereby protecting the cell. To tackle protein's structural dynamics by innovative biophysical methods, Mb has been the best prototype; laser flash technology made it possible to obtain molecular movies by time-resolved Laue crystallography (with ps resolution). This approach unveiled the complexity of the energy landscape and the structural basis of the stretched interconversion between conformational substates of a protein.

Characterization of ligand migration mechanisms inside hemoglobins from the analysis of geminate rebinding kinetics

The presence of internal hydrophobic cavities and packing defects has been demonstrated for several small globular proteins, including hemoglobins. The reduced thermodynamic stability appears to be compensated for by the capability of controlling ligand diffusion through the protein matrix to the active site, possibly by stocking more than one reactant molecule in selected sites. Photolysis of carbon monoxide complexes of hemoglobins encapsulated in silica gels leads to multiphasic geminate rebinding kinetics at room temperature, reflecting rebinding also from different temporary docking sites inside the protein matrix. A careful analysis of the ligand rebinding kinetics allows the determination of the microscopic rates for the underlying reactions, including those governing the migration to and from the docking sites. This chapter describes the experimental approach used to characterize the ligand rebinding kinetics for heme proteins in silica gels after nanosecond laser flash photolysis and the computational methods necessary to retrieve the kinetic parameters.

Allostery and cooperativity revisited

Although phenomenlogical models that account for cooperativity in allosteric systems date back to the early and mid-60's (e.g., the KNF and MWC models), there is resurgent interest in the topic due to the recent experimental and computational studies that attempted to reveal, at an atomistic level, how allostery actually works. In this review, using systems for which atomistic simulations have been carried out in our groups as examples, we describe the current understanding of allostery, how the mechanisms go beyond the classical MWC/Pauling-KNF descriptions, and point out that the "new view" of allostery, emphasizing "population shifts," is, in fact, an "old view." The presentation offers not only an up-to-date description of allostery from a theoretical/computational perspective, but also helps to resolve several outstanding issues concerning allostery.

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.

Time-resolved methods in Biophysics. 2. Monitoring haem proteins at work with nanosecond laser flash photolysis

Haem proteins have long been the most studied proteins in biophysics, and have become paradigms for the characterization of fundamental biomolecular processes as ligand binding and regulatory conformational transitions. The presence of the haem prosthetic group, the absorbance spectrum of which has a ligation sensitive region conveniently located in the UV-visible range, has offered a powerful and sensitive tool for the investigation of molecular functions. The central Fe atom is capable of reversibly binding diatomic ligands, including O(2), CO, and NO. The Fe-ligand bond is photolabile, and a reactive unligated state can be transiently generated with a pulsed laser. The photodissociated ligands quickly rebind to the haem and the process can be monitored by transient absorbance methods. The ligand rebinding kinetics reflects protein dynamics and ligand migration within the protein inner cavities. The characterization of these processes was done in the past mainly by low temperature experiments. The use of silica gels to trap proteins allows the characterization of internal ligand dynamics at room temperature. In order to show the potential of the laser flash photolysis techniques, combined with modern numerical analysis methods, we report experiments conducted on two non-symbiotic haemoglobins from Arabidopsis thaliana. The comparison between time courses recorded on haemoglobins in solution and encapsulated in silica gels allows for the highlighting of different interplays between protein dynamics and ligand migration.

Circular dichroism spectroscopy of tertiary and quaternary conformations of human hemoglobin entrapped in wet silica gels

The relative contributions to changes in visible and near UV circular dichroism spectra of hemoglobin of heme ligation and tertiary and quaternary conformational transitions were separated by exploiting the slowing down of structural relaxations for proteins encapsulated in wet, nanoporous silica gels. Spectral signatures, previously assumed to be characteristic of T and R quaternary states, were demonstrated to be specific to different tertiary conformations. The results support the view that ligation and allosteric effectors can modulate the structural and functional properties of hemoglobin by regulating the equilibrium between the same tertiary species within both quaternary states.

Dynamics of Allostery in Hemoglobin: Roles of the Penultimate Tyrosine H bonds

The tyrosine residues adjacent to the C termini of the hemoglobin (Hb) subunits, alphaY140 and betaY145, are expected to play important structural roles, because the C termini are the loci of T-state quaternary salt-bridges, and because the tyrosine side-chains bridge the H and F helices via H bonds to the alphaV93 and betaV98 carbonyl groups. These roles have been investigated via measurements of oxygen binding, (1)H NMR spectra, resonance Raman (RR) spectra, and time-resolved resonance Raman (TR(3)) spectra on site mutants in which the Hcdots, three dots, centeredF H bonds are eliminated by replacing the tyrosine residues with phenylalanine. The TR(3) spectra confirm the hypothesis, based on TR(3) studies of wild-type Hb, that the Hcdots, three dots, centeredF H bonds break and then re-form during the sub-microsecond phase of the R-T quaternary transition. The TR(3) spectra support the inference from other mutational studies that the alphabeta dimers act as single dynamic units in this early phase, motions of the E and F helices being coupled tightly across the dimer interface. Formation of T quaternary contacts occurs at about the same rate in the mutants as in HbA. However, these contacts are weakened substantially by the Y/F substitutions. Equilibrium perturbations are apparent also, especially for the alpha-subunits, in which relaxation of the Fe-His bond, strengthening of the Acdots, three dots, centeredE interhelical H bond, and weakening of the "switch" quaternary contact in deoxyHb are all apparent. Structural effects are less marked for the beta-chain Y/F replacement, but the Bohr effect is reduced by 25%, indicating that the salt-bridge and H bond interactions of the adjacent C terminus are loosened. The alpha-chain replacement reduces the Bohr effect much more, consistent with the global perturbations detected by the structure probes.

Spectroscopic markers of the T<-->R quaternary transition in human hemoglobin

In this work, we use a sol-gel protocol to trap and compare the R and T quaternary states of both the deoxygenated (deoxyHb) and carbonmonoxide (HbCO) derivatives of human hemoglobin. The near infrared optical absorption band III and the infrared CO stretching band are used to detect the effect of quaternary structure on the spectral properties of deoxyHb and HbCO; comparison with myoglobin allows for an assessment of tertiary and quaternary contributions to the measured band shifts. The R<-->T transition is shown to cause a blue shift of the band III by approximately 35 cm(-1) for deoxyHb and a red shift of the CO stretching band by only approximately 0.3 cm(-1) for HbCO. This clearly shows that quaternary structure changes are transmitted to the heme pocket and that effects on deoxyHb are much larger than on HbCO, at least as far as the band energies are concerned. Experiments performed in the ample temperature interval of 300-10K show that the above quaternary structure effects are "static" and do not influence the dynamic properties of the heme pocket, at least as probed by the temperature dependence of band III and of the CO stretching band. The availability of quaternary structure sensitive spectroscopic markers and the quantitative measurement of the quaternary structure contribution to band shifts will be of considerable help in the analysis of flash-photolysis experiments on hemoglobin. Moreover, it will enable one to characterize the dynamic properties of functionally relevant hemoglobin intermediates and to study the kinetics of both the T-->R and R-->T quaternary transitions through time-resolved spectroscopy.

A speed limit for conformational change of an allosteric membrane protein

Neuromuscular acetylcholine receptors are synaptic ion channels that open and close with rate constants of approximately equal to 48,000 s(-1) and approximately equal to 1,700 s(-1), respectively (in adult mouse, at 24 degrees C, -100 mV membrane potential). Perturbations of many different sites in the protein can change these rate constants, with those in the extracellular domain mainly affecting channel-opening and many of those in the membrane and intracellular domains mainly affecting channel-closing. We used single-channel recordings to measure the total open time per activation (tau(b)) elicited by a low concentration of the natural transmitter, acetylcholine. tau(b) increased in constructs with mutations that increased the gating equilibrium constant by either increasing the opening or decreasing the closing rate constant. However, tau(b) did not approach the same asymptote in fast-opening and slow-closing constructs. The maximum value for the slow closers was about twice that for the fast openers. One interpretation of this difference is that there is an upper limit to the channel-opening rate constant, which we estimate to be approximately 0.86 mus(-1). One possibility is that this limit is the rate of conformational change in the absence of an overall activation barrier and thus reflects the kinetic prefactor for the acetylcholine receptor opening isomerization.

Hemoglobin site-mutants reveal dynamical role of interhelical H-bonds in the allosteric pathway: time-resolved UV resonance Raman evidence for intra-dimer coupling

The dynamical effect of eliminating specific tertiary H-bonds in the hemoglobin (Hb) tetramer has been investigated by site-directed mutagenesis and time-resolved absorption and ultraviolet resonance Raman (UVRR) spectroscopy. The Trp alpha 14...Thr alpha 67 and Trp beta 15...Ser beta 72 H-bonds connect the A and E helices in the alpha and beta chains, and are proposed to break in the earliest protein intermediate (Rdeoxy) following photo-deligation of HbCO, along with a second pair of H-bonds involving tyrosine residues. Mutation of the acceptor residues Thr alpha 67 and Ser beta 72 to Val and Ala eliminates the A-E H-bonds, but has been shown to have no significant effect on ligand-binding affinity or cooperativity, or on spectroscopic markers of the T-state quaternary interactions. However, the mutations have profound and unexpected effects on the character of the Rdeoxy intermediate, and on the dynamics of the subsequent steps leading to the T state. Formation of the initial quaternary contact (RT intermediate) is accelerated, by an order of magnitude, but the locking-in of the T state is delayed by a factor of 2. These rate effects are essentially the same for either mutation, or for the double mutation, suggesting that the alpha beta dimer behaves as a mechanically coupled dynamical unit. Further evidence for intra-dimer coupling is provided by the Rdeoxy UVRR spectrum, in which either or both mutations eliminate the tyrosine difference intensity, although only tryptophan H-bonds are directly affected. A possible mechanism for mechanical coupling is outlined, involving transmission of forces through the alpha(1)beta(1) (and alpha(2)beta(2)) interface. The present observations establish that quaternary motions can occur on the approximately 100 ns time-scale. They show also that a full complement of interhelical H-bonds actually slows the initial quaternary motion in Hb, but accelerates the locking in of the T-contacts.

Time-resolved Absorption and UV Resonance Raman Spectra Reveal Stepwise Formation of T Quaternary Contacts in the Allosteric Pathway of Hemoglobin

Hemoglobin undergoes a series of molecular changes on the nanosecond and microsecond time-scale following photodissociation of CO ligands. We have monitored these processes with a combination of transient absorption and resonance Raman (RR) spectroscopy. The latter have been acquired at higher data rates than previously available, thanks to kilohertz Ti:sapphire laser technology, with frequency-quadrupling into the ultraviolet. As a result of improved resolution of the UVRR time-course, a new intermediate has been identified in the pathway from the R (HbCO) to the T (deoxyHb) state. This intermediate is not detected via absorption transients, since the change in heme absorption is insignificant, but its lifetime agrees with a reported magnetic circular dichroism transient, which has been attributed to a quaternary tryptophan interaction. The new UVRR data allow elaboration of the allosteric pathway by establishing that the T-state quaternary contacts are formed in two well-separated steps, with time constants of 2.9 micros and 21 micros, instead of a single 20 micros process. The first step involves the "hinge" region contacts, as monitored by the Trp beta 37...Asp alpha 94 H-bond, while the second involves the "switch" region, as monitored by the Tyr alpha 42...Asp beta 99 H-bond. A working model for the allosteric pathway is presented.

Laser-induced time-resolved photoacoustic calorimetry study on photo-dissociation of human and bovine oxyhemoglobin

The dynamics of the enthalpy and volume changes related to the photo-dissociation of oxygen from human and bovine oxyhemoglobin are investigated by nanosecond time-resolved photoacoustic calorimetry (PAC). The values of enthalpy and volume change associated with the above process are deltaH = 37.8 +/- 3 kcal/mol, deltaV = 5.0 +/- 1 ml/mol for human HbO(2); and deltaH = 35.7 +/- 3.5 kcal/mol, deltaV = 4.8 +/- 1 ml/mol for bovine HbO(2), respectively. A possible explanation for the similar values between both human and bovine oxyhemoglobin is proposed. In addition, the PAC results for human HbO(2) and HbCO are compared and discussed.

A tertiary two-state allosteric model for hemoglobin

The two-state allosteric model of Monod, Wyman, and Changeux (MWC) provides an excellent description of homotropic effects in a vast array of equilibrium and kinetic measurements on cooperative ligand binding by hemoglobin. However, in contrast to experimental observations, the model does not allow for alteration of the ligand affinity of the T quaternary structure by allosteric effectors. This failure to explain heterotropic effects has been appreciated for over 30 years, and it has been generally assumed to result from tertiary conformational changes in the absence of a quaternary change. Here we explore a model that preserves the essential MWC idea that binding without a quaternary conformational change is non-cooperative, but where tertiary conformations of individual subunits play the primary role instead of the quaternary conformations. In this model, which we call the 'tertiary two-state (TTS) model', the two affinity states correspond to two tertiary conformations of individual subunits rather than the two quaternary conformations of the MWC two-state allosteric model. Ligation and the R quaternary structure bias the subunit population toward the high affinity tertiary conformation, while deligation and the T quaternary structure favor the low affinity tertiary conformation. We show that the model is successful in quantitatively explaining a demanding set of kinetic data from nanosecond carbon monoxide photodissociation experiments at times longer than approximately 1 micros. Better agreement between the model and the submicrosecond kinetic data may result from detailed considerations of the distribution and dynamics of conformational substates of the two tertiary conformations. The model is consistent with the results of solution, gel, and single crystal oxygen binding studies, but underestimates the population of doubly-liganded molecules determined in low-temperature electrophoresis experiments.

Kinetics of nitric oxide binding to R-state hemoglobin

Despite earlier work indicating otherwise, some recent reports have suggested that nitric oxide (NO) binds to hemoglobin cooperatively. In particular, it has been suggested that, under physiological conditions, NO binds to the high-affinity R-state hemoglobin as much as 100 times faster than to the low-affinity T-state hemoglobin. This rapid NO binding could provide a means of preserving NO bioactivity. However, using a flash-flow photolysis technique, we have determined that the rate of NO binding to normal adult R-state hemoglobin is (2.1 +/- 0.1) x 10(7) (s(-1) M(-1)), which is essentially the same as that reported for T-state NO binding. (c)2002 Elsevier Science (USA).

Control of nitric oxide dynamics by guanylate cyclase in its activated state

Soluble guanylate cyclase (sGC) is the target of nitric oxide (NO) released by nitric-oxide synthase in endothelial cells, inducing an increase of cGMP synthesis in response. This heterodimeric protein possesses a regulatory subunit carrying a heme where NO binding occurs, while the second subunit harbors the catalytic site. The binding of NO and the subsequent breaking of the bond between the proximal histidine and the heme-Fe(2+) are assumed to induce conformational changes, which are the origin of the catalytic activation. At the molecular level, the activation and deactivation mechanisms are unknown, as is the dynamics of NO once in the heme pocket. Using ultrafast time-resolved absorption spectroscopy, we measured the kinetics of NO rebinding to sGC after photodissociation. The main spectral transient in the Soret band does not match the equilibrium difference spectrum of NO-liganded minus unliganded sGC, and the geminate rebinding was found to be monoexponential and ultrafast (tau = 7.5 ps), with a relative amplitude close to unity (0.97). These characteristics, so far not observed in other hemoproteins, indicate that NO encounters a high energy barrier for escaping from the heme pocket once the His-Fe(2+) bond has been cleaved; this bond does not reform before NO recombination. The deactivation of isolated sGC cannot occur by only simple diffusion of NO from the heme; therefore, several allosteric states may be inferred, including a desensitized one, to induce NO release. Thus, besides the structural change leading to activation, a consequence of the decoupling of the proximal histidine may also be to induce a change of the heme pocket distal geometry, which raises the energy barrier for NO escape, optimizing the efficiency of NO trapping. The non-single exponential character of the NO picosecond rebinding coexists only with the presence of the protein structure surrounding the heme, and the single exponential rate observed in sGC is very likely to be due to a closed conformation of the heme pocket. Our results emphasize the physiological importance of NO geminate recombination in hemoproteins like nitric-oxide synthase and sGC and show that the protein structure controls NO dynamics in a manner adapted to their function. This control of ligand dynamics provides a regulation at molecular level in the function of these enzymes.

Multiple geminate ligand recombinations in human hemoglobin

The geminate ligand recombination reactions of photolyzed carbonmonoxyhemoglobin were studied in a nanosecond double-excitation-pulse time-resolved absorption experiment. The second laser pulse, delayed by intervals as long as 400 ns after the first, provided a measure of the geminate kinetics by rephotolyzing ligands that have recombined during the delay time. The peak-to-trough magnitude of the Soret band photolysis difference spectrum measured as a function of the delay between excitation pulses showed that the room temperature kinetics of geminate recombination in adult human hemoglobin are best described by two exponential processes, with lifetimes of 36 and 162 ns. The relative amounts of bimolecular recombination to T- and R-state hemoglobins and the temperature dependence of the submicrosecond kinetics between 283 and 323 K are also consistent with biexponential kinetics for geminate recombination. These results are discussed in terms of two models: geminate recombination kinetics modulated by concurrent protein relaxation and heterogeneous kinetics arising from alpha and beta chain differences.

Nanosecond time-resolved spectroscopy of biomolecular processes

Over the past two decades, nanosecond absorption and vibrational spectroscopies have developed into powerful tools for monitoring the secondary, tertiary, and quaternary structural relaxations of biological macromolecules under near-physiological conditions of solvent and temperature. Observed through such methods, the dynamic response of a biomolecule to photoinitiated excursions from equilibrium can reveal valuable information about the structure-function relationship, information beyond that obtained from the static structures provided by X-ray crystallography, nuclear magnetic resonance spectroscopy, and other steady-state methods. Most recently, the development of ultra-sensitive polarization techniques for absorption spectroscopy has greatly enhanced the amount of time-resolved structural information that can be obtained from the broadened electronic spectra of biomolecules. This review examines nanosecond absorption, vibrational, and polarized absorption methods, and their applications to protein function and folding, emphasizing the complementary nature of information obtained from electronic and vibrational spectra measured on the nanosecond time scale.

Nanosecond time-resolved absorption studies of human oxyhemoglobin photolysis intermediates

The time-resolved spectra of photoproducts from ligand photodissociation of oxyhemoglobin are measured in the Soret spectral region for times from 10 ns to 320 microseconds after laser photolysis. Four processes are detected at a heme concentration of 80 microM: a 38-ns geminate recombination, a 137-ns tertiary relaxation, and two bimolecular processes for rebinding of molecular oxygen. The pseudo-first-order rate constants for rebinding to the alpha and beta subunits of hemoglobin are 3.2 x 10(4) s-1 (31 microseconds lifetime) and 9.4 x 10(4) s-1 (11 microseconds lifetime), respectively. The significance of kinetic measurements made at different heme concentrations is discussed in terms of the equilibrium compositions of hemoglobin tetramer and dimer mixtures. The rebinding rate constants for alpha and beta chains are observed to be about two times slower in the dimer than in the tetramer, a finding that appears to support the observation of quaternary enhancement in equilibrium ligand binding by hemoglobin tetramers.