Extended subnanosecond structural dynamics of myoglobin revealed by Laue crystallography

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.

Using photocaging for fast time-resolved structural biology studies

Careful selection of photocaging approaches is critical to achieve fast and well synchronized reaction initiation and perform successful time-resolved structural biology experiments. This review summarizes the best characterized and most relevant photocaging groups previously described in the literature. It also provides a walkthrough of the essential factors to consider in designing a suitable photocaged molecule to address specific biological questions, focusing on photocaging groups with well characterized spectroscopic properties. The relationships between decay rates (k in s-1), quantum yields (ϕ) and molar extinction coefficients (ϵmax in M-1 cm-1) are highlighted for different groups. The effects of the nature of the photocaged group on these properties is also discussed. Four main photocaging scaffolds are presented in detail, o-nitrobenzyls, p-hydroxyphenyls, coumarinyls and nitrodibenzofuranyls, along with three examples of the use of this technology. Furthermore, a subset of specialty photocages are highlighted: photoacids, molecular photoswitches and metal-containing photocages. These extend the range of photocaging approaches by, for example, controlling pH or generating conformationally locked molecules.

Kinetic mechanisms for O2 binding to myoglobins and hemoglobins

Antonini and Brunori's 1971 book "Hemoglobin and Myoglobin in Their Reactions with Ligands" was a truly remarkable publication that summarized almost 100 years of research on O₂ binding to these globins. Over the ensuing 50 years, ultra-fast laser photolysis techniques, high-resolution and time resolved X-ray crystallography, molecular dynamics simulations, and libraries of recombinant myoglobin (Mb) and hemoglobin (Hb) variants have provided structural interpretations of O₂ binding to these proteins. The resultant mechanisms provide quantitative descriptions of the stereochemical factors that govern overall affinity, including proximal and distal steric restrictions that affect iron reactivity and favorable positive electrostatic interactions that preferentially stabilize bound O₂. The pathway for O₂ uptake and release by Mb and subunits of Hb has been mapped by screening libraries of site-directed mutants in laser photolysis experiments. O₂ enters mammalian Mb and the α and β subunits of human HbA through a channel created by upward and outward rotation of the distal His at the E7 helical position, is non-covalently captured in the interior of the distal cavity, and then internally forms a bond with the heme Fe(II) atom. O₂ dissociation is governed by disruption of hydrogen bonding interactions with His (E7), breakage of the Fe(II)-O₂ bond, and then competition between rebinding and escape through the E7-gate. The structural features that govern the rates of both the individual steps and overall reactions have been determined and provide the framework for: (1) defining the physiological functions of specific globins and their evolution; (2) understanding the clinical features of hemoglobinopathies; and (3) designing safer and more efficient acellular hemoglobin-based oxygen carriers (HBOCs) for transfusion therapy, organ preservation, and other commercially relevant O₂ transport and storage processes.

Femtosecond X-ray spectroscopy of haem proteins

We discuss our recently reported femtosecond (fs) X-ray emission spectroscopy results on the ligand dissociation and recombination in nitrosylmyoglobin (MbNO) in the context of previous studies on ferrous haem proteins. We also present a preliminary account of femtosecond X-ray absorption studies on MbNO, pointing to the presence of more than one species formed upon photolysis.

A guide to time-resolved structural analysis of light-activated proteins

Dynamical changes in protein structures are essential for protein function and occur over femtoseconds to seconds timescales. X-ray free electron lasers have facilitated investigations of structural dynamics in proteins with unprecedented temporal and spatial resolution. Light-activated proteins are attractive targets for time-resolved structural studies, as the reaction chemistry and associated protein structural changes can be triggered by short laser pulses. Proteins with different light-absorbing centres have evolved to detect light and harness photon energy to bring about downstream chemical and biological output responses. Following light absorption, rapid chemical/small-scale structural changes are typically localised around the chromophore. These localised changes are followed by larger structural changes propagated throughout the photoreceptor/photocatalyst that enables the desired chemical and/or biological output response. Time-resolved serial femtosecond crystallography (SFX) and solution scattering techniques enable direct visualisation of early chemical change in light-activated proteins on timescales previously inaccessible, whereas scattering gives access to slower timescales associated with more global structural change. Here, we review how advances in time-resolved SFX and solution scattering techniques have uncovered mechanisms of photochemistry and its coupling to output responses. We also provide a prospective on how these time-resolved structural approaches might impact on other photoreceptors/photoenzymes that have not yet been studied by these methods.

