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

Nitric Oxide Ameliorates the Effects of Hypoxia in Mice by Regulating Oxygen Transport by Hemoglobin

Xiaoying Zhou, Wenting Su, Quanwei Bao, Yu Cui, Xiaoxu Li, Yidong Yang, Chengzhong Yang, Chengyuan Wang, Li Jiao, Dewei Chen, and Jian Huang. Nitric oxide ameliorates the effects of hypoxia in mice by regulating oxygen transport by hemoglobin. High Alt Med Biol. 00:00-00, 2024.-Hypoxia is a common pathological and physiological phenomenon in ischemia, cancer, and strenuous exercise. Nitric oxide (NO) acts as an endothelium-derived relaxing factor in hypoxic vasodilation and serves as an allosteric regulator of hemoglobin (Hb). However, the ultimate effects of NO on the hematological system in vivo remain unknown, especially in extreme environmental hypoxia. Whether NO regulation of the structure of Hb improves oxygen transport remains unclear. Hence, we examined whether NO altered the oxygen affinity of Hb (Hb-O2 affinity) to protect extremely hypoxic mice. Mice were exposed to severe hypoxia with various concentrations of NO, and the survival time, exercise capacity, and other physical indexes were recorded. The survival time was prolonged in the 5 ppm NO (6.09 ± 1.29 minutes) and 10 ppm NO (6.39 ± 1.58 minutes) groups compared with the 0 ppm group (4.98 ± 1.23 minutes). Hypoxia of the brain was relieved, and the exercise exhaustion time was prolonged when mice inhaled 20 ppm NO (24.70 ± 6.87 minutes vs. 20.23 ± 6.51 minutes). In addition, the differences in arterial oxygen saturation (SO2%) (49.64 ± 7.29% vs. 42.90 ± 4.30%) and arteriovenous SO2% difference (25.14 ± 8.95% vs. 18.10 ± 6.90%) obviously increased. In ex vivo experiments, the oxygen equilibrium curve (OEC) left shifted as P50 decreased from 43.77 ± 2.49 mmHg (0 ppm NO) to 40.97 ± 1.40 mmHg (100 ppm NO) and 38.36 ± 2.78 mmHg (200 ppm NO). Furthermore, the Bohr effect of Hb was enhanced by the introduction of 200 ppm NO (-0.72 ± 0.062 vs.-0.65 ± 0.051), possibly allowing Hb to more easily offload oxygen in tissue at lower pH. The crystal structure reveals a greater distance between Asp94β-His146β in nitrosyl -Hb(NO-Hb), NO-HbβCSO93, and S-NitrosoHb(SNO-Hb) compared to tense Hb(T-Hb, 3.7 Å, 4.3 Å, and 5.8 Å respectively, versus 3.5 Å for T-Hb). Moreover, hydrogen bonds were less likely to form, representing a key limitation of relaxed Hb (R-Hb). Upon NO interaction with Hb, hydrogen bonds and salt bridges were less favored, facilitating relaxation. We speculated that NO ameliorated the effects of hypoxia in mice by promoting erythrocyte oxygen loading in the lung and offloading in tissues.

Identification of histidine residues that affect the T/R-state conformations of human hemoglobin using constant pH molecular dynamics simulations

Human hemoglobin (Hb) is a tetrameric protein consisting of two α and two β subunits that can adopt a low-affinity T- and high-affinity R-state conformations. Under physiological pH conditions, histidine (His) residues are the main sites for proton binding or release, and their protonation states can affect the T/R-state conformation of Hb. However, it remains unclear which His residues can effectively affect the Hb conformation. Herein, the impact of the 38 His residues of Hb on its T/R-state conformations was evaluated using constant-pH molecular dynamics (CpHMD) simulations at physiological pH while focusing on the His protonation states. Overall, the protonation states of some His residues were found to be correlated with the Hb conformation state. These residues were mainly located in the proximity of the heme (α87 and β92), and at the α1β2 and α2β1 interfaces (α89 and β97). This correlation may be partly explained by how easily hydrogen bonds can be formed, which depends on the protonation states of the His residues. Taken together, these CpHMD-based findings provide new insights into the identification of titratable His residues α87, α89, β92, and β97 that can affect Hb conformational switching under physiological pH conditions.

