Relating the Diffusion of Small Ligands in Human Neuroglobin to Its Structural and Mechanical Properties

Neuroglobin Facilitates Neuronal Oxygenation through Tropic Migration under Hypoxia or Anemia in Rat: How Does the Brain Breathe?

The discovery of neuroglobin (Ngb), a brain- or neuron-specific member of the hemoglobin family, has revolutionized our understanding of brain oxygen metabolism. Currently, how Ngb plays such a role remains far from clear. Here, we report a novel mechanism by which Ngb might facilitate neuronal oxygenation upon hypoxia or anemia. We found that Ngb was present in, co-localized to, and co-migrated with mitochondria in the cell body and neurites of neurons. Hypoxia induced a sudden and prominent migration of Ngb towards the cytoplasmic membrane (CM) or cell surface in living neurons, and this was accompanied by the mitochondria. In vivo, hypotonic and anemic hypoxia induced a reversible Ngb migration toward the CM in cerebral cortical neurons in rat brains but did not alter the expression level of Ngb or its cytoplasm/mitochondria ratio. Knock-down of Ngb by RNA interference significantly diminished respiratory succinate dehydrogenase (SDH) and ATPase activity in neuronal N2a cells. Over-expression of Ngb enhanced SDH activity in N2a cells upon hypoxia. Mutation of Ngb at its oxygen-binding site (His64) significantly increased SDH activity and reduced ATPase activity in N2a cells. Taken together, Ngb was physically and functionally linked to mitochondria. In response to an insufficient oxygen supply, Ngb migrated towards the source of oxygen to facilitate neuronal oxygenation. This novel mechanism of neuronal respiration provides new insights into the understanding and treatment of neurological diseases such as stroke and Alzheimer's disease and diseases that cause hypoxia in the brain such as anemia.

Hydrogen peroxide induces heme degradation and protein aggregation in human neuroglobin: roles of the disulfide bridge and hydrogen-bonding in the distal heme cavity

In the present study, human neuroglobin (hNgb) was found to undergo H₂ O₂ -induced breakdown of the heme center at a much slower rate than other globins, namely in the timescale of hours against minutes. We investigated how the rate of the process is affected by the Cys46/Cys55 disulfide bond and the network of non-covalent interactions in the distal heme side involving Tyr44, Lys67, the His64 heme iron axial ligand and the heme propionate-7. The rate is increased by the Tyr44 to Ala and Phe mutations; however the rate is lowered by Lys67 to Ala swapping. The absence of the disulfide bridge slows down the reaction further. Therefore, the disulfide bond-controlled accessibility of the heme site and the residues at position 44 and 67 affect the activation barrier of the reaction. Wild-type and mutated species form β-amyloid aggregates in the presence of H₂ O₂ producing globular structures. Furthermore, the C46A/C55A, Y44A, Y44F and Y44F/C46A/C55A variants yield potentially harmful fibrils. Finally, the nucleation and growth kinetics for the aggregation of the amyloid structures can be successfully described by the Finke-Watzky model.

Structural and (Pseudo-)Enzymatic Properties of Neuroglobin: Its Possible Role in Neuroprotection

Neuroglobin (Ngb), the third member of the globin family, was discovered in human and murine brains in 2000. This monomeric globin is structurally similar to myoglobin (Mb) and hemoglobin (Hb) α and β subunits, but it hosts a bis-histidyl six-coordinated heme-Fe atom. Therefore, the heme-based reactivity of Ngb is modulated by the dissociation of the distal HisE7-heme-Fe bond, which reflects in turn the redox state of the cell. The high Ngb levels (~100–200 μM) present in the retinal ganglion cell layer and in the optic nerve facilitate the O₂ buffer and delivery. In contrast, the very low levels of Ngb (~1 μM) in most tissues and organs support (pseudo-)enzymatic properties including NO/O₂ metabolism, peroxynitrite and free radical scavenging, nitrite, hydroxylamine, hydrogen sulfide reduction, and the nitration of aromatic compounds. Here, structural and (pseudo-)enzymatic properties of Ngb, which are at the root of tissue and organ protection, are reviewed, envisaging a possible role in the protection from neuronal degeneration of the retina and the optic nerve.

