Ligand binding and conformational motions in myoglobin

CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

Fast collective motions of backbone in transmembrane α helices are critical to water transfer of aquaporin

Fast collective motions are widely present in biomolecules, but their functional relevance remains unclear. Herein, we reveal that fast collective motions of backbone are critical to the water transfer of aquaporin Z (AqpZ) by using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. A total of 212 residue site-specific dipolar order parameters and 158 15N spin relaxation rates of the backbone are measured by combining the 13C- and 1H-detected multidimensional ssNMR spectra. Analysis of these experimental data by theoretic models suggests that the small-amplitude (~10°) collective motions of the transmembrane α helices on the nanosecond-to-microsecond timescales are dominant for the dynamics of AqpZ. The MD simulations demonstrate that these collective motions are critical to the water transfer efficiency of AqpZ by facilitating the opening of the channel and accelerating the water-residue hydrogen bonds renewing in the selectivity filter region.

Probing the functional hotspots inside protein hydrophobic pockets by in situ photochemical trifluoromethylation and mass spectrometry

Due to the complex high-order structures and interactions of proteins within an aqueous solution, a majority of chemical functionalizations happen on the hydrophilic sites of protein external surfaces which are naturally exposed to the solution. However, the hydrophobic pockets inside proteins are crucial for ligand binding and function as catalytic centers and transporting tunnels. Herein, we describe a reagent pre-organization and in situ photochemical trifluoromethylation strategy to profile the functional sites inside the hydrophobic pockets of native proteins. Unbiased mass spectrometry profiling was applied for the characterization of trifluoromethylated sites with high sensitivity. Native proteins including myoglobin, trypsin, haloalkane dehalogenase, and human serum albumin have been engaged in this mild photochemical process and substantial hydrophobic site-specific and structure-selective trifluoromethylation substitutes are obtained without significant interference to their bioactivity and structures. Sodium triflinate is the only reagent required to functionalize the unprotected proteins with wide pH-range tolerance and high biocompatibility. This "in-pocket" activation model provides a general strategy to modify the potential binding pockets and gain essential structural insights into the functional hotspots inside protein hydrophobic pockets.

Coenzyme-binding pathway on glutamate dehydrogenase suggested from multiple-binding sites visualized by cryo-electron microscopy

The structure of hexameric glutamate dehydrogenase (GDH) in the presence of the coenzyme nicotinamide adenine dinucleotide phosphate (NADP) was visualized using cryogenic transmission electron microscopy to investigate the ligand-binding pathways to the active site of the enzyme. Each subunit of GDH comprises one hexamer-forming core domain and one nucleotide-binding domain (NAD domain), which spontaneously opens and closes the active-site cleft situated between the two domains. In the presence of NADP, the potential map of GDH hexamer, assuming D3 symmetry, was determined at a resolution of 2.4 Å, but the NAD domain was blurred due to the conformational variety. After focused classification with respect to the NAD domain, the potential maps interpreted as NADP molecules appeared at five different sites in the active-site cleft. The subunits associated with NADP molecules were close to one of the four metastable conformations in the unliganded state. Three of the five binding sites suggested a pathway of NADP molecules to approach the active-site cleft for initiating the enzymatic reaction. The other two binding modes may rarely appear in the presence of glutamate, as demonstrated by the reaction kinetics. Based on the visualized structures and the results from the enzymatic kinetics, we discussed the binding modes of NADP to GDH in the absence and presence of glutamate.

Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

A Novel Rabbit Anti-Myoglobin Monoclonal Antibody's Potential Application in Rhabdomyolysis Associated Acute Kidney Injury

Myoglobin (Mb) is the main constituent of vertebrate skeletal muscle and myocardium and plays an essential role in oxygen binding, storage, transport, and earliest disease diagnosis. This study focuses on preparing the novel recombinant rabbit anti-Mb monoclonal antibody and applying it to a diagnosis of Mb deposition in rhabdomyolysis-associated acute kidney injury (RM-AKI). The full-length coding sequence of rat Mb was cloned and expressed, and the high-quality and titer rabbit anti-Mb polyclonal antibodies were produced by the immunogen His-Mb fusion protein. A new hybridoma cell was obtained by hybridoma screening technology. With the help of DNA sequencing and a molecular clonal, anti-Mb monoclonal antibody heavy and light chains expression plasmid was constructed. Finally, the recombinant rabbit anti-Mb monoclonal antibody with extraordinarily high affinity (K_D = 1.21 pM) was obtained. Meanwhile, it had broad species reactivity (mouse, rat, human, and horse) and good tissue specificity (skeletal muscle and myocardium). It also had a very good performance in western blotting, immunohistochemistry, and immunofluorescence assay to detect the Mb level in the kidney, myocardium, and skeletal muscle of RM-AKI. This study will be significantly helpful for Mb-associated disease diagnosis, and pathogenesis exploration, and further may act as a neutralizing antibody for disease treatment.

The "violin model": Looking at community networks for dynamic allostery

Allosteric regulation of proteins continues to be an engaging research topic for the scientific community. Models describing allosteric communication have evolved from focusing on conformation-based descriptors of protein structural changes to appreciating the role of internal protein dynamics as a mediator of allostery. Here, we explain a "violin model" for allostery as a contemporary method for approaching the Cooper-Dryden model based on redistribution of protein thermal fluctuations. Based on graph theory, the violin model makes use of community network analysis to functionally cluster correlated protein motions obtained from molecular dynamics simulations. This Review provides the theory and workflow of the methodology and explains the application of violin model to unravel the workings of protein kinase A.

Ruthenium hydrides encapsulated in sol–gel glasses exhibit new ultrafast vibrational dynamics

Vibrational dynamics were measured by IR pump–probe spectroscopy and two-dimensional IR spectroscopy for triruthenium dodecacarbonyl and the undecacarbonyl hydride that forms when it is encapsulated in an alumina sol–gel glass. For comparison, a triruthenium undecacarbonyl hydride salt was also synthesized and studied in neat solution to identify the potential influence of the confined solvent environment on the dynamics experienced by carbon monoxide ligands. The vibrational lifetime was found to be significantly decreased for both hydride species relative to the dodecacarbonyl compound. Conversely, spectral diffusion of the CO vibrations was measured to be faster for the parent compound. The most significant dynamic changes occurred upon transformation from the starting compound to the hydride, while only minor differences were observed between the dynamics of the freely dissolved and sol–gel encapsulated hydrides. The results suggest that the structural change to the hydride has the largest impact on the dynamics and that its improved catalytic properties likely do not originate from confined solvent effects.

