Missense variants in human ACE2 strongly affect binding to SARS-CoV-2 Spike providing a mechanism for ACE2 mediated genetic risk in Covid-19: A case study in affinity predictions of interface variants

SARS-CoV-2 Spike (Spike) binds to human angiotensin-converting enzyme 2 (ACE2) and the strength of this interaction could influence parameters relating to virulence. To explore whether population variants in ACE2 influence Spike binding and hence infection, we selected 10 ACE2 variants based on affinity predictions and prevalence in gnomAD and measured their affinities and kinetics for Spike receptor binding domain through surface plasmon resonance (SPR) at 37°C. We discovered variants that reduce and enhance binding, including three ACE2 variants that strongly inhibited (p.Glu37Lys, ΔΔG = –1.33 ± 0.15 kcal mol^-1 and p.Gly352Val, predicted ΔΔG = –1.17 kcal mol^-1) or abolished (p.Asp355Asn) binding. We also identified two variants with distinct population distributions that enhanced affinity for Spike. ACE2 p.Ser19Pro (ΔΔG = 0.59 ± 0.08 kcal mol^-1) is predominant in the gnomAD African cohort (AF = 0.003) whilst p.Lys26Arg (ΔΔG = 0.26 ± 0.09 kcal mol^-1) is predominant in the Ashkenazi Jewish (AF = 0.01) and European non-Finnish (AF = 0.006) cohorts. We compared ACE2 variant affinities to published SARS-CoV-2 pseudotype infectivity data and confirmed that ACE2 variants with reduced affinity for Spike can protect cells from infection. The effect of variants with enhanced Spike affinity remains unclear, but we propose a mechanism whereby these alleles could cause greater viral spreading across tissues and cell types, as is consistent with emerging understanding regarding the interplay between receptor affinity and cell-surface abundance. Finally, we compared mCSM-PPI2 ΔΔG predictions against our SPR data to assess the utility of predictions in this system. We found that predictions of decreased binding were well-correlated with experiment and could be improved by calibration, but disappointingly, predictions of highly enhanced binding were unreliable. Recalibrated predictions for all possible ACE2 missense variants at the Spike interface were calculated and used to estimate the overall burden of ACE2 variants on Covid-19.     

Regioselectivity mechanism of the Thunbergia alata Δ6-16:0-acyl carrier protein desaturase

Plant plastidial acyl–acyl carrier protein (ACP) desaturases are a soluble class of diiron-containing enzymes that are distinct from the diiron-containing integral membrane desaturases found in plants and other organisms. The archetype of this class is the stearoyl-ACP desaturase which converts stearoyl-ACP into oleoyl (18:1Δ⁹cis)-ACP. Several variants expressing distinct regioselectivity have been described including a Δ⁶-16:0-ACP desaturase from black-eyed Susan vine (Thunbergia alata). We solved a crystal structure of the T. alata desaturase at 2.05 Å resolution. Using molecular dynamics (MD) simulations, we identified a low-energy complex between 16:0-ACP and the desaturase that would position C6 and C7 of the acyl chain adjacent to the diiron active site. The model complex was used to identify mutant variants that could convert the T. alata Δ⁶ desaturase to Δ⁹ regioselectivity. Additional modeling between ACP and the mutant variants confirmed the predicted regioselectivity. To validate the in-silico predictions, we synthesized two variants of the T. alata desaturase and analyzed their reaction products using gas chromatography-coupled mass spectrometry. Assay results confirmed that mutants designed to convert T. alata Δ⁶ to Δ⁹ selectivity exhibited the predicted changes. In complementary experiments, variants of the castor desaturase designed to convert Δ⁹ to Δ⁶ selectivity lost some of their Δ9 desaturation ability and gained the ability to desaturate at the Δ⁶ position. The computational workflow for revealing the mechanistic understanding of regioselectivity presented herein lays a foundation for designing acyl-ACP desaturases with novel selectivities to increase the diversity of monoenes available for bioproduct applications.

