Catalytic Mechanism of Mycobacterium tuberculosis Methionine Sulfoxide Reductase A.

The oxidation of Met to methionine sulfoxide (MetSO) by oxidants such as hydrogen peroxide, hypochlorite, or peroxynitrite has profound effects on protein function. This modification can be reversed by methionine sulfoxide reductases (msr). In the context of pathogen infection, the reduction of oxidized proteins gains significance due to microbial oxidative damage generated by the immune system. For example, Mycobacterium tuberculosis (Mt) utilizes msrs (MtmsrA and MtmsrB) as part of the repair response to the host-induced oxidative stress. The absence of these enzymes makes Mycobacteria prone to increased susceptibility to cell death, pointing them out as potential therapeutic targets. This study provides a detailed characterization of the catalytic mechanism of MtmsrA using a comprehensive approach, including experimental techniques and theoretical methodologies. Confirming a ping-pong type enzymatic mechanism, we elucidate the catalytic parameters for sulfoxide and thioredoxin substrates (kcat/KM = 2656 ± 525 M-1 s-1 and 1.7 ± 0.8 × 106 M-1 s-1, respectively). Notably, the entropic nature of the activation process thermodynamics, representing ∼85% of the activation free energy at room temperature, is underscored. Furthermore, the current study questions the plausibility of a sulfurane intermediate, which may be a transition-state-like structure, suggesting the involvement of a conserved histidine residue as an acid-base catalyst in the MetSO reduction mechanism. This mechanistic insight not only advances our understanding of Mt antioxidant enzymes but also holds implications for future drug discovery and biotechnological applications.

Multiple oxidative post-translational modifications of human glutamine synthetase mediate peroxynitrite-dependent enzyme inactivation and aggregation.

Glutamine synthetase (GS), which catalyzes the ATP-dependent synthesis of L-glutamine from L-glutamate and ammonia, is a ubiquitous and conserved enzyme that plays a pivotal role in nitrogen metabolism across all life domains. In vertebrates, GS is highly expressed in astrocytes, where its activity sustains the glutamate-glutamine cycle at glutamatergic synapses and is thus essential for maintaining brain homeostasis. In fact, decreased GS levels or activity have been associated with neurodegenerative diseases, with these alterations attributed to oxidative post-translational modifications of the protein, in particular tyrosine nitration. In this study, we expressed and purified human GS (HsGS) and performed an in-depth analysis of its oxidative inactivation by peroxynitrite (ONOO-) in vitro. We found that ONOO- exposure led to a dose-dependent loss of HsGS activity, the oxidation of cysteine, methionine, and tyrosine residues and also the nitration of tryptophan and tyrosine residues. Peptide mapping by LC-MS/MS through combined H216O/H218O trypsin digestion identified up to 10 tyrosine nitration sites and five types of dityrosine cross-links; these modifications were further scrutinized by structural analysis. Tyrosine residues 171, 185, 269, 283, and 336 were the main nitration targets; however, tyrosine-to-phenylalanine HsGS mutants revealed that their sole nitration was not responsible for enzyme inactivation. In addition, we observed that ONOO- induced HsGS aggregation and activity loss. Thiol oxidation was a key modification to elicit aggregation, as it was also induced by hydrogen peroxide treatment. Taken together, our results indicate that multiple oxidative events at various sites are responsible for the inactivation and aggregation of human GS.

Mechanism of Peroxynitrite Interaction with Ferric M. tuberculosis Nitrobindin: A Computational Study.