Ligand-Binding Calculations with Metadynamics

All-atom molecular dynamics simulations can capture the dynamic degrees of freedom that characterize molecular recognition, the knowledge of which constitutes the cornerstone of rational approaches to drug design and optimization. In particular, enhanced sampling algorithms, such as metadynamics, are powerful tools to dramatically reduce the computational cost required for a mechanistic description of the binding process. Here, we describe the essential details characterizing these simulation strategies, focusing on the critical step of identifying suitable reaction coordinates, as well as on the different analysis algorithms to estimate binding affinity and residence times. We conclude with a survey of published applications that provides explicit examples of successful simulations for several targets.

Mapping Hydrophobic Tunnels and Cavities in Neuroglobin with Noble Gas under Pressure

Internal cavities are crucial for conformational flexibility of proteins and can be mapped through noble gas diffusion and docking. Here we investigate the hydrophobic cavities and tunnel network in neuroglobin (Ngb), a hexacoordinated heme protein likely to be involved in neuroprotection, using crystallography under noble gas pressure, mostly at room temperature. In murine Ngb, a large internal cavity is involved in the heme sliding mechanism to achieve binding of gaseous ligands through coordination to the heme iron. In this study, we report that noble gases are hosted by two major sites within the internal cavity. We propose that these cavities could store oxygen and allow its relay in the heme proximity, which could correspond to NO location in the nitrite-reductase function of Ngb. Thanks to a recently designed pressurization cell using krypton at high pressure, a new gas binding site has been characterized that reveals an alternate pathway for gaseous ligands. A new gas binding site on the proximal side of the heme has also been characterized, using xenon pressure on a Ngb mutant (V140W) that binds CO with a similar rate and affinity to the wild-type, despite a reshaping of the internal cavity. Moreover, this study, to our knowledge, provides new insights into the determinants of the heme sliding mechanism, suggesting that the shift at the beginning of helix G precedes and drives this process.

ID3 regulates the MDC1-mediated DNA damage response in order to maintain genome stability

MDC1 plays a critical role in the DNA damage response (DDR) by interacting directly with several factors including γ-H2AX. However, the mechanism by which MDC1 is recruited to damaged sites remains elusive. Here, we show that MDC1 interacts with a helix–loop–helix (HLH)-containing protein called inhibitor of DNA-binding 3 (ID3). In response to double-strand breaks (DSBs) in the genome, ATM phosphorylates ID3 at serine 65 within the HLH motif, and this modification allows a direct interaction with MDC1. Moreover, depletion of ID3 results in impaired formation of ionizing radiation (IR)-induced MDC1 foci, suppression of γ-H2AX-bound MDC1, impaired DSB repair, cellular hypersensitivity to IR, and genomic instability. Disruption of the MDC1–ID3 interaction prevents accumulation of MDC1 at sites of DSBs and suppresses DSB repair. Thus, our study uncovers an ID3-dependent mechanism of recruitment of MDC1 to DNA damage sites and suggests that the ID3–MDC1 interaction is crucial for DDR.

MDC1 is a key component of the DNA damage response and interacts with several factors such as γ-H2AX. Here the authors show that MDC1 interacts with ID3, facilitating MDC1 recruitment to sites of damage and repair of breaks.

Computational investigation of O2 diffusion through an intra-molecular tunnel in AlkB; influence of polarization on O2 transport

