Evolutionary and Functional Relationships in the Truncated Hemoglobin Family

Heme d formation in a Shewanella benthica hemoglobin

In our continued investigations of microbial globins, we solved the structure of a truncated hemoglobin from Shewanella benthica, an obligate psychropiezophilic bacterium. The distal side of the heme active site is lined mostly with hydrophobic residues, with the exception of a tyrosine, Tyr34 (CD1) and a histidine, His24 (B13). We found that purified SbHbN, when crystallized in the ferric form with polyethylene glycol as precipitant, turned into a green color over weeks. The electron density obtained from the green crystals accommodated a trans heme d, a chlorin-type derivative featuring a γ-spirolactone and a vicinal hydroxyl group on a pyrroline ring. In solution, exposure of the protein to one equivalent of hydrogen peroxide resulted in a similar green color change, but caused by the formation of multiple products. These were oxidation species released on protein denaturation, likely including heme d, and a species with heme covalently attached to the polypeptide. The Tyr34Phe replacement prevented the formation of both heme d and the covalent linkage. The ready modification of heme b by SbHbN expands the range of chemistries supported by the globin fold and offers a route to a novel heme cofactor.

Globin phylogeny, evolution and function, the newest update

Our globin census update allows us to refine our vision of globin origin, evolution, and structure to function relationship in the context of the currently accepted tree of life. The modern globin domain originates as a single domain, three-over-three α-helical folded structure before the diversification of the kingdoms of life (Bacteria, Archaea, Eukarya). Together with the diversification of prokaryotes, three monophyletic globin families (M, S, and T) emerged, most likely in Proteobacteria and Actinobacteria, displaying specific sequence and structural features, and spread by vertical and horizontal gene transfer, most probably already present in the last universal common ancestor (LUCA). Non-globin domains were added, and eventually lost again, creating multi-domain structures in key branches of M- (FHb and Adgb) and the vast majority of S globins, which with their coevolved multi-domain architectures, have predominantly "sensor" functions. Single domain T-family globins diverged into four major groups and most likely display functions related to reactive nitrogen and oxygen species (RNOS) chemistry, as well as oxygen storage/transport which drives the evolution of its major branches with their characteristic key distal residues (B10, E11, E7, and G8). M-family evolution also lead to distinctive major types (FHb and Fgb, Ngb, Adgb, GbX vertebrate Gbs), and shows the shift from high oxygen affinity controlled by TyrB10-Gln/AsnE11 likely related to RNOS chemistry in microorganisms, to a moderate oxygen affinity storage/transport function controlled by hydrophobic B10/E11-HisE7 in multicellular animals.

The globins of cyanobacteria and green algae: An update

The globin superfamily of proteins is ancient and diverse. Regular assessments based on the increasing number of available genome sequences have elaborated on a complex evolutionary history. In this review, we present a summary of a decade of advances in characterising the globins of cyanobacteria and green algae. The focus is on haem-containing globins with an emphasis on recent experimental developments, which reinforce links to nitrogen metabolism and nitrosative stress response in addition to dioxygen management. Mention is made of globins that do not bind haem to provide an encompassing view of the superfamily and perspective on the field. It is reiterated that an effort toward phenotypical and in-vivo characterisation is needed to elucidate the many roles that these versatile proteins fulfil in oxygenic photosynthetic microbes. It is also proposed that globins from oxygenic organisms are promising proteins for applications in the biotechnology arena.

Bacterial cytochrome P450s: a bioinformatics odyssey of substrate discovery

Bacterial P450 cytochromes (BacCYPs) are versatile heme-containing proteins responsible for oxidation reactions on a wide range of substrates, contributing to the production of valuable natural products with limitless biotechnological potential. While the sequencing of microbial genomes has provided a wealth of BacCYP sequences, functional characterization lags behind, hindering our understanding of their roles. This study employs a comprehensive approach to predict BacCYP substrate specificity, bridging the gap between sequence and function. We employed an integrated approach combining sequence and functional data analysis, genomic context exploration, 3D structural modeling with molecular docking, and phylogenetic clustering. The research begins with an in-depth analysis of BacCYP sequence diversity and structural characteristics, revealing conserved motifs and recurrent residues in the active site. Phylogenetic analysis identifies distinct groups within the BacCYP family based on sequence similarity. However, our study reveals that sequence alone does not consistently predict substrate specificity, necessitating additional perspectives. The study delves into the genetic context of BacCYPs, utilizing neighboring gene information to infer potential substrates, a method proven very effective in many cases. Molecular docking is employed to assess BacCYP-substrate interactions, confirming potential substrates and providing insights into selectivity. Finally, a comprehensive strategy is proposed for predicting BacCYP substrates, involving all the evaluated approaches. The effectiveness of this strategy is demonstrated with two case studies, highlighting its potential for substrate discovery.