Experimental demonstration of the novel "van-Hove integral method (vHI)" for measuring diffusive dynamics by elastic neutron scattering

Quasi-elastic neutron scattering (QENS)-based on the seminal work of Nobel Laureate Brockhouse-has been one of the major methods for studying pico-second to nano-second diffusive dynamics over the past 70 years. This is regarded as an "inelastic" method for dynamics. In contrast, we recently proposed a new neutron-scattering method for dynamics, which uses the elastic line of the scattering to access system dynamics directly in the time domain (Benedetto and Kearley in Sci Rep 9:11284, 2019). This new method has been denoted "vHI" that stands for "van Hove Integral". The reason is that, under certain conditions, the measured elastic intensity corresponds to the running-time integral of the intermediate scattering function, [Formula: see text], up to a time that is inversely proportional to the energy band-width incident on the sample. As a result, [Formula: see text] is accessed from the time derivative of the measured vHI profile. vHI has been supported by numerical and Monte-Carlo simulations, but has been difficult to validate experimentally due to the lack of a suitable instrument. Here we show that vHI works in practice, which we achieved by using a simple modification to the standard QENS backscattering spectrometer methodology. Basically, we varied the neutron-energy band-widths incident at the sample via a step-wise variation of the frequency of the monochromator Doppler-drive. This provides a measurement of the vHI profile at the detectors. The same instrument and sample were also used in standard QENS mode for comparison. The intermediate scattering functions, [Formula: see text], obtained by the two methods-vHI and QENS-are strikingly similar providing a direct experimental validation of the vHI method. Perhaps surprisingly, the counting statistics of the two methods are comparable even though the instrument used was expressly designed for QENS. This shows that the methodology modification adopted here can be used in practice to access vHI profiles at many of the backscattering spectrometers worldwide. We also show that partial integrations of the measured QENS spectrum cannot provide the vHI profile, which clarifies a common misconception. At the same time, we show a novel approach which does access [Formula: see text] from QENS spectra.

A Minireview on Temperature Dependent Protein Conformational Sampling

In this minireview we discuss the role of the more subtle conformational change-protein conformational sampling and connect it to the classic relationship of protein structure and function. The theory of pre-existing functional states of protein are discussed in context of alternate protein conformational sampling. Last, we discuss how temperature, ligand binding and mutations affect the protein conformational sampling mode which is linked to the protein function regulation. The review includes several protein systems that showed temperature dependent protein conformational sampling. We also specifically included two enzyme systems, thermophilic alcohol dehydrogenase (ht-ADH) and thermolysin which we previously studied when discussing temperature dependent protein conformational sampling.

Moving beyond static snapshots: Protein dynamics and the Protein Data Bank

Proteins are the molecular machines of living systems. Their dynamics are an intrinsic part of their evolutionary selection in carrying out their biological functions. Although the dynamics are more difficult to observe than a static, average structure, we are beginning to observe these dynamics and form sound mechanistic connections between structure, dynamics, and function. This progress is highlighted in case studies from myoglobin and adenylate kinase to the ribosome and molecular motors where these molecules are being probed with a multitude of techniques across many timescales. New approaches to time-resolved crystallography are allowing simple “movies” to be taken of proteins in action, and new methods of mapping the variations in cryo-electron microscopy are emerging to reveal a more complete description of life’s machines. The results of these new methods are aided in their dissemination by continual improvements in curation and distribution by the Protein Data Bank and their partners around the world.

Conformational Plasticity in Human Heme-Based Dioxygenases

Human indoleamine 2,3-dioxygenase 1 (hIDO1) and human tryptophan dioxygenase (hTDO) are two important heme proteins that degrade the essential amino acid, l-tryptophan (Trp), along the kynurenine pathway. The two enzymes share a similar active site structure and an analogous catalytic mechanism, but they exhibit a variety of distinct functional properties. Here we used carbon monoxide (CO) as a structural probe to interrogate how the functionalities of the two enzymes are encoded in their structures. With X-ray crystallography, we detected an unexpected photochemical intermediate trapped in a crystal of the hIDO1-CO-Trp complex, where CO is photolyzed from the heme iron by X-rays at cryogenic temperatures (100 K). The CO photolysis triggers a large-scale migration of the substrate Trp, as well as the photolyzed CO, from the active site to a temporary binding site, Sa*. It is accompanied by a large conformational change to an active site loop, JK-Loop^C, despite the severely restricted protein motion under the frozen conditions, which highlights the remarkable conformational plasticity of the hIDO1 protein. Comparative studies of a crystal of the hTDO-CO-Trp complex show that CO and Trp remain bound in the active site under comparable X-ray illumination, indicating a much more rigid protein architecture. The data offer important new insights into the structure and function relationships of the heme-based dioxygenases and provide new guidelines for structure-based design of inhibitors targeting them.

The Impact of Electron Correlation on Describing QM/MM Interactions in the Attendant Molecular Dynamics Simulations of CO in Myoglobin

The impact of the dispersion and electron correlation effects on describing quantum mechanics/molecular mechanics (QM/MM) interactions in QM/MM molecular dynamics (MD) simulations was explored by performing a series of up to 2 ns QM/MM MD simulations on the B states of the myoglobin-carbon monoxide (MbCO) system. The results indicate that both dispersion and electron correlations play significant roles in the simulation of the ratios of two B states (B₁/B₂), which suggests that the inclusion of the electron correlation effects is essential for accurately modeling the interactions between QM and MM subsystems. We found that the QM/MM interaction energies between the CO and the surroundings statistically present a linear correlation with the electric fields along the CO bond. This indicates that QM/MM interactions can be described by a simple physical model of a dipole with constant moment under the action of the electric fields. The treatment provides us with an accurate and effective approach to account for the electron correlation effects in QM/MM MD simulations.