Predictions regarding the mechanism of Δ⁶ regioselectivity of the Thunbergia alata desaturase based on X-ray crystallography and molecular dynamics simulations are confirmed by experiment.     

Gene therapy of Gaucher's and Fabry's diseases: current status and prospects

Gaucher disease and Fabry disease are lysosomal storage disorders characterized by the accumulation of sphingolipids. In both cases, the goal of gene therapy is to permanently provide tissues with enzyme levels allowing to avoid storage of the undigested substrates. Different gene therapy strategies must however be designed as Gaucher disease is due to a deficiency in the membrane-associated enzyme glucocerebrosidase, whereas Fabry disease is caused by a deficiency in the soluble enzyme alpha-galactosidase. Indeed, a soluble enzyme can be provided to tissues is trans by gene-modified cells whereas gene transfer has to target the most affected cells in the case of membrane-bound enzymes. Thus, in non-neurological Gaucher disease (type 1), the hematopoietic tissue has to be targeted as the deficiency affects the monocyte/macrophage lineage. Following promising preclinical studies, clinical protocols have been initiated to explore the feasibility and safety of retroviral transfer of the glucocerebrosidase gene into CD34+ cells from patients with type 1 Gaucher disease. Although gene-marked cells were detected in vivo, the level of corrected cells was very low, a finding indicating that improved vectors along with partial myeloablation may be necessary. Here, lentiviral vectors should enable more gene transduction into the hematopoietic target cells. As concerns the diffuse neurological lesions in types 2 and 3 of Gaucher disease, they will probably be especially difficult to target by gene therapy because of the non soluble nature of glucocerebrosidase. Finally, over the last few years, Fabry disease has become a compelling target for gene therapy as an etiology-based treatment strategy. Indeed, several recent studies aiming at creating a large in vivo source of alpha-galactosidase have yielded positive long-term results in the Fabry knock-out mouse model.     

Plasma Globotriaosylsphingosine and α-Galactosidase A Activity as a Combined Screening Biomarker for Fabry Disease in a Large Japanese Cohort

Fabry disease is an X-linked disorder of α-galactosidase A (GLA) deficiency. Our previous interim analysis (1 July 2014 to 31 December 2015) revealed plasma globotriaosylsphingosine as a promising primary screening biomarker for Fabry disease probands. Herein, we report the final results, including patients enrolled from 1 January to 31 December 2016 for evaluating the potential of plasma globotriaosylsphingosine and GLA activity as a combined screening marker. We screened 5691 patients (3439 males) referred from 237 Japanese specialty clinics based on clinical findings suggestive of Fabry disease using plasma globotriaosylsphingosine and GLA activity as primary screening markers, and GLA variant status as a secondary screening marker. Of the 14 males who tested positive in the globotriaosylsphingosine screen (≥2.0 ng/mL), 11 with low GLA activity (<4.0 nmol/h/mL) displayed GLA variants (four classic, seven late-onset) and one with normal GLA activity and no pathogenic variant displayed lamellar bodies in affected organs, indicating late-onset biopsy-proven Fabry disease. Of the 19 females who tested positive in the globotriaosylsphingosine screen, eight with low GLA activity displayed GLA variants (six classic, two late-onset) and five with normal GLA activity displayed a GLA variant (one classic) and no pathogenic variant (four late-onset biopsy-proven). The combination of plasma globotriaosylsphingosine and GLA activity can be a primary screening biomarker for classic, late-onset, and late-onset biopsy-proven Fabry disease probands.     

Therapy of Fabry disease with pharmacological chaperones: from in silico predictions to in vitro tests

Background

Fabry disease is a rare disorder caused by a large variety of mutations in the gene encoding lysosomal alpha-galactosidase. Many of these mutations are unique to individual families. Fabry disease can be treated with enzyme replacement therapy, but a promising novel strategy relies on small molecules, so called "pharmacological chaperones", which can be administered orally. Unfortunately only 42% of genotypes respond to pharmacological chaperones.

