Deconstructing SARS-CoV-2 neutralization: A modular molecular framework for computational design and comparison of antibodies and nanobodies targeting the spike RBD.

Since 2020 the COVID-19 pandemic has led scientists to search for strategies to predict the transmissibility and virulence of new severe acute respiratory syndrome coronavirus 2 variants based on the estimation of the affinity of the spike receptor binding domain (RBD) for the human angiotensin-converting enzyme 2 (ACE2) receptor and/or neutralizing antibodies. In this context, our lab developed a computational pipeline to quickly quantify the free energy of interaction at the spike RBD/ACE2 protein-protein interface, reflecting the incidence trend observed in the transmissibility/virulence of the investigated variants. In this new study, we used our pipeline to estimate the free energy of interaction between the RBD from 10 variants, and 14 antibodies (ab), or 5 nanobodies (nb), highlighting the RBD regions preferentially targeted by the investigated ab/nb. Our structural comparative analysis and interaction energy calculations allowed us to propose the most promising RBD regions to be targeted by future ab/nb to be designed by site-directed mutagenesis of existing high-affinity ab/nb, to increase their affinity for the target RBD region, for preventing spike-RBD/ACE2 interactions and virus entry in host cells. Furthermore, we evaluated the ability of the investigated ab/nb to simultaneously interact with the three RBD located on the surface of the trimeric spike protein, which can alternatively be in up- or down- (all-3-up-, all-3-down-, 1-up-/2-down-, 2-up-/1-down-) conformations.

N-acetylaspartate release by glutaminolytic ovarian cancer cells sustains protumoral macrophages.

Glutaminolysis is known to correlate with ovarian cancer aggressiveness and invasion. However, how this affects the tumor microenvironment is elusive. Here, we show that ovarian cancer cells become addicted to extracellular glutamine when silenced for glutamine synthetase (GS), similar to naturally occurring GS-low, glutaminolysis-high ovarian cancer cells. Glutamine addiction elicits a crosstalk mechanism whereby cancer cells release N-acetylaspartate (NAA) which, through the inhibition of the NMDA receptor, and synergistically with IL-10, enforces GS expression in macrophages. In turn, GS-high macrophages acquire M2-like, tumorigenic features. Supporting this in␣vitro model, in silico data and the analysis of ascitic fluid isolated from ovarian cancer patients prove that an M2-like macrophage phenotype, IL-10 release, and NAA levels positively correlate with disease stage. Our study uncovers the unprecedented role of glutamine metabolism in modulating macrophage polarization in highly invasive ovarian cancer and highlights the anti-inflammatory, protumoral function of NAA.

CRAT missense variants cause abnormal carnitine acetyltransferase function in an early-onset case of Leigh syndrome.

Leigh syndrome, or subacute necrotizing encephalomyelopathy, is one of the most severe pediatric disorders of the mitochondrial energy metabolism. By performing whole-exome sequencing in a girl affected by Leigh syndrome and her parents, we identified two heterozygous missense variants (p.Tyr110Cys and p.Val569Met) in the carnitine acetyltransferase (CRAT) gene, encoding an enzyme involved in the control of mitochondrial short-chain acyl-CoA concentrations. Biochemical assays revealed carnitine acetyltransferase deficiency in the proband-derived fibroblasts. Functional analyses of recombinant-purified CRAT proteins demonstrated that both missense variants, located in the acyl-group binding site of the enzyme, severely impair its catalytic function toward acetyl-CoA, and the p.Val569Met variant also toward propionyl-CoA and octanoyl-CoA. Although a single recessive variant in CRAT has been recently associated with neurodegeneration with brain iron accumulation (NBIA), this study reports the first kinetic analysis of naturally occurring CRAT variants and demonstrates the genetic basis of carnitine acetyltransferase deficiency in a case of mitochondrial encephalopathy.

SLC25A10 biallelic mutations in intractable epileptic encephalopathy with complex I deficiency.

Mitochondrial diseases are a plethora of inherited neuromuscular disorders sharing defects in mitochondrial respiration, but largely different from one another for genetic basis and pathogenic mechanism. Whole exome sequencing was performed in a familiar trio (trio-WES) with a child affected by severe epileptic encephalopathy associated with respiratory complex I deficiency and mitochondrial DNA depletion in skeletal muscle. By trio-WES we identified biallelic mutations in SLC25A10, a nuclear gene encoding a member of the mitochondrial carrier family. Genetic and functional analyses conducted on patient fibroblasts showed that SLC25A10 mutations are associated with reduction in RNA quantity and aberrant RNA splicing, and to absence of SLC25A10 protein and its transporting function. The yeast SLC25A10 ortholog knockout strain showed defects in mitochondrial respiration and mitochondrial DNA content, similarly to what observed in the patient skeletal muscle, and growth susceptibility to oxidative stress. Albeit patient fibroblasts were depleted in the main antioxidant molecules NADPH and glutathione, transport assays demonstrated that SLC25A10 is unable to transport glutathione. Here, we report the first recessive mutations of SLC25A10 associated to an inherited severe mitochondrial neurodegenerative disorder. We propose that SLC25A10 loss-of-function causes pathological disarrangements in respiratory-demanding conditions and oxidative stress vulnerability.