A Quantitative Comparison of the Counting Significance of van Hove Integral Spectroscopy and Quasielastic Neutron Scattering

We have recently proposed a new method to access system dynamics via neutron scattering based on measuring the elastic scattered intensity: By varying the energy band-width that impinges on the sample (also known as instrumental energy resolution), the purely elastic-scattering from this variation is the running time-integral of the intermediate scattering function (I(t)) [Benedetto and Kearley, Sci. Rep. 9, 11284, 2019]. In this correspondence we denote our method "vHI", which stands for "van Hove Integral". The method is now widely accepted as "valid" and here we focus on the efficiency of the vHI method compared with the standard quasi-elastic neutron scattering (QENS) method. We use a numerical Monte-Carlo simulation of an instrument that is equally capable of measuring QENS and vHI under identical conditions. For an "experiment" in which the same number of neutrons enter the instrument, we present comparisons between QENS and vHI at three levels of data-reduction. Firstly, at the raw-data level vHI achieves 100 times more neutrons at the detector than QENS. Secondly, vHI has a factor of 2 less statistical error, which would translate to an overall gain of 4 for vHI in counting-time. Lastly, we compare the distortions caused in obtaining the final I(t) via time-Fourier transform (QENS) and polynomial time-derivative (vHI). Here, the statistical error is 10 times smaller for vHI. This last comparison is the most important result where the 10 times smaller residual for vHI gives a net gain in counting time of 100 better than QENS to obtain the same underlying dynamics of the system under study.

Effect of Occluded Ligand Migration on the Kinetics and Structural Dynamics of Homodimeric Hemoglobin

Small molecules such as molecular oxygen, nitric oxide, and carbon monoxide play important roles in life, and many proteins require the transport of small molecules to and from the bulk solvent for their function. Ligand migration within a protein molecule is expected to be closely related to the overall structural changes of the protein, but the detailed and quantitative connection remains elusive. For example, despite numerous studies, how occluded ligand migration affects the kinetics and structural dynamics of the R-T transition remains unclear. To shed light on this issue, we chose homodimeric hemoglobin (HbI) with the I114F mutation (I114F), which is known to interfere with ligand migration between the primary and secondary docking sites, and studied its kinetics and structural dynamics using time-resolved X-ray solution scattering. The kinetic analysis shows that I114F has three structurally distinct intermediates (I₁, I₂, and I₃) as in the wild type (WT), but its geminate CO recombination occurs directly from I₁ without the path via I₂ observed in WT. Moreover, the structural transitions, which involve ligand migration (the transitions from I₁ to I₂ and from I₃ to the initial state), are decelerated compared to WT. The structural analysis revealed that I114F involves generally smaller structural changes in all three intermediates compared to WT.

Direct observation of ligand migration within human hemoglobin at work

Human hemoglobin is the textbook example of the stereochemistry of an allosteric protein and of the exquisite control that a protein can exert over ligand binding. However, the fundamental basis by which the protein facilitates the ligand movement remains unknown. In this study, we used cryogenic X-ray crystallography and a high-repetition pulsed laser irradiation technique to elucidate the atomic details of ligand migration processes in hemoglobin after photolysis of the bound CO. Our data clarify the distinct CO migration pathways in the individual subunits of hemoglobin and unravel the functional roles of the internal cavities and neighboring amino acid residues in ligand exit and entry. Our results also demonstrate the high gas permeability and porosity of hemoglobin, facilitating O₂ delivery.

Hemoglobin is one of the best-characterized proteins with respect to structure and function, but the internal ligand diffusion pathways remain obscure and controversial. Here we captured the CO migration processes in the tense (T), relaxed (R), and second relaxed (R2) quaternary structures of human hemoglobin by crystallography using a high-repetition pulsed laser technique at cryogenic temperatures. We found that in each quaternary structure, the photodissociated CO molecules migrate along distinct pathways in the α and β subunits by hopping between the internal cavities with correlated side chain motions of large nonpolar residues, such as α14Trp(A12), α105Leu(G12), β15Trp(A12), and β71Phe(E15). We also observe electron density evidence for the distal histidine [α58/β63His(E7)] swing-out motion regardless of the quaternary structure, although less evident in α subunits than in β subunits, suggesting that some CO molecules have escaped directly through the E7 gate. Remarkably, in T-state Fe(II)-Ni(II) hybrid hemoglobins in which either the α or β subunits contain Ni(II) heme that cannot bind CO, the photodissociated CO molecules not only dock at the cavities in the original Fe(II) subunit, but also escape from the protein matrix and enter the cavities in the adjacent Ni(II) subunit even at 95 K, demonstrating the high gas permeability and porosity of the hemoglobin molecule. Our results provide a comprehensive picture of ligand movements in hemoglobin and highlight the relevance of cavities, nonpolar residues, and distal histidines in facilitating the ligand migration.

Conformational sampling and kinetics changes across a non-Arrhenius break point in the enzyme thermolysin

Numerous studies have suggested a significant role that protein dynamics play in optimizing enzyme catalysis, and changes in conformational sampling offer a window to explore this role. Thermolysin from Bacillus thermoproteolyticus rokko, which is a heat-stable zinc metalloproteinase, serves here as a model system to study changes of protein function and conformational sampling across a temperature range of 16-36 °C. The temperature dependence of kinetics of thermolysin showed a biphasic transition at 26 °C that points to potential conformational and dynamic differences across this temperature. The non-Arrhenius behavior observed resembled results from previous studies of a thermophilic alcohol dehydrogenase enzyme, which also indicated a biphasic transition at ambient temperatures. To explore the non-Arrhenius behavior of thermolysin, room temperature crystallography was applied to characterize structural changes in a temperature range across the biphasic transition temperature. The alternate conformation of side chain fitting to electron density of a group of residues showed a higher variability in the temperature range from 26 to 29 °C, which indicated a change in conformational sampling that correlated with the non-Arrhenius break point.

Role of hydration water in the onset of protein structural dynamics

