Can machines learn the mutation signatures of SARS-CoV-2 and enable viral-genotype guided predictive prognosis?

MotivationContinuous emergence of new variants through appearance/accumulation/disappearance of mutations is a hallmark of many viral diseases. SARS-CoV-2 variants have particularly exerted tremendous pressure on global healthcare system owing to their life threatening and debilitating implications. The sheer plurality of variants and huge scale of genomic data have added to the challenges of tracing the mutations/variants and their relationship to infection severity (if any).ResultsWe explored the suitability of virus-genotype guided machine-learning in infection prognosis and identification of features/mutations-of-interest. Total 199,519 outcome-traced genomes, representing 45,625 nucleotide-mutations, were employed. Among these, post data-cleaning, Low and High severity genomes were classified using an integrated model (employing virus genotype, epitopic-influence and patient-age) with consistently high ROC-AUC (Asia:0.97 ± 0.01, Europe:0.94 ± 0.01, N.America:0.92 ± 0.02, Africa:0.94 ± 0.07, S.America:0.93 ± 03). Although virus-genotype alone could enable high predictivity (0.97 ± 0.01, 0.89 ± 0.02, 0.86 ± 0.04, 0.95 ± 0.06, 0.9 ± 0.04), the performance was not found to be consistent and the models for a few geographies displayed significant improvement in predictivity when the influence of age and/or epitope was incorporated with virus-genotype (Wilcoxon p_BH < 0.05). Neither age or epitopic-influence or clade information could out-perform the integrated features. A sparse model (6 features), developed using patient-age and epitopic-influence of the mutations, performed reasonably well (>0.87 ± 0.03, 0.91 ± 0.01, 0.87 ± 0.03, 0.84 ± 0.08, 0.89 ± 0.05). High-performance models were employed for inferring the important mutations-of-interest using Shapley Additive exPlanations (SHAP). The changes in HLA interactions of the mutated epitopes of reference SARS-CoV-2 were then subsequently probed. Notably, we also describe the significance of a ‘temporal-modeling approach’ to benchmark the models linked with continuously evolving pathogens. We conclude that while machine learning can play a vital role in identifying relevant mutations and factors driving the severity, caution should be exercised in using the genotypic signatures for predictive prognosis.     

Plant synthetic biology: from inspiration to augmentation

1. Download : Download high-res image (240KB)
2. Download : Download full-size image

     

Design and synthesis of boron-containing diphenylpyrimidines as potent BTK and JAK3 dual inhibitors

Bruton's tyrosine kinase (BTK) and Janus kinase 3 (JAK3) are very promising targets for hematological malignancies and autoimmune diseases. In recent years, a few compounds have been approved as a marketed medicine, and several are undergoing clinical trials. By recombining the dominant backbone of known active compounds, constructing a foused library, and screening a broad panel of kinases, we found a class of compounds with dual activities of anti-BTK and anti-JAK3. Some of the compounds have shown 10-folds more active in the enzyme and cell-based assays than a known active compound. Furthermore, liver microsome stability experiments show that these compounds have better stability than ibrutinib. These explorations offered new clues to discover benzoxaborole fragment and pyrimidine scaffold as more effective BTK and JAK3 dual inhibitors.     

Molecular evolution and functional divergence of UDP-hexose 4-epimerases

UDP-glucose 4-epimerase (GalE) catalyzes the interconversion of UDP-glucose (UDP-Glc) and UDP-galactose (UDP-Gal) and/or the interconversion of UDP-N-acetylglucosamine (UDP-GlcNAc) and UDP-N-acetylgalactosamine (UDP-GalNAc) in sugar metabolism. GalEs belong to the short-chain dehydrogenase/reductase superfamily, use a conserved ‘transient keto intermediate’ mechanism and have variable substrate specificity. GalEs have been classified into three groups based on substrate specificity: group 1 prefers UDP-Glc/Gal, group 3 prefers UDP-GlcNAc/GalNAc, and group 2 has comparable activities for both types of the substrates. The phylogenetic relationship and structural basis for the specificities of GalEs revealed possible molecular evolution of UDP-hexose 4-epimerases in various organisms. Based on the recent advances in studies on GalEs and related enzymes, an updated view of their evolutional diversification is presented.     

