Coupling between Protein Stability and Catalytic Activity Determines Pathogenicity of G6PD Variants

Analytical and operational considerations of measuring glucose 6-phosphate dehydrogenase (G6PD) activity using a fully automated assay

OBJECTIVES

This study determined the performance characteristics of a quantitative glucose-6-phosphate dehydrogenase (G6PD) assay with automated lysis and evaluated the robustness of the operational workflow following implementation in a hospital laboratory.

METHODS

The G6PD activity was measured in whole blood using an enzymatic quantitative test on a Roche cobas c501 analyzer with onboard lysis configuration and normalized to hemoglobin (Hb). The performance characteristics of the method and stability of G6PD in whole blood collected in EDTA-containing tubes were evaluated, and the reference interval was established on a population of healthy individuals (n = 279). The robustness of this automated workflow for sample lysis was evaluated during validation and after implementation for routine clinical use for 18 months and in 2,181 patients.

RESULTS

The G6PD assay was linear from 0.7 to 16.5 U/g Hb. Inter- and intra-assay precision using control and patient samples was below 12%. The G6PD results correlated well with a reference laboratory method (r = 0.96, y = 0.9615x - 1.222). The reference interval in our population was 9.8 to 15.5 U/g Hb. There were no interferences by lipemia and icteria, although grossly hemolyzed specimens may be affected. The testing workflow requires analyzing samples within minutes from mixing and loading into the instrument to avoid sample sedimentation. Measures to repeat samples with Hb 8.0 g/dL or less identified sedimented samples. In our patient population, 10.6% and 5.8% of the total males and females tested were G6PD deficient, respectively.

CONCLUSIONS

The G6PD assay with automated lysis is acceptable for patient testing. Several measures ensured the robustness of this workflow in a hospital laboratory.

Screening and the analysis of genotypic and phenotypic characterization of glucose-6-phosphate dehydrogenase (G6PD) deficiency in Fujian province, China

Introduction

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common X-linked hereditary disorder in southern China. However, the incidence rate of G6PD deficiency and the frequency of the most common G6PD gene variants vary widely. The purpose of this study was to investigate the prevalence, genotype, and phenotypic features of G6PD deficiency in neonates in Fujian province, southeastern China.

Methods

This retrospective cohort study enrolled 2,789,002 newborns (1,521,431 males and 1,267,571 females) based on the newborn screening program for G6PD deficiency in Fujian Province between January 2010 and December 2021.

Results

Of the 2,789,002 newborns enrolled, 26,437 cases were diagnosed (22,939 males and 3,498 females), and the estimated prevalence of G6PD deficiency in Fujian province was 0.95%. The prevalence was significantly higher among males (1.51%) than in females (0.28%) (p < 0.00001). Among the 3,198 patients with G6PD deficiency, 3,092 cases (2,145 males and 947 females) were detected to have G6PD gene variants. The top six prevalent genotypes identified represented 90.84% (2095/3,198) of the total and included c.1376G > T (44.93%), c.1388G > A (18.42%), c.1024C > T (9.32%), c.95A > G (8.69%), c.392G > T (5.25%), and c.871G > A (4.22%). The frequency of genotypes with c.1388G > A, c.1024C > T, and c.871G > A was higher in males in the Fujian province than in females, while the frequency of genotypes with c.1376G > T was lower. Furthermore, when comparing the enzyme activities of the top six prevalent genotypes, there were significant differences in the enzyme activities among the genotypes of male hemizygotes and female heterozygotes. According to the new classification of G6PD variants proposed by the World Health Organization (WHO), the variants with c.1376G > T, c.95A > G, and c.871G > A were recognized as Class A, while the c.392G > T, c.1388G > A, and c.1024C > T were recognized as Class B.

Discussion

To the best of our knowledge, this study is the first to systematically describe the overview of epidemiological characteristics of newborn G6PD deficiency in Fujian province, China, including the screening rate, incidence rate, and variant spectrum. Additionally, we elucidated the relationship between the distribution of enzyme activity with specific mutations and their WHO classification patterns. Our results could provide strategies for screening, diagnosis, and genetic counseling of G6PD deficiency in this area.

Acetylation-dependent coupling between G6PD activity and apoptotic signaling

Lysine acetylation has been discovered in thousands of non-histone human proteins, including most metabolic enzymes. Deciphering the functions of acetylation is key to understanding how metabolic cues mediate metabolic enzyme regulation and cellular signaling. Glucose-6-phosphate dehydrogenase (G6PD), the rate-limiting enzyme in the pentose phosphate pathway, is acetylated on multiple lysine residues. Using site-specifically acetylated G6PD, we show that acetylation can activate (AcK89) and inhibit (AcK403) G6PD. Acetylation-dependent inactivation is explained by structural studies showing distortion of the dimeric structure and active site of G6PD. We provide evidence for acetylation-dependent K95/97 ubiquitylation of G6PD and Y503 phosphorylation, as well as interaction with p53 and induction of early apoptotic events. Notably, we found that the acetylation of a single lysine residue coordinates diverse acetylation-dependent processes. Our data provide an example of the complex roles of acetylation as a posttranslational modification that orchestrates the regulation of enzymatic activity, posttranslational modifications, and apoptotic signaling.

The pentose phosphate pathway in health and disease