Proteins are the molecular workhorses in a living organism. Their 3D structures are animated by a multitude of equilibrium fluctuations and specific out-of-equilibrium motions that are required for proteins to be biologically active. When studied as a function of temperature, functionally relevant dynamics are observed at and above the so-called protein dynamical transition (~240 K) in hydrated, but not in dry proteins. In this review we present and discuss the main experimental and computational results that provided evidence for the dynamical transition, with a focus on the role of hydration water dynamics in sustaining functional protein dynamics. The coupling and mutual influence of hydration water dynamics and protein dynamics are discussed and the hypotheses illustrated that have been put forward to explain the physical origin of their onsets.

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.

Dynamic allostery-based molecular workings of kinase:peptide complexes

A dense interplay between structure and dynamics underlies the working of proteins, especially enzymes. Protein kinases are molecular switches that are optimized for their regulation rather than catalytic turnover rates. Using long-simulations dynamic allostery analysis, this study describes an exploration of the dynamic kinase:peptide complex. We have used protein kinase A (PKA) as a model system as a generic prototype of the protein kinase superfamily of signaling enzymes. Our results explain the role of dynamic coupling of active-site residues that must work in coherence to provide for a successful activation or inhibition response from the kinase. Amino acid networks-based community analysis allows us to ponder the conformational entropy of the kinase:nucleotide:peptide ternary complex. We use a combination of 7 peptides that include 3 types of PKA-binding partners: Substrates, products, and inhibitors. The substrate peptides provide for dynamic insights into the enzyme:substrate complex, while the product phospho-peptide allows for accessing modes of enzyme:product release. Mapping of allosteric communities onto the PKA structure allows us to locate the more unvarying and flexible dynamic regions of the kinase. These distributions, when correlated with the structural elements of the kinase core, allow for a detailed exploration of key dynamics-based signatures that could affect peptide recognition and binding at the kinase active site. These studies provide a unique dynamic allostery-based perspective to kinase:peptide complexes that have previously been explored only in a structural or thermodynamic context.

Time-Resolved Macromolecular Crystallography at Pulsed X-ray Sources

The focus of structural biology is shifting from the determination of static structures to the investigation of dynamical aspects of macromolecular function. With time-resolved macromolecular crystallography (TRX), intermediates that form and decay during the macromolecular reaction can be investigated, as well as their reaction dynamics. Time-resolved crystallographic methods were initially developed at synchrotrons. However, about a decade ago, extremely brilliant, femtosecond-pulsed X-ray sources, the free electron lasers for hard X-rays, became available to a wider community. TRX is now possible with femtosecond temporal resolution. This review provides an overview of methodological aspects of TRX, and at the same time, aims to outline the frontiers of this method at modern pulsed X-ray sources.

Effect of Single-Point Mutations on Nitric Oxide Rebinding and the Thermodynamic Stability of Myoglobin

The effect of single amino acid mutations on the rebinding dynamics of nitrogen monoxide (NO) to myoglobin is investigated using reactive molecular dynamics simulations. In particular, mutations of residues surrounding the heme-active site (Leu29, His64, Val68) were considered. Consistent with experiments, all mutations studied here have a significant effect on the kinetics of the NO-rebinding process, which consists of a rapid (several 10 ps) and a slow (100s of ps) time scale. For all modifications considered, the time scales and rebinding fractions agree to within a few percents with results from experiments by adjusting one single, physically meaningful, conformationally averaged quantity: the asymptotic energy separation between the NO-bound (²A) and photodissociated (⁴A) states. It is furthermore shown that the thermodynamic stability of wild-type versus mutant Mb for the ligand-free and ligand-bound variants of the protein can be described by the same computational model. Therefore, ligand kinetics and thermodynamics are related in a direct fashion akin to Φ-value analysis, which establishes a relationship between protein folding rates and thermal stability of proteins.

Kinetics of polymer tumbling in shear flow: A coarse-grained description

We study the classic problem of dynamical evolution of a polymer in a shear flow. Interestingly, the polymer goes through several distinctly identifiable conformations during its passage from coiled to stretched states back and forth. We identify these conformations assumed by the polymer while tumbling and study the kinetics of the process in terms of the residence and recurrence times of individual conformations. The distribution of residence times exhibits exponentially decaying tails which helps us build an effective Markovian picture of the truly non-Markovian problem. We present the explicit W matrix for the coarse-grained evolution via a master equation and study its elements as a function of the Weissenberg number. We show that the timescales of decay of the autocorrelation function for the full Langevin dynamics compare quite well with the approximate results from the master equation approach.

Kinetic Isotope Effects as Probes for Hydrogen Tunneling, Coupled Motion and Dynamics Contributions to Enzyme Catalysis

Since the early days of enzymology attempts have been made to deconvolute the various contributions of physical phenomena to enzyme catalysis. Here we present experimental and theoretical studies that examine the possible role of hydrogen tunneling, coupled motion, and enzyme dynamics in catalysis. In this review, we first introduce basic concepts of enzyme catalysis from a physical chemistry point of view. Then, we present several recent developments in the application of experimental tools that can probe tunneling, coupled motion, dynamic effects and other possible physical phenomena that may contribute to catalysis. These tools include kinetic isotope effects (KIEs), their temperature dependency and H/D/T mutual relations (the Swain–Schaad relationship). Several theories and models that assist in understanding those phenomena are also described. The possibility that these models invoke a direct role for the enzyme's dynamics (environmental fluctuations and rearrangements) is discussed. Finally, the need to compare the enzymatic reaction to the uncatalyzed one while investigating contributions to catalysis is emphasised.

Kinetic Study of Ligand Binding and Conformational Changes in Inducible Nitric Oxide Synthase