Selective Inhibition of Cysteine-Dependent Enzymes by Bioorthogonal Tethering

A general approach for the rapid and selective inhibition of enzymes in cells using a common tool compound would be of great value for research and therapeutic development. We previously reported a chemogenetic strategy that addresses this challenge for kinases, relying on bioorthogonal tethering of a pan inhibitor to a target kinase through a genetically encoded non-canonical amino acid. However, pan inhibitors are not available for many enzyme classes. Here, we expand the scope of the chemogenetic strategy to cysteine-dependent enzymes by bioorthogonal tethering of electrophilic warheads. For proof of concept, selective inhibition of two E2 ubiquitin-conjugating enzymes, UBE2L3 and UBE2D1, was demonstrated in biochemical assays. Further development and optimization of this approach should enable its use in cells as well as with other cysteine-dependent enzymes, facilitating the investigation of their cellular function and validation as therapeutic targets.     

Impact on S. aureus and E. coli Membranes of Treatment with Chlorhexidine and Alcohol Solutions: Insights from Molecular Simulations and Nuclear Magnetic Resonance

Membranes form the first line of defence of bacteria against potentially harmful molecules in the surrounding environment. Understanding the protective properties of these membranes represents an important step towards development of targeted anti-bacterial agents such as sanitizers. Use of propanol, isopropanol and chlorhexidine can significantly decrease the threat imposed by bacteria in the face of growing anti-bacterial resistance via mechanisms that include membrane disruption. Here we have employed molecular dynamics simulations and nuclear magnetic resonance to explore the impact of chlorhexidine and alcohol on the S. aureus cell membrane, as well as the E. coli inner and outer membranes. We identify how sanitizer components partition into these bacterial membranes, and show that chlorhexidine is instrumental in this process.     

Characteristic fragmentation of polyunsaturated fatty acids with allylic vicinal diols in positive-ion LC/ESI-MS/MS

A characteristic fragmentation was observed for PUFAs that contain allylic vicinal diol groups (resolvin D1, D2, D4, E3, lipoxin A4, B4, and maresin 2), which were derivatized with N,N-dimethylethylenediamine (DMED), in positive-ion ESI-MS/MS. The findings indicate that when these compounds contain an allylic hydroxyl group that is located distal to the terminal DMED moiety in the case of resolvin D1, D4, and lipoxin A4, an aldehyde (-CH=O) is predominately formed, which arises from the breakdown in between vicinal diols, whereas, in the case of an allylic hydroxyl group that is located proximal to the DMED moiety, as in resolvin D2, E3, lipoxin B4, and maresin 2, an allylic carbene (-CH=CH-CH:) is formed. These specific fragmentations could be used as diagnostic ions for characterizing the above seven PUFAs. As a result, it was possible to detect resolvin D1, D2, E3, lipoxin A4, and B4 in sera (20 μl) obtained from healthy volunteers by multiple-reaction monitoring using LC/ESI-MS/MS.     

Synthesis of novel 2-pyrrolidinone and pyrrolidine derivatives and study of their inhibitory activity against autotaxin enzyme

Autotaxin (ATX), a glycoprotein (~125 kDa) isolated as an autocrine motility factor from melanoma cells, belongs to a seven-membered family of ectonucleotide pyrophosphatase/phosphodiesterase (ENPP), and exhibits lysophospholipase D activity. ATX is responsible for the hydrolysis of lysophosphatidylcholine (LPC) to produce the bioactive lipid lysophosphatidic acid (LPA), which is upregulated in a variety of pathological inflammatory conditions, including fibrosis, cancer, liver toxicity and thrombosis. Given its role in human disease, the ATX-LPA axis is an interesting target for therapy, and the development of novel potent ATX inhibitors is of great importance. In the present work a novel class of ATX inhibitors, optically active derivatives of 2-pyrrolidinone and pyrrolidine heterocycles were synthesized. Some of them exhibited interesting in vitro activity, namely the hydroxamic acid 16 (IC₅₀ 700 nM) and the carboxylic acid 40b (IC₅₀ 800 nM), while the boronic acid derivatives 3k (IC₅₀ 50 nM), 3l (IC₅₀ 120 nM), 3 m (IC₅₀ 180 nM) and 21 (IC₅₀ 35 nM) were found to be potent inhibitors of ATX.     