The pentose phosphate pathway (PPP) is a glucose-oxidizing pathway that runs in parallel to upper glycolysis to produce ribose 5-phosphate and nicotinamide adenine dinucleotide phosphate (NADPH). Ribose 5-phosphate is used for nucleotide synthesis, while NADPH is involved in redox homoeostasis as well as in promoting biosynthetic processes, such as the synthesis of tetrahydrofolate, deoxyribonucleotides, proline, fatty acids and cholesterol. Through NADPH, the PPP plays a critical role in suppressing oxidative stress, including in certain cancers, in which PPP inhibition may be therapeutically useful. Conversely, PPP-derived NADPH also supports purposeful cellular generation of reactive oxygen species (ROS) and reactive nitrogen species (RNS) for signalling and pathogen killing. Genetic deficiencies in the PPP occur relatively commonly in the committed pathway enzyme glucose-6-phosphate dehydrogenase (G6PD). G6PD deficiency typically manifests as haemolytic anaemia due to red cell oxidative damage but, in severe cases, also results in infections due to lack of leucocyte oxidative burst, highlighting the dual redox roles of the pathway in free radical production and detoxification. This Review discusses the PPP in mammals, covering its roles in biochemistry, physiology and disease.

A mitotic NADPH upsurge promotes chromosome segregation and tumour progression in aneuploid cancer cells

Redox metabolites have been observed to fluctuate through the cell cycle in cancer cells, but the functional impacts of such metabolic oscillations remain unknown. Here, we uncover a mitosis-specific nicotinamide adenine dinucleotide phosphate (NADPH) upsurge that is essential for tumour progression. Specifically, NADPH is produced by glucose 6-phosphate dehydrogenase (G6PD) upon mitotic entry, which neutralizes elevated reactive oxygen species (ROS) and prevents ROS-mediated inactivation of mitotic kinases and chromosome missegregation. Mitotic activation of G6PD depends on the phosphorylation of its co-chaperone protein BAG3 at threonine 285, which results in dissociation of inhibitory BAG3. Blocking BAG3T285 phosphorylation induces tumour suppression. A mitotic NADPH upsurge is present in aneuploid cancer cells with high levels of ROS, while nearly unobservable in near-diploid cancer cells. High BAG3T285 phosphorylation is associated with worse prognosis in a cohort of patients with microsatellite-stable colorectal cancer. Our study reveals that aneuploid cancer cells with high levels of ROS depend on a G6PD-mediated NADPH upsurge in mitosis to protect them from ROS-induced chromosome missegregation.

Functional characterization of MLH1 missense variants unveils mechanisms of pathogenicity and clarifies role in cancer

Lynch syndrome is a heritable condition caused by a heterozygous germline inactivating mutation of the DNA mismatch repair (MMR) genes, most commonly the MLH1 gene. However, one third of the identified alterations are missense variants, for which the clinical significance is unclear in many cases. We have identified three MLH1 missense alterations (p.(Glu736Lys), p.(Pro640Thr) and p.(Leu73Pro)) in six individuals from large Tunisian families. For none of these alterations, a classification of pathogenicity was available, consequently diagnosis, predictive testing and targeted surveillance in affected families was impossible. We therefore performed functional laboratory testing using a system testing stability as well as catalytic activity that includes clinically validated reference variants. Both p.(Leu73Pro) and p.(Pro640Thr) were found to be non-functional due to severe defects in protein stability and catalytic activity. In contrast, p.(Glu736Lys) was comparable to the wildtype protein and therefore considered a neutral substitution. Analysis of residue conservation and of the structural roles of the substituted residues corroborated these findings. In conjunction with the available clinical data, two variants fulfil classification criteria for class 4 "likely pathogenic". The findings of this work clarify the mechanism of pathogenicity of two unclear MLH1 variants and enables predictive testing and targeted surveillance in members of carrier families worldwide.

Targeting Apollo-NADP+ to Image NADPH Generation in Pancreatic Beta-Cell Organelles

NADPH/NADP⁺ redox state supports numerous reactions related to cell growth and survival; yet the full impact is difficult to appreciate due to organelle compartmentalization of NADPH and NADP⁺. To study glucose-stimulated NADPH production in pancreatic beta-cell organelles, we targeted the Apollo-NADP⁺ sensor by first selecting the most pH-stable version of the single-color sensor. We subsequently targeted mTurquoise2-Apollo-NADP⁺ to various organelles and confirmed activity in the cytoplasm, mitochondrial matrix, nucleus, and peroxisome. Finally, we measured the glucose- and glutamine-stimulated NADPH responses by single- and dual-color imaging of the targeted sensors. Overall, we developed multiple organelle-targeted Apollo-NADP⁺ sensors to reveal the prominent role of beta-cell mitochondria in determining NADPH production in the cytoplasm, nucleus, and peroxisome.

Stabilization of glucose-6-phosphate dehydrogenase oligomers enhances catalytic activity and stability of clinical variants

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a genetic trait that can cause hemolytic anemia. To date, over 150 nonsynonymous mutations have been identified in G6PD, with pathogenic mutations clustering near the dimer and/or tetramer interface and the allosteric NADP+-binding site. Recently, our lab identified a small molecule that activates G6PD variants by stabilizing the allosteric NADP+ and dimer complex, suggesting therapeutics that target these regions may improve structural defects. Here, we elucidated the connection between allosteric NADP+ binding, oligomerization, and pathogenicity to determine whether oligomer stabilization can be used as a therapeutic strategy for G6PD deficiency (G6PDdef). We first solved the crystal structure for G6PDK403Q, a mutant that mimics the physiological acetylation of wild-type G6PD in erythrocytes and demonstrated that loss of allosteric NADP+ binding induces conformational changes in the dimer. These structural changes prevent tetramerization, are unique to Class I variants (the most severe form of G6PDdef), and cause the deactivation and destabilization of G6PD. We also introduced nonnative cysteines at the oligomer interfaces and found that the tetramer complex is more catalytically active and stable than the dimer. Furthermore, stabilizing the dimer and tetramer improved protein stability in clinical variants, regardless of clinical classification, with tetramerization also improving the activity of G6PDK403Q and Class I variants. These findings were validated using enzyme activity and thermostability assays, analytical size-exclusion chromatography (SEC), and SEC coupled with small-angle X-ray scattering (SEC-SAXS). Taken together, our findings suggest a potential therapeutic strategy for G6PDdef and provide a foundation for future drug discovery efforts.

