Myoglobin: the hydrogen atom of biology and a paradigm of complexity

Hidden water's influence on rhodopsin activation

Structural biology relies on several powerful techniques, but these tend to be limited in their ability to characterize protein fluctuations and mobility. Overreliance on structural approaches can lead to omission of critical information regarding biological function. Currently there is a need for complementary biophysical methods to visualize these mobile aspects of protein function. Here, we review hydrostatic and osmotic pressure-based techniques to address this shortcoming for the paradigm of rhodopsin. Hydrostatic and osmotic pressure data contribute important examples, which are interpreted in terms of an energy landscape for hydration-mediated protein dynamics. We find that perturbations of rhodopsin conformational equilibria by force-based methods are not unrelated phenomena; rather they probe various hydration states involving functional proton reactions. Hydrostatic pressure acts on small numbers of strongly interacting structural or solvent-shell water molecules with relatively high energies, while osmotic pressure acts on large numbers of weakly interacting bulk-like water molecules with low energies. Local solvent fluctuations due to the hydration shell and collective water interactions affect hydrogen-bonded networks and domain motions that are explained by a hierarchical energy landscape model for protein dynamics.

Osmotic stress studies of G-protein-coupled receptor rhodopsin activation

We summarize and critically review osmotic stress studies of the G-protein-coupled receptor rhodopsin. Although small amounts of structural water are present in these receptors, the effect of bulk water on their function remains uncertain. Studies of the influences of osmotic stress on the GPCR archetype rhodopsin have given insights into the functional role of water in receptor activation. Experimental work has discovered that osmolytes shift the metarhodopsin equilibrium after photoactivation, either to the active or inactive conformations according to their molar mass. At least 80 water molecules are found to enter rhodopsin in the transition to the photoreceptor active state. We infer that this movement of water is both necessary and sufficient for receptor activation. If the water influx is prevented, e.g., by large polymer osmolytes or by dehydration, then the receptor functional transition is back shifted. These findings imply a new paradigm in which rhodopsin becomes solvent swollen in the activation mechanism. Water thus acts as an allosteric modulator of function for rhodopsin-like receptors in lipid membranes.

Nitrobindin versus myoglobin: A comparative structural and functional study

Most hemoproteins display an all-α-helical fold, showing the classical three on three (3/3) globin structural arrangement characterized by seven or eight α-helical segments that form a sandwich around the heme. Over the last decade, a completely distinct class of heme-proteins called nitrobindins (Nbs), which display an all-β-barrel fold, has been identified and characterized from both structural and functional perspectives. Nbs are ten-stranded anti-parallel all-β-barrel heme-proteins found across the evolutionary ladder, from bacteria to Homo sapiens. Myoglobin (Mb), commonly regarded as the prototype of monomeric all-α-helical globins, is involved along with the oligomeric hemoglobin (Hb) in diatomic gas transport, storage, and sensing, as well as in the detoxification of reactive nitrogen and oxygen species. On the other hand, the function(s) of Nbs is still obscure, even though it has been postulated that they might participate to O2/NO signaling and metabolism. This function might be of the utmost importance in poorly oxygenated tissues, such as the eye's retina, where a delicate balance between oxygenation and blood flow (regulated by NO) is crucial. Dysfunction in this balance is associated with several pathological conditions, such as glaucoma and diabetic retinopathy. Here a detailed comparison of the structural, spectroscopic, and functional properties of Mb and Nbs is reported to shed light on the similarities and differences between all-α-helical and all-β-barrel heme-proteins.

Detection of a Geminate Photoproduct of Bovine Cytochrome c Oxidase by Time-Resolved Serial Femtosecond Crystallography

Cytochrome c oxidase (CcO) is a large membrane-bound hemeprotein that catalyzes the reduction of dioxygen to water. Unlike classical dioxygen binding hemeproteins with a heme b group in their active sites, CcO has a unique binuclear center (BNC) composed of a copper atom (CuB) and a heme a3 iron, where O2 binds and is reduced to water. CO is a versatile O2 surrogate in ligand binding and escape reactions. Previous time-resolved spectroscopic studies of the CO complexes of bovine CcO (bCcO) revealed that photolyzing CO from the heme a3 iron leads to a metastable intermediate (CuB-CO), where CO is bound to CuB, before it escapes out of the BNC. Here, with a pump-probe based time-resolved serial femtosecond X-ray crystallography, we detected a geminate photoproduct of the bCcO-CO complex, where CO is dissociated from the heme a3 iron and moved to a temporary binding site midway between the CuB and the heme a3 iron, while the locations of the two metal centers and the conformation of Helix-X, housing the proximal histidine ligand of the heme a3 iron, remain in the CO complex state. This new structure, combined with other reported structures of bCcO, allows for a clearer definition of the ligand dissociation trajectory as well as the associated protein dynamics.