Nitrobindins (Nbs) are all-ß-barrel heme proteins present along the evolutionary ladder. They display a highly solvent-exposed ferric heme group with the iron atom being coordinated by the proximal His residue and a water molecule at the distal position. Ferric nitrobindins (Nb(III)) play a role in the conversion of toxic peroxynitrite (ONOO-) to harmless nitrate, with the value of the second-order rate constant being similar to those of most heme proteins. The value of the second-order rate constant of Nbs increases as the pH decreases; this suggests that Nb(III) preferentially reacts with peroxynitrous acid (ONOOH), although ONOO- is more nucleophilic. In this work, we shed light on the molecular basis of the ONOO- and ONOOH reactivity of ferric Mycobacterium tuberculosis Nb (Mt-Nb(III)) by dissecting the ligand migration toward the active site, the water molecule release, and the ligand binding process by computer simulations. Classical molecular dynamics simulations were performed by employing a steered molecular dynamics approach and the Jarzynski equality to obtain ligand migration free energy profiles for both ONOO- and ONOOH. Our results indicate that ONOO- and ONOOH migration is almost unhindered, consistent with the exposed metal center of Mt-Nb(III). To further analyze the ligand binding process, we computed potential energy profiles for the displacement of the Fe(III)-coordinated water molecule using a hybrid QM/MM scheme at the DFT level and a nudged elastic band approach. These results indicate that ONOO- exhibits a much larger barrier for ligand displacement than ONOOH, suggesting that water displacement is assisted by protonation of the leaving group by the incoming ONOOH.

Assessment of Embedding Schemes in a Hybrid Machine Learning/Classical Potentials (ML/MM) Approach.

Machine learning (ML) methods have reached high accuracy levels for the prediction of in vacuo molecular properties. However, the simulation of large systems solely through ML methods (such as those based on neural network potentials) is still a challenge. In this context, one of the most promising frameworks for integrating ML schemes in the simulation of complex molecular systems are the so-called ML/MM methods. These multiscale approaches combine ML methods with classical force fields (MM), in the same spirit as the successful hybrid quantum mechanics-molecular mechanics methods (QM/MM). The key issue for such ML/MM methods is an adequate description of the coupling between the region of the system described by ML and the region described at the MM level. In the context of QM/MM schemes, the main ingredient of the interaction is electrostatic, and the state of the art is the so-called electrostatic-embedding. In this study, we analyze the quality of simpler mechanical embedding-based approaches, specifically focusing on their application within a ML/MM framework utilizing atomic partial charges derived in vacuo. Taking as reference electrostatic embedding calculations performed at a QM(DFT)/MM level, we explore different atomic charges schemes, as well as a polarization correction computed using atomic polarizabilites. Our benchmark data set comprises a set of about 80k small organic structures from the ANI-1x and ANI-2x databases, solvated in water. The results suggest that the minimal basis iterative stockholder (MBIS) atomic charges yield the best agreement with the reference coupling energy. Remarkable enhancements are achieved by including a simple polarization correction.

Insights into the product release mechanism of dengue virus NS3 helicase.

The non-structural protein 3 helicase (NS3h) is a multifunctional protein that is critical in RNA replication and other stages in the flavivirus life cycle. NS3h uses energy from ATP hydrolysis to translocate along single stranded nucleic acid and to unwind double stranded RNA. Here we present a detailed mechanistic analysis of the product release stage in the catalytic cycle of the dengue virus (DENV) NS3h. This study is based on a combined experimental and computational approach of product-inhibition studies and free energy calculations. Our results support a model in which the catalytic cycle of ATP hydrolysis proceeds through an ordered sequential mechanism that includes a ternary complex intermediate (NS3h-Pi-ADP), which evolves releasing the first product, phosphate (Pi), and subsequently ADP. Our results indicate that in the product release stage of the DENV NS3h a novel open-loop conformation plays an important role that may be conserved in NS3 proteins of other flaviviruses as well.

Autocatalytic Mechanism in the Anaerobic Reduction of Metmyoglobin by Sulfide Species.