Combined effects of double mutations on catalytic activity and structural stability contribute to clinical manifestations of glucose-6-phosphate dehydrogenase deficiency

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy in humans, affecting ~ 500 million worldwide. A detailed study of the structural stability and catalytic activity of G6PD variants is required to understand how different mutations cause varying degrees of enzyme deficiency, reflecting the response of G6PD variants to oxidative stress. Furthermore, for G6PD double variants, investigating how two mutations jointly cause severe enzyme deficiency is important. Here, we characterized the functional and structural properties of nine G6PD variants: G6PD Gaohe, G6PD Mahidol, G6PD Shoklo, G6PD Canton, G6PD Kaiping, G6PD Gaohe + Kaiping, G6PD Mahidol + Canton, G6PD Mahidol + Kaiping and G6PD Canton + Kaiping. All variants were less catalytically active and structurally stable than the wild type enzyme, with G6PD double mutations having a greater impact than single mutations. G6PD Shoklo and G6PD Canton + Kaiping were the least catalytically active single and double variants, respectively. The combined effects of two mutations were observed, with the Canton mutation reducing structural stability and the Kaiping mutation increasing it in the double mutations. Severe enzyme deficiency in the double mutants was mainly determined by the trade-off between protein stability and catalytic activity. Additionally, it was demonstrated that AG1, a G6PD activator, only marginally increased G6PD enzymatic activity and stability.

Treatment strategies for glucose-6-phosphate dehydrogenase deficiency: past and future perspectives

Glucose-6-phosphate dehydrogenase (G6PD) maintains redox balance in a variety of cell types and is essential for erythrocyte resistance to oxidative stress. G6PD deficiency, caused by mutations in the G6PD gene, is present in ~400 million people worldwide, and can cause acute hemolytic anemia. Currently, there are no therapeutics for G6PD deficiency. We discuss the role of G6PD in hemolytic and nonhemolytic disorders, treatment strategies attempted over the years, and potential reasons for their failure. We also discuss potential pharmacological pathways, including glutathione (GSH) metabolism, compensatory NADPH production routes, transcriptional upregulation of the G6PD gene, highlighting potential drug targets. The needs and opportunities described here may motivate the development of a therapeutic for hematological and other chronic diseases associated with G6PD deficiency.

Long-range structural defects by pathogenic mutations in most severe glucose-6-phosphate dehydrogenase deficiency

Mechanism of the loss of activity of the most severe patient-derived mutants of glucose-6-phosphate dehydrogenase (G6PD) deficiency has remained elusive despite the availability of the G6PD structures for decades. Structural and biophysical investigations have revealed a common mechanism and dynamics of how these mutations hinder the substrate-binding site, reducing enzymatic activity. These are triggered by a long-distance propagation of structural defects at the dimer interface and the binding site of the noncatalytic cofactor. These structural distortions are found among all of the class I mutants investigated, providing critical clues for drug design to address G6PD deficiency by correcting the structural defects.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common blood disorder, presenting multiple symptoms, including hemolytic anemia. It affects 400 million people worldwide, with more than 160 single mutations reported in G6PD. The most severe mutations (about 70) are classified as class I, leading to more than 90% loss of activity of the wild-type G6PD. The crystal structure of G6PD reveals these mutations are located away from the active site, concentrating around the noncatalytic NADP⁺-binding site and the dimer interface. However, the molecular mechanisms of class I mutant dysfunction have remained elusive, hindering the development of efficient therapies. To resolve this, we performed integral structural characterization of five G6PD mutants, including four class I mutants, associated with the noncatalytic NADP⁺ and dimerization, using crystallography, small-angle X-ray scattering (SAXS), cryogenic electron microscopy (cryo-EM), and biophysical analyses. Comparisons with the structure and properties of the wild-type enzyme, together with molecular dynamics simulations, bring forward a universal mechanism for this severe G6PD deficiency due to the class I mutations. We highlight the role of the noncatalytic NADP⁺-binding site that is crucial for stabilization and ordering two β-strands in the dimer interface, which together communicate these distant structural aberrations to the active site through a network of additional interactions. This understanding elucidates potential paths for drug development targeting G6PD deficiency.

Prevalence and Molecular Characterization of Glucose-6-Phosphate Dehydrogenase (G6PD) Deficiency in Females from Previously Malaria Endemic Regions in Northeastern Thailand and Identification of a Novel G6PD Variant

Introduction

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common X-linked enzymopathy, highly prevalent in the areas where malaria is or has been endemic. Prevalence of G6PD deficiency and characterization of G6PD variants in females from previously malaria-endemic areas of northeastern Thailand remain unstudied.

Methods

Prevalence of G6PD deficiency was determined by a fluorescent spot test (FST), quantitative G6PD activity assay, and multiplex allele-specific (AS)- and restriction fragment length polymorphic (RFLP)-PCR developed for detection of common G6PD variants in the Thai population.

Results

Prevalence of G6PD deficiency in female samples (n = 355) was 18% by FST, 29.6% by quantitation of G6PD activity, and 28.1% by PCR-based genotyping. The most common variant was G6PD Viangchan (54%), followed by G6PD Canton (11%) and G6PD Union (11%); in addition, a novel heterozygous variant, G6PD Khon Kaen (c.305T>C, p.F102S), was identified. The majority of heterozygotes expressed G6PD activity within the intermediate deficiency range (30–70% median of normal enzyme activity).

Conclusion

High prevalence of G6PD deficiency was present in females from northeastern Thailand, the majority being due to heterozygosity of G6PD variants. The findings will have a bearing on the inclusion of primaquine in antimalarial-based policies for malaria elimination in populations with a high prevalence of G6PD deficiency.

