A Pharmacophore Model for SARS-CoV-2 3CLpro Small Molecule Inhibitors and in Vitro Experimental Validation of Computationally Screened Inhibitors

in Journal of Chemical Information and Modeling (2021), 61(8), 4082-4096

Among the biomedical efforts in response to the current coronavirus (COVID-19) pandemic, pharmacological strategies to reduce viral load in patients with severe forms of the disease are being studied ... [more ▼]

Among the biomedical efforts in response to the current coronavirus (COVID-19) pandemic, pharmacological strategies to reduce viral load in patients with severe forms of the disease are being studied intensively. One of the main drug target proteins proposed so far is the SARS-CoV-2 viral protease 3CLpro (also called Mpro), an essential component for viral replication. Ongoing ligand- and receptor-based computational screening efforts would be facilitated by an improved understanding of the electrostatic, hydrophobic and steric features that characterize small molecule inhibitors binding stably to 3CLpro, as well as by an extended collection of known binders. Here, we present combined virtual screening, molecular dynamics simulation, machine learning and in vitro experimental validation analyses which have led to the identification of small molecule inhibitors of 3CLpro with micromolar activity, and to a pharmacophore model that describes functional chemical groups associated with the molecular recognition of ligands by the 3CLpro binding pocket. Experimentally validated inhibitors using a ligand activity assay include natural compounds with available prior knowledge on safety and bioavailability properties, such as the natural compound rottlerin (IC50 = 37 µM), and synthetic compounds previously not characterized (e.g. compound CID 46897844, IC50 = 31 µM). In combination with the developed pharmacophore model, these and other confirmed 3CLpro inhibitors may provide a basis for further similarity-based screening in independent compound databases and structural design optimization efforts, to identify 3CLpro ligands with improved potency and selectivity. Overall, this study suggests that the integration of virtual screening, molecular dynamics simulations and machine learning can facilitate 3CLpro-targeted small molecule screening investigations. Different receptor-, ligand- and machine learning-based screening strategies provided complementary information, helping to increase the number and diversity of identified active compounds. Finally, the resulting pharmacophore model and experimentally validated small molecule inhibitors for 3CLpro provide resources to support follow-up computational screening efforts for this drug target. [less ▲]

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms.

in Molecular systems biology (2021), 17(10), 10387

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets ... [more ▼]

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective. [less ▲]

Biomarker discovery studies for patient stratification using machine learning analysis of omics data: a scoping review

in BMJ Open (2021), 11(12), 053674

Objective: To review biomarker discovery studies using omics data for patient stratification which led to clinically validated FDA-cleared tests or laboratory developed tests, in order to identify common ... [more ▼]

Objective: To review biomarker discovery studies using omics data for patient stratification which led to clinically validated FDA-cleared tests or laboratory developed tests, in order to identify common characteristics and derive recommendations for future biomarker projects. Design: Scoping review. Methods: We searched PubMed, EMBASE and Web of Science to obtain a comprehensive list of articles from the biomedical literature published between January 2000 and July 2021, describing clinically validated biomarker signatures for patient stratification, derived using statistical learning approaches. All documents were screened to retain only peer-reviewed research articles, review articles or opinion articles, covering supervised and unsupervised machine learning applications for omics-based patient stratification. Two reviewers independently confirmed the eligibility. Disagreements were solved by consensus. We focused the final analysis on omics-based biomarkers which achieved the highest level of validation, that is, clinical approval of the developed molecular signature as a laboratory developed test or FDA approved tests. Results: Overall, 352 articles fulfilled the eligibility criteria. The analysis of validated biomarker signatures identified multiple common methodological and practical features that may explain the successful test development and guide future biomarker projects. These include study design choices to ensure sufficient statistical power for model building and external testing, suitable combinations of non-targeted and targeted measurement technologies, the integration of prior biological knowledge, strict filtering and inclusion/exclusion criteria, and the adequacy of statistical and machine learning methods for discovery and validation. Conclusions: While most clinically validated biomarker models derived from omics data have been developed for personalised oncology, first applications for non-cancer diseases show the potential of multivariate omics biomarker design for other complex disorders. Distinctive characteristics of prior success stories, such as early filtering and robust discovery approaches, continuous improvements in assay design and experimental measurement technology, and rigorous multicohort validation approaches, enable the derivation of specific recommendations for future studies. [less ▲]