The mechanism of the metal centered reduction of metmyoglobin (MbFeIII) by sulfide species (H2S/HS-) under an argon atmosphere has been studied by a combination of spectroscopic, kinetic, and computational methods. Asymmetric S-shaped time-traces for the formation of MbFeII at varying ratios of excess sulfide were observed at pH 5.3 < pH < 8.0 and 25 °C, suggesting an autocatalytic reaction mechanism. An increased rate at more alkaline pHs points to HS- as relevant reactive species for the reduction. The formation of the sulfanyl radical (HS•) in the slow initial phase was assessed using the spin-trap phenyl N-tert-butyl nitrone. This radical initiates the formation of S-S reactive species as disulfanuidyl/ disulfanudi-idyl radical anions and disulfide (HSSH•-/HSS•2- and HSS-, respectively). The autocatalysis has been ascribed to HSS-, formed after HSSH•-/HSS•2- disproportionation, which behaves as a fast reductant toward the intermediate complex MbFeIII(HS-). We propose a reaction mechanism for the sulfide-mediated reduction of metmyoglobin where only ferric heme iron initiates the oxidation of sulfide species. Beside the chemical interest, this insight into the MbFeIII/sulfide reaction under an argon atmosphere is relevant for the interpretation of biochemical aspects of ectopic myoglobins found on hypoxic tissues toward reactive sulfur species.

Novel Lennard-Jones Parameters for Cysteine and Selenocysteine in the AMBER Force Field.

Cysteine is a common amino acid with a thiol group that plays a pivotal role in a variety of scenarios in redox biochemistry. In contrast, selenocysteine, the 21st amino acid, is only present in 25 human proteins. Classical force-field parameters for cysteine and selenocysteine are still scarce. In this context, we present a methodology to obtain Lennard-Jones parameters for cysteine and selenocysteine in different physiologically relevant oxidation and protonation states. The new force field parameters obtained in this work are available at https://github.com/MALBECC/AMBER-parameters-database. The parameters were adjusted to reproduce water radial distribution functions obtained by density functional theory ab initio molecular dynamics. We validated the results by evaluating the impact of the choice of parameters on the structure and dynamics in classical molecular dynamics simulations of representative proteins containing catalytic cysteine/selenocysteine residues. There are significant changes in protein structure and dynamics depending on the parameters choice, specifically affecting the residues close to the catalytic sites.

A Unique Mode of Coenzyme A Binding to the Nucleotide Binding Pocket of Human Metastasis Suppressor NME1.

Coenzyme A (CoA) is a key cellular metabolite which participates in diverse metabolic pathways, regulation of gene expression and the antioxidant defense mechanism. Human NME1 (hNME1), which is a moonlighting protein, was identified as a major CoA-binding protein. Biochemical studies showed that hNME1 is regulated by CoA through both covalent and non-covalent binding, which leads to a decrease in the hNME1 nucleoside diphosphate kinase (NDPK) activity. In this study, we expanded the knowledge on previous findings by focusing on the non-covalent mode of CoA binding to the hNME1. With X-ray crystallography, we solved the CoA bound structure of hNME1 (hNME1-CoA) and determined the stabilization interactions CoA forms within the nucleotide-binding site of hNME1. A hydrophobic patch stabilizing the CoA adenine ring, while salt bridges and hydrogen bonds stabilizing the phosphate groups of CoA were observed. With molecular dynamics studies, we extended our structural analysis by characterizing the hNME1-CoA structure and elucidating possible orientations of the pantetheine tail, which is absent in the X-ray structure due to its flexibility. Crystallographic studies suggested the involvement of arginine 58 and threonine 94 in mediating specific interactions with CoA. Site-directed mutagenesis and CoA-based affinity purifications showed that arginine 58 mutation to glutamate (R58E) and threonine 94 mutation to aspartate (T94D) prevent hNME1 from binding to CoA. Overall, our results reveal a unique mode by which hNME1 binds CoA, which differs significantly from that of ADP binding: the α- and ß-phosphates of CoA are oriented away from the nucleotide-binding site, while 3'-phosphate faces catalytic histidine 118 (H118). The interactions formed by the CoA adenine ring and phosphate groups contribute to the specific mode of CoA binding to hNME1.

Acidity and nucleophilic reactivity of glutathione persulfide.