Sphaerotilus natans hemoglobins have an NADH oxidation activity and promote the yield of limonene in an engineered E. coli strain

Bacterial hemoglobins play important roles inside the cell. Phylogenetically, they belong to three different families: the single domain hemoglobin, flavohemoglobin and truncated hemoglobin. Vitreoscilla hemoglobin (VHb) is the first characterized bacterial hemoglobin, and belongs to the single domain hemoglobin family. Heterologous expression of VHb promotes the growth of host cells under microaerobic conditions, and enhances the yield of products during fermentation. Although VHb has been widely applied in the biotechnology field, other bacterial hemoglobins have not demonstrated similar applications. In this study, we identified four bacterial hemoglobins from the microaerobic growing bacterium Sphaerotilus natans, including one flavohemoglobins (FHB) and three truncated hemoglobins (THB1, THB2 and THB3). Absorption spectrum studies validate the existent of the Soret peak and Q-band characteristic to heme and suggest heme groups in FHB and THB1 are hexa- or penta-coordinated, respectively. Our studies demonstrate that FHB and all three truncated hemoglobins have NADH oxidation and radical production activities, which is surprising since truncated hemoglobins do not have a reductase domain that could bind NADH. However, the M. tuberculosis HbN does not show these activities, indicating they are not universal among truncated hemoglobins. Docking studies suggest the nicotinamide ring of NADH may bind to the distal heme pocket of THB1, suggesting the direct electron transfer from NADH to heme might be possible. Our truncated hemoglobins also show peroxidase activities that in THB2 and THB3 could be inhibited by FdR, indicating possible interactions between FdR and truncate hemoglobins. Expression of FHB and THB1 in E. coli could promote cell growth. THB1 also enhances the production of limonene in an engineered E. coli strain, while VHb does not have this effect, which suggests that studies on truncated hemoglobins may lead to the discovery of new and more powerful tools that could have profound impact on biotechnology.

Architectural digest: Thermodynamic stability and domain structure of a consensus monomeric globin

Artificial proteins representing the consensus of a set of homologous sequences have attracted attention for their increased thermodynamic stability and conserved activity. Here, we applied the consensus approach to a b-type heme-binding protein to inspect the contribution of a dissociable cofactor to enhanced stability and the chemical consequences of creating a generic heme environment. We targeted the group 1 truncated hemoglobin (TrHb1) subfamily of proteins for their small size (∼120 residues) and ease of characterization. The primary structure, derived from a curated set of ∼300 representative sequences, yielded a highly soluble consensus globin (cGlbN) enriched in acidic residues. Optical and NMR spectroscopies revealed high-affinity heme binding in the expected site and in two orientations. At neutral pH, proximal and distal iron coordination was achieved with a pair of histidine residues, as observed in some natural TrHb1s, and with labile ligation on the distal side. As opposed to studied TrHb1s, which undergo additional folding upon heme binding, cGlbN displayed the same extent of secondary structure whether the heme was associated with the protein or not. Denaturation required guanidine hydrochloride and showed that apo- and holoprotein unfolded in two transitions-the first (occurring with a midpoint of ∼2 M) was shifted to higher denaturant concentration in the holoprotein (∼3.7 M) and reflected stabilization due to heme binding, while the second transition (∼6.2 M) was common to both forms. Thus, the consensus sequence stabilized the protein but exposed the existence of two separately cooperative subdomains within the globin architecture, masked as one single domain in TrHb1s with typical stabilities. The results suggested ways in which specific chemical or thermodynamic features may be controlled in artificial heme proteins.

Best Practices on QM/MM Simulations of Biological Systems

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

A novel single sensor hemoglobin domain from the thermophilic cyanobacteria Thermosynechococcus elongatus BP-1 exhibits higher pH but lower thermal stability compared to globins from mesophilic organisms