Globular Proteins and Where to Find Them within a Polymer Brush-A Case Study

Protein adsorption by polymerized surfaces is an interdisciplinary topic that has been approached in many ways, leading to a plethora of theoretical, numerical and experimental insight. There is a wide variety of models trying to accurately capture the essence of adsorption and its effect on the conformations of proteins and polymers. However, atomistic simulations are case-specific and computationally demanding. Here, we explore universal aspects of the dynamics of protein adsorption through a coarse-grained (CG) model, that allows us to explore the effects of various design parameters. To this end, we adopt the hydrophobic-polar (HP) model for proteins, place them uniformly at the upper bound of a CG polymer brush whose multibead-spring chains are tethered to a solid implicit wall. We find that the most crucial factor affecting the adsorption efficiency appears to be the polymer grafting density, while the size of the protein and its hydrophobicity ratio come also into play. We discuss the roles of ligands and attractive tethering surfaces to the primary adsorption as well as secondary and ternary adsorption in the presence of attractive (towards the hydrophilic part of the protein) beads along varying spots of the backbone of the polymer chains. The percentage and rate of adsorption, density profiles and the shapes of the proteins, alongside with the respective potential of mean force are recorded to compare the various scenarios during protein adsorption.

Detection of a geminate photoproduct of bovine cytochrome c oxidase by time-resolved serial femtosecond crystallography

Cytochrome c oxidase (C c O) is a large membrane-bound hemeprotein that catalyzes the reduction of dioxygen to water. Unlike classical dioxygen binding hemeproteins with a heme b group in their active sites, C c O has a unique binuclear center (BNC) comprised of a copper atom (Cu B ) and a heme a 3 iron, where O 2 binds and is reduced to water. CO is a versatile O 2 surrogate in ligand binding and escape reactions. Previous time-resolved spectroscopic studies of the CO complexes of bovine C c O (bC c O) revealed that photolyzing CO from the heme a 3 iron leads to a metastable intermediate (Cu B -CO), where CO is bound to Cu B , before it escapes out of the BNC. Here, with a time-resolved serial femtosecond X-ray crystallography-based pump-probe method, we detected a geminate photoproduct of the bC c O-CO complex, where CO is dissociated from the heme a 3 iron and moved to a temporary binding site midway between the Cu B and the heme a 3 iron, while the locations of the two metal centers and the conformation of the Helix-X, housing the proximal histidine ligand of the heme a 3 iron, remain in the CO complex state. This new structure, combined with other reported structures of bC c O, allows the full definition of the ligand dissociation trajectory, as well as the associated protein dynamics.

Ligand-Based Regulation of Dynamics and Reactivity of Hemoproteins

Hemoproteins include several heme-binding proteins with distinct structure and function. The presence of the heme group confers specific reactivity and spectroscopic properties to hemoproteins. In this review, we provide an overview of five families of hemoproteins in terms of dynamics and reactivity. First, we describe how ligands modulate cooperativity and reactivity in globins, such as myoglobin and hemoglobin. Second, we move on to another family of hemoproteins devoted to electron transport, such as cytochromes. Later, we consider heme-based reactivity in hemopexin, the main heme-scavenging protein. Then, we focus on heme-albumin, a chronosteric hemoprotein with peculiar spectroscopic and enzymatic properties. Eventually, we analyze the reactivity and dynamics of the most recently discovered family of hemoproteins, i.e., nitrobindins.

Perspective on the potential of tandem-ion mobility /mass spectrometry methods for structural proteomics applications

Cellular processes are usually carried out collectively by the entirety of all proteins present in a biological cell, i.e. the proteome. Mass spectrometry-based methods have proven particularly successful in identifying and quantifying the constituent proteins of proteomes, including different molecular forms of a protein. Nevertheless, protein sequences alone do not reveal the function or dysfunction of the identified proteins. A straightforward way to assign function or dysfunction to proteins is characterization of their structures and dynamics. However, a method capable to characterize detailed structures of proteins and protein complexes in a large-scale, systematic manner within the context of cellular processes does not yet exist. Here, we discuss the potential of tandem-ion mobility / mass spectrometry (tandem-IM/MS) methods to provide such ability. We highlight the capability of these methods using two case studies on the protein systems ubiquitin and avidin using the tandem-TIMS/MS technology developed in our laboratory and discuss these results in the context of other developments in the broader field of tandem-IM/MS.

Some further classes of pseudo-differential operators in the p-adic context and their applications

The purpose of this article is to study new non-Archimedean pseudo-differential operators whose symbols are determined from the behavior of two functions defined on the p-adic numbers. Thanks to the characteristics of our symbols, we can find connections between these operators and new types of non-homogeneous differential equations, Feller semigroups, contraction semigroups and strong Markov processes.

Origin of Protein Quake: Energy Waves Conducted by a Precise Mechanical Machine