Regional anesthesia and sickle cell crisis in pediatric patients: An educational-focused review

INTRODUCTION

Sickle cell disease (SCD) is the most common inherited hemoglobinopathy, affecting approximately 100 000 patients in United States and millions worldwide. Although the mainstay of pain management for VOC remains systemic opioids, given the potential for adverse effects including respiratory depression and hypoxemia, there remains interest in the use of regional anesthetic techniques (neuraxial or peripheral nerve blockade).

METHODS

A systematic search of pubMed, Scopus, and Google Scholar was conducted using the terms sickle cell disease, sickle cell crisis, pain crisis, vaso-occlusive crisis, regional anesthesia, peripheral nerve blockade, and neuraxial anesthesia.

RESULTS

We identified 7 publications, all of which were retrospective case series or single case reports, outlining the use of neuraxial anesthesia in a total of 26 patients with SCD. Additionally, we identified 4 publications, including one retrospective case series and 3 single case reports, entailing the use of peripheral blockade in patients with VOC and SCD.

DISCUSSION

The available literature, albeit all retrospective or anecdotal, suggests the potential utility of regional anesthesia to treat pain in patients with SCD. Additional benefits have included avoidance of the potential deleterious physiologic effects of systemic opioids and in one case series, an improvement in respiratory function as judged by pulse oximetry. The anecdotal and retrospective nature of the available reports with an absence of prospective trials limits the evidence based medicine available from which to develop to guidlines for the optimal local anesthetic agent to use, its concentration, the rate of infusion, and the choice of adjunctive agents.

Evaluating sheep hemoglobins with MD simulations as an animal model for sickle cell disease

Sickle cell disease (SCD) affects millions worldwide, yet there are few therapeutic options. To develop effective treatments, preclinical models that recapitulate human physiology and SCD pathophysiology are needed. SCD arises from a single Glu-to-Val substitution at position 6 in the β subunit of hemoglobin (Hb), promoting Hb polymerization and subsequent disease. Sheep share important physiological and developmental characteristics with humans, including the same developmental pattern of fetal to adult Hb switching. Herein, we investigated whether introducing the SCD mutation into the sheep β-globin locus would recapitulate SCD's complex pathophysiology by generating high quality SWISS-MODEL sheep Hb structures and performing MD simulations of normal/sickle human (huHbA/huHbS) and sheep (shHbB/shHbS) Hb, establishing how accurately shHbS mimics huHbS behavior. shHbS, like huHbS, remained stable with low RMSD, while huHbA and shHbB had higher and fluctuating RMSD. shHbB and shHbS also behaved identically to huHbA and huHbS with respect to β2-Glu6 and β1-Asp73 (β1-Asn72 in sheep) solvent interactions. These data demonstrate that introducing the single SCD-causing Glu-to-Val substitution into sheep β-globin causes alterations consistent with the Hb polymerization that drives RBC sickling, supporting the development of a SCD sheep model to pave the way for alternative cures for this debilitating, globally impactful disease.

Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

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 retrospective on statistical mechanical models for hemoglobin allostery

Understanding allosteric interactions in proteins has become one of the major research areas in protein science. The original aim of the famous theoretical model of Monod, Wyman, and Changeux (MWC) was to explain the regulation of enzymatic activity in biochemical pathways. However, its first successful quantitative application was to explain cooperative oxygen binding by hemoglobin, often called the "hydrogen molecule of biology." The combination of its original application and the enormous amount of research on hemoglobin has made it the paradigm for studies of allostery, especially for multi-subunit proteins, and for the development of statistical mechanical models to describe how structure determines function. This article is a historical account of the development of statistical mechanical models for hemoglobin to explain both the cooperative binding of oxygen (called homotropic effects by MWC) and how oxygen binding is affected by ligands that bind distant from the heme oxygen binding site (called heterotropic allosteric effects by MWC). This account makes clear the many remaining challenges for describing the relationship of structure to function for hemoglobin in terms of a satisfactory statistical mechanical model.

Carbon Dioxide and the Carbamate Post-Translational Modification