Thermosynechococcus elongatus-BP1 belongs to the class of photoautotrophic cyanobacterial organisms. The presence of chlorophyll a, carotenoids, and phycocyanobilin are the characteristics that categorize T. elongatus as a photosynthetic organism. Here, we report the structural and spectroscopic characteristics of novel hemoglobin (Hb) Synel Hb from T.elongatus, synonymous with Thermosynechococcus vestitus BP-1. The X-ray crystal structure (2.15 Å) of Synel Hb suggests the presence of a globin domain with a pre-A helix similar to the sensor domain (S) family of Hbs. The rich hydrophobic core accommodates heme in a penta-coordinated state and readily binds an extraneous ligand(imidazole). The absorption and circular dichroic spectral analysis of Synel Hb reiteratedthat the heme is in Fe^III+ state with a predominantly α-helical structure similar to myoglobin. Synel Hb displays higher resistance to structural perturbations induced via external stresses like pH and guanidium hydrochloride, which is comparable to Synechocystis Hb. However, Synel Hb exhibited lower thermal stability compared to mesophilic hemoglobins. Overall, the data is suggestive of the structural sturdiness of Synel Hb, which probably corroborates its origin in extreme thermophilic conditions. The stable globin provides scope for further investigation and may lead to new insights with scope for engineering stability in hemoglobin-based oxygen carriers.

Group II truncated haemoglobin YjbI prevents reactive oxygen species-induced protein aggregation in Bacillus subtilis

Oxidative stress-mediated formation of protein hydroperoxides can induce irreversible fragmentation of the peptide backbone and accumulation of cross-linked protein aggregates, leading to cellular toxicity, dysfunction, and death. However, how bacteria protect themselves from damages caused by protein hydroperoxidation is unknown. Here, we show that YjbI, a group II truncated haemoglobin from Bacillus subtilis, prevents oxidative aggregation of cell-surface proteins by its protein hydroperoxide peroxidase-like activity, which removes hydroperoxide groups from oxidised proteins. Disruption of the yjbI gene in B. subtilis lowered biofilm water repellence, which associated with the cross-linked aggregation of the biofilm matrix protein TasA. YjbI was localised to the cell surface or the biofilm matrix, and the sensitivity of planktonically grown cells to generators of reactive oxygen species was significantly increased upon yjbI disruption, suggesting that YjbI pleiotropically protects labile cell-surface proteins from oxidative damage. YjbI removed hydroperoxide residues from the model oxidised protein substrate bovine serum albumin and biofilm component TasA, preventing oxidative aggregation in vitro. Furthermore, the replacement of Tyr⁶³ near the haem of YjbI with phenylalanine resulted in the loss of its protein peroxidase-like activity, and the mutant gene failed to rescue biofilm water repellency and resistance to oxidative stress induced by hypochlorous acid in the yjbI-deficient strain. These findings provide new insights into the role of truncated haemoglobin and the importance of hydroperoxide removal from proteins in the survival of aerobic bacteria.

Transcript analysis and expression of the glbO gene, encoding truncated hemoglobin,O, of M. smegmatis implicate its role under hypoxia and oxidative stress

Although truncated hemoglobin O, (trHbO), is ubiquitous among mycobacteria, its physiological function is not very obvious and may be diverse. In an attempt to understand role of trHbO in cellular metabolism of a non-pathogenic mycobacterium, we analysed expression profile of the glbO gene, encoding trHbO, in M. smegmatis and studied implications of its overexpression on physiology of its host under different environmental conditions. Quantitative RT-PCR indicated that transcript level of the glbO gene remains low at a basal level under aerobic growth cycle of M. smegmatis but its level gets induced significantly during low oxygen, oxidative stress and macrophage infection. Overexpression of the glbO gene enhanced growth of M. smegmatis under hypoxia, promoted pellicle biofilm formation and provided resistance towards oxidative stress. Additionally, glbO gene overexpressing M. smegmatis exhibited enhanced cell survival over isogenic control cells and altered the level of pro- and anti- inflammatory cytokines during intracellular infection. These results suggested important role of trHbO, in supporting the cellular metabolism and survival of M, smegmatis both under low oxygen and oxidative stress.

Hydroxylamine-induced oxidation of ferrous nitrobindins