l-carnitine: Nutrition,pathology, and health benefits

Carnitine is a medically needful nutrient that contributes in the production of energy and the metabolism of fatty acids. Bioavailability is higher in vegetarians than in people who eat meat. Deficits in carnitine transporters occur as a result of genetic mutations or in combination with other illnesses such like hepatic or renal disease. Carnitine deficit can arise in diseases such endocrine maladies, cardiomyopathy, diabetes, malnutrition, aging, sepsis, and cirrhosis due to abnormalities in carnitine regulation. The exogenously provided molecule is obviously useful in people with primary carnitine deficits, which can be life-threatening, and also some secondary deficiencies, including such organic acidurias: by eradicating hypotonia, muscle weakness, motor skills, and wasting are all improved l-carnitine (LC) have reported to improve myocardial functionality and metabolism in ischemic heart disease patients, as well as athletic performance in individuals with angina pectoris. Furthermore, although some intriguing data indicates that LC could be useful in a variety of conditions, including carnitine deficiency caused by long-term total parenteral supplementation or chronic hemodialysis, hyperlipidemias, and the prevention of anthracyclines and valproate-induced toxicity, such findings must be viewed with caution.     

Structural Basis for the Ribonuclease Activity of a Thermostable CRISPR-Cas13a from Thermoclostridium caenicola

The RNA-targeting type VI CRISPR-Cas effector complexes are widely used in biotechnology applications such as gene knockdown, RNA editing, and molecular diagnostics. Compared with Cas13a from mesophilic organisms, a newly discovered Cas13a from thermophilic bacteria Thermoclostridium caenicola (TccCas13a) shows low sequence similarity, high thermostability, and lacks pre-crRNA processing activity. The thermostability of TccCas13a has been harnessed to make a sensitive and robust tool for nucleic acid detection. Here we present the structures of TccCas13a-crRNA binary complex at 2.8 Å, and TccCas13a at 3.5 Å. Although TccCas13a shares a similarly bilobed architecture with other mesophilic organism-derived Cas13a proteins, TccCas13a displayed distinct structure features. Specifically, it holds a long crRNA 5′-flank, forming extensive polar contacts with Helical-1 and HEPN2 domains. The detailed analysis of the interaction between crRNA 5′-flank and TccCas13a suggested lack of suitable nucleophile to attack the 2′-OH of crRNA 5′-flank may explain why TccCas13a fails to cleave pre-crRNA. The stem-loop segment of crRNA spacer toggles between double-stranded and single-stranded conformational states, suggesting a potential safeguard mechanism for target recognition. Superimposition of the structures of TccCas13a and TccCas13a-crRNA revealed several conformational changes required for crRNA loading, including dramatic movement of Helical-2 domain. Collectively, these structural insights expand our understanding into type VI CRISPR-Cas effectors, and would facilitate the development of TccCas13a-based applications.     

A Novel Spectral Annotation Strategy Streamlines Reporting of Mono-ADP-ribosylated Peptides Derived from Mouse Liver and Spleen in Response to IFN-γ

Mass-spectrometry-enabled ADP-ribosylation workflows are developing rapidly, providing researchers a variety of ADP-ribosylome enrichment strategies and mass spectrometric acquisition options. Despite the growth spurt in upstream technologies, systematic ADP-ribosyl (ADPr) peptide mass spectral annotation methods are lacking. HCD-dependent ADP-ribosylome studies are common, but the resulting MS2 spectra are complex, owing to a mixture of b/y-ions and the m/p-ion peaks representing one or more dissociation events of the ADPr moiety (m-ion) and peptide (p-ion). In particular, p-ions that dissociate further into one or more fragment ions can dominate HCD spectra but are not recognized by standard spectral annotation workflows. As a result, annotation strategies that are solely reliant upon the b/y-ions result in lower spectral scores that in turn reduce the number of reportable ADPr peptides. To improve the confidence of spectral assignments, we implemented an ADPr peptide annotation and scoring strategy. All MS2 spectra are scored for the ADPr m-ions, but once spectra are assigned as an ADPr peptide, they are further annotated and scored for the p-ions. We implemented this novel workflow to ADPr peptides enriched from the liver and spleen isolated from mice post 4 h exposure to systemic IFN-γ. HCD collision energy experiments were first performed on the Orbitrap Fusion Lumos and the Q Exactive, with notable ADPr peptide dissociation properties verified with CID (Lumos). The m-ion and p-ion series score distributions revealed that ADPr peptide dissociation properties vary markedly between instruments and within instrument collision energy settings, with consequences on ADPr peptide reporting and amino acid localization. Consequentially, we increased the number of reportable ADPr peptides by 25% (liver) and 17% (spleen) by validation and the inclusion of lower confidence ADPr peptide spectra. This systematic annotation strategy will streamline future reporting of ADPr peptides that have been sequenced using any HCD/CID-based method.     