Nitric oxide synthases (NOSs) are heme enzymes that generate highly reactive nitric oxide from l-arginine (l-Arg) in a complex mechanism that is still only partially understood. We have studied carbon monoxide (CO) binding to the oxygenase domain of murine inducible NOS (iNOS) by using flash photolysis. The P420 and P450 forms of the enzyme, assigned to a protonated and unprotonated proximal cysteine, through which the heme is anchored to the protein, show markedly different CO rebinding properties. The data suggest that P420 has a widely open distal pocket that admits water. CO rebinding to the P450 form strongly depends on the presence of the substrate l-Arg, the intermediate N^ω-hydroxy-l-arginine, and the cofactor tetrahydrobiopterin. After adding these small molecules to the enzyme solution, the CO kinetics change slowly over the hours. This process can be described as a relaxation from a fast rebinding, metastable species to a slowly rebinding, thermodynamically stable species, which we associate with the enzymatically active form. Our results allow us to determine kinetic parameters of l-Arg binding to the ferrous deoxy iNOS protein for the first time and also provide clues regarding the nature of structural differences between the two interconverting species.

Dynamical comparison between myoglobin and hemoglobin

Myoglobin and hemoglobin are globular hemeproteins, when the former is a monomer and the latter a heterotetramer. Despite the structural similarity of myoglobin to α and β subunits of hemoglobin, there is a functional difference between the two proteins, owing to the quaternary structure of hemoglobin. The effect of the quaternary structure of hemoglobin on the intrinsic dynamics of its subunits is explored by dynamical comparison of the two proteins. Anisotropic Network Model modes of motion were calculated for hemoglobin and myoglobin. Dynamical comparison between the proteins was performed using global and local Anisotropic Network Model mode alignment algorithms based on the algorithms of Smith-Waterman and Needleman-Wunsch for sequence comparison. The results indicate that the quaternary structure of hemoglobin substantially alters the intrinsic dynamics of its subunits, an effect that may contribute to the functional difference between the two proteins. Local dynamics similarity between the proteins is still observed at the major exit route of the ligand.

Selective Permeability of Carboxysome Shell Pores to Anionic Molecules

Carboxysomes are closed polyhedral cellular microcompartments that increase the efficiency of carbon fixation in autotrophic bacteria. Carboxysome shells consist of small proteins that form hexameric units with semipermeable central pores containing binding sites for anions. This feature is thought to selectively allow access to RuBisCO enzymes inside the carboxysome by HCO3- (the dominant form of CO2 in the aqueous solution at pH 7.4) but not O2, which leads to a nonproductive reaction. To test this hypothesis, here we use molecular dynamics simulations to characterize the energetics and permeability of CO2, O2, and HCO3- through the central pores of two different shell proteins, namely, CsoS1A of α-carboxysome and CcmK4 of β-carboxysome shells. We find that the central pores are in fact selectively permeable to anions such as HCO3-, as predicted by the model.

Cold-Adaptation Signatures in the Ligand Rebinding Kinetics to the Truncated Hemoglobin of the Antarctic Bacterium Pseudoalteromonas haloplanktis TAC125

Cold-adapted organisms have evolved proteins endowed with higher flexibility and lower stability in comparison to their thermophilic homologues, resulting in enhanced reaction rates at low temperatures. In this context, protein-bound water molecules were suggested to play a major role, and their weaker interactions at protein active sites have been associated with cold adaptation. In this work, we tested this hypothesis on truncated hemoglobins (a family of microbial heme-proteins of yet-unclear function) applying molecular dynamics simulations and ligand-rebinding kinetics on a protein from the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125 in comparison with its thermophilic Thermobifida fusca homologue. The CO rebinding kinetics of the former highlight several geminate phases, with an unusually long-lived geminate intermediate. An articulated tunnel with at least two distinct docking sites was identified by analysis of molecular dynamics simulations and was suggested to be at the origin of the unusual geminate rebinding phase. Water molecules are present in the distal pocket, but their stabilization by TrpG8, TyrB10, and HisCD1 is much weaker than in thermophilic Thermobifida fusca truncated hemoglobin, resulting in a faster geminate rebinding. Our results support the hypothesis that weaker water-molecule interactions at the reaction site are associated with cold adaptation.

Structure and energy based quantitative missense variant effect analysis provides insights into drug resistance mechanisms of anaplastic lymphoma kinase mutations

Anaplastic lymphoma kinase (ALK) is considered as a validated molecular target in multiple malignancies, such as non-small cell lung cancer (NSCLC). However, the effectiveness of molecularly targeted therapies using ALK inhibitors is almost universally limited by drug resistance. Drug resistance to molecularly targeted therapies has now become a major obstacle to effective cancer treatment and personalized medicine. It is of particular importance to provide an improved understanding on the mechanisms of resistance of ALK inhibitors, thus rational new therapeutic strategies can be developed to combat resistance. We used state-of-the-art computational approaches to systematically explore the mutational effects of ALK mutations on drug resistance properties. We found the activation of ALK was increased by substitution with destabilizing mutations, creating the capacity to confer drug resistance to inhibitors. In addition, results implied that evolutionary constraints might affect the drug resistance properties. Moreover, an extensive profile of drugs against ALK mutations was constructed to give better understanding of the mechanism of drug resistance based on structural transitions and energetic variation. Our work hopes to provide an up-to-date mechanistic framework for understanding the mechanisms of drug resistance induced by ALK mutations, thus tailor treatment decisions after the emergence of resistance in ALK-dependent diseases.

Slow Dynamics around a Protein and Its Coupling to Solvent

Solvent is essential for protein dynamics and function, but its role in regulating the dynamics remains debated. Here, we employ saturation transfer electron spin resonance (ST-ESR) to explore the issue and characterize the dynamics on a longer (from μs to s) time scale than has been extensively studied. We first demonstrate the reliability of ST-ESR by showing that the dynamical changeovers revealed in the spectra agree to liquid-liquid transition (LLT) in the state diagram of the glycerol/water system. Then, we utilize ST-ESR with four different probes to systematically map out the variation in local (site-specific) dynamics around a protein surface at subfreezing temperatures (180-240 K) in 10 mol % glycerol/water mixtures. At highly exposed sites, protein and solvent dynamics are coupled, whereas they deviate from each other when temperature is greater than LLT temperature (∼190 K) of the solvent. At less exposed sites, protein however exhibits a dynamic, which is distinct from the bulk solvent, throughout the temperature range studied. Dominant dynamic components are thus revealed, showing that (from low to high temperatures) the overall structural fluctuation, rotamer dynamics, and internal side-chain dynamics, in turn, dominate the temperature dependence of spin-label motions. The structural fluctuation component is relatively slow, collective, and independent of protein structural segments, which is thus inferred to a fundamental dynamic component intrinsic to protein. This study corroborates that bulk solvent plasticizes protein and facilitates rather than slaves protein dynamics.