Results

A procedure to predict which genotypes responsive to pharmacological chaperones in Fabry disease has been recently proposed. The method uses a position-specific substitution matrix to score the mutations. Using this method, we have screened public databases for predictable responsive cases and selected nine representative mutations as yet untested with pharmacological chaperones. Mutant lysosomal alpha galactosidases were produced by site directed mutagenesis and expressed in mammalian cells. Seven out of nine mutations responded to pharmacological chaperones. Nineteen other mutations that were tested with pharmacological chaperones, but were not included in the training of the predictive method, were gathered from literature and analyzed in silico. In this set all five mutations predicted to be positive were responsive to pharmacological chaperones, bringing the percentage of responsive mutations among those predicted to be positive and not used to train the classifier to 86% (12/14). This figure differs significantly from the percentage of responsive cases observed among all the Fabry mutants tested so far.

Conclusions

In this paper we provide experimental support to an "in silico" method designed to predict missense mutations in the gene encoding lysosomal alpha galactosidase responsive to pharmacological chaperones. We demonstrated that responsive mutations can be predicted with a low percentage of false positive cases. Most of the mutations tested to validate the method were described in the literature as associated to classic or mild classic phenotype. The analysis can provide a guideline for the therapy with pharmacological chaperones supported by experimental results obtained in vitro. We are aware that our results were obtained in vitro and cannot be translated straightforwardly into benefit for patients, but need to be validated by clinical trials.

     

Acute-Phase CD8 T Cell Responses That Select for Escape Variants Are Needed to Control Live Attenuated Simian Immunodeficiency Virus

The overall CD8 T cell response to human/simian immunodeficiency virus (HIV/SIV) targets a collection of discrete epitope specificities. Some of these epitope-specific CD8 T cells emerge in the weeks and months following infection and rapidly select for sequence variants, whereas other CD8 T cell responses develop during the chronic infection phase and rarely select for sequence variants. In this study, we tested the hypothesis that acute-phase CD8 T cell responses that do not rapidly select for escape variants are unable to control viral replication in vivo as well as those that do rapidly select for escape variants. We created a derivative of live attenuated SIV (SIVmac239Δnef) in which we ablated five epitopes that elicit early CD8 T cell responses and rapidly accumulate sequence variants in SIVmac239-infected Mauritian cynomolgus macaques (MCMs) that are homozygous for the M3 major histocompatibility complex (MHC) haplotype. This live attenuated SIV variant was called m3KOΔnef. Viremia was significantly higher in M3 homozygous MCMs infected with m3KOΔnef than in either MHC-mismatched MCMs infected with m3KOΔnef or MCMs infected with SIVmac239Δnef. Three CD8 T cell responses, including two that do not rapidly select for escape variants, predominated during early m3KOΔnef infection in the M3 homozygous MCMs, but these animals were unable to control viral replication. These results provide evidence that acute-phase CD8 T cell responses that have the potential to rapidly select for escape variants in the early phase of infection are needed to establish viral control in vivo.     

Enzyme replacement therapy: conception,chaos and culmination

Soon after the enzymatic defects in Gaucher disease and in Niemann-Pick disease were discovered, enzyme replacement or enzyme supplementation was proposed as specific treatment for patients with these and related metabolic storage disorders. While relatively straightforward in concept, successful implementation of this approach required many years of intensive effort to bring it to fruition. Procedures were eventually developed to produce sufficient quantities of the requisite enzymes for clinical trials and to target therapeutic enzymes to lipid-storing cells. These achievements led to the development of effective enzyme replacement therapy for patients with Gaucher disease and for Fabry disease. These demonstrations provide strong incentive for the application of this strategy for the treatment of many human disorders of metabolism.     

Microbiota-Independent Ameliorative Effects of Antibiotics on Spontaneous Th2-Associated Pathology of the Small Intestine