A long-standing challenge in protein biophysics is to understand protein quake in myoglobin─the structural dynamics responsible for redistributing the excess heme energy after photolysis. Despite extensive efforts, the molecular mechanism of this process remains elusive. Using the energy flow theory, we uncovered a fundamental new phenomenon: the heme energy is redistributed by sinusoidal waves with a ubiquitous fundamental frequency and two overtones. The energy waves emanate from the heme into the myoglobin backbone via a conduit of five consecutive dihedrals of the proximal histidine and then travel quickly along the backbone to reach sidechains across the protein. This mechanism is far more effective than the diffusion-based mechanism from previous studies because waves are systematic while diffusion is random. To propagate energy waves, coordinates must cooperate, resulting in collective modes that are singular vectors of the generalized work functional. These modes show task partitioning: a handful of high-energy modes generate large-scale breathing motion, which loosens up the protein matrix to enable hundreds of low-energy vibrational modes for energy transduction.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Stationary and Time-Dependent Carbon Monoxide Stretching Mode Features in Carboxy Myoglobin: A Theoretical-Computational Reappraisal

The stationary and time-dependent infrared spectrum (IR) of the CO stretching mode (ν_CO) in carboxymyoglobin (MbCO), a longstanding problem of biophysical chemistry, has been modeled through a theoretical-computational method specifically designed for simulating quantum observables in complex atomic-molecular systems and based on a combined application of long time scale molecular dynamics simulations and quantum-chemical calculations. This study is basically focused on two aspects: (i) the origin of the stationary IR substates (termed as A₀, A₁, and A₃) and (ii) the modeling and the interpretation of the ν_CO energy relaxation. The results, strengthened by a more than satisfactory agreement with the experimental data, concisely indicate that (i) the conformational His64-FeCO relevant substates, i.e., characterized by the formation-disruption of the H-bond between the above moieties, are the main responsible of the presence of two distinct and well separated (A₀ and A₁/A₃) spectroscopic regions; (ii) the characteristic bimodal shape of the A₁/A₃ spectral region, according to our model, is the result of the fluctuation of the electric field pattern as provided by the protein-solvent framework perturbing the bound His64-CO-Heme complex; and (iii) the electric field pattern, in conjunction with the relatively high density of MbCO vibrational states, is also the main determinant of the ν_CO energy relaxation, characterizing its kinetic efficiency.

Structural dynamics in the evolution of a bilobed protein scaffold

Proteins conduct numerous complex biological functions by use of tailored structural dynamics. The molecular details of how these emerged from ancestral peptides remains mysterious. How does nature utilize the same repertoire of folds to diversify function? To shed light on this, we analyzed bilobed proteins with a common structural core, which is spread throughout the tree of life and is involved in diverse biological functions such as transcription, enzymatic catalysis, membrane transport, and signaling. We show here that the structural dynamics of the structural core differentiate predominantly via terminal additions during a long-period evolution. This diversifies substrate specificity and, ultimately, biological function.

Novel biophysical tools allow the structural dynamics of proteins and the regulation of such dynamics by binding partners to be explored in unprecedented detail. Although this has provided critical insights into protein function, the means by which structural dynamics direct protein evolution remain poorly understood. Here, we investigated how proteins with a bilobed structure, composed of two related domains from the periplasmic-binding protein–like II domain family, have undergone divergent evolution, leading to adaptation of their structural dynamics. We performed a structural analysis on ∼600 bilobed proteins with a common primordial structural core, which we complemented with biophysical studies to explore the structural dynamics of selected examples by single-molecule Förster resonance energy transfer and Hydrogen–Deuterium exchange mass spectrometry. We show that evolutionary modifications of the structural core, largely at its termini, enable distinct structural dynamics, allowing the diversification of these proteins into transcription factors, enzymes, and extracytoplasmic transport-related proteins. Structural embellishments of the core created interdomain interactions that stabilized structural states, reshaping the active site geometry, and ultimately altered substrate specificity. Our findings reveal an as-yet-unrecognized mechanism for the emergence of functional promiscuity during long periods of evolution and are applicable to a large number of domain architectures.

Kinetic mechanisms for O2 binding to myoglobins and hemoglobins

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

Understanding the close encounter of heme proteins with carboxylated multiwalled carbon nanotubes: a case study of contradictory stability trend for hemoglobin and myoglobin