Molecular evolution of myoglobin in the Tibetan Plateau endemic schizothoracine fish (Cyprinidae, Teleostei) and tissue-specific expression changes under hypoxia

Myoglobin (Mb) is an oxygen-binding hemoprotein that was once thought to be exclusively expressed in oxidative myocytes of skeletal and cardiac muscle where it serves in oxygen storage and facilitates intracellular oxygen diffusion. In this study, we cloned the coding sequence of the Mb gene from four species, representing three groups, of the schizothoracine fish endemic to the Qinghai-Tibetan Plateau (QTP), then conducted molecular evolution analyses. We also investigated tissue expression patterns of Mb and the expression response to moderate and severe hypoxia at the mRNA and protein levels in a representative of the highly specialized schizothoracine fish species, Schizopygopsis pylzovi. Molecular evolution analyses showed that Mb from the highly specialized schizothoracine fish have undergone positive selection and one positively selected residue (81L) was identified, which is located in the F helix, close to or in contact with the heme. We present tentative evidence that the Mb duplication event occurred in the ancestor of the schizothoracine and Cyprininae fish (common carp and goldfish), and that the Mb2 paralog was subsequently lost in the schizothoracine fish. In S. pylzovi, Mb mRNA is expressed in various tissues with the exception of the intestine and gill, but all such tissues, including the liver, muscle, kidney, brain, eye, and skin, expressed very low levels of Mb mRNA (< 8.0%) relative to that of the heart. The trace levels of Mb expression in non-muscle tissues are perhaps the major reason why non-muscle Mb remained undiscovered for so long. The expression response of the Mb gene to hypoxia at the mRNA and protein levels was strikingly different in S. pylzovi compared to that found in the common carp, medaka, zebrafish, and goldfish, suggesting that the hypoxia response of Mb in fish may be species and tissue-specific. Notably, severe hypoxia induced significant expression of Mb at the mRNA and protein levels in the S. pylzovi heart, which suggests Mb has a major role in the supply of oxygen to the heart of Tibetan Plateau fish.

A machine learning approach towards the prediction of protein–ligand binding affinity based on fundamental molecular properties

There is an exigency of transformation of the enormous amount of biological data available in various forms into some significant knowledge. We have tried to implement Machine Learning (ML) algorithm models on the protein–ligand binding affinity data already available to predict the binding affinity of the unknown. ML methods are appreciably faster and cheaper as compared to traditional experimental methods or computational scoring approaches. The prerequisites of this prediction are sufficient and unbiased features of training data and a prediction model which can fit the data well. In our study, we have applied Random forest and Gaussian process regression algorithms from the Weka package on protein–ligand binding affinity, which encompasses protein and ligand binding information from PdbBind database. The models are trained on the basis of selective fundamental information of both proteins and ligand, which can be effortlessly fetched from online databases or can be calculated with the availability of structure. The assessment of the models was made on the basis of correlation coefficient (R²) and root mean square error (RMSE). The Random forest model gave R² and RMSE of 0.76 and 1.31 respectively. We have also used our features and prediction models on the dataset used by others and found that our model with our features outperformed the existing ones.

There is an exigency of transformation of the enormous amount of biological data available in various forms into some significant knowledge.

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.

Ultrafast CO Kinetics in Heme Proteins: Adiabatic Ligand Binding and Heavy Atom Tunneling

The ultrafast kinetics of CO rebinding to carbon monoxide oxidation activator protein (ChCooA) are measured over a wide temperature range and compared with the kinetics of CO binding in other heme systems such as myoglobin (Mb) and hemoglobin (Hb). The Arrhenius prefactor for CO binding to ChCooA and protoheme (∼10¹¹ s^-1) is similar to what is found for spin-allowed NO binding to heme proteins and is several orders of magnitude larger than the prefactor of Mb and Hb (∼10⁹ s^-1). This indicates that the CO binding reaction is adiabatic, in contrast to the commonly held view that it is nonadiabatic due to spin-forbidden (ΔS = 2) selection rules. Under the adiabatic condition, entropic factors, rather than spin-selection rules, are the source of the reduced Arrhenius prefactors associated with CO binding in Mb and Hb. The kinetic response of ChCooA-CO is nonexponential at all temperatures, including 298 K, and is described quantitatively using a distribution of enthalpic rebinding barriers associated with heterogeneity in the heme doming conformation. Above the solvent glass transition (T_g ∼ 180 K), the rebinding progress slows as temperature increases, and this is ascribed to an evolution of the distribution toward increased heme doming and larger enthalpic barriers. Between T_g and ∼60 K, the nonexponential rebinding slows down as the temperature is lowered and the survival fraction follows the predictions expected for a quenched barrier distribution. Below ∼60 K the rebinding kinetics do not follow these predictions unless quantum mechanical tunneling along the heme doming coordinate is also included as an active channel for CO binding.

Probing the Role of the Heme Distal and Proximal Environment in Ligand Dynamics in the Signal Transducer Protein HemAT by Time-Resolved Step-Scan FTIR and Resonance Raman Spectroscopy

HemAT is a heme-containing oxygen sensor protein that controls aerotaxis. Time-resolved step-scan FTIR studies were performed on the isolated sensor domain and full-length HemAT proteins as well as on the Y70F (B-helix), L92A (E-helix), T95A (E-helix), and Y133F (G-helix) mutants to elucidate the effect of the site-specific mutations on the ligand dynamics subsequent to CO photolysis. The mutations aimed to perturb H-bonding and electrostatic interactions near the heme Fe-bound gaseous ligand (CO) and the heme proximal environment. Rebinding of CO to the heme Fe is biphasic in the sensor domain and full-length HemAT as well as in the mutants, with the exception of the Y133F mutant protein. The monophasic rebinding of CO in Y133F suggests that in the absence of the H-bond between Y133 and the heme proximal H123 residue the ligand rebinding process is significantly affected. The role of the proximal environment is also probed by resonance Raman photodissociation experiments, in which the Fe-His mode of the photoproduct of sensor domain HemAT-CO is detected at a frequency higher than that of the deoxy form in the difference resonance Raman spectra. The role of the conformational changes of Y133 (G-helix) and the role of the distal L92 and T95 residues (E-helix) in regulating ligand dynamics in the heme pocket are discussed.