CNS anti-depressant,anxiolytic and analgesic effects of Ganoderma applanatum (mushroom) along with ligand-receptor binding screening provide new insights: Multi-disciplinary approaches

This research was designed to evaluate the CNS depressant, anxiolytic, and analgesic action of aqueous and ethanol extract of Ganoderma applanatum, a valuable medicinal fungus used in multiple disorders belongs to Ganodermataceae family. Two extracts of G. applanatum were prepared using distilled water and ethanol as solvents and named AEGA and EEGA. Open field method, rotarod method, tail suspension method, and hole cross method were utilized for the CNS depressant action. In contrast, elevated plus-maze test and hole board method were utilized for the anxiolytic action. For determining the analgesic potential, acetic acid-induced writhing test, hot plate method, and tail immersion test were used. Besides, molecular docking has been implemented by using Discovery studio 2020, UCSF Chimera and PyRx autodock vina. At both doses (200 and 400 mg/kg) of AEGA and EEGA showed significant CNS depressant effect (p < 0.05 to 0.001) against all four tests used for CNS depressant activity. Both doses of AEGA and EEGA exhibited important anxiolytic activity effect (p < 0.05 to 0.001)against the EPM and hole board test. Both doses of AEGA and EEGA also exhibited a potential analgesic effect (p < 0.05 to 0.001) against all three tests used for analgesic action. In addition, in the molecular docking the compounds obtained the scores of ?5.2 to ?12.8 kcal/mol. Ganoapplanin, sphaeropsidin D and cytosporone C showed the best binding affinity to the selected recptors. It can be concluded that AEGA and EEGA have potential CNS depressant, anxiolytic, and analgesic action, which can be used as a natural antidepressant, anxiolytic, and analgesic source.     

Crystal structure of a nucleotide-binding domain of fatty acid kinase FakA from Thermus thermophilus HB8

Fatty acid kinase is necessary for the incorporation of exogenous fatty acids into membrane phospholipids. Fatty acid kinase consists of two components: a kinase component, FakA, that phosphorylates a fatty acid bound to a fatty acid-binding component, FakB. However, the molecular details underlying the phosphotransfer reaction remain to be resolved. We determined the crystal structure of the N-terminal domain of FakA bound to ADP from Thermus thermophilus HB8. The overall structure of this domain showed that the helical barrel fold is similar to the nucleotide-binding component of dihydroxyacetone kinase. The structure of the nucleotide-binding site revealed the roles of the conserved residues in recognition of ADP and Mg²⁺, but the N-terminal domain of FakA lacked the ADP-capping loop found in the dihydroxyacetone kinase component. Based on the structural similarity to the two subunits of dihydroxyacetone kinase complex, we constructed a model of the complex of T. thermophilus FakB and the N-terminal domain of FakA. In this model, the invariant Arg residue of FakB occupied a position that was spatially similar to that of the catalytically important Arg residue of dihydroxyacetone kinase, which predicted a composite active site in the Fatty acid kinase complex.     

SP3-Enabled Rapid and High Coverage Chemoproteomic Identification of Cell-State–Dependent Redox-Sensitive Cysteines