Hemoglobin and myoglobin are generally taken as molecular models of all-α-helical heme-proteins. On the other hand, nitrophorins and nitrobindins (Nb), which are arranged in 8 and 10 β-strands, respectively, represent the molecular models of all-β-barrel heme-proteins. Here, kinetics of the hydroxylamine- (HA-) mediated oxidation of ferrous Mycobacterium tuberculosis, Arabidopsis thaliana, and Homo sapiens nitrobindins (Mt-Nb(II), At-Nb(II), and Hs-Nb(II), respectively), at pH 7.0 and 20.0 °C, are reported. Of note, HA displays antibacterial properties and is a good candidate for the treatment and/or prevention of reactive nitrogen species- (RNS-) linked aging-related pathologies, such as macular degeneration. Under anaerobic conditions, mixing the Mt-Nb(II), At-Nb(II), and Hs-Nb(II) solutions with the HA solutions brings about absorbance spectral changes reflecting the formation of the ferric derivative (i.e., Mt-Nb(III), At-Nb(III), and Hs-Nb(III), respectively). Values of the second order rate constant for the HA-mediated oxidation of Mt-Nb(II), At-Nb(II), and Hs-Nb(II) are 1.1 × 10⁴ M^-1 s^-1, 6.5 × 10⁴ M^-1 s^-1, and 2.2 × 10⁴ M^-1 s^-1, respectively. Moreover, the HA:Nb(II) stoichiometry is 1:2 as reported for ferrous deoxygenated and carbonylated all-α-helical heme-proteins. A comparative look of the HA reduction kinetics by several ferrous heme-proteins suggests that an important role might be played by residues (such as His or Tyr) in the proximity of the heme-Fe atom either coordinating it or not. In this respect, Nbs seem to exploit somewhat different structural aspects, indicating that redox mechanisms for the heme-Fe(II)-to-heme-Fe(III) conversion might differ between all-α-helical and all-β-barrel heme-proteins.

Functional Classification and Characterization of the Fungal Glycoside Hydrolase 28 Protein Family

Pectin is a major constituent of the plant cell wall, comprising compounds with important industrial applications such as homogalacturonan, rhamnogalacturonan and xylogalacturonan. A large array of enzymes is involved in the degradation of this amorphous substrate. The Glycoside Hydrolase 28 (GH28) family includes polygalacturonases (PG), rhamnogalacturonases (RG) and xylogalacturonases (XG) that share a structure of three to four pleated β-sheets that form a rod with the catalytic site amidst a long, narrow groove. Although these enzymes have been studied for many years, there has been no systematic analysis. We have collected a comprehensive set of GH28 encoding sequences to study their evolution in fungi, directed at obtaining a functional classification, as well as at the identification of substrate specificity as functional constraint. Computational tools such as Alphafold, Consurf and MEME were used to identify the subfamilies’ characteristics. A hierarchic classification defines the major classes of endoPG, endoRG and endoXG as well as three exoPG classes. Ascomycete endoPGs are further classified in two subclasses whereas we identify four exoRG subclasses. Diversification towards exomode is explained by loops that appear inserted in a number of turns. Substrate-driven diversification can be identified by various specificity determining positions that appear to surround the binding groove.

Mycobacterium avium subsp. paratuberculosis Virulence: A Review

To propose a solution for control of Mycobacterium avium subsp. paratuberculosis (MAP) infections in animals as well as in humans, and develop effective prevention, diagnostic and treatment strategies, it is essential to understand the molecular mechanisms of MAP pathogenesis. In the present review, we discuss the mechanisms utilised by MAP to overcome the host defense system to achieve the virulence status. Putative MAP virulence genes are mentioned and their probable roles in view of other mycobacteria are discussed. This review provides information on MAP strain diversity, putative MAP virulence factors and highlights the knowledge gaps regarding MAP virulence mechanisms that may be important in control and prevention of paratuberculosis.

Control of distal lysine coordination in a monomeric hemoglobin: A role for heme peripheral interactions

THB1 is a monomeric truncated hemoglobin (TrHb) found in the cytoplasm of the green alga Chlamydomonas reinhardtii. The canonical heme coordination scheme in hemoglobins is a proximal histidine ligand and an open distal site. In THB1, the latter site is occupied by Lys53, which is likely to facilitate Fe(II)/Fe(III) redox cycling but hinders dioxygen binding, two features inherent to the NO dioxygenase activity of the protein. TrHb surveys show that a lysine at a position aligning with Lys53 is an insufficient determinant of coordination, and in this study, we sought to identify factors controlling lysine affinity for the heme iron. We solved the "Lys-off" X-ray structure of THB1, represented by the cyanide adduct of the Fe(III) protein, and hypothesized that interactions that differ between the known "Lys-on" structure and the Lys-off structure participate in the control of Lys53 affinity for the heme iron. We applied an experimental approach (site-directed mutagenesis, heme modification, pH titrations in the Fe(III) and Fe(II) states) and a computational approach (MD simulations in the Fe(II) state) to assess the role of heme propionate-protein interactions, distal helix capping, and the composition of the distal pocket. All THB1 modifications resulted in a weakening of lysine affinity and affected the coupling between Lys53 proton binding and heme redox potential. The results supported the importance of specific heme peripheral interactions for the pH stability of iron coordination and the ability of the protein to undergo redox reactions.