Carbon nanotubes (CNTs) are one of the unique and promising nanomaterials that possess plenty of applications, such as biosensors, advanced drug delivery systems and biotechnology. CNTs bind rapidly with proteins, which result in the formation of a protein coating layer known as a "protein corona" around the surface of the nanomaterial. This hinders their applications as a drug carrier and influences the properties of biological macromolecules. The present work focuses on studying the thermal stability and molecular level interactions of two heme proteins, hemoglobin (Hb) and myoglobin (Mb), in the presence of carboxylated functionalized multi-walled CNTs (CA-MWCNTs). Through the current study, the following steps have been taken to distinguish the biocompatibility of the hydrophilic surface CA-MWCNTs for heme proteins via a series of spectroscopic techniques and differential scanning calorimetry (DSC). UV-Visible and steady-state fluorescence spectroscopy were used to reveal changes in the aromatic amino acid residues of heme proteins upon the addition of CA-MWCNTs. Circular dichroism spectroscopy (CD) shows the alteration in the native structure of proteins in the presence of the nanomaterial. A tremendous increase in the size of the protein CA-MWCNTs system is observed in dynamic light scattering (DLS), which clearly manifests the protein corona formation. Unexpectedly, both proteins interact differently with CA-MWCNTs, which is observed in CD spectroscopy and DSC. In the presence of CA-MWCNTs, an increase in the transition temperature (T_m) was observed for Hb, while the T_m value decreases for Mb. Different interactions with proteins at the molecular scale may be the reason for this unexpected behavior. Henceforth, the present results can help in the design of the next-generation drug carrier nanomaterials with the idea of the heme protein corona formation prior to development.

Cellular NADH and NADPH Conformation as a Real-Time Fluorescence-Based Metabolic Indicator under Pressurized Conditions