Low-Temperature Decoupling of Water and Protein Dynamics Measured by Neutron Scattering

Water plays a major role in biosystems, greatly contributing to determine their structure, stability, and function. It is well known, for instance, that proteins require a minimum amount of water to be fully functional. Despite many years of intensive research, however, the detailed nature of protein-hydration water interactions is still partly unknown. The widely accepted "protein dynamical transition" scenario is based on perfect coupling between the dynamics of proteins and that of their hydration water, which has never been probed in depth experimentally. I present here high-resolution elastic neutron scattering measurements of the atomistic dynamics of lysozyme in water. The results show for the first time that the dynamics of proteins and of their hydration water are actually decoupled at low temperatures. This important result challenges the "protein dynamical transition" scenario and requires a new model to link protein dynamics to the dynamics of its hydration water.

Enhanced vibrational solvatochromism and spectral diffusion by electron rich substituents on small molecule silanes

Fourier transform infrared and two-dimensional IR (2D-IR) spectroscopies were applied to two different silanes in three different solvents. The selected solutes exhibit different degrees of vibrational solvatochromism for the Si-H vibration. Density functional theory calculations confirm that this difference in sensitivity is the result of higher mode polarization with more electron withdrawing ligands. This mode sensitivity also affects the extent of spectral diffusion experienced by the silane vibration, offering a potential route to simultaneously optimize the sensitivity of vibrational probes in both steady-state and time-resolved measurements. Frequency-frequency correlation functions obtained by 2D-IR show that both solutes experience dynamics on similar time scales and are consistent with a picture in which weakly interacting solvents produce faster, more homogeneous fluctuations. Molecular dynamics simulations confirm that the frequency-frequency correlation function obtained by 2D-IR is sensitive to the presence of hydrogen bonding dynamics in the surrounding solvation shell.

Size and Shape Controlled Crystallization of Hemoglobin for Advanced Crystallography

While high-throughput screening for protein crystallization conditions have rapidly evolved in the last few decades, it is also becoming increasingly necessary for the control of crystal size and shape as increasing diversity of protein crystallographic experiments. For example, X-ray crystallography (XRC) combined with photoexcitation and/or spectrophotometry requires optically thin but well diffracting crystals. By contrast, large-volume crystals are needed for weak signal experiments, such as neutron crystallography (NC) or recently developed X-ray fluorescent holography (XFH). In this article, we present, using hemoglobin as an example protein, some techniques for obtaining the crystals of controlled size, shape, and adequate quality. Furthermore, we describe a few case studies of applications of the optimized hemoglobin crystals for implementing the above mentioned crystallographic experiments, providing some hints and tips for the further progress of advanced protein crystallography.

Substrate Binding Primes Human Tryptophan 2,3-Dioxygenase for Ligand Binding

The human heme enzyme tryptophan 2,3-dioxygenase (hTDO) catalyzes the insertion of dioxygen into its cognate substrate, l-tryptophan (l-Trp). Its active site structure is highly dynamic, and the mechanism of enzyme-substrate-ligand complex formation and the ensuing enzymatic reaction is not yet understood. Here we have studied complex formation in hTDO by using time-resolved optical and infrared spectroscopy with carbon monoxide (CO) as a ligand. We have observed that both substrate-free and substrate-bound hTDO coexist in two discrete conformations with greatly different ligand binding rates. In the fast rebinding hTDO conformation, there is facile ligand access to the heme iron, but it is greatly hindered in the slowly rebinding conformation. Spectroscopic evidence implicates active site solvation as playing a crucial role for the observed kinetic differences. Substrate binding shifts the conformational equilibrium markedly toward the fast species and thus primes the active site for subsequent ligand binding, ensuring that formation of the ternary complex occurs predominantly by first binding l-Trp and then the ligand. Consequently, the efficiency of catalysis is enhanced because O₂ binding prior to substrate binding, resulting in nonproductive oxidation of the heme iron, is greatly suppressed.

Migration of small ligands in globins: Xe diffusion in truncated hemoglobin N

In heme proteins, the efficient transport of ligands such as NO or O₂ to the binding site is achieved via ligand migration networks. A quantitative assessment of ligand diffusion in these networks is thus essential for a better understanding of the function of these proteins. For this, Xe migration in truncated hemoglobin N (trHbN) of Mycobacterium Tuberculosis was studied using molecular dynamics simulations. Transitions between pockets of the migration network and intra-pocket relaxation occur on similar time scales (10 ps and 20 ps), consistent with low free energy barriers (1-2 kcal/mol). Depending on the pocket from where Xe enters a particular transition, the conformation of the side chains lining the transition region differs which highlights the coupling between ligand and protein degrees of freedom. Furthermore, comparison of transition probabilities shows that Xe migration in trHbN is a non-Markovian process. Memory effects arise due to protein rearrangements and coupled dynamics as Xe moves through it.

Binding and transport of ligands in proteins is essential, in particular in globular proteins which often exhibit internal cavities. In truncated Hemoglobin N (trHbN) these cavities are arranged as a network with particular connectivities. The present work supports the notion that ligand diffusion in trHbN is an active process and coupled to the protein dynamics. Furthermore, transition probabilities between neighboring pockets depend on the location from where the ligand entered the transition, which is typical for non-Markovian processes. Hence, ligand migration in trHbN exhibits memory effects due to dynamical coupling between the protein and ligand motion.

On the indirect relationship between protein dynamics and enzyme activity