Persulfides (RSSH/RSS-) participate in sulfur trafficking and metabolic processes, and are proposed to mediate the signaling effects of hydrogen sulfide (H2S). Despite their growing relevance, their chemical properties are poorly understood. Herein, we studied experimentally and computationally the formation, acidity, and nucleophilicity of glutathione persulfide (GSSH/GSS-), the derivative of the abundant cellular thiol glutathione (GSH). We characterized the kinetics and equilibrium of GSSH formation from glutathione disulfide and H2S. A pKa of 5.45 for GSSH was determined, which is 3.49 units below that of GSH. The reactions of GSSH with the physiologically relevant electrophiles peroxynitrite and hydrogen peroxide, and with the probe monobromobimane, were studied and compared with those of thiols. These reactions occurred through SN2 mechanisms. At neutral pH, GSSH reacted faster than GSH because of increased availability of the anion and, depending on the electrophile, increased reactivity. In addition, GSS- presented higher nucleophilicity with respect to a thiolate with similar basicity. This can be interpreted in terms of the so-called α effect, i.e. the increased reactivity of a nucleophile when the atom adjacent to the nucleophilic atom has high electron density. The magnitude of the α effect correlated with the Brønsted nucleophilic factor, ßnuc, for the reactions with thiolates and with the ability of the leaving group. Our study constitutes the first determination of the pKa of a biological persulfide and the first examination of the α effect in sulfur nucleophiles, and sheds light on the chemical basis of the biological properties of persulfides.

Catalysis of Peroxide Reduction by Fast Reacting Protein Thiols.

Life on Earth evolved in the presence of hydrogen peroxide, and other peroxides also emerged before and with the rise of aerobic metabolism. They were considered only as toxic byproducts for many years. Nowadays, peroxides are also regarded as metabolic products that play essential physiological cellular roles. Organisms have developed efficient mechanisms to metabolize peroxides, mostly based on two kinds of redox chemistry, catalases/peroxidases that depend on the heme prosthetic group to afford peroxide reduction and thiol-based peroxidases that support their redox activities on specialized fast reacting cysteine/selenocysteine (Cys/Sec) residues. Among the last group, glutathione peroxidases (GPxs) and peroxiredoxins (Prxs) are the most widespread and abundant families, and they are the leitmotif of this review. After presenting the properties and roles of different peroxides in biology, we discuss the chemical mechanisms of peroxide reduction by low molecular weight thiols, Prxs, GPxs, and other thiol-based peroxidases. Special attention is paid to the catalytic properties of Prxs and also to the importance and comparative outlook of the properties of Sec and its role in GPxs. To finish, we describe and discuss the current views on the activities of thiol-based peroxidases in peroxide-mediated redox signaling processes.

Structure and function of crocodilian hemoglobins and allosteric regulation by chloride, ATP, and CO2.

Hemoglobins (Hbs) of crocodilians are reportedly characterized by unique mechanisms of allosteric regulatory control, but there are conflicting reports regarding the importance of different effectors, such as chloride ions, organic phosphates, and CO2. Progress in understanding the unusual properties of crocodilian Hbs has also been hindered by a dearth of structural information. Here, we present the first comparative analysis of blood properties and Hb structure and function in a phylogenetically diverse set of crocodilian species. We examine mechanisms of allosteric regulation in the Hbs of 13 crocodilian species belonging to the families Crocodylidae and Alligatoridae. We also report new amino acid sequences for the α- and ß-globins of these taxa, which, in combination with structural analyses, provide insights into molecular mechanisms of allosteric regulation. All crocodilian Hbs exhibited a remarkably strong sensitivity to CO2, which would permit effective O2 unloading to tissues in response to an increase in metabolism during intense activity and diving. Although the Hbs of all crocodilians exhibit similar intrinsic O2-affinities, there is considerable variation in sensitivity to Cl- ions and ATP, which appears to be at least partly attributable to variation in the extent of NH2-terminal acetylation. Whereas chloride appears to be a potent allosteric effector of all crocodile Hbs, ATP has a strong, chloride-independent effect on Hb-O2 affinity only in caimans. Modeling suggests that allosteric ATP binding has a somewhat different structural basis in crocodilian and mammalian Hbs.