Glucose-6-phosphate dehydrogenase deficiency

Glucose 6-phosphate dehydrogenase (G6PD) deficiency is 1 of the commonest human enzymopathies, caused by inherited mutations of the X-linked gene G6PD. G6PD deficiency makes red cells highly vulnerable to oxidative damage, and therefore susceptible to hemolysis. Over 200 G6PD mutations are known: approximately one-half are polymorphic and therefore common in various populations. Some 500 million persons with any of these mutations are mostly asymptomatic throughout their lifetime; however, any of them may develop acute and sometimes very severe hemolytic anemia when triggered by ingestion of fava beans, by any of a number of drugs (for example, primaquine, rasburicase), or, more rarely, by infection. Approximately one-half of the G6PD mutations are instead sporadic: rare patients with these mutations present with chronic nonspherocytic hemolytic anemia. Almost all G6PD mutations are missense mutations, causing amino acid replacements that entail deficiency of G6PD enzyme activity: they compromise the stability of the protein, the catalytic activity is decreased, or a combination of both mechanisms occurs. Thus, genotype-phenotype correlations have been reasonably well clarified in many cases. G6PD deficiency correlates remarkably, in its geographic distribution, with past/present malaria endemicity: indeed, it is a unique example of an X-linked human polymorphism balanced through protection of heterozygotes from malaria mortality. Acute hemolytic anemia can be managed effectively provided it is promptly diagnosed. Reliable diagnostic procedures are available, with point-of-care tests becoming increasingly important where primaquine and its recently introduced analog tafenoquine are required for the elimination of malaria.

Genetic spectrum and clinical early natural history of glucose-6-phosphate dehydrogenase deficiency in Mexican children detected through newborn screening

BACKGROUND

Glucose-6-phosphate dehydrogenase deficiency (G6PDd) newborn screening is still a matter of debate due to its highly heterogeneous birth prevalence and clinical expression, as well as, the lack of enough knowledge on its natural history. Herein, we describe the early natural clinical course and the underlying GDPD genotypes in infants with G6PDd detected by newborn screening and later studied in a single follow-up center. G6PDd newborns were categorized into three groups: group 1: hospitalized with or without neonatal jaundice (NNJ); group 2: non-hospitalized with NNJ; and group 3: asymptomatic. Frequencies of homozygous UGT1A1*28 (rs34983651) genotypes among G6PDd patients with or without NNJ were also explored.

RESULTS

A total of 81 newborns (80 males, one female) were included. Most individuals (46.9%) had NNJ without other symptoms, followed by asymptomatic (42.0%) and hospitalized (11.1%) patients, although the hospitalization of only 3 of these patients was related to G6PDd, including NNJ or acute hemolytic anemia (AHA). Nine different G6PDd genotypes were found; the G6PD A-202A/376G genotype was the most frequent (60.5%), followed by the G6PD A-376G/968C (22.2%) and the Union-Maewo (rs398123546, 7.4%) genotypes. These genotypes produce a wide range of clinical and biochemical phenotypes with significant overlapping residual enzymatic activity values among class I, II or III variants. Some G6PD A-202A/376G individuals had enzymatic values that were close to the cutoff value (5.3 U/g Hb, 4.6 and 4.8 U/g Hb in the groups with and without NNJ, respectively), while others showed extremely low enzymatic values (1.1 U/g Hb and 1.4 U/g Hb in the groups with and without NNJ, respectively). Homozygosity for UGT1A1*28 among G6PDd patients with (11.9%, N = 5/42) or without (10.3%, N = 4/39) NNJ did not shown significant statistical difference (p = 0.611).

CONCLUSION

Wide variability in residual enzymatic activity was noted in G6PDd individuals with the same G6PD genotype. This feature, along with a documented heterogeneous mutational spectrum, makes it difficult to categorize G6PD variants according to current WHO classification and precludes the prediction of complications such as AHA, which can occur even with > 10% of residual enzymatic activity and/or be associated with the common and mild G6PD A-376G/968C and G6PD A-202A/376G haplotypes.

Integrative identification of the pathogenic role of a novel G6PD missense mutation c.697G>C

Background

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a hereditary disease caused by pathogenic mutations of G6PD. While most of the pathogenic variants of G6PD have been annotated, hemolysis of unknown etiology but analogous to that in G6PD deficiency persists, implying the existence of undocumented pathogenic variants. In our previous study, we reported four novel G6PD variants in China, for which the pathogenicity remains to be verified.

Methods

The variants were verified by exogenous expression in HEK-293 cells, and their functions were predicted by PolyPhen-2 and SIFT. The CRISPR/Cas9 system was exploited to edit the G6PD c.697G>C variant in HEK-293 cells and K562 cells. The expression of G6PD was detected by quantitative PCR (qPCR) and western blotting. The cell growth capacity was detected by the CCK-8 assay and crystal violet staining. The G6PD enzyme activity was reflected by the G6P/6PG ratio test. The apoptosis of cells was detected by Annexin V-APC/7-AAD staining. The secondary and crystallographic structures were denoted according to the literature and PyMOL software. The G6PD protein was purified from lysis of transformed Escherichia coli (E. coli) cell with Ni-charged Resin Column. The enzymatic activity was detected at different temperatures.

Results

The G6PD activity of exogenous G6PD c.697G>C in HEK-293 cells was significantly lower than that of wild type (WT) G6PD, a finding that was consistent with the observation in clinical samples. The functional predictions conducted by different algorithms indicated the damage role of the G6PD c.697G>C variant in its enzymatic activity. We recapitulated the G6PD c.697G>C variant both in HEK-293 cells and K562 cells by adapting the CRISPR/Cas9 strategy. Using distinct cell lines expressing the G6PD c.697G>C variant endogenously, we confirmed the deteriorative role of the G6PD c.697G>C variant in its enzymatic activity. Regarding the secondary and crystallographic structure, we found a mutated amino acid approaching the structural NADP⁺ binding site. Finally, we demonstrated the c.697G>C variant compromised the thermal stability of G6PD protein.