E. Coli AlkB catalyzes the direct dealkylation of various alkylated bases in damaged DNA. The diffusion of molecular oxygen to the active site in AlkB is an essential step for the oxidative dealkylation activity. Despite detailed studies on the stepwise oxidation mechanism of AlkB, there is no conclusive picture of how O2 molecules reach the active site of the protein. Yu et al. (Nature, 439, 879) proposed the existence of an intra-molecular tunnel based on their initial crystal structures of AlkB. We have employed computational simulations to investigate possible migration pathways inside AlkB for O2 molecules. Extensive molecular dynamics (MD) simulations, including explicit ligand sampling and potential of mean force (PMF) calculations, have been performed to provide a microscopic description of the O2 delivery pathway in AlkB. Analysis of intra-molecular tunnels using the CAVER software indicates two possible pathways for O2 to diffuse into the AlkB active site. Explicit ligand sampling simulations suggests that only one of these tunnels provides a viable route. The free energy path for an oxygen molecule to travel along each of these tunnels has been determined with AMBER and AMOEBA. Both PMFs indicate passive transport of O2 from the surface of the protein. However, the inclusion of explicit polarization shows a very large barrier for diffusion of the co-substrate out of the active site, compared with the non-polarizable potential. In addition, our results suggest that the mutation of a conserved residue along the tunnel, Y178, has dramatic effects on the dynamics of AlkB and on the transport of O2 along the tunnel.

Determinants of neuroglobin plasticity highlighted by joint coarse-grained simulations and high pressure crystallography

Investigating the effect of pressure sheds light on the dynamics and plasticity of proteins, intrinsically correlated to functional efficiency. Here we detail the structural response to pressure of neuroglobin (Ngb), a hexacoordinate globin likely to be involved in neuroprotection. In murine Ngb, reversible coordination is achieved by repositioning the heme more deeply into a large internal cavity, the “heme sliding mechanism”. Combining high pressure crystallography and coarse-grain simulations on wild type Ngb as well as two mutants, one (V101F) with unaffected and another (F106W) with decreased affinity for CO, we show that Ngb hinges around a rigid mechanical nucleus of five hydrophobic residues (V68, I72, V109, L113, Y137) during its conformational transition induced by gaseous ligand, that the intrinsic flexibility of the F-G loop appears essential to drive the heme sliding mechanism, and that residue Val 101 may act as a sensor of the interaction disruption between the heme and the distal histidine.

Structural Plasticity in Globins: Role of Protein Dynamics in Defining Ligand Migration Pathways

Globins are a family of proteins characterized by the presence of the heme prosthetic group and involved in variety of biological functions in the cell. Due to their biological relevance and widespread distribution in all kingdoms of life, intense research efforts have been devoted to disclosing the relationships between structural features, protein dynamics, and function. Particular attention has been paid to the impact of differences in amino acid sequence on the topological features of docking sites and cavities and to the influence of conformational flexibility in facilitating the migration of small ligands through these cavities. Often, tunnels are carved in the interior of globins, and ligand exchange is regulated by gating residues. Understanding the subtle intricacies that relate the differences in sequence with the structural and dynamical features of globins with the ultimate aim of rationalizing the thermodynamics and kinetics of ligand binding continues to be a major challenge in the field. Due to the evolution of computational techniques, significant advances into our understanding of these questions have been made. In this review we focus our attention on the analysis of the ligand migration pathways as well as the function of the structural cavities and tunnels in a series of representative globins, emphasizing the synergy between experimental and theoretical approaches to gain a comprehensive knowledge into the molecular mechanisms of this diverse family of proteins.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

(1)H, (15)N and (13)C backbone resonance assignments of murine met-neuroglobin, free and in complex with cyanide

Neuroglobin is a globin present in the brain and retina of mammals. This hexacoordinated hemoprotein binds small diatomic molecules, albeit with lower affinity compared with other globins. We report here the resonance assignment of murine met-Neuroglobine, free and in complex with cyanide.

A possible mechanism for redox control of human neuroglobin activity

Neuroglobin (Ngb) promotes neuron survival under hypoxic/ischemic conditions. In vivo and in vitro assays provide evidence for redox-regulated functioning of Ngb. On the basis of X-ray crystal structures and our MD simulations, a mechanism for redox control of human Ngb (hNgb) activity via the influence of the CD loop on the active site is proposed. We provide evidence that the CD loop undergoes a strand-to-helix transition when the external environment becomes sufficiently oxidizing, and that this CD loop conformational transition causes critical restructuring of the active site. We postulate that the strand-to-helix mechanics of the CD loop allows hNgb to utilize the lability of Cys46/Cys55 disulfide bonding and of the Tyr44/His64/heme propionate interaction network for redox-controlled functioning of hNgb.

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.