The benefits of low COVID-19 incidence in Europe

in The Lancet (2021), 398(10303), 838-839

How will the coronavirus disease 2019 (COVID-19) pandemic develop in the coming months and years? Based on an expert survey, we examine key aspects that are likely to influence the COVID-19 pandemic in ... [more ▼]

How will the coronavirus disease 2019 (COVID-19) pandemic develop in the coming months and years? Based on an expert survey, we examine key aspects that are likely to influence the COVID-19 pandemic in Europe. The challenges and developments will strongly depend on the progress of national and global vaccination programs, the emergence and spread of variants of concern (VOCs), and public responses to non-pharmaceutical interventions (NPIs). In the short term, many people remain unvaccinated, VOCs continue to emerge and spread, and mobility and population mixing are expected to increase. Therefore, lifting restrictions too much and too early risk another damaging wave. This challenge remains despite the reduced opportunities for transmission given vaccination progress and reduced indoor mixing in summer 2021. In autumn 2021, increased indoor activity might accelerate the spread again, whilst a necessary reintroduction of NPIs might be too slow. The incidence may strongly rise again, possibly filling intensive care units, if vaccination levels are not high enough. A moderate, adaptive level of NPIs will thus remain necessary. These epidemiological aspects combined with economic, social, and health-related consequences provide a more holistic perspective on the future of the COVID-19 pandemic. [less ▲]

Altered sphingolipid function in Alzheimer's disease; a gene regulatory network approach

in Neurobiology of Aging (2021), in press(in press),

Sphingolipids (SLs) are bioactive lipids involved in various important physiological functions. The SL pathway has been shown to be affected in several brain-related disorders, including Alzheimer’s ... [more ▼]

Sphingolipids (SLs) are bioactive lipids involved in various important physiological functions. The SL pathway has been shown to be affected in several brain-related disorders, including Alzheimer’s disease (AD). Recent evidence suggests that epigenetic dysregulation plays an important role in the pathogenesis of AD as well. Here, we use an integrative approach to better understand the relationship between epigenetic and transcriptomic processes in regulating SL function in the middle temporal gyrus of AD patients. Transcriptomic analysis of 252 SL-related genes, selected based on GO term annotations, from 46 AD patients and 32 healthy age-matched controls, revealed 103 differentially expressed SL-related genes in AD patients. Additionally, methylomic analysis of the same subjects revealed parallel hydroxymethylation changes in PTGIS, GBA, and ITGB2 in AD. Subsequent gene regulatory network-based analysis identified three candidate genes, i.e. SELPLG, SPHK1 and CAV1 whose alteration holds the potential to revert the gene expression program from a diseased towards a healthy state. Together, this epigenomic and transcriptomic approach highlights the importance of SL-related genes in AD, and may provide novel biomarkers and therapeutic alternatives to traditionally investigated biological pathways in AD. [less ▲]

Characterization of DNA Methylomic signatures in induced pluripotent stem cells during neuronal differentiation

in Frontiers in Cell and Developmental Biology (2021)

In development, differentiation from a pluripotent state results in global epigenetic changes, although the extent to which this occurs in induced pluripotent stem cell based neuronal models has not been ... [more ▼]

In development, differentiation from a pluripotent state results in global epigenetic changes, although the extent to which this occurs in induced pluripotent stem cell based neuronal models has not been extensively characterized. In the present study, induced pluripotent stem cell colonies (33Qn1 line) were differentiated and collected at four time-points, with DNA methylation assessed using the Illumina Infinium Human Methylation EPIC BeadChip array. Dynamic changes in DNA methylation occurring during differentiation were investigated using a data-driven trajectory inference method. We identified a large number of Bonferroni-significant loci that showed progressive alterations in DNA methylation during neuronal differentiation. A gene-gene interaction network analysis identified 60 densely connected genes that were influential in the differentiation of neurons, with STAT3 being the gene with the highest connectivity. [less ▲]