Conclusions

Our data delineated the pathogenic role of G6PD c.697G>C variant for G6PD deficiency, implying the wide usage of CRISPR/Cas9 for genetic disease research.

Glucose-6-Phosphate Dehydrogenase Is Not Essential for K-Ras-Driven Tumor Growth or Metastasis

The enzyme glucose-6-phosphate dehydrogenase (G6PD) is a major contributor to NADPH production and redox homeostasis and its expression is upregulated and correlated with negative patient outcomes in multiple human cancer types. Despite these associations, whether G6PD is essential for tumor initiation, growth, or metastasis remains unclear. Here, we employ modern genetic tools to evaluate the role of G6PD in lung, breast, and colon cancer driven by oncogenic K-Ras. Human HCT116 colorectal cancer cells lacking G6PD exhibited metabolic indicators of oxidative stress, but developed into subcutaneous xenografts with growth comparable with that of wild-type controls. In a genetically engineered mouse model of non-small cell lung cancer driven by K-Ras G12D and p53 deficiency, G6PD knockout did not block formation or proliferation of primary lung tumors. In MDA-MB-231-derived human triple-negative breast cancer cells implanted as orthotopic xenografts, loss of G6PD modestly decreased primary site growth without ablating spontaneous metastasis to the lung and moderately impaired the ability of breast cancer cells to colonize the lung when delivered via tail vein injection. Thus, in the studied K-Ras tumor models, G6PD was not strictly essential for tumorigenesis and at most modestly promoted disease progression. SIGNIFICANCE: K-Ras-driven tumors can grow and metastasize even in the absence of the oxidative pentose pathway, a main NADPH production route.

A small molecule G6PD inhibitor reveals immune dependence on pentose phosphate pathway

Glucose is catabolized by two fundamental pathways, glycolysis to make ATP and the oxidative pentose phosphate pathway to make NADPH. The first step of the oxidative pentose phosphate pathway is catalyzed by the enzyme glucose-6-phosphate dehydrogenase (G6PD). Here we develop metabolite reporter and deuterium tracer assays to monitor cellular G6PD activity. Using these, we show that the most widely cited G6PD antagonist, dehydroepiandosterone (DHEA), does not robustly inhibit G6PD in cells. We then identify a small molecule (G6PDi-1) that more effectively inhibits G6PD. Across a range of cultured cells, G6PDi-1 depletes NADPH most strongly in lymphocytes. In T cells but not macrophages, G6PDi-1 markedly decreases inflammatory cytokine production. In neutrophils, it suppresses respiratory burst. Thus, we provide a cell-active small molecule tool for oxidative pentose phosphate pathway inhibition, and use it to identify G6PD as a pharmacological target for modulating immune response.

Quantification of glucose-6-phosphate dehydrogenase activity by spectrophotometry: A systematic review and meta-analysis

BACKGROUND

The radical cure of Plasmodium vivax and P. ovale requires treatment with primaquine or tafenoquine to clear dormant liver stages. Either drug can induce haemolysis in individuals with glucose-6-phosphate dehydrogenase (G6PD) deficiency, necessitating screening. The reference diagnostic method for G6PD activity is ultraviolet (UV) spectrophotometry; however, a universal G6PD activity threshold above which these drugs can be safely administered is not yet defined. Our study aimed to quantify assay-based variation in G6PD spectrophotometry and to explore the diagnostic implications of applying a universal threshold.

METHODS AND FINDINGS

Individual-level data were pooled from studies that used G6PD spectrophotometry. Studies were identified via PubMed search (25 April 2018) and unpublished contributions from contacted authors (PROSPERO: CRD42019121414). Studies were excluded if they assessed only individuals with known haematological conditions, were family studies, or had insufficient details. Studies of malaria patients were included but analysed separately. Included studies were assessed for risk of bias using an adapted form of the Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool. Repeatability and intra- and interlaboratory variability in G6PD activity measurements were compared between studies and pooled across the dataset. A universal threshold for G6PD deficiency was derived, and its diagnostic performance was compared to site-specific thresholds. Study participants (n = 15,811) were aged between 0 and 86 years, and 44.4% (7,083) were women. Median (range) activity of G6PD normal (G6PDn) control samples was 10.0 U/g Hb (6.3-14.0) for the Trinity assay and 8.3 U/g Hb (6.8-15.6) for the Randox assay. G6PD activity distributions varied significantly between studies. For the 13 studies that used the Trinity assay, the adjusted male median (AMM; a standardised metric of 100% G6PD activity) varied from 5.7 to 12.6 U/g Hb (p < 0.001). Assay precision varied between laboratories, as assessed by variance in control measurements (from 0.1 to 1.5 U/g Hb; p < 0.001) and study-wise mean coefficient of variation (CV) of replicate measures (from 1.6% to 14.9%; p < 0.001). A universal threshold of 100% G6PD activity was defined as 9.4 U/g Hb, yielding diagnostic thresholds of 6.6 U/g Hb (70% activity) and 2.8 U/g Hb (30% activity). These thresholds diagnosed individuals with less than 30% G6PD activity with study-wise sensitivity from 89% (95% CI: 81%-94%) to 100% (95% CI: 96%-100%) and specificity from 96% (95% CI: 89%-99%) to 100% (100%-100%). However, when considering intermediate deficiency (<70% G6PD activity), sensitivity fell to a minimum of 64% (95% CI: 52%-75%) and specificity to 35% (95% CI: 24%-46%). Our ability to identify underlying factors associated with study-level heterogeneity was limited by the lack of availability of covariate data and diverse study contexts and methodologies.