In Silico Insight into the Reductive Nitrosylation of Ferric Hemeproteins.

A combination of in silico methods was used to extend the experimental description of the reductive nitrosylation mechanism in ferric hemeproteins with the molecular details of the role of surrounding amino acids. The computational strategy consisted in the estimation of potential energy profiles for the transition process associated with the interactions of the coordinated N(NO) moiety with O(H2O) or O(OH-) as nucleophiles, and with distal amino acids as proton acceptors or affecting the stability of transition states. We inspected the reductive nitrosylation in three representative hemeproteins -sperm whale metmyoglobin, α subunit of human methemoglobin and nitrophorin 4 of Rhodnius prolixus. For each case, classical molecular dynamics simulations were performed in order to obtain relevant reactive conformations, and a potential energy profile for the reactive step was obtained using adiabatic mapping or nudged elastic band approaches at the QM/MM level. Specifically, we report the role of a charged Arg45 of myoglobin in destabilizing the transition state when H2O acts as nucleophile, differently to the neutral Pro43 of the hemoglobin subunit. The case of the nitrophorin is unique in that the access of the required water molecules is scarce, thus, preventing the reaction.

Spectroscopy in Complex Environments from QM-MM Simulations.

The applications of multiscale quantum-classical (QM-MM) approaches have shown an extraordinary expansion and diversification in the last couple of decades. A great proportion of these efforts have been devoted to interpreting and reproducing spectroscopic experiments in a variety of complex environments such as solutions, interfaces, and biological systems. Today, QM-MM-based computational spectroscopy methods constitute accomplished tools with refined predictive power. The present review summarizes the advances that have been made in QM-MM approaches to UV-visible, Raman, IR, NMR, electron paramagnetic resonance, and Mössbauer spectroscopies, providing in every case an introductory discussion of the corresponding methodological background. A representative number of applications are presented to illustrate the historical evolution and the state of the art of this field, highlighting the advantages and limitations of the available methodologies. Finally, we present our view of the perspectives and open challenges in the field.

Kinetics, subcellular localization, and contribution to parasite virulence of a Trypanosoma cruzi hybrid type A heme peroxidase (TcAPx-CcP).

The Trypanosoma cruzi ascorbate peroxidase is, by sequence analysis, a hybrid type A member of class I heme peroxidases [TcAPx-cytochrome c peroxidase (CcP)], suggesting both ascorbate (Asc) and cytochrome c (Cc) peroxidase activity. Here, we show that the enzyme reacts fast with H2O2 (k = 2.9 × 107 M-1⋅s-1) and catalytically decomposes H2O2 using Cc as the reducing substrate with higher efficiency than Asc (kcat/Km = 2.1 × 105 versus 3.5 × 104 M-1⋅s-1, respectively). Visible-absorption spectra of purified recombinant TcAPx-CcP after H2O2 reaction denote the formation of a compound I-like product, characteristic of the generation of a tryptophanyl radical-cation (Trp233•+). Mutation of Trp233 to phenylalanine (W233F) completely abolishes the Cc-dependent peroxidase activity. In addition to Trp233•+, a Cys222-derived radical was identified by electron paramagnetic resonance spin trapping, immunospin trapping, and MS analysis after equimolar H2O2 addition, supporting an alternative electron transfer (ET) pathway from the heme. Molecular dynamics studies revealed that ET between Trp233 and Cys222 is possible and likely to participate in the catalytic cycle. Recognizing the ability of TcAPx-CcP to use alternative reducing substrates, we searched for its subcellular localization in the infective parasite stages (intracellular amastigotes and extracellular trypomastigotes). TcAPx-CcP was found closely associated with mitochondrial membranes and, most interestingly, with the outer leaflet of the plasma membrane, suggesting a role at the host-parasite interface. TcAPx-CcP overexpressers were significantly more infective to macrophages and cardiomyocytes, as well as in the mouse model of Chagas disease, supporting the involvement of TcAPx-CcP in pathogen virulence as part of the parasite antioxidant armamentarium.

