Frontier Residues Lining Globin Internal Cavities Present Specific Mechanical Properties

Solid-State Photochemical Cascade Process Boosting Smart Ultralong Room-Temperature Phosphorescence in Bismuth Halides

Molecular ultralong room-temperature phosphorescence (RTP), exhibiting multiple stimuli-responsive characteristics, has garnered considerable attention due to its potential applications in light-emitting devices, sensors, and information safety. This work proposes the utilization of photochemical cascade processes (PCCPs) in molecular crystals to design a stepwise smart RTP switch. By harnessing the sequential dynamics of photo-burst movement (induced by [2+2] photocycloaddition) and photochromism (induced by photogenerated radicals) in a bismuth (Bi)-based metal-organic halide (MOH), a continuous and photo-responsive ultralong RTP can be achieved. Furthermore, utilizing the same Bi-based MOH, diverse application demonstrations, such as multi-mode anti-counterfeiting and information encryption, can be easily implemented. This work thus not only serves as a proof-of-concept for the development of solid-state PCCPs that integrate photosalient effect and photochromism with light-chemical-mechanical energy conversion, but also lays the groundwork for designing new Bi-based MOHs with dynamically responsive ultralong RTP.

Exploring Binding Pockets in the Conformational States of the SARS-CoV-2 Spike Trimers for the Screening of Allosteric Inhibitors Using Molecular Simulations and Ensemble-Based Ligand Docking

Understanding mechanisms of allosteric regulation remains elusive for the SARS-CoV-2 spike protein, despite the increasing interest and effort in discovering allosteric inhibitors of the viral activity and interactions with the host receptor ACE2. The challenges of discovering allosteric modulators of the SARS-CoV-2 spike proteins are associated with the diversity of cryptic allosteric sites and complex molecular mechanisms that can be employed by allosteric ligands, including the alteration of the conformational equilibrium of spike protein and preferential stabilization of specific functional states. In the current study, we combine conformational dynamics analysis of distinct forms of the full-length spike protein trimers and machine-learning-based binding pocket detection with the ensemble-based ligand docking and binding free energy analysis to characterize the potential allosteric binding sites and determine structural and energetic determinants of allosteric inhibition for a series of experimentally validated allosteric molecules. The results demonstrate a good agreement between computational and experimental binding affinities, providing support to the predicted binding modes and suggesting key interactions formed by the allosteric ligands to elicit the experimentally observed inhibition. We establish structural and energetic determinants of allosteric binding for the experimentally known allosteric molecules, indicating a potential mechanism of allosteric modulation by targeting the hinges of the inter-protomer movements and blocking conformational changes between the closed and open spike trimer forms. The results of this study demonstrate that combining ensemble-based ligand docking with conformational states of spike protein and rigorous binding energy analysis enables robust characterization of the ligand binding modes, the identification of allosteric binding hotspots, and the prediction of binding affinities for validated allosteric modulators, which is consistent with the experimental data. This study suggested that the conformational adaptability of the protein allosteric sites and the diversity of ligand bound conformations are both in play to enable efficient targeting of allosteric binding sites and interfere with the conformational changes.

Coarse-Grained Molecular Simulations and Ensemble-Based Mutational Profiling of Protein Stability in the Different Functional Forms of the SARS-CoV-2 Spike Trimers: Balancing Stability and Adaptability in BA.1, BA.2 and BA.2.75 Variants