We have previously generated a mouse model of spontaneous Th2-associated disease of the small intestine called TRAF6ΔDC, in which dendritic cell (DC)-intrinsic expression of the signaling mediator TRAF6 is ablated. Interestingly, broad-spectrum antibiotic treatment ameliorates TRAF6ΔDC disease, implying a role for commensal microbiota in disease development. However, the relationship between the drug effects and commensal microbiota status remains to be formally demonstrated. To directly assess this relationship, we have now generated TRAF6ΔDC bone marrow chimera mice under germ-free (GF) conditions lacking commensal microbiota, and found, unexpectedly, that Th2-associated disease is actually exacerbated in GF TRAF6ΔDC mice compared to specific pathogen-free (SPF) TRAF6ΔDC mice. At the same time, broad-spectrum antibiotic treatment of GF TRAF6ΔDC mice has an ameliorative effect similar to that observed in antibiotics-treated SPF TRAF6ΔDC mice, implying a commensal microbiota-independent effect of broad-spectrum antibiotic treatment. We further found that treatment of GF TRAF6ΔDC mice with broad-spectrum antibiotics increases Foxp3⁺ Treg populations in lymphoid organs and the small intestine, pointing to a possible mechanism by which treatment may directly exert an immunomodulatory effect. To investigate links between the exacerbated phenotype of the small intestines of GF TRAF6ΔDC mice and local microbiota, we performed microbiotic profiling of the luminal contents specifically within the small intestines of diseased TRAF6ΔDC mice, and, when compared to co-housed control mice, found significantly increased total bacterial content characterized by specific increases in Firmicutes Lactobacillus species. These data suggest a protective effect of Firmicutes Lactobacillus against the spontaneous Th2-related inflammation of the small intestine of the TRAF6ΔDC model, and may represent a potential mechanism for related disease phenotypes.     

Evidence for new C-terminally truncated variants of α- and β-tubulins

Cellular α-tubulin can bear various carboxy-terminal sequences: full-length tubulin arising from gene neosynthesis is tyrosinated, and two truncated variants, corresponding to detyrosinated and Δ2 α‑tubulin, result from the sequential cleavage of one or two C-terminal residues, respectively. Here, by using a novel antibody named 3EG that is highly specific to the –EEEG C-terminal sequence, we demonstrate the occurrence in neuronal tissues of a new αΔ3‑tubulin variant corresponding to α1A/B‑tubulin deleted of its last three residues (EEY). αΔ3‑tubulin has a specific distribution pattern: its quantity in the brain is similar to that of αΔ2-tubulin around birth but is much lower in adult tissue. This truncated α1A/B-tubulin variant can be generated from αΔ2-tubulin by the deglutamylases CCP1, CCP4, CCP5, and CCP6 but not by CCP2 and CCP3. Moreover, using 3EG antibody, we identify a C‑terminally truncated β-tubulin form with the same –EEEG C-terminal sequence. Using mass spectrometry, we demonstrate that β2A/B-tubulin is modified by truncation of the four C-terminal residues (EDEA). We show that this newly identified βΔ4-tubulin is ubiquitously present in cells and tissues and that its level is constant throughout the cell cycle. These new C-terminally truncated α- and β-tubulin variants, both ending with –EEEG sequence, are expected to regulate microtubule physiology. Of interest, the αΔ3-tubulin seems to be related to dynamic microtubules, resembling tyrosinated-tubulin rather than the other truncated variants, and may have critical function(s) in neuronal development.     

Functionally Relevant Domains of the Prion Protein Identified In Vivo