CONCLUSIONS

Our findings indicate that there is substantial variation in G6PD measurements by spectrophotometry between sites. This is likely due to variability in laboratory methods, with possible contribution of unmeasured population factors. While an assay-specific, universal quantitative threshold offers robust diagnosis at the 30% level, inter-study variability impedes performance of universal thresholds at the 70% level. Caution is advised in comparing findings based on absolute G6PD activity measurements across studies. Novel handheld quantitative G6PD diagnostics may allow greater standardisation in the future.

Loss of angiogenin function is related to earlier ALS onset and a paradoxical increase in ALS duration

0.5-1% of ALS (Amyotrophic Lateral Sclerosis) and Parkinson's disease (PD) are associated with mutations in the angiogenin (ANG). These mutations are thought to cause disease through a loss of ANG function, but this hypothesis has not been evaluated statistically. In addition, the potential for ANG to promote disease has not been considered. With the goal of better defining the etiology of ANG-ALS, we assembled all clinical onset and disease duration data and determined if these were correlated with biochemical properties of ANG variants. Loss of ANG stability and ribonuclease activity were found to correlate with early ALS onset, confirming an aspect of the prevailing model of ANG-ALS. Conversely, loss of ANG stability and ribonuclease activity correlated with longer survival following diagnosis, which is inconsistent with the prevailing model. These results indicate that functional ANG appears to decrease the risk of developing ALS but exacerbate ALS once in progress. These findings are rationalized in terms of studies demonstrating that distinct mechanisms contribute to ALS onset and progression and propose that ANG replacement or stabilization would benefit pre-symptomatic ANG-ALS patients. However, this study challenges the prevailing hypothesis that augmenting ANG will benefit symptomatic ANG-ALS patients. Instead, our results suggest that silencing of ANG activity may be beneficial for symptomatic ALS patients. This study will serve as a call-to-arms for neurologists to consistently publish ALS and PD patient's clinical data-if all ANG-ALS patients' data were available our findings could be tested with considerable statistical power.

Small-Molecule Activators of Glucose-6-phosphate Dehydrogenase (G6PD) Bridging the Dimer Interface

We recently identified AG1, a small-molecule activator that functions by promoting oligomerization of glucose-6-phosphate dehydrogenase (G6PD) to the catalytically competent forms. Biochemical experiments indicate that the activation of G6PD by the original hit molecule (AG1) is noncovalent and that one C2 -symmetric region of the G6PD homodimer is important for ligand function. Consequently, the disulfide in AG1 is not required for activation of G6PD, and a number of analogues were prepared without this reactive moiety. Our study supports a mechanism of action whereby AG1 bridges the dimer interface at the structural nicotinamide adenine dinucleotide phosphate (NADP+ ) binding sites of two interacting G6PD monomers. Small molecules that promote G6PD oligomerization have the potential to provide a first-in-class treatment for G6PD deficiency. This general strategy could be applied to other enzyme deficiencies in which control of oligomerization can enhance enzymatic activity and/or stability.

Revealing evolutionary constraints on proteins through sequence analysis

Statistical analysis of alignments of large numbers of protein sequences has revealed "sectors" of collectively coevolving amino acids in several protein families. Here, we show that selection acting on any functional property of a protein, represented by an additive trait, can give rise to such a sector. As an illustration of a selected trait, we consider the elastic energy of an important conformational change within an elastic network model, and we show that selection acting on this energy leads to correlations among residues. For this concrete example and more generally, we demonstrate that the main signature of functional sectors lies in the small-eigenvalue modes of the covariance matrix of the selected sequences. However, secondary signatures of these functional sectors also exist in the extensively-studied large-eigenvalue modes. Our simple, general model leads us to propose a principled method to identify functional sectors, along with the magnitudes of mutational effects, from sequence data. We further demonstrate the robustness of these functional sectors to various forms of selection, and the robustness of our approach to the identification of multiple selected traits.

Inhibition of Drp1/Fis1 interaction slows progression of amyotrophic lateral sclerosis

Bioenergetic failure and oxidative stress are common pathological hallmarks of amyotrophic lateral sclerosis (ALS), but whether these could be targeted effectively for novel therapeutic intervention needs to be determined. One of the reported contributors to ALS pathology is mitochondrial dysfunction associated with excessive mitochondrial fission and fragmentation, which is predominantly mediated by Drp1 hyperactivation. Here, we determined whether inhibition of excessive fission by inhibiting Drp1/Fis1 interaction affects disease progression. We observed mitochondrial excessive fragmentation and dysfunction in several familial forms of ALS patient‐derived fibroblasts as well as in cultured motor neurons expressing SOD1 mutant. In both cell models, inhibition of Drp1/Fis1 interaction by a selective peptide inhibitor, P110, led to a significant reduction in reactive oxygen species levels, and to improvement in mitochondrial structure and functions. Sustained treatment of mice expressing G93A SOD1 mutation with P110, beginning at the onset of disease symptoms at day 90, produced an improvement in motor performance and survival, suggesting that Drp1 hyperactivation may be an attractive target in the treatment of ALS patients.

Associations of Mitochondrial and Nuclear Mitochondrial Variants and Genes with Seven Metabolic Traits