Ultrafast Structural Changes Decomposed from Serial Crystallographic Data

Direct visualization of electronic and molecular events during biochemical reactions is essential to mechanistic insights. This Letter presents an in-depth analysis of the serial crystallographic data sets collected by Barends and Schlichting et al. ( Science 2015 , 350 , 445 ) that probe the ligand photodissociation in carbonmonoxy myoglobin. This analysis reveals electron density changes caused by the formation of high-spin 3d atomic orbitals of the heme iron upon photolysis and their dynamic behaviors within the first few picoseconds. The heme iron is found popping out of and recoiling back into the heme plane in succession. These findings provide long-awaited visual validations for previous works using ultrafast spectroscopy and molecular dynamics simulations. Electron density variations are also found largely in the solvent during the first period of a low-frequency oscillation. This work demonstrates the importance of the analytical methods in detecting and isolating weak, transient signals of electronic changes arising from chemical reactions.

Ligand pathways in neuroglobin revealed by low-temperature photodissociation and docking experiments

A combined biophysical approach was applied to map gas-docking sites within murine neuroglobin (Ngb), revealing snapshots of events that might govern activity and dynamics in this unique hexacoordinate globin, which is most likely to be involved in gas-sensing in the central nervous system and for which a precise mechanism of action remains to be elucidated. The application of UV-visible microspectroscopy in crystallo, solution X-ray absorption near-edge spectroscopy and X-ray diffraction experiments at 15-40 K provided the structural characterization of an Ngb photolytic intermediate by cryo-trapping and allowed direct observation of the relocation of carbon monoxide within the distal heme pocket after photodissociation. Moreover, X-ray diffraction at 100 K under a high pressure of dioxygen, a physiological ligand of Ngb, unravelled the existence of a storage site for O2 in Ngb which coincides with Xe-III, a previously described docking site for xenon or krypton. Notably, no other secondary sites were observed under our experimental conditions.

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.

A short history of structure based research on the photocycle of photoactive yellow protein

The goals of time-resolved macromolecular crystallography are to extract the molecular structures of the reaction intermediates and the reaction dynamics from time-resolved X-ray data alone. To develop the techniques of time-resolved crystallography, biomolecules with special properties are required. The Photoactive Yellow Protein is the most sparkling of these.

Ultrafast X-Ray Crystallography and Liquidography

Time-resolved X-ray diffraction provides direct information on three-dimensional structures of reacting molecules and thus can be used to elucidate structural dynamics of chemical and biological reactions. In this review, we discuss time-resolved X-ray diffraction on small molecules and proteins with particular emphasis on its application to crystalline (crystallography) and liquid-solution (liquidography) samples. Time-resolved X-ray diffraction has been used to study picosecond and slower dynamics at synchrotrons and can now access even femtosecond dynamics with the recent arrival of X-ray free-electron lasers.

Picosecond Photobiology: Watching a Signaling Protein Function in Real Time via Time-Resolved Small- and Wide-Angle X-ray Scattering

The capacity to respond to environmental changes is crucial to an organism's survival. Halorhodospira halophila is a photosynthetic bacterium that swims away from blue light, presumably in an effort to evade photons energetic enough to be genetically harmful. The protein responsible for this response is believed to be photoactive yellow protein (PYP), whose chromophore photoisomerizes from trans to cis in the presence of blue light. We investigated the complete PYP photocycle by acquiring time-resolved small and wide-angle X-ray scattering patterns (SAXS/WAXS) over 10 decades of time spanning from 100 ps to 1 s. Using a sequential model, global analysis of the time-dependent scattering differences recovered four intermediates (pR0/pR1, pR2, pB0, pB1), the first three of which can be assigned to prior time-resolved crystal structures. The 1.8 ms pB0 to pB1 transition produces the PYP signaling state, whose radius of gyration (Rg = 16.6 Å) is significantly larger than that for the ground state (Rg = 14.7 Å) and is therefore inaccessible to time-resolved protein crystallography. The shape of the signaling state, reconstructed using GASBOR, is highly anisotropic and entails significant elongation of the long axis of the protein. This structural change is consistent with unfolding of the 25 residue N-terminal domain, which exposes the β-scaffold of this sensory protein to a potential binding partner. This mechanistically detailed description of the complete PYP photocycle, made possible by time-resolved crystal and solution studies, provides a framework for understanding signal transduction in proteins and for assessing and validating theoretical/computational approaches in protein biophysics.