Allostery in the nitric oxide dioxygenase mechanism of flavohemoglobin

The substrates O₂ and NO cooperatively activate the NO dioxygenase function of Escherichia coli flavohemoglobin. Steady-state and transient kinetic measurements support a structure-based mechanistic model in which O₂ and NO movements and conserved amino acids at the E11, G8, E2, E7, B10, and F7 positions within the globin domain control activation. In the cooperative and allosteric mechanism, O₂ migrates to the catalytic heme site via a long hydrophobic tunnel and displaces LeuE11 away from the ferric iron, which forces open a short tunnel to the catalytic site gated by the ValG8/IleE15 pair and LeuE11. NO permeates this tunnel and leverages upon the gating side chains triggering the CD loop to furl, which moves the E and F-helices and switches an electron transfer gate formed by LysF7, GlnE7, and water. This allows FADH₂ to reduce the ferric iron, which forms the stable ferric–superoxide–TyrB10/GlnE7 complex. This complex reacts with internalized NO with a bimolecular rate constant of 10¹⁰ M⁻¹ s⁻¹ forming nitrate, which migrates to the CD loop and unfurls the spring-like structure. To restart the cycle, LeuE11 toggles back to the ferric iron. Actuating electron transfer with O₂ and NO movements averts irreversible NO poisoning and reductive inactivation of the enzyme. Together, structure snapshots and kinetic constants provide glimpses of intermediate conformational states, time scales for motion, and associated energies.

Acanthamoeba castellanii as a Screening Tool for Mycobacterium avium Subspecies paratuberculosis Virulence Factors with Relevance in Macrophage Infection

The high prevalence of Johne's disease has driven a continuous effort to more readily understand the pathogenesis of the etiological causative bacterium, Mycobacterium avium subsp. paratuberculosis (MAP), and to develop effective preventative measures for infection spread. In this study, we aimed to create an in vivo MAP infection model employing an environmental protozoan host and used it as a tool for selection of bacterial virulence determinants potentially contributing to MAP survival in mammalian host macrophages. We utilized Acanthamoeba castellanii (amoeba) to explore metabolic consequences of the MAP-host interaction and established a correlation between metabolic changes of this phagocytic host and MAP virulence. Using the library of gene knockout mutants, we identified MAP clones that can either enhance or inhibit amoeba metabolism and we discovered that, for most part, it mirrors the pattern of MAP attenuation or survival during infection of macrophages. It was found that MAP mutants that induced an increase in amoeba metabolism were defective in intracellular growth in macrophages. However, MAP clones that exhibited low metabolic alteration in amoeba were able to survive at a greater rate within mammalian cells, highlighting importance of both category of genes in bacterial pathogenesis. Sequencing of MAP mutants has identified several virulence factors previously shown to have a biological relevance in mycobacterial survival and intracellular growth in phagocytic cells. In addition, we uncovered new genetic determinants potentially contributing to MAP pathogenicity. Results of this study support the use of the amoeba model system as a quick initial screening tool for selection of virulence factors of extremely slow-grower MAP that is challenging to study.

Multi-scale simulation reveals that an amino acid substitution increases photosensitizing reaction inputs in Rhodopsins

Evaluating the availability of molecular oxygen (O₂ ) and energy of excited states in the retinal binding site of rhodopsin is a crucial challenging first step to understand photosensitizing reactions in wild-type (WT) and mutant rhodopsins by absorbing visible light. In the present work, energies of the ground and excited states related to 11-cis-retinal and the O₂ accessibility to the β-ionone ring are evaluated inside WT and human M207R mutant rhodopsins. Putative O₂ pathways within rhodopsins are identified by using molecular dynamics simulations, Voronoi-diagram analysis, and implicit ligand sampling while retinal energetic properties are investigated through density functional theory, and quantum mechanical/molecular mechanical methods. Here, the predictions reveal that an amino acid substitution can lead to enough energy and O₂ accessibility in the core hosting retinal of mutant rhodopsins to favor the photosensitized singlet oxygen generation, which can be useful in understanding retinal degeneration mechanisms and in designing blue-lighting-absorbing proteic photosensitizers.