Evolutionary and functional studies have suggested that the emergence of Omicron variants can be determined by multiple fitness tradeoffs including immune escape, binding affinity, conformational plasticity, protein stability, and allosteric modulation. In this study, we embarked on a systematic comparative analysis of the conformational dynamics, electrostatics, protein stability, and allostery in the different functional states of spike trimers for BA.1, BA.2, and BA.2.75 variants. Using efficient and accurate coarse-grained simulations and atomistic reconstruction of the ensembles, we examined the conformational dynamics of the spike trimers that agree with the recent functional studies, suggesting that BA.2.75 trimers are the most stable among these variants. A systematic mutational scanning of the inter-protomer interfaces in the spike trimers revealed a group of conserved structural stability hotspots that play a key role in the modulation of functional dynamics and are also involved in the inter-protomer couplings through local contacts and interaction networks with the Omicron mutational sites. The results of mutational scanning provided evidence that BA.2.75 trimers are more stable than BA.2 and comparable in stability to the BA.1 variant. Using dynamic network modeling of the S Omicron BA.1, BA.2, and BA.2.75 trimers, we showed that the key network mediators of allosteric interactions are associated with the major stability hotspots that are interconnected along potential communication pathways. The network analysis of the BA.1, BA.2, and BA.2.75 trimers suggested that the increased thermodynamic stability of the BA.2.75 variant may be linked with the organization and modularity of the residue interaction network that allows for allosteric communications between structural stability hotspots and Omicron mutational sites. This study provided a plausible rationale for a mechanism in which Omicron mutations may evolve by targeting vulnerable sites of conformational adaptability to elicit immune escape while maintaining their control on balancing protein stability and functional fitness through robust allosteric communications with the stability hotspots.

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

Frustration-driven allosteric regulation and signal transmission in the SARS-CoV-2 spike omicron trimer structures: a crosstalk of the omicron mutation sites allosterically regulates tradeoffs of protein stability and conformational adaptability

Dissecting the regulatory principles underlying function and activity of the SARS-CoV-2 spike protein at the atomic level is of paramount importance for understanding the mechanisms of virus transmissibility and immune escape. In this work, we introduce a hierarchical computational approach for atomistic modeling of allosteric mechanisms in the SARS-CoV-2 Omicron spike proteins and present evidence of a frustration-based allostery as an important energetic driver of the conformational changes and spike activation. By examining conformational landscapes and the residue interaction networks in the SARS-CoV-2 Omicron spike protein structures, we have shown that the Omicron mutational sites are dynamically coupled and form a central engine of the allosterically regulated spike machinery that regulates the balance and tradeoffs between conformational plasticity, protein stability, and functional adaptability. We have found that the Omicron mutational sites at the inter-protomer regions form regulatory hotspot clusters that control functional transitions between the closed and open states. Through perturbation-based modeling of allosteric interaction networks and diffusion analysis of communications in the closed and open spike states, we have quantified the allosterically regulated activation mechanism and uncover specific regulatory roles of the Omicron mutations. Atomistic reconstruction of allosteric communication pathways and kinetic modeling using Markov transient analysis reveal that the Omicron mutations form the inter-protomer electrostatic bridges that operate as a network of coupled regulatory switches that could control global conformational changes and signal transmission in the spike protein. The results of this study have revealed distinct and yet complementary roles of the Omicron mutation sites as a network of hotspots that enable allosteric modulation of structural stability and conformational changes which are central for spike activation and virus transmissibility.

Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery

Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.

Integrating Conformational Dynamics and Perturbation-Based Network Modeling for Mutational Profiling of Binding and Allostery in the SARS-CoV-2 Spike Variant Complexes with Antibodies: Balancing Local and Global Determinants of Mutational Escape Mechanisms

In this study, we combined all-atom MD simulations, the ensemble-based mutational scanning of protein stability and binding, and perturbation-based network profiling of allosteric interactions in the SARS-CoV-2 spike complexes with a panel of cross-reactive and ultra-potent single antibodies (B1-182.1 and A23-58.1) as well as antibody combinations (A19-61.1/B1-182.1 and A19-46.1/B1-182.1). Using this approach, we quantify the local and global effects of mutations in the complexes, identify protein stability centers, characterize binding energy hotspots, and predict the allosteric control points of long-range interactions and communications. Conformational dynamics and distance fluctuation analysis revealed the antibody-specific signatures of protein stability and flexibility of the spike complexes that can affect the pattern of mutational escape. A network-based perturbation approach for mutational profiling of allosteric residue potentials revealed how antibody binding can modulate allosteric interactions and identified allosteric control points that can form vulnerable sites for mutational escape. The results show that the protein stability and binding energetics of the SARS-CoV-2 spike complexes with the panel of ultrapotent antibodies are tolerant to the effect of Omicron mutations, which may be related to their neutralization efficiency. By employing an integrated analysis of conformational dynamics, binding energetics, and allosteric interactions, we found that the antibodies that neutralize the Omicron spike variant mediate the dominant binding energy hotpots in the conserved stability centers and allosteric control points in which mutations may be restricted by the requirements of the protein folding stability and binding to the host receptor. This study suggested a mechanism in which the patterns of escape mutants for the ultrapotent antibodies may not be solely determined by the binding interaction changes but are associated with the balance and tradeoffs of multiple local and global factors, including protein stability, binding affinity, and long-range interactions.