Mapping the Accessible Conformational Landscape of an Insect Carboxylesterase Using Conformational Ensemble Analysis and Kinetic Crystallography

The proper function of enzymes often depends upon their efficient interconversion between particular conformational sub-states on a free-energy landscape. Experimentally characterizing these sub-states is challenging, which has limited our understanding of the role of protein dynamics in many enzymes. Here, we have used a combination of kinetic crystallography and detailed analysis of crystallographic protein ensembles to map the accessible conformational landscape of an insect carboxylesterase (LcαE7) as it traverses all steps in its catalytic cycle. LcαE7 is of special interest because of its evolving role in organophosphate insecticide resistance. Our results reveal that a dynamically coupled network of residues extends from the substrate-binding site to a surface loop. Interestingly, the coupling of this network that is apparent in the apoenzyme appears to be reduced in the phosphorylated enzyme intermediate. Altogether, the results of this work highlight the importance of protein dynamics to enzyme function and the evolution of new activity.

Benchmarking Rapid TLES Simulations of Gas Diffusion in Proteins: Mapping O2 Migration and Escape in Myoglobin as a Case Study

Standard molecular dynamics (MD) simulations of gas diffusion consume considerable computational time and resources even for small proteins. To combat this, temperature-controlled locally enhanced sampling (TLES) examines multiple diffusion trajectories per simulation by accommodating multiple noninteracting copies of a gas molecule that diffuse independently, while the protein and water molecules experience an average interaction from all copies. Furthermore, gas migration within a protein matrix can be accelerated without altering protein dynamics by increasing the effective temperature of the TLES copies. These features of TLES enable rapid simulations of gas diffusion within a protein matrix at significantly reduced (∼98%) computational cost. However, the results of TLES and standard MD simulations have not been systematically compared, which limits the adoption of the TLES approach. We address this drawback here by benchmarking TLES against standard MD in the simulation of O2 diffusion in myoglobin (Mb) as a case study since this model system has been extensively characterized. We find that 2 ns TLES and 108 ns standard simulations map the same network of diffusion tunnels in Mb and uncover the same docking sites, barriers, and escape portals. We further discuss the influence of simulation time as well as the number of independent simulations on the O2 population density within the diffusion tunnels and on the sampling of Mb's conformational space as revealed by principal component analysis. Overall, our comprehensive benchmarking reveals that TLES is an appropriate and robust tool for the rapid mapping of gas diffusion in proteins when the kinetic data provided by standard MD are not required. Furthermore, TLES provides explicit ligand diffusion pathways, unlike most rapid methods.

Understanding the kinetics of ligand binding to globins with molecular dynamics simulations: the necessity of multiple state models

Residue conformational changes and internal cavity migration processes play a key role in regulating the kinetics of ligand migration and binding events in globins. Molecular dynamics simulations have demonstrated their value in the study of these processes in different haemoglobins, but derivation of kinetic data demands the use of more complex techniques like enhanced sampling molecular dynamics methods. This review discusses the different methodologies that are currently applied to study the ligand migration process in globins and highlight those specially developed to derive kinetic data.

Direct observation of ultrafast collective motions in CO myoglobin upon ligand dissociation

The hemoprotein myoglobin is a model system for the study of protein dynamics. We used time-resolved serial femtosecond crystallography at an x-ray free-electron laser to resolve the ultrafast structural changes in the carbonmonoxy myoglobin complex upon photolysis of the Fe-CO bond. Structural changes appear throughout the protein within 500 femtoseconds, with the C, F, and H helices moving away from the heme cofactor and the E and A helices moving toward it. These collective movements are predicted by hybrid quantum mechanics/molecular mechanics simulations. Together with the observed oscillations of residues contacting the heme, our calculations support the prediction that an immediate collective response of the protein occurs upon ligand dissociation, as a result of heme vibrational modes coupling to global modes of the protein.

Identification of Mutational Hot Spots for Substrate Diffusion: Application to Myoglobin