Cellular conformation of reduced pyridine nucleotides NADH and NADPH sensed using autofluorescence spectroscopy is presented as a real-time metabolic indicator under pressurized conditions. The approach provides information on the role of pressure in energy metabolism and antioxidant defense with applications in agriculture and food technologies. Here, we use spectral phasor analysis on UV-excited autofluorescence from Saccharomyces cerevisiae (baker's yeast) to assess the involvement of one or multiple NADH- or NADPH-linked pathways based on the presence of two-component spectral behavior during a metabolic response. To demonstrate metabolic monitoring under pressure, we first present the autofluorescence response to cyanide (a respiratory inhibitor) at 32 MPa. Although ambient and high-pressure responses remain similar, pressure itself also induces a response that is consistent with a change in cellular redox state and ROS production. Next, as an example of an autofluorescence response altered by pressurization, we investigate the response to ethanol at ambient, 12 MPa, and 30 MPa pressure. Ethanol (another respiratory inhibitor) and cyanide induce similar responses at ambient pressure. The onset of non-two-component spectral behavior upon pressurization suggests a change in the mechanism of ethanol action. Overall, results point to new avenues of investigation in piezophysiology by providing a way of visualizing metabolism and mitochondrial function under pressurized conditions.

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.

Interplays of electron and nuclear motions along CO dissociation trajectory in myoglobin revealed by ultrafast X-rays and quantum dynamics calculations

Ultrafast structural dynamics with different spatial and temporal scales were investigated during photodissociation of carbon monoxide (CO) from iron(II)-heme in bovine myoglobin during the first 3 ps following laser excitation. We used simultaneous X-ray transient absorption (XTA) spectroscopy and X-ray transient solution scattering (XSS) at an X-ray free electron laser source with a time resolution of 80 fs. Kinetic traces at different characteristic X-ray energies were collected to give a global picture of the multistep pathway in the photodissociation of CO from heme. In order to extract the reaction coordinates along different directions of the CO departure, XTA data were collected with parallel and perpendicular relative polarizations of the laser pump and X-ray probe pulse to isolate the contributions of electronic spin state transition, bond breaking, and heme macrocycle nuclear relaxation. The time evolution of the iron K-edge X-ray absorption near edge structure (XANES) features along the two major photochemical reaction coordinates, i.e., the iron(II)-CO bond elongation and the heme macrocycle doming relaxation were modeled by time-dependent density functional theory calculations. Combined results from the experiments and computations reveal insight into interplays between the nuclear and electronic structural dynamics along the CO photodissociation trajectory. Time-resolved small-angle X-ray scattering data during the same process are also simultaneously collected, which show that the local CO dissociation causes a protein quake propagating on different spatial and temporal scales. These studies are important for understanding gas transport and protein deligation processes and shed light on the interplay of active site conformational changes and large-scale protein reorganization.

Mycobacterial and Human Ferrous Nitrobindins: Spectroscopic and Reactivity Properties

Structural and functional properties of ferrous Mycobacterium tuberculosis (Mt-Nb) and human (Hs-Nb) nitrobindins (Nbs) were investigated. At pH 7.0 and 25.0 °C, the unliganded Fe(II) species is penta-coordinated and unlike most other hemoproteins no pH-dependence of its coordination was detected over the pH range between 2.2 and 7.0. Further, despite a very open distal side of the heme pocket (as also indicated by the vanishingly small geminate recombination of CO for both Nbs), which exposes the heme pocket to the bulk solvent, their reactivity toward ligands, such as CO and NO, is significantly slower than in most hemoproteins, envisaging either a proximal barrier for ligand binding and/or crowding of H₂O molecules in the distal side of the heme pocket which impairs ligand binding to the heme Fe-atom. On the other hand, liganded species display already at pH 7.0 and 25 °C a severe weakening (in the case of CO) and a cleavage (in the case of NO) of the proximal Fe-His bond, suggesting that the ligand-linked movement of the Fe(II) atom onto the heme plane brings about a marked lengthening of the proximal Fe-imidazole bond, eventually leading to its rupture. This structural evidence is accompanied by a marked enhancement of both ligands dissociation rate constants. As a whole, these data clearly indicate that structural–functional relationships in Nbs strongly differ from what observed in mammalian and truncated hemoproteins, suggesting that Nbs play a functional role clearly distinct from other eukaryotic and prokaryotic hemoproteins.

Control of Oxygen Affinity in Mammalian Hemoglobins: Implications for a System Biology Description of the Respiratory Properties of the Red Blood Cell

Hemoglobin and myoglobin have been considered for a long time the paradigmatic model systems for protein function, to the point of being defined the "hydrogen atom[s] of biology". Given this privileged position and the huge amount of quantitative information available on these proteins, the red blood cell might appear as the model system and"hydrogen atom" of system biology. Indeed, since the red cell's main function is O2 transport by hemoglobin, the gap between the protein and the cell may appear quite small. Yet, a surprisingly large amount of detailed biochemical information is required for the modelization of the respiratory properties of the erythrocyte. This problem is compounded if modelization aims at uncovering or explaining evolutionarily selected functional properties of hemoglobin. The foremost difficulty lies in the fact that hemoglobins having different intrinsic properties and relatively ancient evolutionary divergence may behave similarly in the complex milieu of blood, whereas very similar hemoglobins sharing a substantial sequence similarity may present important functional differences because of the mutation of a few key residues. Thus, the functional properties of hemoglobin and blood may reflect more closely the recent environmental challenges than the remote evolutionary history of the animal. We summarize in this review the case of hemoglobins from mammals, in an attempt to provide a reasoned summary of their complexity that, we hope, may be of help to scientists interested in the quantitative exploration of the evolutionary physiology of respiration. Indeed the basis of a meaningful modelization of the red cell requires a large amount of information collected in painstaking and often forgotten studies of the biochemical properties of hemoglobin carried out over more than a century.

p-Adic mathematics and theoretical biology

The principal objective of this article is a brief overview of the main parts of p-adic mathematics, which have already had valuable applications and may have a significant impact in the near future on the further development of some fields of theoretical and mathematical biology. In particular, we present the basics of ultrametrics, p-adic numbers and p-adic analysis, as well as insight into their applications for modeling some cognitive processes, genetic code and protein dynamics. We also argue that ultrametric concepts and p-adic mathematics are natural tools for the viable description of biological systems and phenomena with a hierarchical structure.

Dioxygen Binding and Sensing Proteins

Oxygen binding proteins (O₂BIP) have been actively investigated for the past five decades due to their rich redox chemistry and function as O₂ carriers in blood cells, as well as their function as gasotransmitters and sensors that modulate cellular signaling. A series of meetings on the periodic advances in the knowledge gained in the field of globin structure and function are conducted typically on a biannual basis. In the fall of 2018, the XXth International Conference was conducted, and very important articles with breakthrough discoveries were presented and very enthusiastically discussed. This was yet another highly successful meeting in the series. Select articles from this meeting were recently reviewed, updated, and published over several issues of Antioxidants and Redox Signaling, as Forum articles communicating the latest advances in this important area of redox biology. This Forum editorial introduces these articles and highlights their scientific significance in advancing the field. Each of these articles grew out of lectures presented in the meeting, and appears either as an original contribution or a comprehensive review in the journal. Overall, the articles published in the Forum provide in-depth details on the recent developments in the field as well as point the way to future directions. These Forum articles thus serve as an important summary of progress and the ongoing direction of this field, and serve to highlight recent advances in our understanding of O₂BIP.

Human nitrobindin: the first example of an all-β-barrel ferric heme-protein that catalyzes peroxynitrite detoxification

Nitrobindins (Nbs), constituting a heme‐protein family spanning from bacteria to Homo sapiens, display an all‐β‐barrel structural organization. Human Nb has been described as a domain of the nuclear protein named THAP4, whose physiological function is still unknown. We report the first evidence of the heme‐Fe(III)‐based detoxification of peroxynitrite by the all‐β‐barrel C‐terminal Nb‐like domain of THAP4. Ferric human Nb (Nb(III)) catalyzes the conversion of peroxynitrite to NO3− and impairs the nitration of free l‐tyrosine. The rate of human Nb(III)‐mediated scavenging of peroxynitrite is similar to those of all‐α‐helical horse heart and sperm whale myoglobin and human hemoglobin, generally taken as the prototypes of all‐α‐helical heme‐proteins. The heme‐Fe(III) reactivity of all‐β‐barrel human Nb(III) and all‐α‐helical prototypical heme‐proteins possibly reflects the out‐to‐in‐plane transition of the heme‐Fe(III)‐atom preceding peroxynitrite binding. Human Nb(III) not only catalyzes the detoxification of peroxynitrite but also binds NO, possibly representing a target of reactive nitrogen species.

Dynamic multiple-scattering treatment of X-ray absorption: Parameterization of a new molecular dynamics force field for myoglobin

We present a detailed analysis of the X-ray absorption near-edge structure (XANES) data on the Fe K-edge of CO Myoglobin based on a combined procedure of Molecular Dynamics (MD) calculations and MXAN (Minuit XANes) data analysis that we call D-MXAN. The ability of performing quantitative XANES data analysis allows us to refine classical force field MD parameters, thus obtaining a reliable tool for the atomic investigation of this important model system for biological macromolecules. The iterative procedure here applied corrects the greatest part of the structural discrepancy between classical MD sampling and experimental determinations. Our procedure, moreover, is able to discriminate between different heme conformational basins visited during the MD simulation, thus demonstrating the necessity of a sampling on the order of tens of nanoseconds, even for an application such X-ray absorption spectroscopy data analysis.

Nitrophorins and nitrobindins: structure and function

Classical all α-helical globins are present in all living organisms and are ordered in three lineages: (i) flavohemoglobins and single domain globins, (ii) protoglobins and globin coupled sensors and (iii) truncated hemoglobins, displaying the 3/3 or the 2/2 all α-helical fold. However, over the last two decades, all β-barrel and mixed α-helical-β-barrel heme-proteins displaying heme-based functional properties (e.g. ligand binding, transport and sensing) closely similar to those of all α-helical globins have been reported. Monomeric nitrophorins (NPs) and α1-microglobulin (α1-m), belonging to the lipocalin superfamily and nitrobindins (Nbs) represent prototypical heme-proteins displaying the all β-barrel and mixed α-helical-β-barrel folds. NPs are confined to the Reduviidae and Cimicidae families of Heteroptera, whereas α1-m and Nbs constitute heme-protein families spanning bacteria to Homo sapiens. The structural organization and the reactivity of the stable ferric solvent-exposed heme-Fe atom suggest that NPs and Nbs are devoted to NO transport, storage and sensing, whereas Hs-α1-m participates in heme metabolism. Here, the structural and functional properties of NPs and Nbs are reviewed in parallel with those of sperm whale myoglobin, which is generally taken as the prototype of monomeric globins.

Quantum chemical calculations of tryptophan → heme electron and excitation energy transfer rates in myoglobin

The development of optical multidimensional spectroscopic techniques has opened up new possibilities for the study of biological processes. Recently, ultrafast two-dimensional ultraviolet spectroscopy experiments have determined the rates of tryptophan → heme electron transfer and excitation energy transfer for the two tryptophan residues in myoglobin (Consani et al., Science, 2013, 339, 1586). Here, we show that accurate prediction of these rates can be achieved using Marcus theory in conjunction with time-dependent density functional theory. Key intermediate residues between the donor and acceptor are identified, and in particular the residues Val68 and Ile75 play a critical role in calculations of the electron coupling matrix elements. Our calculations demonstrate how small changes in structure can have a large effect on the rates, and show that the different rates of electron transfer are dictated by the distance between the heme and tryptophan residues, while for excitation energy transfer the orientation of the tryptophan residues relative to the heme is important. © 2017 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Inferring repeat-protein energetics from evolutionary information

Natural protein sequences contain a record of their history. A common constraint in a given protein family is the ability to fold to specific structures, and it has been shown possible to infer the main native ensemble by analyzing covariations in extant sequences. Still, many natural proteins that fold into the same structural topology show different stabilization energies, and these are often related to their physiological behavior. We propose a description for the energetic variation given by sequence modifications in repeat proteins, systems for which the overall problem is simplified by their inherent symmetry. We explicitly account for single amino acid and pair-wise interactions and treat higher order correlations with a single term. We show that the resulting evolutionary field can be interpreted with structural detail. We trace the variations in the energetic scores of natural proteins and relate them to their experimental characterization. The resulting energetic evolutionary field allows the prediction of the folding free energy change for several mutants, and can be used to generate synthetic sequences that are statistically indistinguishable from the natural counterparts.

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.

Ultrafast Excited State Relaxation of a Metalloporphyrin Revealed by Femtosecond X-ray Absorption Spectroscopy

Photoexcited Nickel(II) tetramesitylporphyrin (NiTMP), like many open-shell metalloporphyrins, relaxes rapidly through multiple electronic states following an initial porphyrin-based excitation, some involving metal centered electronic configuration changes that could be harnessed catalytically before excited state relaxation. While a NiTMP excited state present at 100 ps was previously identified by X-ray transient absorption (XTA) spectroscopy at a synchrotron source as a relaxed (d,d) state, the lowest energy excited state (J. Am. Chem. Soc., 2007, 129, 9616 and Chem. Sci., 2010, 1, 642), structural dynamics before thermalization were not resolved due to the ∼100 ps duration of the available X-ray probe pulse. Using the femtosecond (fs) X-ray pulses of the Linac Coherent Light Source (LCLS), the Ni center electronic configuration from the initial excited state to the relaxed (d,d) state has been obtained via ultrafast Ni K-edge XANES (X-ray absorption near edge structure) on a time scale from hundreds of femtoseconds to 100 ps. This enabled the identification of a short-lived Ni(I) species aided by time-dependent density functional theory (TDDFT) methods. Computed electronic and nuclear structure for critical excited electronic states in the relaxation pathway characterize the dependence of the complex's geometry on the electron occupation of the 3d orbitals. Calculated XANES transitions for these excited states assign a short-lived transient signal to the spectroscopic signature of the Ni(I) species, resulting from intramolecular charge transfer on a time scale that has eluded previous synchrotron studies. These combined results enable us to examine the excited state structural dynamics of NiTMP prior to thermal relaxation and to capture intermediates of potential photocatalytic significance.

A molecule for all seasons: The heme

If life without heme-Fe were at all possible, it would definitely be different. Indeed this complex and versatile iron-porphyrin macrocycle upon binding to different “globins” yields hemeproteins crucial to sustain a variety of vital functions, generally classified, for convenience, in a limited number of functional families. Over-and-above the array of functions briefly outlined below, the spectacular progress in molecular genetics seen over the last 30 years led to the discovery of many hitherto unknown novel hemeproteins in prokaryotes and eukaryotes. Here, we highlight a few basic aspects of the chemistry of the hemeprotein universe, in particular those that are relevant to the control of heme-Fe reactivity and specialization, as sculpted by a variety of interactions with the protein moiety.

Nitrobindin: An Ubiquitous Family of All β-Barrel Heme-proteins

Rhodnius prolixus nitrophorins (Rp-NPs), Arabidopsis thaliana nitrobindin (At-Nb), and Homo sapiens THAP4 (Hs-THAP4) are the unique known proteins that use a β-barrel fold to bind ferric heme, which is devoted to NO transport and/or catalysis. The eight-stranded antiparallel β-barrel Rp-NPs, which represent the only heme-binding lipocalins, are devoted to deliver NO into the blood vessel of the host and to scavenge histamine during blood sucking. Regarding Nbs, crystallographic data suggest the ability of At-Nb and Hs-THAP4 to bind ferric heme; however, no data are available with respect to these functions in the natural host. Here, a bioinformatics investigation based on the amino acid sequences and three-dimensional structures of At-Nb and Hs-THAP4 suggests a conservation of the 10-stranded antiparallel β-barrel Nb structural module in all life kingdoms of the evolutionary ladder. In particular, amino acid residues involved in the heme recognition and in the structure stabilization of the Nb structural module are highly conserved (identity > 29%; homology > 83%). Moreover, molecular models of putative Nbs from different organisms match very well with each other and known three-dimensional structures of Nbs. Furthermore, phylogenetic tree reconstruction indicates that NPs and Nbs group in distinct clades. These data indicate that 10-stranded β-barrel Nbs constitute a new ubiquitous heme protein family spanning from bacteria to Homo sapiens. © 2016 IUBMB Life, 68(6):423-428, 2016.

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.

Picosecond to nanosecond dynamics provide a source of conformational entropy for protein folding

Myoglobin can be trapped in fully folded structures, partially folded molten globules, and unfolded states under stable equilibrium conditions.

Investigating the Selectivity of Metalloenzyme Inhibitors in the Presence of Competing Metalloproteins

Metalloprotein inhibitors (MPi) are an important class of therapeutics for the treatment of a variety of diseases, including hypertension, cancer, and HIV/AIDS. However, despite their clinical success, there is an apprehension that MPi may be less selective than other small-molecule therapeutics and more prone to inhibit off-target metalloenzymes. We examined the issue of MPi specificity by investigating the selectivity of a variety of MPi against a representative panel of metalloenzymes in the presence of competing metalloproteins (metallothionein, myoglobin, carbonic anhydrase, and transferrin). Our findings reveal that a wide variety of MPi do not exhibit a decrease in inhibitory activity in the presence of large excesses of competing metalloproteins, suggesting that the competing proteins do not titrate the MPi away from its intended target. This study represents a rudimentary but important means to mimic the biological milieu, which contains other metalloproteins that could compete the MPi away from its target. The strategy used in this study may be a useful approach to examine the selectivity of other MPi in development.

Observing heme doming in myoglobin with femtosecond X-ray absorption spectroscopy

We report time-resolved X-ray absorption measurements after photolysis of carbonmonoxy myoglobin performed at the LCLS X-ray free electron laser with nearly 100 fs (FWHM) time resolution. Data at the Fe K-edge reveal that the photoinduced structural changes at the heme occur in two steps, with a faster (∼70 fs) relaxation preceding a slower (∼400 fs) one. We tentatively attribute the first relaxation to a structural rearrangement induced by photolysis involving essentially only the heme chromophore and the second relaxation to a residual Fe motion out of the heme plane that is coupled to the displacement of myoglobin F-helix.

Ultrafast myoglobin structural dynamics observed with an X-ray free-electron laser

Light absorption can trigger biologically relevant protein conformational changes. The light-induced structural rearrangement at the level of a photoexcited chromophore is known to occur in the femtosecond timescale and is expected to propagate through the protein as a quake-like intramolecular motion. Here we report direct experimental evidence of such 'proteinquake' observed in myoglobin through femtosecond X-ray solution scattering measurements performed at the Linac Coherent Light Source X-ray free-electron laser. An ultrafast increase of myoglobin radius of gyration occurs within 1 picosecond and is followed by a delayed protein expansion. As the system approaches equilibrium it undergoes damped oscillations with a ~3.6-picosecond time period. Our results unambiguously show how initially localized chemical changes can propagate at the level of the global protein conformation in the picosecond timescale.

Molecular dynamics study of naturally existing cavity couplings in proteins

Couplings between protein sub-structures are a common property of protein dynamics. Some of these couplings are especially interesting since they relate to function and its regulation. In this article we have studied the case of cavity couplings because cavities can host functional sites, allosteric sites, and are the locus of interactions with the cell milieu. We have divided this problem into two parts. In the first part, we have explored the presence of cavity couplings in the natural dynamics of 75 proteins, using 20 ns molecular dynamics simulations. For each of these proteins, we have obtained two trajectories around their native state. After applying a stringent filtering procedure, we found significant cavity correlations in 60% of the proteins. We analyze and discuss the structure origins of these correlations, including neighbourhood, cavity distance, etc. In the second part of our study, we have used longer simulations (≥100 ns) from the MoDEL project, to obtain a broader view of cavity couplings, particularly about their dependence on time. Using moving window computations we explored the fluctuations of cavity couplings along time, finding that these couplings could fluctuate substantially during the trajectory, reaching in several cases correlations above 0.25/0.5. In summary, we describe the structural origin and the variations with time of cavity couplings. We complete our work with a brief discussion of the biological implications of these results.

Biophysical highlights from 54 years of macromolecular crystallography

The United Nations has declared 2014 the International Year of Crystallography, and in commemoration, this review features a selection of 54 notable macromolecular crystal structures that have illuminated the field of biophysics in the 54 years since the first excitement of the myoglobin and hemoglobin structures in 1960. Chronological by publication of the earliest solved structure, each illustrated entry briefly describes key concepts or methods new at the time and key later work leveraged by knowledge of the three-dimensional atomic structure.

Frustration in biomolecules

Biomolecules are the prime information processing elements of living matter. Most of these inanimate systems are polymers that compute their own structures and dynamics using as input seemingly random character strings of their sequence, following which they coalesce and perform integrated cellular functions. In large computational systems with finite interaction-codes, the appearance of conflicting goals is inevitable. Simple conflicting forces can lead to quite complex structures and behaviors, leading to the concept of frustration in condensed matter. We present here some basic ideas about frustration in biomolecules and how the frustration concept leads to a better appreciation of many aspects of the architecture of biomolecules, and especially how biomolecular structure connects to function by means of localized frustration. These ideas are simultaneously both seductively simple and perilously subtle to grasp completely. The energy landscape theory of protein folding provides a framework for quantifying frustration in large systems and has been implemented at many levels of description. We first review the notion of frustration from the areas of abstract logic and its uses in simple condensed matter systems. We discuss then how the frustration concept applies specifically to heteropolymers, testing folding landscape theory in computer simulations of protein models and in experimentally accessible systems. Studying the aspects of frustration averaged over many proteins provides ways to infer energy functions useful for reliable structure prediction. We discuss how frustration affects folding mechanisms. We review here how the biological functions of proteins are related to subtle local physical frustration effects and how frustration influences the appearance of metastable states, the nature of binding processes, catalysis and allosteric transitions. In this review, we also emphasize that frustration, far from being always a bad thing, is an essential feature of biomolecules that allows dynamics to be harnessed for function. In this way, we hope to illustrate how Frustration is a fundamental concept in molecular biology.

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

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