Proteinaceous cysteine residues act as privileged sensors of oxidative stress. As reactive oxygen and nitrogen species have been implicated in numerous pathophysiological processes, deciphering which cysteines are sensitive to oxidative modification and the specific nature of these modifications is essential to understanding protein and cellular function in health and disease. While established mass spectrometry-based proteomic platforms have improved our understanding of the redox proteome, the widespread adoption of these methods is often hindered by complex sample preparation workflows, prohibitive cost of isotopic labeling reagents, and requirements for custom data analysis workflows. Here, we present the SP3-Rox redox proteomics method that combines tailored low cost isotopically labeled capture reagents with SP3 sample cleanup to achieve high throughput and high coverage proteome-wide identification of redox-sensitive cysteines. By implementing a customized workflow in the free FragPipe computational pipeline, we achieve accurate MS1-based quantitation, including for peptides containing multiple cysteine residues. Application of the SP3-Rox method to cellular proteomes identified cysteines sensitive to the oxidative stressor GSNO and cysteine oxidation state changes that occur during T cell activation.     

Temporal Analysis of Protein Ubiquitylation and Phosphorylation During Parkin-Dependent Mitophagy

Mitophagy, the selective degradation of mitochondria by autophagy, affects defective mitochondria following damage or stress. At the onset of mitophagy, parkin ubiquitylates proteins on the mitochondrial outer membrane. While the role of parkin at the onset of mitophagy is well understood, less is known about its activity during later stages in the process. Here, we used HeLa cells expressing catalytically active or inactive parkin to perform temporal analysis of the proteome, ubiquitylome, and phosphoproteome during 18 h after induction of mitophagy by mitochondrial uncoupler carbonyl cyanide m-chlorophenyl hydrazine. Abundance profiles of proteins downregulated in parkin-dependent manner revealed a stepwise and “outside–in” directed degradation of mitochondrial subcompartments. While ubiquitylation of mitochondrial outer membrane proteins was enriched among early parkin-dependent targets, numerous mitochondrial inner membrane, matrix, and cytosolic proteins were also found ubiquitylated at later stages of mitophagy. Phosphoproteome analysis revealed a possible crosstalk between phosphorylation and ubiquitylation during mitophagy on key parkin targets, such as voltage-dependent anion channel 2.     

Methylglyoxal attenuates isoproterenol-induced increase in uncoupling protein 1 expression through activation of JNK signaling pathway in beige adipocytes

Methylglyoxal (MG) is a metabolite derived from glycolysis whose levels in the blood and tissues of patients with diabetes are higher than those of healthy individuals, suggesting that MG is associated with the development of diabetic complications. However, it remains unknown whether high levels of MG are a cause or consequence of diabetes. Here, we show that MG negatively affects the expression of uncoupling protein 1 (UCP1), which is involved in thermogenesis and the regulation of systemic metabolism. Decreased Ucp1 expression is associated with obesity and type 2 diabetes. We found that MG attenuated the increase in Ucp1 expression following treatment with isoproterenol in beige adipocytes. However, MG did not affect protein kinase A signaling, the core coordinator of isoproterenol-induced Ucp1 expression. Instead, MG activated c-Jun N-terminal kinase (JNK) and p38 mitogen-activated protein kinases. We found that JNK inhibition, but not p38, recovered isoproterenol-stimulated Ucp1 expression under MG treatment. Altogether, these results suggest an inhibitory role of MG on the thermogenic function of beige adipocytes through the JNK signaling pathway.     

Streptomyces poriferorum sp. nov., a novel marine sponge-derived Actinobacteria species expressing anti-MRSA activity