The pathways by which small molecules (substrates or inhibitors) access active sites are a key aspect of the function of enzymes and other proteins. A key problem in designing or altering such proteins is to identify sites for mutation that will have the desired effect on the substrate transport properties. While specific access channels have been invoked in the past, molecular simulations suggest that multiple routes are possible, complicating the analysis. This complexity, however, can be captured by a Markov State Model (MSM) of the ligand diffusion process. We have developed a sensitivity analysis of the resulting rate matrix, which identifies the locations where mutations should have the largest effect on the diffusive on rate. We apply this method to myoglobin, which is the best characterized example both from experiment and simulation. We validate the approach by translating the sensitivity parameter obtained from this method into the CO binding rates in myoglobin upon mutation, resulting in a semi-quantitative correlation with experiments. The model is further validated against an explicit simulation for one of the experimental mutants.

Full kinetics of CO entry, internal diffusion, and exit in myoglobin from transition-path theory simulations

We use Markovian milestoning molecular dynamics (MD) simulations on a tessellation of the collective variable space for CO localization in myoglobin to estimate the kinetics of entry, exit, and internal site-hopping. The tessellation is determined by analysis of the free-energy surface in that space using transition-path theory (TPT), which provides criteria for defining optimal milestones, allowing short, independent, cell-constrained MD simulations to provide properly weighted kinetic data. We coarse grain the resulting kinetic model at two levels: first, using crystallographically relevant internal cavities and their predicted interconnections and solvent portals; and second, as a three-state side-path scheme inspired by similar models developed from geminate recombination experiments. We show semiquantitative agreement with experiment on entry and exit rates and in the identification of the so-called "histidine gate" at position 64 through which ≈90% of flux between solvent and the distal pocket passes. We also show with six-dimensional calculations that the minimum free-energy pathway of escape through the histidine gate is a "knock-on" mechanism in which motion of the ligand and the gate are sequential and interdependent. In total, these results suggest that such TPT simulations are indeed a promising approach to overcome the practical time-scale limitations of MD to allow reliable estimation of transition mechanisms and rates among metastable states.

Quantitative delineation of how breathing motions open ligand migration channels in myoglobin and its mutants

Ligand migration and binding are central to the biological functions of many proteins such as myoglobin (Mb) and it is widely thought that protein breathing motions open up ligand channels dynamically. However, how a protein exerts its control over the opening and closing of these channels through its intrinsic dynamics is not fully understood. Specifically, a quantitative delineation of the breathing motions that are needed to open ligand channels is lacking. In this work, we present and apply a novel normal mode-based method to quantitatively delineate what and how breathing motions open ligand migration channels in Mb and its mutants. The motivation behind this work springs from the observation that normal mode motions are closely linked to the breathing motions that are thought to open ligand migration channels. In addition, the method provides a direct and detailed depiction of the motions of each and every residue that lines a channel and can identify key residues that play a dominating role in regulating the channel. The all-atom model and the full force-field employed in the method provide a realistic energetics on the work cost required to open a channel, and as a result, the method can be used to efficiently study the effects of mutations on ligand migration channels and on ligand entry rates. Our results on Mb and its mutants are in excellent agreement with MD simulation results and experimentally determined ligand entry rates.

Time-resolved serial crystallography captures high-resolution intermediates of photoactive yellow protein

Serial femtosecond crystallography using ultrashort pulses from x-ray free electron lasers (XFELs) enables studies of the light-triggered dynamics of biomolecules. We used microcrystals of photoactive yellow protein (a bacterial blue light photoreceptor) as a model system and obtained high-resolution, time-resolved difference electron density maps of excellent quality with strong features; these allowed the determination of structures of reaction intermediates to a resolution of 1.6 angstroms. Our results open the way to the study of reversible and nonreversible biological reactions on time scales as short as femtoseconds under conditions that maximize the extent of reaction initiation throughout the crystal.

Probing the electronic and geometric structure of ferric and ferrous myoglobins in physiological solutions by Fe K-edge absorption spectroscopy