Mitochondria (MT), the major site of cellular energy production, are under dual genetic control by 37 mitochondrial DNA (mtDNA) genes and numerous nuclear genes (MT-nDNA). In the CHARGEmtDNA+ Consortium, we studied genetic associations of mtDNA and MT-nDNA associations with body mass index (BMI), waist-hip-ratio (WHR), glucose, insulin, HOMA-B, HOMA-IR, and HbA1c. This 45-cohort collaboration comprised 70,775 (insulin) to 170,202 (BMI) pan-ancestry individuals. Validation and imputation of mtDNA variants was followed by single-variant and gene-based association testing. We report two significant common variants, one in MT-ATP6 associated (p ≤ 5E-04) with WHR and one in the D-loop with glucose. Five rare variants in MT-ATP6, MT-ND5, and MT-ND6 associated with BMI, WHR, or insulin. Gene-based meta-analysis identified MT-ND3 associated with BMI (p ≤ 1E-03). We considered 2,282 MT-nDNA candidate gene associations compiled from online summary results for our traits (20 unique studies with 31 dataset consortia's genome-wide associations [GWASs]). Of these, 109 genes associated (p ≤ 1E-06) with at least 1 of our 7 traits. We assessed regulatory features of variants in the 109 genes, cis- and trans-gene expression regulation, and performed enrichment and protein-protein interactions analyses. Of the identified mtDNA and MT-nDNA genes, 79 associated with adipose measures, 49 with glucose/insulin, 13 with risk for type 2 diabetes, and 18 with cardiovascular disease, indicating for pleiotropic effects with health implications. Additionally, 21 genes related to cholesterol, suggesting additional important roles for the genes identified. Our results suggest that mtDNA and MT-nDNA genes and variants reported make important contributions to glucose and insulin metabolism, adipocyte regulation, diabetes, and cardiovascular disease.

Inferring interaction partners from protein sequences using mutual information

Functional protein-protein interactions are crucial in most cellular processes. They enable multi-protein complexes to assemble and to remain stable, and they allow signal transduction in various pathways. Functional interactions between proteins result in coevolution between the interacting partners, and thus in correlations between their sequences. Pairwise maximum-entropy based models have enabled successful inference of pairs of amino-acid residues that are in contact in the three-dimensional structure of multi-protein complexes, starting from the correlations in the sequence data of known interaction partners. Recently, algorithms inspired by these methods have been developed to identify which proteins are functional interaction partners among the paralogous proteins of two families, starting from sequence data alone. Here, we demonstrate that a slightly higher performance for partner identification can be reached by an approximate maximization of the mutual information between the sequence alignments of the two protein families. Our mutual information-based method also provides signatures of the existence of interactions between protein families. These results stand in contrast with structure prediction of proteins and of multi-protein complexes from sequence data, where pairwise maximum-entropy based global statistical models substantially improve performance compared to mutual information. Our findings entail that the statistical dependences allowing interaction partner prediction from sequence data are not restricted to the residue pairs that are in direct contact at the interface between the partner proteins.

Correcting glucose-6-phosphate dehydrogenase deficiency with a small-molecule activator

Glucose-6-phosphate dehydrogenase (G6PD) deficiency, one of the most common human genetic enzymopathies, is caused by over 160 different point mutations and contributes to the severity of many acute and chronic diseases associated with oxidative stress, including hemolytic anemia and bilirubin-induced neurological damage particularly in newborns. As no medications are available to treat G6PD deficiency, here we seek to identify a small molecule that corrects it. Crystallographic study and mutagenesis analysis identify the structural and functional defect of one common mutant (Canton, R459L). Using high-throughput screening, we subsequently identify AG1, a small molecule that increases the activity of the wild-type, the Canton mutant and several other common G6PD mutants. AG1 reduces oxidative stress in cells and zebrafish. Furthermore, AG1 decreases chloroquine- or diamide-induced oxidative stress in human erythrocytes. Our study suggests that a pharmacological agent, of which AG1 may be a lead, will likely alleviate the challenges associated with G6PD deficiency.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency provides insufficient protection from oxidative stress, contributing to diverse human pathologies. Here, the authors identify a small molecule that increases the activity and/or stability of mutant G6PD and show that it reduces oxidative stress in zebrafish and hemolysis in isolated human erythrocytes.

Extracting Complementary Insights from Molecular Phenotypes for Prioritization of Disease-Associated Mutations

Rapid advances in next-generation sequencing technology have resulted in an explosion of whole-exome/genome sequencing data, providing an unprecedented opportunity to identify disease- and trait-associated variants in humans on a large scale. To date, the long-standing paradigm has leveraged fitness-based approximations to translate this ever-expanding sequencing data into causal insights in disease. However, while this approach robustly identifies variants under evolutionary constraint, it fails to provide molecular insights. Moreover, complex disease phenomena often violate standard assumptions of a direct organismal phenotype to overall fitness effect relationship. Here we discuss the potential of a molecular phenotype-oriented paradigm to uniquely identify candidate disease-causing mutations from the human genetic background. By providing a direct connection between single nucleotide mutations and observable organismal and cellular phenotypes associated with disease, we suggest that molecular phenotypes can readily incorporate alongside established fitness-based methodologies to provide complementary insights to the functional impact of human mutations. Lastly, we discuss how integrated approaches between molecular phenotypes and fitness-based perspectives facilitate new insights into the molecular mechanisms underlying disease-associated mutations while also providing a platform for improved interpretation of epistasis in human disease.

Validation of the quantitative point-of-care CareStart biosensor for assessment of G6PD activity in venous blood

INTRODUCTION

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common enzymopathy in the human population affecting an estimated 8% of the world population, especially those living in areas of past and present malaria endemicity. Decreased G6PD enzymatic activity is associated with drug-induced hemolysis and increased risk of severe neonatal hyperbilirubinemia leading to brain damage. The G6PD gene is on the X chromosome therefore mutations cause enzymatic deficiency in hemizygote males and homozygote females while the majority of heterozygous females have an intermediate activity (between 30-80% of normal) with a large distribution into the range of deficiency and normality. Current G6PD qualitative tests are unable to diagnose G6PD intermediate activities which could hinder wide use of 8-aminoquinolines for Plasmodium vivax elimination. The aim of the study was to assess the diagnostic performances of the new Carestart G6PD quantitative biosensor.

METHODS

A total of 150 samples of venous blood with G6PD deficient, intermediate and normal phenotypes were collected among healthy volunteers living along the north-western Thailand-Myanmar border. Samples were analyzed by complete blood count, by gold standard spectrophotometric assay using Trinity kits and by the latest model of Carestart G6PD biosensor which analyzes both G6PD and hemoglobin.

RESULTS

Bland-Altman comparison of the CareStart normalized G6PD values to that of the gold standard assay showed a strong bias in values resulting in poor area under-the-curve values for both 30% and 80% thresholds. Performing a receiver operator curve identified threshold values for the CareStart product equivalent to the 30% and 80% gold standard values with good sensitivity and specificity values, 100% and 92% (for 30% G6PD activity) and 92% and 94% (for 80% activity) respectively.

CONCLUSION

The Carestart G6PD biosensor represents a significant improvement for quantitative diagnosis of G6PD deficiency over previous versions. Further improvements and validation studies are required to assess its utility for informing radical cure decisions in malaria endemic settings.

Prevalence, genetic variants and clinical implications of G-6-PD deficiency in Burkina Faso: a systematic review

BACKGROUND

It is now well-known that some antimalarials such as primaquine may induce severe hemolytic anemia in people with G-6-PD deficiency. Antimalarial drug prescriptions must, therefore take into account the patient's G-6-PD status in malaria endemic areas such as Burkina Faso, where the prevalence of this genetic abnormality is relatively high. Although great clinical heterogeneity is observed depending on the molecular nature of the deficiency and the residual enzyme activity in the red blood cell, there is very poor data on the prevalence of G-6-PD deficiency and the distribution of involved genetic variants in Burkina Faso. In this systematic review, we present a synthesis of the various studies carried out on the G-6-PD deficiency in Burkina Faso in order to determine its prevalence, probable distribution of the genetic variants involved and their clinical implications for a national systematic screening policy among the groups most vulnerable to malaria.

METHODS

A systematic review was carried out to analyze available published data on the prevalence, phenotypes and mutations responsible for G-6-PD deficiency in Burkina Faso. The key words used were "G-6-PD deficiency AND Burkina Faso" or "Déficit en G-6-PD AND Burkina Faso" in French. To identify the relevant articles, two independent reviewers reviewed the titles, abstracts and the full text of the selected papers.

RESULTS

An average prevalence of 16.6% (183/1100; CI 95%: 0.145-0.190) and 6.5% (69/1066; CI 95%: 0.051-0.081) of G-6-PD deficiency was found respectively in men and women in this systematic review. Although the predominance (99.8% of G-6-PD deficient cases) of 202A/376G G-6-PD A- variant, the Santamaria and Betica Selma variants were identified in Burkina Faso. Independently of the method used, the enzymatic deficiency was significantly higher in males (2.5-20.5%) compared to females (3.3-12.3%).

CONCLUSION

This systematic review suggests that despite the ubiquity of the 202A/376G G-6-PD A- variant in Burkina Faso, it will be necessary to consider the Santamaria and Betica Selma variants although their frequencies remain to be specified. A systematic screening of the G-6-PD deficiency is also needed to prevent the occurrence of iatrogenic hemolytic accidents.

A trade off between catalytic activity and protein stability determines the clinical manifestations of glucose-6-phosphate dehydrogenase (G6PD) deficiency

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common polymorphism and enzymopathy in humans, affecting approximately 400 million people worldwide. It is responsible for various clinical manifestations, including favism, hemolytic anemia, chronic non-spherocytic hemolytic anemia, spontaneous abortion, and neonatal hyperbilirubinemia. Understanding the molecular mechanisms underlying the severity of G6PD deficiency is of great importance but that of many G6PD variants are still unknown. In this study, we report the construction, expression, purification, and biochemical characterization in terms of kinetic properties and stability of five clinical G6PD variants-G6PD Bangkok, G6PD Bangkok noi, G6PD Songklanagarind, G6PD Canton+Bangkok noi, and G6PD Union+Viangchan. G6PD Bangkok and G6PD Canton+Bangkok noi showed a complete loss of catalytic activity and moderate reduction in thermal stability when compared with the native G6PD. G6PD Bangkok noi and G6PD Union+Viangchan showed a significant reduction in catalytic efficiency, whereas G6PD Songklanagarind showed a catalytic activity comparable to the wild-type enzyme. The Union+Viangchan mutation showed a remarkable effect on the global stability of the enzyme. In addition, our results indicate that the location of mutations in G6PD variants affects their catalytic activity, stability, and structure. Hence, our results provide a molecular explanation for clinical manifestations observed in individuals with G6PD deficiency.

Structural analysis of clinically relevant pathogenic G6PD variants reveals the importance of tetramerization for G6PD activity

Over 220 different amino acid variants have been identified in human glucose-6-phosphate dehydrogenase (G6PD), covering over 30% of the protein sequence. Many of these variants are pathogenic, causing varying degrees of G6PD deficiency with symptoms ranging from severe chronic anemia (class I) to milder triggered hemolytic episodes (classes II and III). The phenotypic effects of most G6PD variants have been reported, providing an opportunity to correlate phenotypic and structural information. In particular, we sought to investigate the tetramer interface of G6PD in relation to pathogenic variation, as there are conflicting reports indicating the importance of tetramerization for G6PD activity. Using a three-dimensional spatial scan statistic, hotspots of structural enrichment were identified for each class of pathogenic G6PD variants. Class I variants, the most phenotypically severe, were enriched at the dimer interface, consistent with previous evidence that dimerization is essential for G6PD activity. Class II variants were enriched near the tetramer interface, suggesting that tetramerization is also important for G6PD activity. This analysis explains why these two classes, both yielding 10% or less G6PD activity as compared to normal, lead to different clinical outcomes.