Predictive and interpretable models via the stacked elastic net

in Bioinformatics (2021), 37(14), 20122016

Motivation: Machine learning in the biomedical sciences should ideally provide predictive and interpretable models. When predicting outcomes from clinical or molecular features, applied researchers often ... [more ▼]

Motivation: Machine learning in the biomedical sciences should ideally provide predictive and interpretable models. When predicting outcomes from clinical or molecular features, applied researchers often want to know which features have effects, whether these effects are positive or negative, and how strong these effects are. Regression analysis includes this information in the coefficients but typically renders less predictive models than more advanced machine learning techniques. Results: Here we propose an interpretable meta-learning approach for high-dimensional regression. The elastic net provides a compromise between estimating weak effects for many features and strong effects for some features. It has a mixing parameter to weight between ridge and lasso regularisation. Instead of selecting one weighting by tuning, we combine multiple weightings by stacking. We do this in a way that increases predictivity without sacrificing interpretability. Availability and Implementation: The R package starnet is available on GitHub: https://github.com/rauschenberger/starnet. [less ▲]

Scoping review on machine learning methods for stratification

Presentation (2020, December 02)

A rare loss-of function variant of ADAM17 is associated with late-onset familial Alzheimer disease

in Molecular Psychiatry (2020), 25(3), 629-639

Common variants of about 20 genes contributing to AD risk have so far been identified through genome-wide association studies (GWAS). However, there is still a large proportion of heritability that might ... [more ▼]

Common variants of about 20 genes contributing to AD risk have so far been identified through genome-wide association studies (GWAS). However, there is still a large proportion of heritability that might be explained by rare but functionally important variants. One of the so far identified genes with rare AD causing variants is ADAM10. Using whole-genome sequencing we now identified a single rare nonsynonymous variant (SNV) rs142946965 [p.R215I] in ADAM17 co-segregating with an autosomal-dominant pattern of late-onset AD in one family. Subsequent genotyping and analysis of available whole-exome sequencing data of additional case/control samples from Germany, the UK and the USA identified five variant carriers among AD patients only. The mutation inhibits pro-protein cleavage and the formation of the active enzyme, thus leading to loss-of-function of ADAM17 α-secretase. Further, we identified a strong negative correlation between ADAM17 and APP gene expression in human brain and present in vitro evidence that ADAM17 negatively controls the expression of APP. As a consequence, p.R215I mutation of ADAM17 leads to elevated Aß formation in vitro. Together our data supports a causative association of the identified ADAM17 variant in the pathogenesis of AD. [less ▲]

Gene selection for optimal prediction of cell position in tissues from single-cell transcriptomics

in Life Science Alliance (2020), 3(11), 202000867

Single-cell RNA-seq (scRNAseq) technologies are rapidly evolving. While very informative, in standard scRNAseq experiments the spatial organization of the cells in the tissue of origin is lost. Conversely ... [more ▼]

Single-cell RNA-seq (scRNAseq) technologies are rapidly evolving. While very informative, in standard scRNAseq experiments the spatial organization of the cells in the tissue of origin is lost. Conversely, spatial RNA-seq technologies designed to maintain cell localization have limited throughput and gene coverage. Mapping scRNAseq to genes with spatial information increases coverage while providing spatial location. However, methods to perform such mapping have not yet been benchmarked. To fill this gap, we organized the DREAM Single-Cell Transcriptomics challenge focused on the spatial reconstruction of cells from the Drosophila embryo from scRNAseq data, leveraging as silver standard, genes with in situ hybridization data from the Berkeley Drosophila Transcription Network Project reference atlas. The 34 participating teams used diverse algorithms for gene selection and location prediction, while being able to correctly localize clusters of cells. Selection of predictor genes was essential for this task. Predictor genes showed a relatively high expression entropy, high spatial clustering and included prominent developmental genes such as gap and pair-rule genes and tissue markers. Application of the Top-10 methods to a zebrafish embryo dataset yielded similar performance and statistical properties of the selected genes than in the Drosophila data. This suggests that methods developed in this challenge are able to extract generalizable properties of genes that are useful to accurately reconstruct the spatial arrangement of cells in tissues. [less ▲]