We present an iron K-edge X-ray absorption study of carboxymyoglobin (MbCO), nitrosylmyoglobin (MbNO), oxymyoglobin (MbO2), cyanomyoglobin (MbCN), aquomet myoglobin (metMb) and unligated myoglobin (deoxyMb) in physiological media. The analysis of the XANES region is performed using the full-multiple scattering formalism, implemented within the MXAN package. This reveals trends within the heme structure, absent from previous crystallographic and X-ray absorption analysis. In particular, the iron-nitrogen bond lengths in the porphyrin ring converge to a common value of about 2 Å, except for deoxyMb whose bigger value is due to the doming of the heme. The trends of the Fe-Nε (His93) bond length is found to be consistent with the effect of ligand binding to the iron, with the exception of MbNO, which is explained in terms of the repulsive trans effect. We derive a high resolution description of the relative geometry of the ligands with respect to the heme and quantify the magnitude of the heme doming in the deoxyMb form. Finally, time-dependent density functional theory is used to simulate the pre-edge spectra and is found to be in good agreement with the experiment. The XAS spectra typically exhibit one pre-edge feature which arises from transitions into the unoccupied dσ and dπ - πligand* orbitals. 1s → dπ transitions contribute weakly for MbO2, metMb and deoxyMb. However, despite this strong Fe d contribution these transitions are found to be dominated by the dipole (1s → 4p) moment due to the low symmetry of the heme environment.

The crystal structure of wild-type human brain neuroglobin reveals flexibility of the disulfide bond that regulates oxygen affinity

Neuroglobin plays an important function in the supply of oxygen in nervous tissues. In human neuroglobin, a cysteine at position 46 in the loop connecting the C and D helices of the globin fold is presumed to form an intramolecular disulfide bond with Cys55. Rupture of this disulfide bridge stabilizes bi-histidyl haem hexacoordination, causing an overall decrease in the affinity for oxygen. Here, the first X-ray structure of wild-type human neuroglobin is reported at 1.74 Å resolution. This structure provides a direct observation of two distinct conformations of the CD region containing the intramolecular disulfide link and highlights internal cavities that could be involved in ligand migration and/or are necessary to enable the conformational transition between the low and high oxygen-affinity states following S-S bond formation.

Femtosecond crystallography with ultrabright electrons and x-rays: capturing chemistry in action

With the recent advances in ultrabright electron and x-ray sources, it is now possible to extend crystallography to the femtosecond time domain to literally light up atomic motions involved in the primary processes governing structural transitions. This review chronicles the development of brighter and brighter electron and x-ray sources that have enabled atomic resolution to structural dynamics for increasingly complex systems. The primary focus is on achieving sufficient brightness using pump-probe protocols to resolve the far-from-equilibrium motions directing chemical processes that in general lead to irreversible changes in samples. Given the central importance of structural transitions to conceptualizing chemistry, this emerging field has the potential to significantly improve our understanding of chemistry and its connection to driving biological processes.

Enzyme transient state kinetics in crystal and solution from the perspective of a time-resolved crystallographer

With recent technological advances at synchrotrons [Graber et al., J. Synchrotron Radiat. 18, 658-670 (2011)], it is feasible to rapidly collect time-resolved crystallographic data at multiple temperature settings [Schmidt et al., Acta Crystallogr. D 69, 2534-2542 (2013)], from which barriers of activation can be extracted. With the advent of fourth generation X-ray sources, new opportunities emerge to investigate structure and dynamics of biological macromolecules in real time [M. Schmidt, Adv. Condens. Matter Phys. 2013, 1-10] in crystals and potentially from single molecules in random orientation in solution [Poon et al., Adv. Condens. Matter Phys. 2013, 750371]. Kinetic data from time-resolved experiments on short time-scales must be interpreted in terms of chemical kinetics [Steinfeld et al., Chemical Kinetics and Dynamics, 2nd ed. (Prentience Hall, 1985)] and tied to existing time-resolved experiments on longer time-scales [Schmidt et al., Acta Crystallogr. D 69, 2534-2542 (2013); Jung et al., Nat. Chem. 5, 212-220 (2013)]. With this article, we will review and outline steps that are required to routinely determine the energetics of reactions in biomolecules in crystal and solution with newest X-ray sources. In eight sections, we aim to describe concepts and experimental details that may help to inspire new approaches to collect and interpret these data.

Long-range electron transfer with myoglobin immobilized at Au/mixed-SAM junctions: mechanistic impact of the strong protein confinement