Between Two Walls: Modeling the Adsorption Behavior of β-Glucosidase A on Bare and SAM-Functionalized Gold Surfaces

The efficient immobilization of enzymes on surfaces remains a complex but central issue in the biomaterials field, which requires us to understand this process at the atomic level. Using a multiscale approach combining all-atom molecular dynamics and coarse-grain Brownian dynamics simulations, we investigated the adsorption behavior of β-glucosidase A (βGA) on bare and self-assembled monolayer (SAM)-functionalized gold surfaces. We monitored the enzyme position and orientation during the molecular dynamics (MD) trajectories and measured the contacts it forms with both surfaces. While the adsorption process has little impact on the protein conformation, it can nonetheless perturb its mechanical properties and catalytic activity. Our results show that compared to the SAM-functionalized surface, the adsorption of βGA on bare gold is more stable, but less specific, and more likely to disrupt the enzyme's function. This observation emphasizes the fact that the structural organization of proteins at the solid interface is a key point when designing devices based on enzyme immobilization, as one must find an acceptable stability-activity trade-off.

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.

Loss of Dihydroxyacid Dehydratase Induces Auxotrophy in Bacillus anthracis

Anthrax disease is caused by infection with the bacteria Bacillus anthracis which, if left untreated, can result in fatal bacteremia and toxemia. Current treatment for infection requires prolonged administration of antibiotics. Despite this, inhalational and gastrointestinal anthrax still result in lethal disease. By identifying key metabolic steps that B. anthracis uses to grow in host-like environments, new targets for antibacterial strategies can be identified. Here, we report that the ilvD gene, which encodes dihydroxyacid dehydratase in the putative pathway for synthesizing branched chain amino acids, is necessary for B. anthracis to synthesize isoleucine de novo in an otherwise limiting microenvironment. We observed that ΔilvD B. anthracis cannot grow in media lacking isoleucine, but growth is restored when exogenous isoleucine is added. In addition, ΔilvD bacilli are unable to utilize human hemoglobin or serum albumin to overcome isoleucine auxotrophy, but can when provided with the murine forms. This species-specific effect is due to the lack of isoleucine in human hemoglobin. Furthermore, even when supplemented with physiological levels of human serum albumin, apotransferrin, fibrinogen, and IgG, the ilvD knockout strain grew poorly relative to nonsupplemented wild type. In addition, comparisons upon infecting humanized mice suggest that murine hemoglobin is a key source of isoleucine for both WT and ΔilvD bacilli. Further growth comparisons in murine and human blood show that the auxotrophy is detrimental for growth in human blood, not murine. This report identifies ilvD as necessary for isoleucine production in B. anthracis, and that it plays a key role in allowing the bacilli to effectively grow in isoleucine poor hosts. IMPORTANCE Anthrax disease, caused by B. anthracis, can cause lethal bacteremia and toxemia, even following treatment with antibiotics. This report identifies the ilvD gene, which encodes a dihydroxyacid dehydratase, as necessary for B. anthracis to synthesize the amino acid isoleucine in a nutrient-limiting environment, such as its mammalian host. The use of this strain further demonstrated a unique species-dependent utilization of hemoglobin as an exogenous source of extracellular isoleucine. By identifying mechanisms that B. anthracis uses to grow in host-like environments, new targets for therapeutic intervention are revealed.

Coarse-grain simulations on NMR conformational ensembles highlight functional residues in proteins

