Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis

Discovery of potent DNMT1 inhibitors against sickle cell disease using structural-based virtual screening, MM-GBSA and molecular dynamics simulation-based approaches

Sickle cell disease (SCD) is an autosomal recessive genetic disorder affecting millions of people worldwide. A reversible and selective DNMT1 inhibitor, GSK3482364, has been known to decrease the overall methylation activity of DNMT1, resulting in the increase of HbF levels and percentage of HbF-expressing erythrocytes in an in vitro and in vivo model. In this study, a structure-based virtual screening was done with GSK3685032, a co-crystalized ligand of DNMT1 (PDB ID: 6X9K) with an IC50 value of 0.036 μM and identified 3988 compounds from three databases (ChEMBL, PubChem and Drug Bank). Using this screening method, we identified around 15 compounds with XP docking scores greater than -8 kcal/mol. Further, prime MM-GBSA calculations have been performed and found compound SCHEMBL19716714 with the highest binding free energy of -83.31 kcal/mol. Finally, four compounds were identified based on glide energy and ΔG bind scores that have the most binding with DG7, DG19, DG20 bases and Lys1535, His1507, Trp1510, Ser1230, which were required for the target enzyme inhibition. Furthermore, molecular dynamics simulation studies of top ligands validate the stability of the docked complexes by examining root mean square deviations, root mean square fluctuations, solvent accessible surface area, and radius of gyration graphs from simulation trajectories. These findings suggest that the top four hit compounds may be capable of inhibiting DNMT1 and that additional in vitro and in vivo studies will be essential to prove the clinical effectiveness of the selected lead compounds.Communicated by Ramaswamy H. Sarma.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Metabolite genome-wide association studies (mGWAS) have advanced our understanding of the genetic control of metabolite levels. However, interpreting these associations remains challenging due to a lack of tools to annotate gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we introduce the shortest reactional distance (SRD) metric, drawing from the comprehensive KEGG database, to enhance the biological interpretation of mGWAS results. We applied this approach to three independent mGWAS, including a case study on sickle cell disease patients. Our analysis reveals an enrichment of small SRD values in reported mGWAS pairs, with SRD values significantly correlating with mGWAS p values, even beyond the standard conservative thresholds. We demonstrate the utility of SRD annotation in identifying potential false negatives and inaccuracies within current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs, suitable to integrate statistical evidence to biological networks.

Metabolic signatures of cardiorenal dysfunction in plasma from sickle cell patients as a function of therapeutic transfusion and hydroxyurea treatment

Metabolomics studies in sickle cell disease (SCD) have been so far limited to tens of samples, owing to technical and experimental limitations. To overcome these limitations, we performed plasma metabolomics analyses on 596 samples from patients with SCD enrolled in the WALK-PHaSST study (clinicaltrials gov. Identifier: NCT00492531). Clinical covariates informed the biological interpretation of metabolomics data, including genotypes (hemoglobin [Hb] SS, hemoglobin SC), history of recent transfusion (HbA%), response to hydroxyurea treatment (fetal Hb%). We investigated metabolic correlates to the degree of intravascular hemolysis, cardiorenal function, as determined by tricuspid regurgitation velocity (TRV), estimated glomerular filtration rate (eGFR), and overall hazard ratio (unadjusted or adjusted by age). Recent transfusion events or hydroxyurea treatment were associated with elevation in plasma-free fatty acids and decreases in acyl-carnitines, urate, kynurenine, indoles, carboxylic acids, and glycine- or taurine-conjugated bile acids. High levels of these metabolites, along with low levels of plasma S1P and L-arginine were identified as top markers of hemolysis, cardiorenal function (TRV, eGFR), and overall hazard ratio. We thus uploaded all omics and clinical data on a novel online portal that we used to identify a potential mechanism of dysregulated red cell S1P synthesis and export as a contributor to the more severe clinical manifestations in patients with the SS genotype compared to SC. In conclusion, plasma metabolic signatures - including low S1P, arginine and elevated kynurenine, acyl-carnitines and bile acids - are associated with clinical manifestation and therapeutic efficacy in SCD patients, suggesting new avenues for metabolic interventions in this patient population.