Horse muscle myoglobin (Mb) was tightly immobilized at Au-deposited ~15-Å-thick mixed-type (1:1) alkanethiol SAMs, HS-(CH₂)₁₁-COOH/HS-(CH₂)₁₁-OH, and placed in contact with buffered H₂O or D₂O solutions. Fast-scan cyclic voltammetry (CV) and a Marcus-equation-based analysis were applied to determine unimolecular standard rate constants and reorganization free energies for electron transfer (ET), under variable-temperature (15-55 °C) and -pressure (0.01-150 MPa) conditions. The CV signal was surprisingly stable and reproducible even after multiple temperature and pressure cycles. The data analysis revealed the following values: standard rate constant, 33 s⁻¹ (25 °C, 0.01 MPa, H₂O); reorganization free energy, 0.5 ± 0.1 eV (throughout); activation enthalpy, 12 ± 3 kJ mol⁻¹; activation volume, -3.1 ± 0.2 cm³ mol⁻¹; and pH-dependent solvent kinetic isotope effect (k(H)⁰/k(D)⁰), 0.7-1.4. Furthermore, the values for the rate constant and reorganization free energy are very similar to those previously found for cytochrome c electrostatically immobilized at the monocomponent Au/HS-(CH₂)₁₁-COOH junction. In vivo, Mb apparently forms a natural electrostatic complex with cytochrome b₅ (cyt-b₅) through the "dynamic" (loose) docking pattern, allowing for a slow ET that is intrinsically coupled to the water's removal from the "defective" heme iron (altogether shaping the biological repair mechanism for Mb's "met" form). In contrary, our experiments rather mimic the case of a "simple" (tight) docking of the redesigned (mutant) Mb with cyt-b₅ (Nocek et al. J. Am. Chem. Soc. 2010, 132, 6165-6175). According to our analysis, in this configuration, Mb's distal pocket (linked to the "ligand channel") seems to be arrested within the restricted configuration, allowing the rate-determining reversible ET process to be coupled only to the inner-sphere reorganization (minimal elongation/shortening of an Fe-OH₂ bond) rather than the pronounced detachment (rebinding) of water and, hence, to be much faster.

Real-time tracking of CO migration and binding in the α and β subunits of human hemoglobin via 150-ps time-resolved Laue crystallography

We have developed the method of picosecond Laue crystallography and used this capability to probe ligand dynamics in tetrameric R-state hemoglobin (Hb). Time-resolved, 2 Å-resolution electron density maps of photolyzed HbCO reveal the time-dependent population of CO in the binding (A) and primary docking (B) sites of both α and β subunits from 100 ps to 10 μs. The proximity of the B site in the β subunit is about 0.25 Å closer to its A binding site, and its k_BA rebinding rate (~300 μs^-1) is six times faster, suggesting distal control of the rebinding dynamics. Geminate rebinding in the β subunit exhibits both prompt and delayed geminate phases. We developed a microscopic model to quantitatively explain the observed kinetics, with three states for the α subunit and four states for the β subunit. This model provides a consistent framework for interpreting rebinding kinetics reported in prior studies of both HbCO and HbO₂.

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.

A computational study of water and CO migration sites and channels inside myoglobin

Pathways are computed for transport of H₂O and CO in myoglobin (Mb), using the single sweep and zero-temperature string methods in a fully atomistic, explicitly solvated model system. Our predictions of sites and barriers in the pathways for CO transport agree with previous studies. For H₂O, we predict a binding site in the distal pocket (DP), in agreement with crystallographic observations, and another one close to Leu 29 which explains the importance of this residue in controlling the pocket's hydrophobicity, as well as disordered minima in the largely apolar xenon cavities. In particular, H₂O can occupy and transition among the xenon cavities, Xe4, Xe2, and Xe3. Our results support the hypothesis that the thermodynamically most favorable entry/exit portal for H₂O is the so-called histidine gate (HG), the same as for CO. This result, along with the observation of water occupation of both DP and apolar Xe cavities, suggest that water and small gas molecules like CO compete for access to the protein interior, and therefore models of gas molecule transport within proteins should also explicitly consider water transport.

Determination of ligand pathways in globins: apolar tunnels versus polar gates

Although molecular dynamics simulations suggest multiple interior pathways for O(2) entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.

Multiphoton absorption of myoglobin-nitric oxide complex: relaxation by D-NEMD of a stationary state