Plant hemoglobins: a journey from unicellular green algae to vascular plants

Globins (Glbs) are widely distributed in archaea, bacteria and eukaryotes. They can be classified into proteins with 2/2 or 3/3 α-helical folding around the heme cavity. Both types of Glbs occur in green algae, bryophytes and vascular plants. The Glbs of angiosperms have been more intensively studied, and several protein structures have been solved. They can be hexacoordinate or pentacoordinate, depending on whether a histidine is coordinating or not at the sixth position of the iron atom. The 3/3 Glbs of class 1 and the 2/2 Glbs (also called class 3 in plants) are present in all angiosperms, whereas the 3/3 Glbs of class 2 have been only found in early angiosperms and eudicots. The three Glb classes are expected to play different roles. Class 1 Glbs are involved in hypoxia responses and modulate NO concentration, which may explain their roles in plant morphogenesis, hormone signaling, cell fate determination, nutrient deficiency, nitrogen metabolism and plant-microorganism symbioses. Symbiotic Glbs derive from class 1 or class 2 Glbs and transport O₂ in nodules. The physiological roles of class 2 and class 3 Glbs are poorly defined but could involve O₂ and NO transport and/or metabolism, respectively. More research is warranted on these intriguing proteins to determine their non-redundant functions.

Lucina pectinata oxyhemoglobin (II-III) heterodimer pH susceptibility

Lucina pectinata live in high concentrations of hydrogen sulfide (H₂S) and contains one hemoglobin, Hemoglobin I (HbI), transporting H₂S and two hemoglobins, Hemoglobin II (HbII) and Hemoglobin (HbIII), transferring dioxygen to symbionts. HbII and HbIII contain B10 tyrosine (Tyr) and E7 glutamine (Gln) in the heme pocket generating an efficient hydrogen bonding network with the (HbII-HbIII)-O₂ species, leading to very low ligand dissociation rates. The results indicate that the oxy-hemeprotein is susceptible to pH from 4 to 9, at acidic conditions, and as a function of the potassium ferricyanide concentration, 100% of the met-aquo derivative is produced. Without a strong oxidant, pH 5 generates a small concentration of the met-aquo complex. The process is accelerated by the presence of salts, as indicated by the crystallization structures and UV-Vis spectra. The results suggest that acidic pH generates conformational changes associated with B10 and E7 heme pocket amino acids, weakening the (HbII-HbIII)-O₂ hydrogen bond network. The observation is supported by X-ray crystallography, since at pH 4 and 5, the heme-Fe tends to oxidize, while at pH 7, the oxy-heterodimer is present. Conformational changes also are observed at higher pH by the presence of a 605 nm transition associated with the iron heme-Tyr interaction. Therefore, pH is one crucial factor regulating the (HbII-HbIII)-O₂ complex hydrogen-bonding network. Thus, it can be proposed that the hydrogen bonding adjustments between the heme bound O₂ and the Tyr and Gln amino acids contribute to oxygen dissociation from the (HbII-HbIII)-O₂ system.

Conformational Flexibility Drives Cold Adaptation in Pseudoalteromonas haloplanktis TAC125 Globins

Significance: Temperature is one of the most important drivers in shaping protein adaptations. Many biochemical and physiological processes are influenced by temperature. Proteins and enzymes from organisms living at low temperature are less stable in comparison to high-temperature adapted proteins. The lower stability is generally due to greater conformational flexibility. Recent Advances: Adaptive changes in the structure of cold-adapted proteins may occur at subunit interfaces, distant from the active site, thus producing energy changes associated with conformational transitions transmitted to the active site by allosteric modulation, valid also for monomeric proteins in which tertiary structural changes may play an essential role. Critical Issues: Despite efforts, the current experimental and computational methods still fail to produce general principles on protein evolution, since many changes are protein and species dependent. Environmental constraints or other biological cellular signals may override the ancestral information included in the structure of the protein, thus introducing inaccuracy in estimates and predictions on the evolutionary adaptations of proteins in response to cold adaptation. Future Directions: In this review, we describe the studies and approaches used to investigate stability and flexibility in the cold-adapted globins of the Antarctic marine bacterium Pseudoalteromonas haloplanktis TAC125. In fact, future research directions will be prescient on more detailed investigation of cold-adapted proteins and the role of fluctuations between different conformational states.