Carbon dioxide is essential for life. It is at the beginning of every life process as a substrate of photosynthesis. It is at the end of every life process as the product of post-mortem decay. Therefore, it is not surprising that this gas regulates such diverse processes as cellular chemical reactions, transport, maintenance of the cellular environment, and behaviour. Carbon dioxide is a strategically important research target relevant to crop responses to environmental change, insect vector-borne disease and public health. However, we know little of carbon dioxide’s direct interactions with the cell. The carbamate post-translational modification, mediated by the nucleophilic attack by carbon dioxide on N-terminal α-amino groups or the lysine ɛ-amino groups, is one mechanism by which carbon dioxide might alter protein function to form part of a sensing and signalling mechanism. We detail known protein carbamates, including the history of their discovery. Further, we describe recent studies on new techniques to isolate this problematic post-translational modification.

Going Nuclear with Amino Acids and Proteins - Basic Biochemistry and Molecular Biology Primer for the Technologist

In recent years, there has been an influx of new tracers into the field of nuclear medicine and molecular imaging. Most of these tracers that have been FDA approved for clinical imaging exploit various mechanisms of protein biochemistry and molecular biology to bring about their actions, such as amino acid metabolism, protein folding, receptor-ligand interactions, and surface transport mechanisms. In this review, we attempt to paint a clear picture of the basic biochemistry and molecular biology of protein structure, translation, transcription, post-translational modifications, and protein targeting, in the context of the various radiopharmaceuticals currently used clinically, all in an easy-to-understand language for entry level technologists in the field. Tracer characteristics, including indications, dosage, injection-to-imaging time, and the logic behind the normal and pathophysiologic biodistribution of these newer molecular tracers, are also discussed.

Analysis of Fluctuation in the Heme-Binding Pocket and Heme Distortion in Hemoglobin and Myoglobin

Heme is located in the active site of proteins and has diverse and important biological functions, such as electron transfer and oxygen transport and/or storage. The distortion of heme porphyrin is considered an important factor for the diverse functions of heme because it correlates with the physical properties of heme, such as oxygen affinity and redox potential. Therefore, clarification of the relationship between heme distortion and the protein environment is crucial in protein science. Here, we analyzed the fluctuation in heme distortion in the protein environment for hemoglobin and myoglobin using molecular dynamics (MD) simulations and quantum mechanical (QM) calculations as well as statistical analysis of the protein structures of hemoglobin and myoglobin stored in Protein Data Bank. Our computation and statistical analysis showed that the protein environment for hemoglobin and myoglobin prominently affects the doming distortion of heme porphyrin, which correlates with its oxygen affinity, and that the magnitude of distortion is different between hemoglobin and myoglobin. These results suggest that heme distortion is affected by its protein environment and fluctuates around its fitted conformation, leading to physical properties that are appropriate for protein functions.

Dynamics of camel and human hemoglobin revealed by molecular simulations

Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

Sample Processing and Stability for Urine Biomarker Studies

BACKGROUND

Current methods of processing and storing urine samples have not been compared systematically to determine optimal conditions for advancing research on urinary biomarkers. Often, preanalytical handling is nonideal, especially considering the COVID-19 pandemic; consequently, we compared the effects of different short-term storage and processing methods on urinary biomarker measurements.

METHODS

Spot urine samples were collected via a Foley catheter from 20 hospitalized patients from the Yale New Haven Hospital within 48 hours postcardiac surgery. The effects of 3 urine storage and processing methods on biomarkers were tested: (a) 48-hour temporary storage at 4 °C prior to freezing at -80 °C, (b) 48-hour temporary storage at 25 °C prior to freezing at -80 °C, and (c) no centrifugation and immediate storage at -80 °C. Established Meso-Scale Device assay methods were used to measure the urine concentrations of 18 biomarkers: interferon gamma (IFN-ɣ), interleukin (IL)-10, IL-12p70, IL-13, IL-1β, IL-2, IL-4, IL-6, IL-8, IL-18, tumor necrosis factor alpha (TNF-α), epidermal growth factor (EGF), neutrophil gelatinase-associated lipocalin (NGAL), osteopontin (OPN), uromodulin (UMOD), kidney injury molecule-1 (KIM-1), monocyte chemoattractant protein-1 (MCP-1), and chitinase-3-like protein 1 (YKL-40).

RESULTS

