1
|
Kuczynski CE, Porada CD, Atala A, Cho SS, Almeida-Porada G. Evaluating sheep hemoglobins with MD simulations as an animal model for sickle cell disease. Sci Rep 2024; 14:276. [PMID: 38168584 PMCID: PMC10761887 DOI: 10.1038/s41598-023-50707-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 12/23/2023] [Indexed: 01/05/2024] Open
Abstract
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.
Collapse
Affiliation(s)
| | | | - Anthony Atala
- Wake Forest Institute for Regenerative Medicine, Winston-Salem, NC, 27101, USA
| | - Samuel S Cho
- Department of Physics, Wake Forest University, Winston-Salem, NC, 27109, USA.
- Department of Computer Science, Wake Forest University, Winston-Salem, NC, 27109, USA.
| | | |
Collapse
|
2
|
Neto V, Victor BL, Galamba N. Cyclic Peptides as Aggregation Inhibitors for Sickle Cell Disease. J Med Chem 2023; 66:16062-16074. [PMID: 37988411 DOI: 10.1021/acs.jmedchem.3c01484] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2023]
Abstract
Sickle cell disease is a missense genetic disorder characterized by the aggregation of deoxy-HbS into helical fibers that distort erythrocytes into a sickle-like shape. Herein, we investigate, through molecular dynamics, the effect of nine 5-mer cyclic peptides (CPs), tailor-designed to block key lateral contacts of the fibers. Our results show that the CPs bind orthogonally to the main HbS pocket involved in the latter contacts, with some revealing exceedingly long residence times. These CPs display moderate to high specificity, exhibiting molecular recognition events even at a HbS/CP (1:1) ratio. A much lower HbS-CP binding free energy, longer residence times, and higher specificity are also found relative to a previously reported CP with modest in vitro antisickling activity. These results indicate that some of these CPs have the potential to reduce the concentration of aggregation-competent deoxy-HbS, precluding or delaying the formation of lateral contact at the homogeneous nucleation stage.
Collapse
Affiliation(s)
- Vasco Neto
- Biosystems and Integrative Sciences Institute, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande 1749-016, Lisboa, Portugal
| | - Bruno Lourenço Victor
- Biosystems and Integrative Sciences Institute, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande 1749-016, Lisboa, Portugal
| | - Nuno Galamba
- Biosystems and Integrative Sciences Institute, Faculdade de Ciências da Universidade de Lisboa, Edifício C8, Campo Grande 1749-016, Lisboa, Portugal
| |
Collapse
|
3
|
Olagunju MO, Loschwitz J, Olubiyi OO, Strodel B. Multiscale
MD
simulations of wild‐type and sickle hemoglobin aggregation. Proteins 2022; 90:1811-1824. [DOI: 10.1002/prot.26352] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 04/02/2022] [Accepted: 04/20/2022] [Indexed: 11/06/2022]
Affiliation(s)
- Maryam O. Olagunju
- Institute of Biological Information Processing Structural Biochemistry, Forschungszentrum Jülich Jülich Germany
| | - Jennifer Loschwitz
- Institute of Biological Information Processing Structural Biochemistry, Forschungszentrum Jülich Jülich Germany
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf Düsseldorf Germany
| | - Olujide O. Olubiyi
- Institute of Biological Information Processing Structural Biochemistry, Forschungszentrum Jülich Jülich Germany
- Department of Pharmaceutical and Medicinal Chemistry, College of Pharmacy Afe Babalola University Ado‐Ekiti Nigeria
- Institute of Drug Research and Development, Bogoro Centre Afe Babalola University Ado‐Ekiti Nigeria
| | - Birgit Strodel
- Institute of Biological Information Processing Structural Biochemistry, Forschungszentrum Jülich Jülich Germany
- Institute of Theoretical and Computational Chemistry, Heinrich Heine University Düsseldorf Düsseldorf Germany
| |
Collapse
|
4
|
Maity D, Pal D. Molecular Dynamics of Hemoglobin Reveals Structural Alterations and Explains the Interactions Driving Sickle Cell Fibrillation. J Phys Chem B 2021; 125:9921-9933. [PMID: 34459602 DOI: 10.1021/acs.jpcb.1c01684] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
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.
Collapse
Affiliation(s)
- Dibyajyoti Maity
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore 560012, Karnataka, India
| | - Debnath Pal
- Department of Computational and Data Sciences, Indian Institute of Science, Bangalore 560012, Karnataka, India
| |
Collapse
|
5
|
Jakubec D, Vondrášek J. Efficient Estimation of Absolute Binding Free Energy for a Homeodomain-DNA Complex from Nonequilibrium Pulling Simulations. J Chem Theory Comput 2020; 16:2034-2041. [PMID: 32208691 DOI: 10.1021/acs.jctc.0c00006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Estimation of binding free energies is one of the central aims of simulations of biomolecular complexes. We explore the accuracy and efficiency of setups based on nonequilibrium pulling simulations applied to the estimation of binding affinities of DNA-binding proteins. Absolute binding free energies are calculated over a range of temperatures and compared to results obtained previously using an equilibrium method. We show that realistic binding affinities can be obtained with the presented nonequilibrium approach, which also entails lower computational requirements. Errors of the binding free energy estimates are investigated and are shown to be comparable to those observed previously. Bounds are provided on the convergence of the errors with respect to the number of pulling simulations performed and with respect to the applied pull rate.
Collapse
Affiliation(s)
- David Jakubec
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry of the CAS, 166 10 Praha 6, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science, Charles University, 128 43 Praha 2, Czech Republic
| | - Jiří Vondrášek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry of the CAS, 166 10 Praha 6, Czech Republic
| |
Collapse
|
6
|
Buehler PW, Humar R, Schaer DJ. Haptoglobin Therapeutics and Compartmentalization of Cell-Free Hemoglobin Toxicity. Trends Mol Med 2020; 26:683-697. [PMID: 32589936 DOI: 10.1016/j.molmed.2020.02.004] [Citation(s) in RCA: 46] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2019] [Revised: 02/17/2020] [Accepted: 02/18/2020] [Indexed: 02/06/2023]
Abstract
Hemolysis and accumulation of cell-free hemoglobin (Hb) in the circulation or in confined tissue compartments such as the subarachnoid space is an important driver of disease. Haptoglobin is the Hb binding and clearance protein in human plasma and an efficient antagonist of Hb toxicity resulting from physiological red blood cell turnover. However, endogenous concentrations of haptoglobin are insufficient to provide protection against Hb-driven disease processes in conditions such as sickle cell anemia, sepsis, transfusion reactions, medical-device associated hemolysis, or after a subarachnoid hemorrhage. As a result, there is increasing interest in developing haptoglobin therapeutics to target 'toxic' cell-free Hb exposures. Here, we discuss key concepts of Hb toxicity and provide a perspective on the use of haptoglobin as a therapeutic protein.
Collapse
Affiliation(s)
- Paul W Buehler
- Department of Pathology, University of Maryland School of Medicine, Baltimore, MD, USA; Center for Blood Oxygen Transport and Hemostasis, Department of Pediatrics, University of Maryland School of Medicine, Baltimore, MD, USA.
| | - Rok Humar
- Division of Internal Medicine, University Hospital, Zurich, Switzerland
| | - Dominik J Schaer
- Division of Internal Medicine, University Hospital, Zurich, Switzerland.
| |
Collapse
|
7
|
Galamba N. On the Nonaggregation of Normal Adult Hemoglobin and the Aggregation of Sickle Cell Hemoglobin. J Phys Chem B 2019; 123:10735-10745. [PMID: 31747289 DOI: 10.1021/acs.jpcb.9b09727] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Sickle cell disease is a genetic disorder associated with a single mutation (Glu-β6 → Val-β6) in the β chains of hemoglobin, causing the polymerization of deoxygenated sickle cell hemoglobin (deoxy-HbS). The deoxy-HbS binding free energy was recently studied through molecular simulations, and a value of -14 ± 1 kcal mol-1 was found. Here, we studied the binding free energy of normal adult hemoglobin (deoxy-HbA), which does not polymerize at normal physiological conditions, with the aim of elucidating the importance of the presence of Val-β6 and of the absence of Glu-β6 on the aggregation of deoxy-HbS. A binding free energy of -4.4 ± 0.5 kcal mol-1 was found from a one-dimensional potential of mean force. Hydrophobic interactions are shown to represent less than 20% of the interactions in the contact interface, and despite similarly strong hydrogen-bonded ion pairs (i.e., salt bridges) and water bridged electrostatic interactions are found for deoxy-HbA and deoxy-HbS, a large repulsive potential energy is associated with Glu-β6, whereas a mild attractive potential energy is connected with Val-β6. Interestingly, Asp-β73 switches from forming a major electrostatic repulsive pair with Glu-β6 in deoxy-HbA, to forming a major attractive residue pair with Val-β6 in deoxy-HbS, consistent with the view that damping of electrostatic repulsions involving Glu-β6, namely, those associated with Asp-β73, could be responsible for the polymerization of deoxy-HbA at high potassium phosphate concentrations. Solvation analysis shows that functional groups forming salt bridges and water bridged interactions preserve a nearly intact first hydration sphere, avoiding a complete dewetting free energy penalty. These results support the view that the absence of Glu-β6 is more important than the presence of Val-β6, and that although hydrophobic effects, associated with the Val-β6 dehydration and interaction with the hydrophobic pocket in the neighbor tetramer, are important, electrostatic interactions are dominant, opposite to a picture where HbS association is driven by hydrophobic interactions.
Collapse
Affiliation(s)
- N Galamba
- Centre of Chemistry and Biochemistry and Biosystems and Integrative Sciences Institute , Faculdade de Ciências da Universidade de Lisboa , Edifício C8, Campo Grande , 1749-016 Lisboa , Portugal
| |
Collapse
|
8
|
Collective Transformation of Water between Hyperactive Antifreeze Proteins: RiAFPs. CRYSTALS 2019. [DOI: 10.3390/cryst9040188] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
We demonstrate, by molecular dynamics simulations, that water confined between a pair of insect hyperactive antifreeze proteins from the longhorn beetle Rhagium inquisitor is discontinuously expelled as the two proteins approach each other at a certain distance. The extensive striped hydrophobic–hydrophilic pattern on the surface, comprising arrays of threonine residues, enables water to form three independent ice channels through the assistance of hydroxyl groups, even at 300 K. The transformation is reminiscent of a freezing–melting transition rather than a drying transition and governs the stable protein–protein separation in the evaluation of the potential of mean force. The collectivity of water penetration or expulsion and the hysteresis in the time scale of ten nanoseconds predict a potential first-order phase transition at the limit of infinite size and provide a new framework for the water-mediated interaction between solutes.
Collapse
|
9
|
Jakubec D, Vondrášek J. Can All-Atom Molecular Dynamics Simulations Quantitatively Describe Homeodomain-DNA Binding Equilibria? J Chem Theory Comput 2019; 15:2635-2648. [PMID: 30807142 DOI: 10.1021/acs.jctc.8b01144] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
We systematically investigate the applicability of a molecular dynamics-based setup for the calculations of standard binding free energies of biologically relevant protein-DNA complexes. The free energies are extracted from a potential of mean force calculated using umbrella sampling simulations. Two protein-DNA systems derived from a homeodomain transcription factor complex are studied in order to investigate the binding of both disordered and globular proteins. Free energies and trajectories obtained using two modern molecular mechanical force fields are compared to each other and to experimental data. The temperature dependence of the calculated standard binding free energies is investigated by performing all simulations over a range of temperatures. We show that the values of standard binding free energies obtained from these simulations are overestimated compared to experimental results. Significant differences are observed between the two protein-DNA systems and between the two force fields, which are explained by different propensities to form inter- and intramolecular contacts. The number of protein-DNA contacts increases with increasing temperature, in agreement with the experimentally known temperature dependence of enthalpies of binding. However, conclusions about the temperature dependence of the standard binding free energies cannot be made with confidence, as the differences among the values are on the order of statistical uncertainty.
Collapse
Affiliation(s)
- David Jakubec
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic.,Department of Physical and Macromolecular Chemistry, Faculty of Science , Charles University , 128 43 Praha 2, Czech Republic
| | - Jiří Vondrášek
- Bioinformatics Group, Institute of Organic Chemistry and Biochemistry , Czech Academy of Sciences , 166 10 Praha 6, Czech Republic
| |
Collapse
|