Nitrite reduction and cardiovascular protection

Inorganic nitrite, a metabolite of endogenously produced nitric oxide (NO) from NO synthases (NOS), provides the largest endocrine source of directly bioavailable NO. The conversion of nitrite to NO occurs mainly through enzymatic reduction, mediated by a range of proteins, including haem-globins, molybdo-flavoproteins, mitochondrial proteins, cytochrome P450 enzymes, and NOS. Such nitrite reduction is particularly favoured under hypoxia, when endogenous formation of NO from NOS is impaired. Under normoxic conditions, the majority of these nitrite reductases also scavenge NO, or diminish its bioavailability via reactive oxygen species (ROS) production, suggesting an intricate balance. Moreover, nitrite, whether produced endogenously, or derived from exogenous nitrite or nitrate administration (including dietary sources via the Nitrate-Nitrite-NO pathway) beneficially modulates many key cardiovascular pathological processes. In this review, we highlight the landmark studies which revealed nitrite's function in biological systems, and inspect its evolving role in cardiovascular protection. Whilst these effects have mainly been ascribed to the activity of one or more nitrite reductases, we also discuss newly-identified mechanisms, including nitrite anhydration, the involvement of s-nitrosothiols, nitro-fatty acids, and direct nitrite normoxic signalling, involving modification of mitochondrial structure and function, and ROS production. This article is part of a Special Issue entitled "Redox Signalling in the Cardiovascular System".

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.

Small ligand-globin interactions: reviewing lessons derived from computer simulation

In this work we review the application of classical and quantum-mechanical atomistic computer simulation tools to the investigation of small ligand interaction with globins. In the first part, studies of ligand migration, with its connection to kinetic association rate constants (kon), are presented. In the second part, we review studies for a variety of ligands such as O2, NO, CO, HS(-), F(-), and NO2(-) showing how the heme structure, proximal effects, and the interactions with the distal amino acids can modulate protein ligand binding. The review presents mainly results derived from our previous works on the subject, in the context of other theoretical and experimental studies performed by others. The variety and extent of the presented data yield a clear example of how computer simulation tools have, in the last decade, contributed to our deeper understanding of small ligand interactions with globins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Molecular dynamics study of hell's gate globin I (HGbI) from a methanotrophic extremophile: oxygen migration through a large cavity

Hell's gate globin I (HGbI), a heme-containing protein from the extremophile Methylacidiphilum infernorum, has fast oxygen-binding/slow release characteristics due to its distal residues Gln and Tyr. The combination of Gln/Tyr distal iron coordination, adaptation to extreme environmental conditions, and lack of a D helix suggests that ligand migration in HGbI differs from other previously studied globins. Locally enhanced molecular dynamics trajectories of oxygen migration indicate a large internal cavity. This may increase the tendency of oxygen to exit from portals other than the most direct exit from the space near the heme. Oxygen may reside transiently in shallow surface depressions around the exits. Such surface trapping may enhance both oxygen uptake by increasing contact time between molecules, and decrease release by increasing the probability of oxygen reentry from the vicinity of the portal.

Conformational dynamics in human neuroglobin: effect of His64, Val68, and Cys120 on ligand migration