The photodissociation and geminate recombination of nitric oxide in myoglobin, under continuous illumination, is modeled computationally. The relaxation of the photon energy into the protein matrix is also considered in a single simulation scheme that mimics a complete experimental setup. The dynamic approach to non-equilibrium molecular dynamics is used, starting from a steady state, to compute its relaxation to equilibrium. Simulations are conducted for the native form of sperm whale myoglobin and for two other mutants, V68W and L29F, illustrating a fair diversity of spatial and temporal geminate recombination processes. Energy flow to the heme and immediate protein environment provide hints to allostery. In particular, a pathway of energy flow between the heme and the FG loop is illustrated. Although the simulations were conducted for myoglobin only, the thermal fluctuations of the FG corner are in agreement with the large structural shifts of FG during the allosteric transition of tetrameric hemoglobin.

Time-resolved protein nanocrystallography using an X-ray free-electron laser

We demonstrate the use of an X-ray free electron laser synchronized with an optical pump laser to obtain X-ray diffraction snapshots from the photoactivated states of large membrane protein complexes in the form of nanocrystals flowing in a liquid jet. Light-induced changes of Photosystem I-Ferredoxin co-crystals were observed at time delays of 5 to 10 µs after excitation. The result correlates with the microsecond kinetics of electron transfer from Photosystem I to ferredoxin. The undocking process that follows the electron transfer leads to large rearrangements in the crystals that will terminally lead to the disintegration of the crystals. We describe the experimental setup and obtain the first time-resolved femtosecond serial X-ray crystallography results from an irreversible photo-chemical reaction at the Linac Coherent Light Source. This technique opens the door to time-resolved structural studies of reaction dynamics in biological systems.

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

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

BioCARS: a synchrotron resource for time-resolved X-ray science

BioCARS, a NIH-supported national user facility for macromolecular time-resolved X-ray crystallography at the Advanced Photon Source (APS), has recently completed commissioning of an upgraded undulator-based beamline optimized for single-shot laser-pump X-ray-probe measurements with time resolution as short as 100 ps. The source consists of two in-line undulators with periods of 23 and 27 mm that together provide high-flux pink-beam capability at 12 keV as well as first-harmonic coverage from 6.8 to 19 keV. A high-heat-load chopper reduces the average power load on downstream components, thereby preserving the surface figure of a Kirkpatrick-Baez mirror system capable of focusing the X-ray beam to a spot size of 90 µm horizontal by 20 µm vertical. A high-speed chopper isolates single X-ray pulses at 1 kHz in both hybrid and 24-bunch modes of the APS storage ring. In hybrid mode each isolated X-ray pulse delivers up to ~4 × 10(10) photons to the sample, thereby achieving a time-averaged flux approaching that of fourth-generation X-FEL sources. A new high-power picosecond laser system delivers pulses tunable over the wavelength range 450-2000 nm. These pulses are synchronized to the storage-ring RF clock with long-term stability better than 10 ps RMS. Monochromatic experimental capability with Biosafety Level 3 certification has been retained.

Revisiting and computing reaction coordinates with Directional Milestoning

The method of Directional Milestoning is revisited. We start from an exact and more general expression and state the conditions and validity of the memory-loss approximation. An algorithm to compute a reaction coordinate from Directional Milestoning data is presented. The reaction coordinate is calculated as a set of discrete jumps between Milestones that maximizes the flux between two stable states. As an application we consider a conformational transition in solvated adenosine. We compare a long molecular dynamic trajectory with Directional Milestoning and discuss the differences between the maximum flux path and minimum energy coordinates.

Cluster analysis of time-dependent crystallographic data: Direct identification of time-independent structural intermediates