Distinctive structural properties of THB11, a pentacoordinate Chlamydomonas reinhardtii truncated hemoglobin with N- and C-terminal extensions

Hemoglobins (Hbs) utilize heme b as a cofactor and are found in all kingdoms of life. The current knowledge reveals an enormous variability of Hb primary sequences, resulting in topological, biochemical and physiological individuality. As Hbs appear to modulate their reactivities through specific combinations of structural features, predicting the characteristics of a given Hb is still hardly possible. The unicellular green alga Chlamydomonas reinhardtii contains 12 genes encoding diverse Hbs of the truncated lineage, several of which possess extended N- or C-termini of unknown function. Studies on some of the Chlamydomonas Hbs revealed yet unpredictable structural and biochemical variations, which, along with a different expression of their genes, suggest diverse physiological roles. Chlamydomonas thus represents a promising system to analyze the diversification of Hb structure, biochemistry and physiology. Here, we report the crystal structure, resolved to 1.75 Å, of the heme-binding domain of cyanomet THB11 (Cre16.g662750), one of the pentacoordinate algal Hbs, which offer a free Fe-coordination site in the reduced state. The overall fold of THB11 is conserved, but individual features such as a kink in helix E, a tilted heme plane and a clustering of methionine residues at a putative tunnel exit appear to be unique. Both N- and C-termini promote the formation of oligomer mixtures, and the absence of the C terminus results in reduced nitrite reduction rates. This work widens the structural and biochemical knowledge on the 2/2Hb family and suggests that the N- and C-terminal extensions of the Chlamydomonas 2/2Hbs modulate their reactivity by intermolecular interactions.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O₂ and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Contribution of YjbIH to Virulence Factor Expression and Host Colonization in Staphylococcus aureus

To persist within the host and cause disease, Staphylococcus aureus relies on its ability to precisely fine-tune virulence factor expression in response to rapidly changing environments. During an unbiased transposon mutant screen, we observed that disruption of a two-gene operon, yjbIH, resulted in decreased levels of pigmentation and aureolysin (Aur) activity relative to the wild-type strain. Further analyses revealed that YjbH, a predicted thioredoxin-like oxidoreductase, is predominantly responsible for the observed yjbIH mutant phenotypes, though a minor role exists for the putative truncated hemoglobin YjbI. These differences were due to significantly decreased expression of crtOPQMN and aur Previous studies found that YjbH targets the disulfide- and oxidative stress-responsive regulator Spx for degradation by ClpXP. The absence of yjbH or yjbI resulted in altered sensitivities to nitrosative and oxidative stress and iron deprivation. Additionally, aconitase activity was altered in the yjbH and yjbI mutant strains. Decreased levels of pigmentation and aureolysin (Aur) activity in the yjbH mutant were found to be Spx dependent. Lastly, we used a murine sepsis model to determine the effect of the yjbIH deletion on pathogenesis and found that the mutant was better able to colonize the kidneys and spleens during an acute infection than the wild-type strain. These studies identified changes in pigmentation and protease activity in response to YjbIH and are the first to have shown a role for these proteins during infection.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

Structure-function analysis of Sedolisins: evolution of tripeptidyl peptidase and endopeptidase subfamilies in fungi

Background

Sedolisins are acid proteases that are related to the basic subtilisins. They have been identified in all three superkingdoms but are not ubiquitous, although fungi that secrete acids as part of their lifestyle can have up to six paralogs. Both TriPeptidyl Peptidase (TPP) and endopeptidase activity have been identified and it has been suggested that these correspond to separate subfamilies.

Results

We studied eukaryotic sedolisins by computational analysis. A maximum likelihood tree shows one major clade containing non-fungal sequences only and two major as well as two minor clades containing only fungal sequences. One of the major fungal clades contains all known TPPs whereas the other contains characterized endosedolisins. We identified four Cluster Specific Inserts (CSIs) in endosedolisins, of which CSIs 1, 3 and 4 appear as solvent exposed according to structure modeling. Part of CSI2 is exposed but a short stretch forms a novel and partially buried α-helix that induces a conformational change near the binding pocket. We also identified a total of 15 specificity determining positions (SDPs) of which five, identified in two independent analyses, form highly connected SDP sub-networks. Modeling of virtual mutants suggests a key role for the W307A or F307A substitution. The remaining four key SDPs physically interact at the interface of the catalytic domain and the enzyme’s prosegment. Modeling of virtual mutants suggests these SDPs are indeed required to compensate the conformational change induced by CSI2 and the A307. One of the two small fungal clades concerns a subfamily that contains 213 sequences, is mostly similar to the major TPP subfamily but differs, interestingly, in position 307, showing mostly isoleucine and threonine.