Neuroglobin belongs to the family of hexacoordinate hemoglobins and has been implicated in the protection of neuronal tissue under hypoxic and ischemic conditions. Here we present transient absorption and photoacoustic calorimetry studies of CO photodissociation and bimolecular rebinding to neuroglobin focusing on the ligand migration process and the role of distal pocket residues (His64 and Val68) and two Cys residues (Cys55 and Cys120). Our results indicate that His64 has a minor impact on the migration of CO between the distal heme pocket and protein exterior, whereas the Val68 side chain regulates the transition of the photodissociated ligand between the distal pocket and internal hydrophobic cavities, which is evident from the increased geminate quantum yield in this mutated protein (Φ(gem) = 0.32 for WT and His64Gln, and Φ(gem) = 0.85 for Val68Phe). The interface between helix G and the A-B loop provides an escape pathway for the photodissociated ligand, which is evident from a decrease in the reaction enthalpy for the transition between the CO-bound hNgb and five-coordinate hNgb in the Cys120Ser mutant (ΔH = -3 ± 4 kcal mol(-1)) compared to that of the WT protein (ΔH = 20 ± 4 kcal mol(-1)). The extensive electrostatic/hydrogen binding network that includes heme propionate groups, Lys67, His64, and Tyr44 not only restricts the heme binding but also modulates the energetics of binding of CO to the five-coordinate hNgb as substitution of His64 with Gln leads to an endothermic association of CO with the five-coordinate hNgb (ΔH = 6 ± 3 kcal mol(-1)).

Biophysical characterisation of neuroglobin of the icefish, a natural knockout for hemoglobin and myoglobin. Comparison with human neuroglobin

The Antarctic icefish Chaenocephalus aceratus lacks the globins common to most vertebrates, hemoglobin and myoglobin, but has retained neuroglobin in the brain. This conserved globin has been cloned, over-expressed and purified. To highlight similarities and differences, the structural features of the neuroglobin of this colourless-blooded fish were compared with those of the well characterised human neuroglobin as well as with the neuroglobin from the retina of the red blooded, hemoglobin and myoglobin-containing, closely related Antarctic notothenioid Dissostichus mawsoni. A detailed structural and functional analysis of the two Antarctic fish neuroglobins was carried out by UV-visible and Resonance Raman spectroscopies, molecular dynamics simulations and laser-flash photolysis. Similar to the human protein, Antarctic fish neuroglobins can reversibly bind oxygen and CO in the Fe(2+) form, and show six-coordination by distal His in the absence of exogenous ligands. A very large and structured internal cavity, with discrete docking sites, was identified in the modelled three-dimensional structures of the Antarctic neuroglobins. Estimate of the free-energy barriers from laser-flash photolysis and Implicit Ligand Sampling showed that the cavities are accessible from the solvent in both proteins.Comparison of structural and functional properties suggests that the two Antarctic fish neuroglobins most likely preserved and possibly improved the function recently proposed for human neuroglobin in ligand multichemistry. Despite subtle differences, the adaptation of Antarctic fish neuroglobins does not seem to parallel the dramatic adaptation of the oxygen carrying globins, hemoglobin and myoglobin, in the same organisms.

Heme orientation modulates histidine dissociation and ligand binding kinetics in the hexacoordinated human neuroglobin

Neuroglobin (Ngb) is a globin present in the brain and retina of mammals. This hexacoordinated hemoprotein binds small diatomic molecules, albeit with lower affinity compared with other globins. Another distinctive feature of most mammalian Ngb is their ability to form an internal disulfide bridge that increases ligand affinity. As often seen for prosthetic heme b containing proteins, human Ngb exhibits heme heterogeneity with two alternative heme orientations within the heme pocket. To date, no details are available on the impact of heme orientation on the binding properties of human Ngb and its interplay with the cysteine oxidation state. In this work, we used ¹H NMR spectroscopy to probe the cyanide binding properties of different Ngb species in solution, including wild-type Ngb and the single (C120S) and triple (C46G/C55S/C120S) mutants. We demonstrate that in the disulfide-containing wild-type protein cyanide ligation is fivefold faster for one of the two heme orientations (the A isomer) compared with the other isomer, which is attributed to the lower stability of the distal His64–iron bond and reduced steric hindrance at the bottom of the cavity for heme sliding in the A conformer. We also attribute the slower cyanide reactivity in the absence of a disulfide bridge to the tighter histidine–iron bond. More generally, enhanced internal mobility in the CD loop bearing the disulfide bridge hinders access of the ligand to heme iron by stabilizing the histidine–iron bond. The functional impact of heme disorder and cysteine oxidation state on the properties of the Ngb ligand is discussed.