The behaviors of simple thermal systems have been well studied in physical chemistry and the principles obtained from such studies have been applied to complex thermal systems, such as proteins and enzymes. But the simple application of such principles is questionable and may lead to mistakes under some circumstances. In enzymology, the transition state theory of chemical reactions has been accepted as a fundamental theory, but the role of protein dynamics in enzyme catalysis is controversial in the context of transition state theory. By studying behaviors of complex thermal systems, we have revised the Arrhenius equation and transition state theory and our model is validated in enzymology. Formally speaking, the revised Arrhenius equation is apparently similar to a conventional Arrhenius equation, but the physical meanings of its parameters differ from that of traditional forms in principle. Within this model, the role of protein dynamics in enzyme catalysis is well defined and quantified.

Elastic Scattering Spectroscopy (ESS): an Instrument-Concept for Dynamics of Complex (Bio-) Systems From Elastic Neutron Scattering

A new type of neutron-scattering spectroscopy is presented that is designed specifically to measure dynamics in bio-systems that are difficult to obtain in any other way. The temporal information is largely model-free and is analogous to relaxation processes measured with dielectric spectroscopy, but provides additional spacial and geometric aspects of the underlying dynamics. Numerical simulations of the basic instrument design show the neutron beam can be highly focussed, giving efficiency gains that enable the use of small samples. Although we concentrate on continuous neutron sources, the extension to pulsed neutron sources is proposed, both requiring minimal data-treatment and being broadly analogous with dielectric spectroscopy, they will open the study of dynamics to new areas of biophysics.

Nanosecond ligand migration and functional protein relaxation in ba3 oxidoreductase: Structures of the B0, B1 and B2 intermediate states

Nanosecond time-resolved step-scan FTIR spectroscopy (nTRS (2) -FTIR) has been applied to literally probe the active site of the carbon monoxide (CO)-bound thermophilic ba3 heme-copper oxidoreductase as it executes its function. The nTRS (2) - snapshots of the photolysed heme a3 Fe-CO/CuB species captured a "transition state" whose side chains prevent the photolysed CO to enter the docking cavity. There are three sets of ba3 photoproduct bands of docked CO with different orientation exhibiting different kinetics. The trajectories of the "docked" CO at 2122, 2129 and 2137cm(-1) is referred to in the literature as B2, B1 and B0 intermediate states, respectively. The present data provided direct evidence for the role of water in controlling ligand orientation in an intracavity protein environment.

Microscopic relaxations in a protein sustained down to 160K in a non-glass forming organic solvent

BACKGROUND

We have studied microscopic dynamics of a protein in carbon disulfide, a non-glass forming solvent, down to its freezing temperature of ca. 160K.

METHODS

We have utilized quasielastic neutron scattering.

RESULTS

A comparison of lysozyme hydrated with water and dissolved in carbon disulfide reveals a stark difference in the temperature dependence of the protein's microscopic relaxation dynamics induced by the solvent. In the case of hydration water, the common protein glass-forming solvent, the protein relaxation slows down in response to a large increase in the water viscosity on cooling down, exhibiting a well-known protein dynamical transition. The dynamical transition disappears in non-glass forming carbon disulfide, whose viscosity remains a weak function of temperature all the way down to freezing at just below 160K. The microscopic relaxation dynamics of lysozyme dissolved in carbon disulfide is sustained down to the freezing temperature of its solvent at a rate similar to that measured at ambient temperature.

CONCLUSIONS

Our results demonstrate that protein dynamical transition is not merely solvent-assisted, but rather solvent-induced, or, more precisely, is a reflection of the temperature dependence of the solvent's glass-forming dynamics.

GENERAL SIGNIFICANCE

We hypothesize that, if the long debated idea regarding the direct link between the microscopic relaxations and the biological activity in proteins is correct, then not only the microscopic relaxations, but also the activity, could be sustained in proteins all the way down to the freezing temperature of a non-glass forming solvent with a weak temperature dependence of its viscosity. This article is part of a Special Issue entitled "Science for Life" Guest Editor: Dr. Austen Angell, Dr. Salvatore Magazù and Dr. Federica Migliardo.

Effect of Outer-Sphere Side Chain Substitutions on the Fate of the trans Iron-Nitrosyl Dimer in Heme/Nonheme Engineered Myoglobins (Fe(B)Mbs): Insights into the Mechanism of Denitrifying NO Reductases

Denitrifying NO reductases are transmembrane protein complexes that utilize a heme/nonheme diiron center at their active sites to reduce two NO molecules to the innocuous gas N2O. Fe(B)Mb proteins, with their nonheme iron sites engineered into the heme distal pocket of sperm whale myoglobin, are attractive models for studying the molecular details of the NO reduction reaction. Spectroscopic and structural studies of Fe(B)Mb constructs have confirmed that they reproduce the metal coordination spheres observed at the active site of the cytochrome c-dependent NO reductase from Pseudomonas aeruginosa. Exposure of Fe(B)Mb to excess NO, as examined by analytical and spectroscopic techniques, results primarily in the formation of a five-coordinate heme-nitrosyl complex without N2O production. However, substitution of the outer-sphere residue Ile107 with a glutamic acid (i.e., I107E) decreases the formation rate of the five-coordinate heme-nitrosyl complex and allows for the substoichiometric production of N2O. Here, we aim to better characterize the formation of the five-coordinate heme-nitrosyl complex and to explain why the level of N2O production increases with the I107E substitution. We follow the formation of the five-coordinate heme-nitrosyl inhibitory complex through the sequential exposure of Fe(B)Mb to different NO isotopomers using rapid-freeze-quench resonance Raman spectroscopy. The data show that the complex is formed by the displacement of the proximal histidine by a new NO molecule after the weakening of the Fe(II)-His bond in the intermediate six-coordinate low-spin (6cLS) heme-nitrosyl complex. These results lead us to explore diatomic migration within the scaffold of myoglobin and whether substitutions at residue 107 can be sufficient to control access to the proximal heme cavities. Results on a new Fe(B)Mb construct with an I107F substitution (Fe(B)Mb3) show an increased rate for the formation of the five-coordinate low-spin heme-nitrosyl complex without N2O production. Taken together, our results suggest that production of N2O from the [6cLS heme {FeNO}(7)/{Fe(B)NO}(7)] trans iron-nitrosyl dimer intermediate requires a proton transfer event facilitated by an outer-sphere residue such as E107 in Fe(B)Mb2 and E280 in P. aeruginosa cNOR.

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.