The prion consists essentially of PrP^Sc, a misfolded and aggregated conformer of the cellular protein PrP^C. Whereas PrP^C deficient mice are clinically healthy, expression of PrP^C variants lacking its central domain (PrP_ΔCD), or of the PrP-related protein Dpl, induces lethal neurodegenerative syndromes which are repressed by full-length PrP. Here we tested the structural basis of these syndromes by grafting the amino terminus of PrP^C (residues 1–134), or its central domain (residues 90–134), onto Dpl. Further, we constructed a soluble variant of the neurotoxic PrP_ΔCD mutant that lacks its glycosyl phosphatidyl inositol (GPI) membrane anchor. Each of these modifications abrogated the pathogenicity of Dpl and PrP_ΔCD in transgenic mice. The PrP-Dpl chimeric molecules, but not anchorless PrP_ΔCD, ameliorated the disease of mice expressing truncated PrP variants. We conclude that the amino proximal domain of PrP exerts a neurotrophic effect even when grafted onto a distantly related protein, and that GPI-linked membrane anchoring is necessary for both beneficial and deleterious effects of PrP and its variants.     

Defining the Functional Domain of Programmed Cell Death 10 through Its Interactions with Phosphatidylinositol-3,4,5-Trisphosphate

Cerebral cavernous malformations (CCM) are vascular abnormalities of the central nervous system predisposing blood vessels to leakage, leading to hemorrhagic stroke. Three genes, Krit1 (CCM1), OSM (CCM2), and PDCD10 (CCM3) are involved in CCM development. PDCD10 binds specifically to PtdIns(3,4,5)P3 and OSM. Using threading analysis and multi-template modeling, we constructed a three-dimensional model of PDCD10. PDCD10 appears to be a six-helical-bundle protein formed by two heptad-repeat-hairpin structures (α1–3 and α4–6) sharing the closest 3D homology with the bacterial phosphate transporter, PhoU. We identified a stretch of five lysines forming an amphipathic helix, a potential PtdIns(3,4,5)P3 binding site, in the α5 helix. We generated a recombinant wild-type (WT) and three PDCD10 mutants that have two (Δ2KA), three (Δ3KA), and five (Δ5KA) K to A mutations. Δ2KA and Δ3KA mutants hypothetically lack binding residues to PtdIns(3,4,5)P3 at the beginning and the end of predicted helix, while Δ5KA completely lacks all predicted binding residues. The WT, Δ2KA, and Δ3KA mutants maintain their binding to PtdIns(3,4,5)P3. Only the Δ5KA abolishes binding to PtdIns(3,4,5)P3. Both Δ5KA and WT show similar secondary and tertiary structures; however, Δ5KA does not bind to OSM. When WT and Δ5KA are co-expressed with membrane-bound constitutively-active PI3 kinase (p110-CAAX), the majority of the WT is co-localized with p110-CAAX at the plasma membrane where PtdIns(3,4,5)P3 is presumably abundant. In contrast, the Δ5KA remains in the cytoplasm and is not present in the plasma membrane. Combining computational modeling and biological data, we propose that the CCM protein complex functions in the PI3K signaling pathway through the interaction between PDCD10 and PtdIns(3,4,5)P3.     

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient’s RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZ_EI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.     

Targeting autophagy as a novel strategy for facilitating the therapeutic action of potentiators on ΔF508 cystic fibrosis transmembrane conductance regulator