Conclusions

Analysis confirms there are at least two sedolisin subfamilies in fungi: TPPs and endopeptidases, and suggests a third subfamily with unknown characteristics. Sequence and functional diversification was centered around buried SDP307 and resulted in a conformational change of the pocket. Mutual Information network analysis forms a useful instrument in the corroboration of predicted SDPs.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2404-y) contains supplementary material, which is available to authorized users.

Structure and function of haemoglobins

Haemoglobin (Hb) is widely known as the iron-containing protein in blood that is essential for O₂ transport in mammals. Less widely recognised is that erythrocyte Hb belongs to a large family of Hb proteins with members distributed across all three domains of life-bacteria, archaea and eukaryotes. This review, aimed chiefly at researchers new to the field, attempts a broad overview of the diversity, and common features, in Hb structure and function. Topics include structural and functional classification of Hbs; principles of O₂ binding affinity and selectivity between O₂/NO/CO and other small ligands; hexacoordinate (containing bis-imidazole coordinated haem) Hbs; bacterial truncated Hbs; flavohaemoglobins; enzymatic reactions of Hbs with bioactive gases, particularly NO, and protection from nitrosative stress; and, sensor Hbs. A final section sketches the evolution of work on the structural basis for allosteric O₂ binding by mammalian RBC Hb, including the development of newer kinetic models. Where possible, reference to historical works is included, in order to provide context for current advances in Hb research.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -“r” and “t”-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the β subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

Multidomain truncated hemoglobins: New members of the globin family exhibiting tandem repeats of globin units and domain fusion

Truncated hemoglobins (trHbs) are considered the most primitive members of globin superfamily and traditionally exist as a single domain heme protein in three distinct structural organizations, type I (trHb1_N), type II (trHb2_O) and type III (trHb3_P). Our search of microbial and lower eukaryotic genomes revealed a broad array of multidomain organization, representing multiunit and chimeric forms of trHbs, where multiple units of trHbs are joined together and/or integrated with distinct functional domains. Globin motifs of these multidomain trHbs were from all three groups of trHbs and unambiguously assigned to trHb1_N, trHb2_O and trHb3_P. Multiunit and chimeric forms of trHb1_N were identified exclusively in ciliated protozoan parasites, where multiple units of trHb are integrated in tandem and/or fused with another redox active or signalling domain, presenting an interesting example of gene duplication and fusion in lower eukaryotes. In contrast, trHb2_O and trHb3_P trHbs were found only in bacteria in two or multidomain organization, where amino or carboxy terminus of trHb unit is integrated with different redox-active or oxidoreductase domains. The identification of these new multiunit and chimeric trHbs and their specific phyletic distribution presents an interesting and challenging finding to explore and understand complex functionalities of these novel multidomain trHbs. © 2017 IUBMB Life, 69(7):479-488, 2017.

A quantitative model for oxygen uptake and release in a family of hemeproteins

MOTIVATION

Hemeproteins have many diverse functions that largely depend on the rate at which they uptake or release small ligands, like oxygen. These proteins have been extensively studied using either simulations or experiments, albeit only qualitatively and one or two proteins at a time.

RESULTS

We present a physical-chemical model, which uses data obtained exclusively from computer simulations, to describe the uptake and release of oxygen in a family of hemeproteins, called truncated hemoglobins (trHbs). Through a rigorous statistical analysis we demonstrate that our model successfully recaptures all the reported experimental oxygen association and dissociation kinetic rate constants, thus allowing us to establish the key factors that determine the rates at which these hemeproteins uptake and release oxygen. We found that internal tunnels as well as the distal site water molecules control ligand uptake, whereas oxygen stabilization by distal site residues controls ligand release. Because these rates largely determine the functions of these hemeproteins, these approaches will also be important tools in characterizing the trHbs members with unknown functions.

CONTACT

lboechi@ic.fcen.uba.ar

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.