Mitochondrial Membranes of Human SH-SY5Y Neuroblastoma Cells Express Serotonin 5-HT7 Receptor.

Mitochondria in neurons contribute to energy supply, the regulation of synaptic transmission, Ca2+ homeostasis, neuronal excitability, and stress adaptation. In recent years, several studies have highlighted that the neurotransmitter serotonin (5-HT) plays an important role in mitochondrial biogenesis in cortical neurons, and regulates mitochondrial activity and cellular function in cardiomyocytes. 5-HT exerts its diverse actions by binding to cell surface receptors that are classified into seven distinct families (5-HT1 to 5-HT7). Recently, it was shown that 5-HT3 and 5-HT4 receptors are located on the mitochondrial membrane and participate in the regulation of mitochondrial function. Furthermore, it was observed that activation of brain 5-HT7 receptors rescued mitochondrial dysfunction in female mice from two models of Rett syndrome, a rare neurodevelopmental disorder characterized by severe behavioral and physiological symptoms. Our Western blot analyses performed on cell-lysate and purified mitochondria isolated from neuronal cell line SH-SY5Y showed that 5-HT7 receptors are also expressed into mitochondria. Maximal binding capacity (Bmax) obtained by Scatchard analysis on purified mitochondrial membranes was 0.081 pmol/mg of 5-HT7 receptor protein. Lastly, we evaluated the effect of selective 5-HT7 receptor agonist LP-211 and antagonist (inverse agonist) SB-269970 on mitochondrial respiratory chain (MRC) cytochrome c oxidase activity on mitochondria from SH-SY5Y cells. Our findings provide the first evidence that 5-HT7 receptor is also expressed in mitochondria.

Modeling SARS-CoV-2 spike/ACE2 protein-protein interactions for predicting the binding affinity of new spike variants for ACE2, and novel ACE2 structurally related human protein targets, for COVID-19 handling in the 3PM context.

Aims: The rapid spread of new SARS-CoV-2 variants has highlighted the crucial role played in the infection by mutations occurring at the SARS-CoV-2 spike receptor binding domain (RBD) in the interactions with the human ACE2 receptor. In this context, it urgently needs to develop new rapid tools for quickly predicting the affinity of ACE2 for the SARS-CoV-2 spike RBD protein variants to be used with the ongoing SARS-CoV-2 genomic sequencing activities in the clinics, aiming to gain clues about the transmissibility and virulence of new variants, to prevent new outbreaks and to quickly estimate the severity of the disease in the context of the 3PM. Methods: In our study, we used a computational pipeline for calculating the interaction energies at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface for a selected group of characterized infectious variants of concern/interest (VoC/VoI). By using our pipeline, we built 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for the VoC B.1.1.7-United Kingdom (carrying the mutations of concern/interest N501Y, S494P, E484K at the RBD), P.1-Japan/Brazil (RBD mutations: K417T, E484K, N501Y), B.1.351-South Africa (RBD mutations: K417N, E484K, N501Y), B.1.427/B.1.429-California (RBD mutations: L452R), the B.1.141 (RBD mutations: N439K), and the recent B.1.617.1-India (RBD mutations: L452R; E484Q) and the B.1.620 (RBD mutations: S477N; E484K). Then, we used the obtained 3D comparative models of the SARS-CoV-2 spike RBD/ACE2 protein complexes for predicting the interaction energies at the protein-protein interface. Results: Along SARS-CoV-2 mutation database screening and mutation localization analysis, it was ascertained that the most dangerous mutations at VoC/VoI spike proteins are located mainly at three regions of the SARS-CoV-2 spike "boat-shaped" receptor binding motif, on the RBD domain. Notably, the P.1 Japan/Brazil variant present three mutations, K417T, E484K, N501Y, located along the entire receptor binding motif, which apparently determines the highest interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, among those calculated. Conversely, it was also observed that the replacement of a single acidic/hydrophilic residue with a basic residue (E484K or N439K) at the "stern" or "bow" regions, of the boat-shaped receptor binding motif on the RBD, appears to determine an interaction energy with ACE2 receptor higher than that observed with single mutations occurring at the "hull" region or with other multiple mutants. In addition, our pipeline allowed searching for ACE2 structurally related proteins, i.e., THOP1 and NLN, which deserve to be investigated for their possible involvement in interactions with the SARS-CoV-2 spike protein, in those tissues showing a low expression of ACE2, or as a novel receptor for future spike variants. A freely available web-tool for the in silico calculation of the interaction energy at the SARS-CoV-2 spike RBD/ACE2 protein-protein interface, starting from the sequences of the investigated spike and/or ACE2 variants, was made available for the scientific community at: https://www.mitoairm.it/covid19affinities. Conclusion: In the context of the PPPM/3PM, the employment of the described pipeline through the provided webservice, together with the ongoing SARS-CoV-2 genomic sequencing, would help to predict the transmissibility of new variants sequenced from future patients, depending on SARS-CoV-2 genomic sequencing activities and on the specific amino acid replacement and/or on its location on the SARS-CoV-2 spike RBD, to put in play all the possible counteractions for preventing the most deleterious scenarios of new outbreaks, taking into consideration that a greater transmissibility has not to be necessarily related to a more severe manifestation of the disease. Supplementary Information: The online version contains supplementary material available at 10.1007/s13167-021-00267-w.