Channel activators (potentiators) of cystic fibrosis (CF) transmembrane conductance regulator (CFTR), can be used for the treatment of the small subset of CF patients that carry plasma membrane-resident CFTR mutants. However, approximately 90% of CF patients carry the misfolded ΔF508-CFTR and are poorly responsive to potentiators, because ΔF508-CFTR is intrinsically unstable at the plasma membrane (PM) even if rescued by pharmacological correctors. We have demonstrated that human and mouse CF airways are autophagy deficient due to functional sequestration of BECN1 and that the tissue transglutaminase-2 inhibitor, cystamine, or antioxidants restore BECN1-dependent autophagy and reduce SQSTM1/p62 levels, thus favoring ΔF508-CFTR trafficking to the epithelial surface. Here, we investigated whether these treatments could facilitate the beneficial action of potentiators on ΔF508-CFTR homozygous airways. Cystamine or the superoxide dismutase (SOD)/catalase-mimetic EUK-134 stabilized ΔF508-CFTR at the plasma membrane of airway epithelial cells and sustained the expression of CFTR at the epithelial surface well beyond drug withdrawal, overexpressing BECN1 and depleting SQSTM1. This facilitates the beneficial action of potentiators in controlling inflammation in ex vivo ΔF508-CFTR homozygous human nasal biopsies and in vivo in mouse ΔF508-CFTR lungs. Direct depletion of Sqstm1 by shRNAs in vivo in ΔF508-CFTR mice synergized with potentiators in sustaining surface CFTR expression and suppressing inflammation. Cystamine pre-treatment restored ΔF508-CFTR response to the CFTR potentiators genistein, Vrx-532 or Vrx-770 in freshly isolated brushed nasal epithelial cells from ΔF508-CFTR homozygous patients. These findings delineate a novel therapeutic strategy for the treatment of CF patients with the ΔF508-CFTR mutation in which patients are first treated with cystamine and subsequently pulsed with CFTR potentiators.     

Epithelial–Mesenchymal Transition in Kidney Tubular Epithelial Cells Induced by Globotriaosylsphingosine and Globotriaosylceramide

Fabry disease is a lysosomal storage disorder caused by deficiency of alpha-galactosidase A (α-gal A), which results in the deposition of globotriaosylceramide (Gb3) in the vascular endothelium. Globotriaosylsphingosine (lyso-Gb3), a deacylated Gb3, is also increased in the plasma of patients with Fabry disease. Renal fibrosis is a key feature of advanced Fabry disease patients. Therefore, we evaluated the association of Gb3 and lyso-Gb3 accumulation and the epithelial–mesenchymal transition (EMT) on tubular epithelial cells of the kidney. In HK2 cells, exogenous treatments of Gb3 and lyso-Gb3 increased the expression of TGF-β, EMT markers (N-cadherin and α-SMA), and phosphorylation of PI3K/AKT, and decreased the expression of E-cadherin. Lyso-Gb3, rather than Gb3, strongly induced EMT in HK2 cells. In the mouse renal mesangial cell line, SV40 MES 13 cells, Gb3 strongly induced phenotype changes. The EMT induced by Gb3 was inhibited by enzyme α-gal A treatment, but EMT induced by lyso-Gb3 was not abrogated by enzyme treatment. However, TGF-β receptor inhibitor (TRI, SB525334) inhibited the activation of TGF-β and EMT markers in HK2 cells with Gb3 and lyso-Gb3 treatments. This study suggested that increased plasma lyso-Gb3 has a crucial role in the development of renal fibrosis through the cell-specific induction of the EMT in Fabry disease, and that TRI treatment, alongside enzyme replacement therapy, could be a potential therapeutic option for patients with Fabry disease.     

5’-adenosine monophosphate mediated cooling treatment enhances ΔF508-Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) stability in vivo

Background

Gene mutations that produce misprocessed proteins are linked to many human disorders. Interestingly, some misprocessed proteins retained their biological function when stabilized by low temperature treatment of cultured cells in vitro. Here we investigate whether low temperature treatment in vivo can rescue misfolded proteins by applying 5’-AMP mediated whole body cooling to a Cystic Fibrosis (CF) mouse model carrying a mutant cystic fibrosis transmembrane conductance regulator (CFTR) with a deletion of the phenylalanine residue in position 508 (ΔF508-CFTR). Low temperature treatment of cultured cells was previously shown to be able to alleviate the processing defect of ΔF508-CFTR, enhancing its plasma membrane localization and its function in mediating chloride ion transport.

Results

Here, we report that whole body cooling enhanced the retention of ΔF508-CFTR in intestinal epithelial cells. Functional analysis based on β-adrenergic dependent salivary secretion and post-natal mortality rate revealed a moderate but significant improvement in treated compared with untreated CF mice.