Measurements of most biomarkers investigated remained stable after temporary storage at 4 °C. IL-6, IL-8, KIM-1, MCP-1, YKL-40, EGF, and NGAL were stable across all 3 processing conditions. IL-12p70 and IL-4 demonstrated significant differences in all tested conditions compared to the reference standard.

CONCLUSIONS

We identified several notable biomarkers that are robust to variations in preanalytical techniques and can be reliably investigated with nonideal handling conditions.

Molecular Dynamics of Hemoglobin Reveals Structural Alterations and Explains the Interactions Driving Sickle Cell Fibrillation

In sickle cell anemia, deoxyhemoglobin deforms RBCs by forming fibrils inside that disintegrate on oxygenation. We studied 100 ns long all-atom molecular dynamics (MD) for sickle and normal hemoglobin fibril models to understand this process, complemented by multiple 1 μs MD for a single tetramer of sickle and normal hemoglobin in deoxy and oxy states. We find that the presence of hydrophobic residues without a bulky side chain at β-6 in hemoglobin is the reason for the stability of the fibrils. Moreover, the free energy landscapes from MD of hemoglobin starting in the tensed (T) state capture the putative transition from T to relaxed (R) state, associated with oxygen binding. The three conformational wells in the landscapes are characterized by the quaternary changes where one αβ dimer rotates with respect to the other. The conformational changes from the oxygenation of sickle hemoglobin hinder the intermolecular contacts necessary for fibril formation.

Impact of hemoglobin biophysical studies on molecular pathogenesis and drug therapy for sickle cell disease

Basic research on hemoglobin has been essential for understanding the origin and treatment of many hematological disorders due to abnormal hemoglobins. The most important of the hemoglobinopathies is sickle cell disease - Linus Pauling's "molecular disease" that gave birth to molecular medicine. In this review, I will describe the contributions of basic biophysical research on normal and sickle cell hemoglobin (HbS) to understanding the molecular pathogenesis of the disease and providing the conceptual basis for the various approaches to drug therapy that target HbS polymerization. Most prominent among these are the experimental results on the solubility of HbS as a function of oxygen saturation explained by the allosteric model of Monod, Wyman, and Changeux and the Gill-Wyman thermodynamic linkage relation between solubility and oxygen binding, the solubility of mixtures of HbS with normal or fetal hemoglobin explained by Minton's thermodynamic model, and the highly unusual kinetics of HbS polymerization explained by a novel double nucleation mechanism that also accounts for the aggregation kinetics of the Alzheimer's peptide. The HbS polymerization kinetics are of great importance to understanding the pathophysiology and clinical course, as well as guiding drug development for treating this common and severe disease. The article focuses primarily on experimental and theoretical results from my lab, so it is not a comprehensive review of the subject.

Advanced Bioluminescence System for In Vivo Imaging with Brighter and Red-Shifted Light Emission

In vivo bioluminescence imaging (BLI), which is based on luminescence emitted by the luciferase-luciferin reaction, has enabled continuous monitoring of various biochemical processes in living animals. Bright luminescence with a high signal-to-background ratio, ideally red or near-infrared light as the emission maximum, is necessary for in vivo animal experiments. Various attempts have been undertaken to achieve this goal, including genetic engineering of luciferase, chemical modulation of luciferin, and utilization of bioluminescence resonance energy transfer (BRET). In this review, we overview a recent advance in the development of a bioluminescence system for in vivo BLI. We also specifically examine the improvement in bioluminescence intensity by mutagenic or chemical modulation on several beetle and marine luciferase bioluminescence systems. We further describe that intramolecular BRET enhances luminescence emission, with recent attempts for the development of red-shifted bioluminescence system, showing great potency in in vivo BLI. Perspectives for future improvement of bioluminescence systems are discussed.

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 CO₂. 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 CO₂, which would permit effective O₂ unloading to tissues in response to an increase in metabolism during intense activity and diving. Although the Hbs of all crocodilians exhibit similar intrinsic O₂-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 NH₂-terminal acetylation. Whereas chloride appears to be a potent allosteric effector of all crocodile Hbs, ATP has a strong, chloride-independent effect on Hb-O₂ affinity only in caimans. Modeling suggests that allosteric ATP binding has a somewhat different structural basis in crocodilian and mammalian Hbs.

Curating the gnomAD database: Report of novel variants in the globin-coding genes and bioinformatics analysis