In vivo evaluation of the effect of sickle cell hemoglobin S, C and therapeutic transfusion on erythrocyte metabolism and cardiorenal dysfunction

Despite a wealth of exploratory plasma metabolomics studies in sickle cell disease (SCD), no study to date has evaluate a large and well phenotyped cohort to compare the primary erythrocyte metabolome of hemoglobin SS, SC and transfused AA red blood cells (RBCs) in vivo. The current study evaluates the RBC metabolome of 587 subjects with sickle cell sickle cell disease (SCD) from the WALK-PHaSST clinical cohort. The set includes hemoglobin SS, hemoglobin SC SCD patients, with variable levels of HbA related to RBC transfusion events. Here we explore the modulating effects of genotype, age, sex, severity of hemolysis, and transfusion therapy on sickle RBC metabolism. Results show that RBCs from patients with Hb SS genotypes-compared to AA RBCs from recent transfusion events or SC RBCs-are characterized by significant alterations of RBC acylcarnitines, pyruvate, sphingosine 1-phosphate, creatinine, kynurenine and urate metabolism. Surprisingly, the RBC metabolism of SC RBCs is dramatically different from SS, with all glycolytic intermediates significantly elevated in SS RBCs, with the exception of pyruvate. This result suggests a metabolic blockade at the ATP-generating phosphoenolpyruvate to pyruvate step of glycolysis, which is catalyzed by redox-sensitive pyruvate kinase. Metabolomics, clinical and hematological data were collated in a novel online portal. In conclusion, we identified metabolic signatures of HbS RBCs that correlate with the degree of steady state hemolytic anemia, cardiovascular and renal dysfunction and mortality.

Metabolic signatures of cardiorenal dysfunction in plasma from sickle cell patients, as a function of therapeutic transfusion and hydroxyurea treatment

Metabolomics studies in sickle cell disease (SCD) have been so far limited to tens of samples, owing to technical and experimental limitations. To overcome these limitations, we performed plasma metabolomics analyses on 596 samples from patients with sickle cell sickle cell disease (SCD) enrolled in the WALK-PHaSST study. Clinical covariates informed the biological interpretation of metabolomics data, including genotypes (hemoglobin SS, hemoglobin SC), history of recent transfusion (HbA%), response to hydroxyurea treatment (HbF%). We investigated metabolic correlates to the degree of hemolysis, cardiorenal function, as determined by tricuspid regurgitation velocity (TRV), estimated glomerular filtration rate (eGFR), and overall hazard ratio (unadjusted or adjusted by age). Recent transfusion events or hydroxyurea treatment were associated with elevation in plasma free fatty acids and decreases in acyl-carnitines, urate, kynurenine, indoles, carboxylic acids, and glycine- or taurine-conjugated bile acids. High levels of these metabolites, along with low levels of plasma S1P and L-arginine were identified as top markers of hemolysis, cardiorenal function (TRV, eGFR), and overall hazard ratio. We thus uploaded all omics and clinical data on a novel online portal that we used to identify a potential mechanism of dysregulated red cell S1P synthesis and export as a contributor to the more severe clinical manifestations in patients with the SS genotype compared to SC. In conclusion, plasma metabolic signatures - including low S1P, arginine and elevated kynurenine, acyl-carnitines and bile acids - are associated with clinical manifestation and therapeutic efficacy in SCD patients, suggesting new avenues for metabolic interventions in this patient population.

Gene-metabolite annotation with shortest reactional distance enhances metabolite genome-wide association studies results

Studies combining metabolomics and genetics, known as metabolite genome-wide association studies (mGWAS), have provided valuable insights into our understanding of the genetic control of metabolite levels. However, the biological interpretation of these associations remains challenging due to a lack of existing tools to annotate mGWAS gene-metabolite pairs beyond the use of conservative statistical significance threshold. Here, we computed the shortest reactional distance (SRD) based on the curated knowledge of the KEGG database to explore its utility in enhancing the biological interpretation of results from three independent mGWAS, including a case study on sickle cell disease patients. Results show that, in reported mGWAS pairs, there is an excess of small SRD values and that SRD values and p-values significantly correlate, even beyond the standard conservative thresholds. The added-value of SRD annotation is shown for identification of potential false negative hits, exemplified by the finding of gene-metabolite associations with SRD ≤1 that did not reach standard genome-wide significance cut-off. The wider use of this statistic as an mGWAS annotation would prevent the exclusion of biologically relevant associations and can also identify errors or gaps in current metabolic pathway databases. Our findings highlight the SRD metric as an objective, quantitative and easy-to-compute annotation for gene-metabolite pairs that can be used to integrate statistical evidence to biological networks.