Conclusions

Our findings demonstrate that temperature sensitive processing of mutant proteins can be responsive to low temperature treatment in vivo.     

Prediction of Protein Relative Enthalpic Stability from Molecular Dynamics Simulations of the Folded and Unfolded States

For proteins of known structure, the relative enthalpic stability with respect to wild-type, ΔΔH^U, can be estimated by direct computation of the folded and unfolded state energies. We propose a model by which the change in stability upon mutation can be predicted from all-atom molecular dynamics simulations for the folded state and a peptide-based model for the unfolded state. The unfolding enthalpies are expressed in terms of environmental and hydration-solvent reorganization contributions that readily allow a residue-specific analysis of ΔΔH^U. The method is applied to estimate the relative enthalpic stability of variants with buried charged groups in T4 lysozyme. The predicted relative stabilities are in good agreement with experimental data. Environmental factors are observed to contribute more than hydration to the overall ΔΔH^U. The residue-specific analysis finds that the effects of burying charge are both localized and long-range. The enthalpy for hydration-solvent reorganization varies considerably among different amino-acid types, but because the variant folded state structures are similar to those of the wild-type, the hydration-solvent reorganization contribution to ΔΔH^U is localized at the mutation site, in contrast to environmental contributions. Overall, mutation of apolar and polar amino acids to charged amino acids are destabilizing, but the reasons are complex and differ from site to site.     

Gaucher disease and Fabry disease: New markers and insights in pathophysiology for two distinct glycosphingolipidoses

Gaucher disease (GD) and Fabry disease (FD) are two relatively common inherited glycosphingolipidoses caused by deficiencies in the lysosomal glycosidases glucocerebrosidase and alpha-galactosidase A, respectively. For both diseases enzyme supplementation is presently used as therapy. Cells and tissues of GD and FD patients are uniformly deficient in enzyme activity, but the two diseases markedly differ in cell types showing lysosomal accumulation of the glycosphingolipid substrates glucosylceramide and globotriaosylceramide, respectively. The clinical manifestation of Gaucher disease and Fabry disease is consequently entirely different and the response to enzyme therapy is only impressive in the case of GD patients. This review compares both glycosphingolipid storage disorders with respect to similarities and differences. Presented is an update on insights regarding pathophysiological mechanisms as well as recently available biochemical markers and diagnostic tools for both disorders. Special attention is paid to sphingoid bases of the primary storage lipids in both diseases. The value of elevated glucosylsphingosine in Gaucher disease and globotriaosylsphingosine in Fabry disease for diagnosis and monitoring of disease is discussed as well as the possible contribution of the sphingoid bases to (patho)physiology. This article is part of a Special Issue entitled New Frontiers in Sphingolipid Biology.     

ΔFBA—Predicting metabolic flux alterations using genome-scale metabolic models and differential transcriptomic data

Genome-scale metabolic models (GEMs) provide a powerful framework for simulating the entire set of biochemical reactions in a cell using a constraint-based modeling strategy called flux balance analysis (FBA). FBA relies on an assumed metabolic objective for generating metabolic fluxes using GEMs. But, the most appropriate metabolic objective is not always obvious for a given condition and is likely context-specific, which often complicate the estimation of metabolic flux alterations between conditions. Here, we propose a new method, called ΔFBA (deltaFBA), that integrates differential gene expression data to evaluate directly metabolic flux differences between two conditions. Notably, ΔFBA does not require specifying the cellular objective. Rather, ΔFBA seeks to maximize the consistency and minimize inconsistency between the predicted flux differences and differential gene expression. We showcased the performance of ΔFBA through several case studies involving the prediction of metabolic alterations caused by genetic and environmental perturbations in Escherichia coli and caused by Type-2 diabetes in human muscle. Importantly, in comparison to existing methods, ΔFBA gives a more accurate prediction of flux differences.