Massive parallel sequencing technologies are facilitating the faster identification of sequence variants with the consequent capability of untangling the molecular bases of many human genetic syndromes. However, it is not always easy to understand the impact of novel variants, especially for missense changes, which can lead to a spectrum of phenotypes. This study presents a custom-designed multistep methodology to evaluate the impact of novel variants aggregated in the genome aggregation database for the HBB, HBA2, and HBA1 genes, by testing and improving its performance with a dataset of previously described alterations affecting those same genes. This approach scored high sensitivity and specificity values and showed an overall better performance than sequence-derived predictors, highlighting the importance of protein conformation and interaction specific analyses in curating variant databases. This study also describes the strengths and limitations of these structural studies and allows identifying residues in the globin chains more prone to tolerate substitutions.

The importance of the quaternary structure to represent conformational ensembles of the major Mycobacterium tuberculosis drug target

Flexibility is a feature intimately related to protein function, since conformational changes can be used to describe environmental changes, chemical modifications, protein-protein and protein-ligand interactions. In this study, we have investigated the influence of the quaternary structure of 2-trans-enoyl-ACP (CoA) reductase or InhA, from Mycobacterium tuberculosis, to its flexibility. We carried out classical molecular dynamics simulations using monomeric and tetrameric forms to elucidate the enzyme's flexibility. Overall, we observed statistically significant differences between conformational ensembles of tertiary and quaternary structures. In addition, the enzyme's binding site is the most affected region, reinforcing the importance of the quaternary structure to evaluate the binding affinity of small molecules, as well as the effect of single point mutations to InhA protein dynamics.

Investigation of quaternary structure of aggregating 3-ketosteroid dehydrogenase from Sterolibacterium denitrificans: In the pursuit of consensus of various biophysical techniques

In this work we analyzed the quaternary structure of FAD-dependent 3-ketosteroid dehydrogenase (AcmB) from Sterolibacterium denitrificans, the protein that in solution forms massive aggregates (>600 kDa). Using size-excursion chromatography (SEC), dynamic light scattering (DLS), native-PAGE and atomic force microscopy (AFM) we studied the nature of enzyme aggregation. Partial protein de-aggregation was facilitated by the presence of non-ionic detergent such as Tween 20 or by a high degree of protein dilution but not by addition of a reducing agent or an increase of ionic strength. De-aggregating influence of Tween 20 had no impact on either enzyme's specific activity or FAD reconstitution to recombinant AcmB. The joint experimental (DLS, isoelectric focusing) and theoretical investigations demonstrated gradual shift of enzyme's isoelectric point upon aggregation from 8.6 for a monomeric form to even 5.0. The AFM imaging on mica or highly oriented pyrolytic graphite (HOPG) surface enabled observation of individual protein monomers deposited from a highly diluted solution (0.2 μg/ml). Such approach revealed that native AcmB can indeed be monomeric. AFM imaging supported by theoretical random sequential adsorption (RSA) kinetics allowed estimation of distribution enzyme forms in the bulk solution: 5%, monomer, 11.4% dimer and 12% trimer. Finally, based on results of AFM as well as analysis of the surface of AcmB homology models we have observed that aggregation is most probably initiated by hydrophobic forces and then assisted by electrostatic attraction between negatively charged aggregates and positively charged monomers.

Lessons learned about steered molecular dynamics simulations and free energy calculations

The calculation of free energy profiles is central in understanding differential enzymatic activity, for instance, involving chemical reactions that require QM-MM tools, ligand migration, and conformational rearrangements that can be modeled using classical potentials. The use of steered molecular dynamics (sMD) together with the Jarzynski equality is a popular approach in calculating free energy profiles. Here, we first briefly review the application of the Jarzynski equality to sMD simulations, then revisit the so-called stiff-spring approximation and the consequent expectation of Gaussian work distributions and, finally, reiterate the practical utility of the second-order cumulant expansion, as it coincides with the parametric maximum-likelihood estimator in this scenario. We illustrate this procedure using simulations of CO, both in aqueous solution and in a carbon nanotube as a model system for biologically relevant nanoheterogeneous environments. We conclude the use of the second-order cumulant expansion permits the use of faster pulling velocities in sMD simulations, without introducing bias due to large dispersion in the non-equilibrium work distribution.

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.