Metabolic correlates to critical speed in murine models of sickle cell disease

Introduction: Exercise intolerance is a common clinical manifestation in patients with sickle cell disease (SCD), though the mechanisms are incompletely understood. Methods: Here we leverage a murine mouse model of sickle cell disease, the Berkeley mouse, to characterize response to exercise via determination of critical speed (CS), a functional measurement of mouse running speed upon exerting to exhaustion. Results: Upon observing a wide distribution in critical speed phenotypes, we systematically determined metabolic aberrations in plasma and organs-including heart, kidney, liver, lung, and spleen-from mice ranked based on critical speed performances (top vs. bottom 25%). Results indicated clear signatures of systemic and organ-specific alterations in carboxylic acids, sphingosine 1-phosphate and acylcarnitine metabolism. Metabolites in these pathways showed significant correlations with critical speed across all matrices. Findings from murine models were thus further validated in 433 sickle cell disease patients (SS genotype). Metabolomics analyses of plasma from 281 subjects in this cohort (with HbA < 10% to decrease confounding effects of recent transfusion events) were used to identify metabolic correlates to sub-maximal exercise test performances, as measure by 6 min walking test in this clinical cohort. Results confirmed strong correlation between test performances and dysregulated levels of circulating carboxylic acids (especially succinate) and sphingosine 1-phosphate. Discussion: We identified novel circulating metabolic markers of exercise intolerance in mouse models of sickle cell disease and sickle cell patients.

In vivo evaluation of the effect of sickle cell hemoglobin S, C and therapeutic transfusion on erythrocyte metabolism and cardiorenal dysfunction

Despite a wealth of exploratory plasma metabolomics studies in sickle cell disease (SCD), no study to date has evaluate a large and well phenotyped cohort to compare the primary erythrocyte metabolome of hemoglobin SS, SC and transfused AA red blood cells (RBCs) in vivo . The current study evaluates the RBC metabolome of 587 subjects with sickle cell sickle cell disease (SCD) from the WALK-PHaSST clinical cohort. The set includes hemoglobin SS, hemoglobin SC SCD patients, with variable levels of HbA related to RBC transfusion events, and HbF related to hydroxyurea therapy. Here we explore the modulating effects of genotype, age, sex, severity of hemolysis, and hydroxyurea and transfusion therapy on sickle RBC metabolism. Data - collated in an online portal - show that the Hb SS genotype is associated with significant alterations of RBC acylcarnitines, pyruvate, sphingosine 1-phosphate, creatinine, kynurenine and urate metabolism. Surprisingly, the RBC metabolism of SC RBCs is dramatically different from SS, with all glycolytic intermediates significantly elevated in SS RBCs, with the exception of pyruvate. This result suggests a metabolic blockade at the ATP-generating phosphoenolpyruvate to pyruvate step of glycolysis, which is catalyzed by redox-sensitive pyruvate kinase. Increasing in vivo concentrations of HbA improved glycolytic flux and normalized the HbS erythrocyte metabolome. An unexpectedly limited metabolic effect of hydroxyurea and HbF was observed, possibly related to the modest induction of HbF in this cohort. The metabolic signature of HbS RBCs correlated with the degree of steady state hemolytic anemia, cardiovascular and renal dysfunction and mortality.

Key points

In vivo dysregulation of RBC metabolism by HbS is evaluated by metabolic profiling of 587 patients with variable HbA, HbC and HbF levels;RBC acyl-carnitines, urate, pyruvate metabolism, S1P, kynurenine relate to hemolysis and cardiorenal dysfunction, respond to transfusion.