Carbon Dioxide Migration Pathways in Proteins

Some of the most important biological processes, such as carbon fixation, are dependent on protein-gas interactions. The motion of CO2 through the enzyme phosphoenolpyruvate carboxykinase was investigated using extensive all-atom molecular dynamics simulations. Three discrete migration pathways were located, suggesting the protein directs the movement of CO2. The chemical nature of these pathways is discussed, as are their biotechnological ramifications.

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.

Molecular motions in drug design: the coming age of the metadynamics method

Metadynamics is emerging as a useful free energy method in physics, chemistry and biology. Recently, it has been applied also to investigate ligand binding to biomolecules of pharmacological interest. Here, after introducing the basic idea of the method, we review applications to challenging targets for pharmaceutical intervention. We show that this methodology, especially when combined with a variety of other computational approaches such as molecular docking and/or molecular dynamics simulation, may be useful to predict structure and energetics of ligand/target complexes even when the targets lack a deep binding cavity, such as DNA and proteins undergoing fibrillation in neurodegenerative diseases. Furthermore, the method allows investigating the routes of molecular recognition and the associated binding energy profiles, providing a molecular interpretation to experimental data.

Time resolved thermodynamics associated with ligand photorelease in heme peroxidases and globins: Open access channels versus gated ligand release

Heme proteins represent a diverse class of biomolecules responsible for an extremely diverse array of physiological functions including electron transport, monooxygenation, ligand transport and storage, cellular signaling, respiration, etc. An intriguing aspect of these proteins is that such functional diversity is accomplished using a single type of heme macrocycle based upon iron protoporphyrin IX. The functional diversity originates from a delicate balance of inter-molecular interactions within the protein matrix together with well choreographed dynamics that modulate the heme electronic structure as well as ligand entry/exit pathways from the bulk solvent to the active site. Of particular interest are the dynamics and energetics associated with the entry/exit of ligands as this process plays a significant role in regulating the rates of heme protein activity. Time-resolved photoacoustic calorimetry (PAC) has emerged as a powerful tool through which to probe the underlying energetics associated with small molecule dissociation and release to the bulk solvent in heme proteins on time scales from tens of nanoseconds to several microseconds. In this review, the results of PAC studies on various classes of heme proteins are summarized highlighting how different protein structures affect the thermodynamics of ligand migration. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Protein dynamics and ligand migration interplay as studied by computer simulation

Since proteins are dynamic systems in living organisms, the employment of methodologies contemplating this crucial characteristic results fundamental to allow revealing several aspects of their function. In this work, we present results obtained using classical mechanical atomistic simulation tools applied to understand the connection between protein dynamics and ligand migration. Firstly, we will present a review of the different sampling schemes used in the last years to obtain both ligand migration pathways and the thermodynamic information associated with the process. Secondly, we will focus on representative examples in which the schemes previously presented are employed, concerning the following: i) ligand migration, tunnels, and cavities in myoglobin and neuroglobin; ii) ligand migration in truncated hemoglobin members; iii) NO escape and conformational changes in nitrophorins; iv) ligand selectivity in catalase and hydrogenase; and v) larger ligand migration: the P450 and haloalkane dehalogenase cases. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Probing the role of the internal disulfide bond in regulating conformational dynamics in neuroglobin

In this report, we demonstrate that the internal disulfide bridge in human neuroglobin modulates structural changes associated with ligand photo-dissociation from the heme active site. This is evident from time-resolved photothermal studies of CO photo-dissociation, which reveal a 13.4+/-0.9 mL mol(-1) volume expansion upon ligand photo-release from human neuroglobin, whereas the CO dissociation from rat neuroglobin leads to a significantly smaller volume change (DeltaV=4.6+/-0.3 mL mol(-1)). Reduction of the internal disulfide bond in human neuroglobin leads to conformational changes (reflected by DeltaV) nearly identical to those observed for rat Ngb. Our data favor the hypothesis that the disulfide bond between Cys46 and Cys55 modulates the functioning of human neuroglobin.