Marine sponges represent a rich source of uncharacterized microbial diversity, and many are host to microorganisms that produce biologically active specialized metabolites. Here, a polyphasic approach was used to characterize two Actinobacteria strains, P01-B04^T and P01-F02, that were isolated from the marine sponges Geodia barretti (Bowerbank, 1858) and Antho dichotoma (Esper, 1794), respectively. Phylogenetic analysis based on 16S rRNA gene sequences indicated that strains P01-B04^T and P01-F02 are closely related to Streptomyces beijiangensis DSM 41794^T, Streptomyces laculatispora NRRL B-24909^T, and Streptomyces brevispora NRRL B-24910^T. The two strains showed nearly identical 16S rRNA gene sequences (99.93%), and the average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) relatedness values were 99.96% and 99.6%, respectively, suggesting that these strains are affiliated with the same species. Chemotaxonomic and culture characteristics of both strains were also consistent with the genus Streptomyces, while phenotypic properties, genome-based comparisons, and phylogenomic analyses distinguished strains P01-B04^T and P01-F02 from their closest phylogenetic relatives. In silico analysis predicted that the 8.9 Mb genome of P01-B04^T contains at least 41 biosynthetic gene clusters (BGCs) encoding secondary metabolites, indicating that this strain could express diverse bioactive metabolites; in support of this prediction, this strain expressed antibacterial activity against Gram-positive bacteria including a clinical isolate of methicillin-resistant Staphylococcus aureus (MRSA) EAMC30. Based on these results, the marine sponge-associated isolates represent a novel species of the genus Streptomyces, for which the name Streptomyces poriferorum sp. nov. is proposed, with P01-B04^T (=DSM 111306^T = CCM 9048^T) as the type strain.     

Mechanistic aspects of the transamination reactions catalyzed by D-amino acid transaminase from Haliscomenobacter hydrossis

Pyridoxal-5′-phosphate-(PLP-) dependent D-amino acid transaminases (DAATs) catalyze stereoselective reversible transfer of the amino group between D-amino acids and keto acids. In vivo DAATs are commonly known to synthesize D-glutamate for cell wall peptidoglycans. Today DAATs meet increasing attention for application in the synthesis of D-amino acids, whereas little is known about the mechanism of substrate recognition and catalytic steps of the D-amino acids conversion by DAATs. In this work, the pre-steady-state kinetics of the half-reactions of DAAT from Haliscomenobacter hydrossis with D-glutamate, D-alanine, D-leucine, and D-phenylalanine was examined at two wavelengths, 416 and 330 nm, using a stopped-flow technique. Monophasic kinetics was observed with specific substrates D-glutamate and D-alanine, whereas half-reactions with D-leucine and D-phenylalanine exhibited biphasic kinetics. All half-reactions proceeded until the complete conversion of PLP due to the release of the pyridoxamine-5′-phosphate form of cofactor from the holoenzyme . Comparison of kinetic parameters of half-reactions and the overall transamination reactions for D-leucine, D-phenylalanine revealed the increase in the rates of deamination of these substrates in the overall reaction with α-ketoglutarate. In the overall transamination reaction, the catalytic turnover rates for D-leucine and D-phenylalanine increased by 260 and 60 times, correspondingly, comparing with the slowest step rate constants in the half-reactions. We suggested the activating effect by a specific substrate α-ketoglutarate in the overall transamination reaction. The study of half-reactions helped to quantify the specificity of DAAT from H. hydrossis for D-amino acids with different properties. The results obtained are the first detailed analysis of half-reactions catalyzed by DAAT.     

Glyceraldehyde-derived advanced glycation end-products having pyrrolopyridinium-based crosslinks

Reducing sugars and reactive aldehydes, such as glyceraldehyde, non-enzymatically react with amino or guanidino groups of proteins to form advanced glycation end-products (AGEs) by the Maillard reaction that involves Schiff base formation followed by Amadori rearrangement. AGEs are found relatively in abundance in the human eye and to accumulate at a higher rate in diseases that impair vision such as cataract, diabetic retinopathy or age-related macular degeneration. We identified two novel AGEs of pyrrolopyridinium lysine dimer derived from glyceraldehyde, PPG1 and PPG2, in the Maillard reaction of N^α-acetyl-l-lysine with glyceraldehyde under physiological conditions. Having fluorophores similar to that of vesperlysine A, which was isolated from the human lens, PPGs were found to act as photosensitizers producing singlet oxygen in response to blue light irradiation. Moreover, PPG2 interacts with receptor for AGE (RAGE) in vitro with a higher binding affinity than GLAP, a well-known ligand of the receptor. We also proposed a pathway to form PPGs and discussed how they would be formed in vitro. As glyceraldehyde-derived AGEs have been studied extensively in connection with various hyperglycemia-related diseases, further studies will be required to find PPGs in vivo such as in the lens or other tissues.