The effects of glutamine supplementation on markers of apoptosis and autophagy in sickle cell disease peripheral blood mononuclear cells

OBJECTIVES

L-Glutamine was FDA-approved for sickle cell disease (SCD) in 2017, yet the mechanism(s)-of-action are poorly understood. This study investigates the potential activation of autophagy as a previously unexplored mechanism-of-benefit.

DESIGN

Prospective, open-label, 8-week, phase-2 trial of oral L-glutamine (10 g TID) in patients with SCD at risk for pulmonary hypertension identified by Doppler-echocardiography by an elevated tricuspid-regurgitant-jet-velocity (TRV)≥ 2.5 m/s. Peripheral blood mononuclear cells (PBMCs) were isolated from blood samples taken from SCD patients at baseline, two, four, six and eight weeks of glutamine therapy, and from controls at baseline; BAX (pro-apoptotic marker) and LC3-II/LC3-I (autophagy marker) were measured via western blot analysis to assess apoptosis and autophagy respectively.

SETTING

Comprehensive SCD Center in Oakland, California.

RESULTS

Patients with SCD (n = 8) had a mean age of 44 ± 16, 50% were male; 63% Hb-SS, and mean TRV= 3.1 ± 0.7 m/s. Controls' mean age (n = 5) was 32 ± 12% and 57% were male; all were Hb-AA with a mean TRV= 1.8 ± 0.6. At baseline, SCD-PBMCs had 2-times higher levels of BAX and LC3-I versus controls (both p = 0.03). Levels of BAX expression increased by 300% after 8-weeks of glutamine supplementation (p = 0.005); LC3-I protein levels decreased while LC3-II levels increased by 70%, giving a significant increase in the LC3-II/LC3-I ratio (p = 0.02).

CONCLUSION

PBMCs from glutamine-supplemented SCD patients have upregulated apoptotic and autophagy proteins. The parallel increase in BAX and the LC3-II / LC3-I ratio with glutamine supplementation suggest a possible role of autophagic cell death. The increase in apoptotic markers provide insight into a possible mechanism used by peripheral PBMCs during glutamine supplementation in patients with SCD.

Flow adhesion of whole blood to P-selectin: a prognostic biomarker for vaso-occlusive crisis in sickle cell disease

Blood cell adhesion to P-selectin and vascular cell adhesion molecule-1 (VCAM-1) contributes to the pathophysiology of vaso-occlusion crisis (VOC) events in individuals with sickle cell disease (SCD). We evaluated the use of standardized flow adhesion biomarkers in a six-month, 35-subjects longitudinal study (ELIPSIS). Flow adhesion of whole blood on P-selectin (FA-WB-Psel) and VCAM1 (FA-WB-VCAM), and of isolated white blood cells on P-selectin (FA-WBC-Psel) and VCAM-1 (FA-WBC-VCAM) were elevated on VOC days compared with non-VOC days, but only FA-WB-Psel reached statistical significance (P = 0·015). Optimal cut-off values were established with Cox regression models for FA-WB-Psel [46 cells/mm²; hazard ratio (HR): 2·3; 95% confidence interval (CI):1·4-4·0; P = 0·01] and FA-WB-VCAM (408 cells/mm², HR:1·8; 95% CI: 0·9-3·45; P = 0·01). A combined (FA-WB-Psel and FA-WB-VCAM) multimarker risk score was also significantly (P = 0·0006) correlated with VOC risk that was two-fold higher for intermediate and 5·64-fold higher for high score. The concordance (C)-index for the multimarker score was 0·63 in the six-month period (95% CI: 0·56-0·70), indicating a better ability to distinguish patient risk of VOC, compared to individual biomarkers FA-WB-VCAM (C-index: 0·57; 95% CI: 0·49-0·65) or FA-WB-Psel (C-index: 0·58; 95% CI: 0·53-0·62). The presented multimarker score can be used to risk-stratify individuals with SCD during their steady state into low, intermediate, and high-risk strata for self-reported VOCs. Such risk stratification could help focus healthcare resources more efficiently to maintiain health, personalize treatment selection to each patient's individual needs, and potentially reduce healthcare costs.