The initial output of a time-resolved macromolecular crystallography experiment is a time-dependent series of difference electron density maps that displays the time-dependent changes in underlying structure as a reaction progresses. The goal is to interpret such data in terms of a small number of crystallographically refinable, time-independent structures, each associated with a reaction intermediate; to establish the pathways and rate coefficients by which these intermediates interconvert; and thereby to elucidate a chemical kinetic mechanism. One strategy toward achieving this goal is to use cluster analysis, a statistical method that groups objects based on their similarity. If the difference electron density at a particular voxel in the time-dependent difference electron density (TDED) maps is sensitive to the presence of one and only one intermediate, then its temporal evolution will exactly parallel the concentration profile of that intermediate with time. The rationale is therefore to cluster voxels with respect to the shapes of their TDEDs, so that each group or cluster of voxels corresponds to one structural intermediate. Clusters of voxels whose TDEDs reflect the presence of two or more specific intermediates can also be identified. From such groupings one can then infer the number of intermediates, obtain their time-independent difference density characteristics, and refine the structure of each intermediate. We review the principles of cluster analysis and clustering algorithms in a crystallographic context, and describe the application of the method to simulated and experimental time-resolved crystallographic data for the photocycle of photoactive yellow protein.

Quantitative analysis of ligand migration from transition networks

In this work we use transition network analysis for the first time to investigate ligand migration in truncated hemoglobin (trHbN) and obtain kinetic information about the docking-site dynamics in the protein. A comparison with explicit water molecular dynamics simulations (100 ns in total) shows that the rate constants derived from the network analysis are realistic. The transition network analysis provides 1) The time-resolved connectivity network in the protein; 2) The half-lives of the docking sites; 3) The transition timescales between two given docking sites; and 4) The extent of population transfer among different docking sites of the protein as a function of lag time. We investigate the role of the Tyr33 and Gln58 residues in ligand migration by studying ligand migration in four mutants of trHbN. The mutation study suggests that residues Tyr33 and Gln58 stabilize the NO ligand in the Xe2 docking site of trHbN, thus facilitating the efficiency of the NO detoxification reaction.

Kinetics of carbon monoxide migration and binding in solvated neuroglobin as revealed by molecular dynamics simulations and quantum mechanical calculations

Neuroglobin (Ngb) is a globular protein that reversibly binds small ligands at the six coordination position of the heme. With respect to other globins similar to myoglobin, Ngb displays some peculiarities as the topological reorganization of the internal cavities coupled to the sliding of the heme, or the binding of the endogenous distal histidine to the heme in the absence of an exogenous ligand. In this Article, by using multiple (independent) molecular dynamics trajectories (about 500 ns in total), the migration pathways of photolized carbon monoxide (CO) within solvated Ngb were analyzed, and a quantitative description of CO migration and corresponding kinetics was obtained. MD results, combined with quantum mechanical calculations on the CO-heme binding-unbinding reaction step in Ngb, allowed construction of a quantitative model representing the relevant steps of CO migration and rebinding.

‘It's hollow’: the function of pores within myoglobin

Despite a century of research, the cellular function of myoglobin (Mb), the mechanism regulating oxygen (O2) transport in the cell and the structure–function relationship of Mb remain incompletely understood. In particular, the presence and function of pores within Mb have attracted much recent attention. These pores can bind to Xe as well as to other ligands. Indeed, recent cryogenic X-ray crystallographic studies using novel techniques have captured snapshots of carbon monoxide (CO) migrating through these pores. The observed movement of the CO molecule from the heme iron site to the internal cavities and the associated structural changes of the amino acid residues around the cavities confirm the integral role of the pores in forming a ligand migration pathway from the protein surface to the heme. These observations resolve a long-standing controversy – but how these pores affect the physiological function of Mb poses a striking question at the frontier of biology.

Picosecond primary structural transition of the heme is retarded after nitric oxide binding to heme proteins

We investigated the ultrafast structural transitions of the heme induced by nitric oxide (NO) binding for several heme proteins by subpicosecond time-resolved resonance Raman and femtosecond transient absorption spectroscopy. We probed the heme iron motion by the evolution of the iron-histidine Raman band intensity after NO photolysis. Unexpectedly, we found that the heme response and iron motion do not follow the kinetics of NO rebinding. Whereas NO dissociation induces quasi-instantaneous iron motion and heme doming (<0.6 ps), the reverse process results in a much slower picosecond movement of the iron toward the planar heme configuration after NO binding. The time constant for this primary domed-to-planar heme transition varies among proteins (approximately 30 ps for myoglobin and its H64V mutant, approximately 15 ps for hemoglobin, approximately 7 ps for dehaloperoxidase, and approximately 6 ps for cytochrome c) and depends upon constraints exerted by the protein structure on the heme cofactor. This observed phenomenon constitutes the primary structural transition in heme proteins induced by NO binding.