Dynamics are a key feature of protein function, and this is especially true of gating residues, which occupy cavity or tunnel lining positions in the protein structure, and will reversibly switch between open and closed conformations in order to control the diffusion of small molecules within a protein's internal matrix. Earlier work on globins and hydrogenases have shown that these gating residues can be detected using a multiscale scheme combining all-atom classic molecular dynamics simulations and coarse-grain calculations of the resulting conformational ensemble mechanical properties. Here, we show that the structural variations observed in the conformational ensembles produced by NMR spectroscopy experiments are sufficient to induce noticeable mechanical changes in a protein, which in turn can be used to identify residues important for function and forming a mechanical nucleus in the protein core. This new approach, which combines experimental data and rapid coarse-grain calculations and no longer needs to resort to time-consuming all-atom simulations, was successfully applied to five different protein families.

Conformational Stability Adaptation of a Double-Stranded RNA-Binding Domain to Transfer RNA Ligand

The double-stranded RNA-binding domain (dsRBD) is a broadly distributed domain among RNA-maturing enzymes. Although this domain recognizes dsRNA's structures via a conserved canonical structure adopting an α₁-β₁β₂β₃-α₂ topology, several dsRBDs can accommodate discrete structural extensions expanding further their functional repertoire. How these structural elements engage cooperative communications with the canonical structure and how they contribute to the dsRBD's overall folding are poorly understood. Here, we addressed these issues using the dsRBD of human dihydrouridine synthase-2 (hDus2) (hDus2-dsRBD) as a model. This dsRBD harbors N- and C-terminal extensions, the former being directly involved in the recognition of tRNA substrate of hDus2. These extensions engage residues that form a long-range hydrophobic network (LHN) outside the RNA-binding interface. We show by coarse-grain Brownian dynamics that the Nt-extension and its residues F359 and Y364 rigidify the major folding nucleus of the canonical structure via an indirect effect. hDus2-dsRBD unfolds following a two-state cooperative model, whereas both F359A and Y364A mutants, designed to destabilize this LHN, unfold irreversibly. Structural and computational analyses show that these mutants are unstable due to an increase in the dynamics of the two extensions favoring solvent exposure of α2-helix and weakening the main folding nucleus rigidity. This LHN appears essential for maintaining a thermodynamic stability of the overall system and eventually a functional conformation for tRNA recognition. Altogether, our findings suggest that functional adaptability of extended dsRBDs is promoted by a cooperative hydrophobic coupling between the extensions acting as effectors and the folding nucleus of the canonical structure.

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.

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.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

Bridging Enzymatic Structure Function via Mechanics: A Coarse-Grain Approach

Flexibility is a central aspect of protein function, and ligand binding in enzymes involves a wide range of structural changes, ranging from large-scale domain movements to small loop or side-chain rearrangements. In order to understand how the mechanical properties of enzymes, and the mechanical variations that are induced by ligand binding, relate to enzymatic activity, we carried out coarse-grain Brownian dynamics simulations on a set of enzymes whose structures in the unbound and ligand-bound forms are available in the Protein Data Bank. Our results show that enzymes are remarkably heterogeneous objects from a mechanical point of view and that the local rigidity of individual residues is tightly connected to their part in the protein's overall structure and function. The systematic comparison of the rigidity of enzymes in their unbound and bound forms highlights the fact that small conformational changes can induce large mechanical effects, leading to either more or less flexibility depending on the enzyme's architecture and the location of its ligand-biding site. These mechanical variations target a limited number of specific residues that occupy key locations for enzymatic activity, and our approach thus offers a mean to detect perturbation-sensitive sites in enzymes, where the addition or removal of a few interactions will lead to important changes in the proteins internal dynamics.

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

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

Volumetric properties underlying ligand binding in a monomeric hemoglobin: a high-pressure NMR study

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using (1)H-(15)N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

X-ray crystallographic structural characteristics of Arabidopsis hemoglobin I and their functional implications