Genome-Wide Identification and Validation of Gene Expression Biomarkers in the Diagnosis of Ovarian Serous Cystadenocarcinoma.

Ovarian cancer is the second most prevalent gynecologic malignancy, and ovarian serous cystadenocarcinoma (OSCA) is the most common and lethal subtype of ovarian cancer. Current screening methods have strong limits on early detection, and the majority of OSCA patients relapse. In this work, we developed and cross-validated a method for detecting gene expression biomarkers able to discriminate OSCA tissues from healthy ovarian tissues and other cancer types with high accuracy. A preliminary ranking-based approach was applied, resulting in a panel of 41 over-expressed genes in OSCA. The RNA quantity gene expression of the 41 selected genes was then cross-validated by using NanoString nCounter technology. Moreover, we showed that the RNA quantity of eight genes (ADGRG1, EPCAM, ESRP1, MAL2, MYH14, PRSS8, ST14 and WFDC2) discriminates each OSCA sample from each healthy sample in our data set with sensitivity of 100% and specificity of 100%. For the other three genes (MUC16, PAX8 and SOX17) in combination, their RNA quantity may distinguish OSCA from other 29 tumor types.

Differential Expression of ADP/ATP Carriers as a Biomarker of Metabolic Remodeling and Survival in Kidney Cancers.

ADP/ATP carriers (AACs) are mitochondrial transport proteins playing a strategic role in maintaining the respiratory chain activity, fueling the cell with ATP, and also regulating mitochondrial apoptosis. To understand if AACs might represent a new molecular target for cancer treatment, we evaluated AAC expression levels in cancer/normal tissue pairs available on the Tissue Cancer Genome Atlas database (TCGA), observing that AACs are dysregulated in most of the available samples. It was observed that at least two AACs showed a significant differential expression in all the available kidney cancer/normal tissue pairs. Thus, we investigated AAC expression in the corresponding kidney non-cancer (HK2)/cancer (RCC-Shaw and CaKi-1) cell lines, grown in complete medium or serum starvation, for investigating how metabolic alteration induced by different growth conditions might influence AAC expression and resistance to mitochondrial apoptosis initiators, such as "staurosporine" or the AAC highly selective inhibitor "carboxyatractyloside". Our analyses showed that AAC2 and AAC3 transcripts are more expressed than AAC1 in all the investigated kidney cell lines grown in complete medium, whereas serum starvation causes an increase of at least two AAC transcripts in kidney cancer cell lines compared to non-cancer cells. However, the total AAC protein content is decreased in the investigated cancer cell lines, above all in the serum-free medium. The observed decrease in AAC protein content might be responsible for the decrease of OXPHOS activity and for the observed lowered sensitivity to mitochondrial apoptosis induced by staurosporine or carboxyatractyloside. Notably, the cumulative probability of the survival of kidney cancer patients seriously decreases with the decrease of AAC1 expression in KIRC and KIRP tissues making AAC1 a possible new biomarker of metabolic remodeling and survival in kidney cancers.

6-Thioguanine and Its Analogs Promote Apoptosis of Castration-Resistant Prostate Cancer Cells in a BRCA2-Dependent Manner.