Naturally occurring fluorescence in frogs.

Fluorescence, the absorption of short-wavelength electromagnetic radiation reemitted at longer wavelengths, has been suggested to play several biological roles in metazoans. This phenomenon is uncommon in tetrapods, being restricted mostly to parrots and marine turtles. We report fluorescence in amphibians, in the tree frog Hypsiboas punctatus, showing that fluorescence in living frogs is produced by a combination of lymph and glandular emission, with pigmentary cell filtering in the skin. The chemical origin of fluorescence was traced to a class of fluorescent compounds derived from dihydroisoquinolinone, here named hyloins. We show that fluorescence contributes 18-29% of the total emerging light under twilight and nocturnal scenarios, largely enhancing brightness of the individuals and matching the sensitivity of night vision in amphibians. These results introduce an unprecedented source of pigmentation in amphibians and highlight the potential relevance of fluorescence in visual perception in terrestrial environments.

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 O2 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.

Evaluation of quinoxaline compounds as ligands of a site adjacent to S2 (AS2) of cruzain.

The binding of ten quinoxaline compounds (1-10) to a site adjacent to S2 (AS2) of cruzain (CRZ) was evaluated by a protocol that include a first analysis through docking experiments followed by a second analysis using the Molecular Mechanics-Poisson-Boltzmann Surface Area method (MM-PBSA). Through them we demonstrated that quinoxaline compounds bearing substituents of different sizes at positions 3 or 4 of the heterocyclic ring might interact with the AS2, particularly interesting site for drug design. These compounds showed docking scores (ΔGdock) which were similar to those estimated for inhibitors that bind to the enzyme through non-covalent interactions. Nevertheless, the free binding energies (ΔG) values estimated by MM-PBSA indicated that the derivatives 8-10, which bear bulky substituents at position 3 of the heterocycle ring, became detached from the binding site under a dynamic study. Surprisingly, the evaluation of the inhibitory activity of cruzipain (CZ) of some derivatives showed that they increase the enzymatic activity. These results lead us to conclude about the relevance of AS2 as a pocket for compounds binding site, but not necessarily for the design of anti-chagasic compounds.

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.

Mechanism of Sulfide Binding by Ferric Hemeproteins.

The reaction of hydrogen sulfide (H2S) with hemeproteins is a key physiological reaction; still, its mechanism and implications are not completely understood. In this work, we propose a combination of experimental and theoretical tools to shed light on the reaction in model system microperoxidase 11 (MP11-FeIII) and myoglobin (Mb-FeIII), from the estimation of the intrinsic binding constants of the species H2S and hydrosulfide (HS-), and the computational description of the overall binding process. Our results show that H2S and HS- are the main reactive species in Mb-FeIII and MP11-FeIII, respectively, and that the magnitude of their intrinsic binding constants are similar to most of the binding constants reported so far for hemeproteins systems and model compounds. However, while the binding of HS- to Mb-FeIII was negligible, the binding of H2S to MP11-FeIII was significant, providing a frame for a discriminated analysis of both species and revealing differential mechanistic aspects. A joint inspection of the kinetic data and the free energy profiles of the binding processes suggests that a dissociative mechanism with the release of a coordinated water molecule as rate limiting step is operative in the binding of H2S to Mb-FeIII and that the binding of HS- is prevented in the access to the protein matrix. For the MP11-FeIII case, where no access restrictions for the ligands are present, an associative component in the mechanism seems to be operative. Overall, the results suggest that if accessing the active site then both H2S and HS- are capable of binding a ferric heme moiety.

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.