Genome of the model dicot flowering plant, Arabidopsis thaliana, a popular tool for understanding molecular biology of plant physiology, encodes all three classes of plant hemoglobins that differ in their sequence, ligand binding and spectral properties. As such these globins are of considerable attention. Crystal structures of few members of plant class I nonsymbiotic hemoglobin have been described earlier. Here we report the crystal structure of Arabidopsis class I hemoglobin (AHb1) to 2.2Ǻ and compare its key features with the structures of similar nonsymbiotic hemoglobin from other species. Crystal structure of AHb1 is homologous to the related members with similar globin fold and heme pocket architecture. The structure is homodimeric in the asymmetric unit with both distal and proximal histidines coordinating to the heme iron atom. Residues lining the dimeric interface are also conserved in AHb1 with the exception of additional electrostatic interaction between H112 and E113 of each subunit and that involving Y119 through two water molecules. In addition, differences in heme pocket non-covalent interactions, a novel Ser residue at F7 position, Xe binding site variability, internal cavity topology differences, CD loop conformation and stability and other such properties might explain kinetic variability in AHb1. Detailed cavity analysis of AHb1 showed the presence of a novel long tunnel connecting the distal pockets of both the monomers. Presence of such tunnel, along with conformational heterogeneity observed in the two chains, might suggest cooperative ligand binding and support its role in NO scavenging. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Oxygen migration pathways in NO-bound truncated hemoglobin

Atomistic simulations of dioxygen (O(2)) dynamics and migration in nitric oxide-bound truncated Hemoglobin N (trHbN) of Mycobacterium tuberculosis are reported. From more than 100 ns of simulations the connectivity network involving the metastable states for localization of the O(2) ligand is built and analyzed. It is found that channel I is the primary entrance point for O(2) whereas channel II is predominantly an exit path although access to the protein active site is also possible. For O(2) a new site compared to nitric oxide, from which reaction with the heme group can occur, was found. As this site is close to the heme iron, it could play an important role in the dioxygenation mechanism as O(2) can remain there for hundreds of picoseconds after which it can eventually leave the protein, while NO is localized in Xe2. The present study supports recent experimental work which proposed that O(2) docks in alternative pockets than Xe close to the reactive site. Similar to other proteins, a phenylalanine residue (Phe62) plays the role of a gate along the access route in channel I. The most highly connected site is the Xe3 pocket which is a "hub" and free energy barriers between the different metastable states are ≈1.5 kcal mol(-1) which allows facile O(2) migration within the protein.

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.

xPyder: a PyMOL plugin to analyze coupled residues and their networks in protein structures

A versatile method to directly identify and analyze short- or long-range coupled or communicating residues in a protein conformational ensemble is of extreme relevance to achieve a complete understanding of protein dynamics and structural communication routes. Here, we present xPyder, an interface between one of the most employed molecular graphics systems, PyMOL, and the analysis of dynamical cross-correlation matrices (DCCM). The approach can also be extended, in principle, to matrices including other indexes of communication propensity or intensity between protein residues, as well as the persistence of intra- or intermolecular interactions, such as those underlying protein dynamics. The xPyder plugin for PyMOL 1.4 and 1.5 is offered as Open Source software via the GPL v2 license, and it can be found, along with the installation package, the user guide, and examples, at http://linux.btbs.unimib.it/xpyder/.

Thermal fluctuations of haemoglobin from different species: adaptation to temperature via conformational dynamics

Thermodynamic stability, configurational motions and internal forces of haemoglobin (Hb) of three endotherms (platypus, Ornithorhynchus anatinus; domestic chicken, Gallus gallus domesticus and human, Homo sapiens) and an ectotherm (salt water crocodile, Crocodylus porosus) were investigated using circular dichroism, incoherent elastic neutron scattering and coarse-grained Brownian dynamics simulations. The experimental results from Hb solutions revealed a direct correlation between protein resilience, melting temperature and average body temperature of the different species on the 0.1 ns time scale. Molecular forces appeared to be adapted to permit conformational fluctuations with a root mean square displacement close to 1.2 Å at the corresponding average body temperature of the endotherms. Strong forces within crocodile Hb maintain the amplitudes of motion within a narrow limit over the entire temperature range in which the animal lives. In fully hydrated powder samples of human and chicken, Hb mean square displacements and effective force constants on the 1 ns time scale showed no differences over the whole temperature range from 10 to 300 K, in contrast to the solution case. A complementary result of the study, therefore, is that one hydration layer is not sufficient to activate all conformational fluctuations of Hb in the pico- to nanosecond time scale which might be relevant for biological function. Coarse-grained Brownian dynamics simulations permitted to explore residue-specific effects. They indicated that temperature sensing of human and chicken Hb occurs mainly at residues lining internal cavities in the β-subunits.