Background: Mutations in the oncosuppressor gene BReast CAncer susceptibility gene 2 (BRCA2) predispose to aggressive forms of prostate cancer which show poor response to taxane-based therapy, the standard treatment for castration-resistant, aggressive prostate cancer. Herein, we addressed the question whether changes in BRCA2 expression, a potential surrogate marker for BRCA2 activity, may affect the response of castration-resistant prostate cancer cells to 6-thioguanine (6-TG), a thiopurine used in the treatment of haematological malignancies. Methods: Yeast, normal prostate cells and castration-resistant prostate cancer cells were treated with 6-TG or its analogues, in presence or absence of paclitaxel, or with olaparib, a poly-(ADP-ribose) polymerase (PARP) inhibitor currently in clinical trials for treatment of metastatic castration-resistant prostate cancer, and cell proliferation, apoptosis and androgen receptor (AR) levels were measured. Results: 6-TG inhibited cell proliferation in yeast, normal and castration-resistant prostate cancer cells but promoted apoptosis only in cancer cells. Suppression of BRCA2 expression by siRNA or shRNA increased the sensitivity to 6-TG- and olaparib-induced apoptosis but did not affect cancer cell response to taxane. Intriguingly, 6-TG reduced AR expression levels independently on BRCA2 expression. Instead, olaparib decreased AR levels only in BRCA2-knockdown prostate cancer cells. Notably, overexpression of BRCA2 resulted in resistance of castration-resistant prostate cancer cells to 6-TG-, taxane- and olaparib-based treatment but promoted sensitivity to apoptosis induced by 2-amino-6-bromopurine and 2,6-dithiopurine, two 6-TG analogues. Conclusions: Our results provide a pre-clinical rationale for the use of 6-TG in the treatment of BRCA2-deficient castration-resistant prostate cancers, and of certain 6-TG analogues for treatment of BRCA2-proficient prostate cancers.

FAD/NADH Dependent Oxidoreductases: From Different Amino Acid Sequences to Similar Protein Shapes for Playing an Ancient Function.

Flavoprotein oxidoreductases are members of a large protein family of specialized dehydrogenases, which include type II NADH dehydrogenase, pyridine nucleotide-disulphide oxidoreductases, ferredoxin-NAD+ reductases, NADH oxidases, and NADH peroxidases, playing a crucial role in the metabolism of several prokaryotes and eukaryotes. Although several studies have been performed on single members or protein subgroups of flavoprotein oxidoreductases, a comprehensive analysis on structure-function relationships among the different members and subgroups of this great dehydrogenase family is still missing. Here, we present a structural comparative analysis showing that the investigated flavoprotein oxidoreductases have a highly similar overall structure, although the investigated dehydrogenases are quite different in functional annotations and global amino acid composition. The different functional annotation is ascribed to their participation in species-specific metabolic pathways based on the same biochemical reaction, i.e., the oxidation of specific cofactors, like NADH and FADH2. Notably, the performed comparative analysis sheds light on conserved sequence features that reflect very similar oxidation mechanisms, conserved among flavoprotein oxidoreductases belonging to phylogenetically distant species, as the bacterial type II NADH dehydrogenases and the mammalian apoptosis-inducing factor protein, until now retained as unique protein entities in Bacteria/Fungi or Animals, respectively. Furthermore, the presented computational analyses will allow consideration of FAD/NADH oxidoreductases as a possible target of new small molecules to be used as modulators of mitochondrial respiration for patients affected by rare diseases or cancer showing mitochondrial dysfunction, or antibiotics for treating bacterial/fungal/protista infections.

The transcription factors ADR1 or CAT8 are required for RTG pathway activation and evasion from yeast acetic acid-induced programmed cell death in raffinose.

Yeast Saccharomyces cerevisiae grown on glucose undergoes programmed cell death (PCD) induced by acetic acid (AA-PCD), but evades PCD when grown in raffinose. This is due to concomitant relief of carbon catabolite repression (CCR) and activation of mitochondrial retrograde signaling, a mitochondria-to-nucleus communication pathway causing up-regulation of various nuclear target genes, such as CIT2, encoding peroxisomal citrate synthase, dependent on the positive regulator RTG2 in response to mitochondrial dysfunction. CCR down-regulates genes mainly involved in mitochondrial respiratory metabolism. In this work, we investigated the relationships between the RTG and CCR pathways in the modulation of AA-PCD sensitivity under glucose repression or de-repression conditions. Yeast single and double mutants lacking RTG2 and/or certain factors regulating carbon source utilization, including MIG1, HXK2, ADR1, CAT8, and HAP4, have been analyzed for their survival and CIT2 expression after acetic acid treatment. ADR1 and CAT8 were identified as positive regulators of RTG-dependent gene transcription. ADR1 and CAT8 interact with RTG2 and with each other in inducing cell resistance to AA-PCD in raffinose and controlling the nature of cell death. In the absence of ADR1 and CAT8, AA-PCD evasion is acquired through activation of an alternative factor/pathway repressed by RTG2, suggesting that RTG2 may play a function in promoting necrotic cell death in repressing conditions when RTG pathway is inactive. Moreover, our data show that simultaneous mitochondrial retrograde pathway activation and SNF1-dependent relief of CCR have a key role in central carbon metabolism reprogramming which modulates the yeast acetic acid-stress response.