Discovering potential interactions between rare diseases and COVID-19 by combining mechanistic models of viral infection with statistical modeling.
Ontology highlight
ABSTRACT:
SUBMITTER: Lopez-Sanchez M
PROVIDER: S-EPMC9239744 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:PURPOSE:We used human stem and progenitor cells to develop a genetically accurate novel model of MYC-driven Group 3 medulloblastoma. We also developed a new informatics method, Disease-model Signature versus Compound-Variety Enriched Response ("DiSCoVER"), to identify novel therapeutics that target this specific disease subtype. EXPERIMENTAL DESIGN:Human neural stem and progenitor cells derived from the cerebellar anlage were transduced with oncogenic elements associated with aggressive medulloblastoma. An in silico analysis method for screening drug sensitivity databases (DiSCoVER) was used in multiple drug sensitivity datasets. We validated the top hits from this analysis in vitro and in vivo RESULTS:Human neural stem and progenitor cells transformed with c-MYC, dominant-negative p53, constitutively active AKT and hTERT formed tumors in mice that recapitulated Group 3 medulloblastoma in terms of pathology and expression profile. DiSCoVER analysis predicted that aggressive MYC-driven Group 3 medulloblastoma would be sensitive to cyclin-dependent kinase (CDK) inhibitors. The CDK 4/6 inhibitor palbociclib decreased proliferation, increased apoptosis, and significantly extended the survival of mice with orthotopic medulloblastoma xenografts. CONCLUSIONS:We present a new method to generate genetically accurate models of rare tumors, and a companion computational methodology to find therapeutic interventions that target them. We validated our human neural stem cell model of MYC-driven Group 3 medulloblastoma and showed that CDK 4/6 inhibitors are active against this subgroup. Our results suggest that palbociclib is a potential effective treatment for poor prognosis MYC-driven Group 3 medulloblastoma tumors in carefully selected patients. Clin Cancer Res; 22(15); 3903-14. ©2016 AACR.
| S-EPMC5055054 | biostudies-literature
Project description:There is heightened interest in using next-generation sequencing technologies to identify rare variants that influence complex human diseases and traits. Meta-analysis is essential to this endeavor because large sample sizes are required for detecting associations with rare variants. In this article, we provide a comprehensive overview of statistical methods for meta-analysis of sequencing studies for discovering rare-variant associations. Specifically, we discuss the calculation of relevant summary statistics from participating studies, the construction of gene-level association tests, the choice of transformation for quantitative traits, the use of fixed-effects versus random-effects models, and the removal of shadow association signals through conditional analysis. We also show that meta-analysis based on properly calculated summary statistics is as powerful as joint analysis of individual-participant data. In addition, we demonstrate the performance of different meta-analysis methods by using both simulated and empirical data. We then compare four major software packages for meta-analysis of rare-variant associations-MASS, RAREMETAL, MetaSKAT, and seqMeta-in terms of the underlying statistical methodology, analysis pipeline, and software interface. Finally, we present PreMeta, a software interface that integrates the four meta-analysis packages and allows a consortium to combine otherwise incompatible summary statistics.
| S-EPMC4571037 | biostudies-literature
Project description:Knowledge of the time of HIV-1 infection and the multiplicity of viruses that establish HIV-1 infection is crucial for the in-depth analysis of clinical prevention efficacy trial outcomes. Better estimation methods would improve the ability to characterize immunological and genetic sequence correlates of efficacy within preventive efficacy trials of HIV-1 vaccines and monoclonal antibodies. We developed new methods for infection timing and multiplicity estimation using maximum likelihood estimators that shift and scale (calibrate) estimates by fitting true infection times and founder virus multiplicities to a linear regression model with independent variables defined by data on HIV-1 sequences, viral load, diagnostics, and sequence alignment statistics. Using Poisson models of measured mutation counts and phylogenetic trees, we analyzed longitudinal HIV-1 sequence data together with diagnostic and viral load data from the RV217 and CAPRISA 002 acute HIV-1 infection cohort studies. We used leave-one-out cross validation to evaluate the prediction error of these calibrated estimators versus that of existing estimators and found that both infection time and founder multiplicity can be estimated with improved accuracy and precision by calibration. Calibration considerably improved all estimators of time since HIV-1 infection, in terms of reducing bias to near zero and reducing root mean squared error (RMSE) to 5-10 days for sequences collected 1-2 months after infection. The calibration of multiplicity assessments yielded strong improvements with accurate predictions (ROC-AUC above 0.85) in all cases. These results have not yet been validated on external data, and the best-fitting models are likely to be less robust than simpler models to variation in sequencing conditions. For all evaluated models, these results demonstrate the value of calibration for improved estimation of founder multiplicity and of time since HIV-1 infection.
| S-EPMC6669737 | biostudies-literature
Project description:Recent advances in sequencing technologies have made it significantly easier to find the genetic roots of rare hereditary pediatric diseases. These novel methods are not panaceas, however, and they often give ambiguous results, highlighting multiple possible causative mutations in affected patients. Furthermore, even when the mapping results are unambiguous, the affected gene might be of unknown function. In these cases, understanding how a particular genotype can result in a phenotype also needs carefully designed experimental work. Model organism genetics can offer a straightforward experimental setup for hypothesis testing. Containing orthologs for over 80% of the genes involved in human diseases, zebrafish (Danio rerio) has emerged as one of the top disease models over the past decade. A plethora of genetic tools makes it easy to create mutations in almost any gene of the zebrafish genome and these mutant strains can be used in high-throughput preclinical screens for active molecules. As this small vertebrate species offers several other advantages as well, its popularity in biomedical research is bound to increase, with "aquarium to bedside" drug development pipelines taking a more prevalent role in the near future.
| S-EPMC6023479 | biostudies-literature
Project description:A protocol is presented for the global refinement of homology models of proteins. It combines the advantages of temperature-based replica-exchange molecular dynamics (REMD) for conformational sampling and the use of statistical potentials for model selection. The protocol was tested using 21 models. Of these 14 were models of 10 small proteins for which high-resolution crystal structures were available, the remainder were targets of the recent CASPR exercise. It was found that REMD in combination with currently available force fields could sample near-native conformational states starting from high-quality homology models. Conformations in which the backbone RMSD of secondary structure elements (SSE-RMSD) was lower than the starting value by 0.5-1.0 A were found for 15 out of the 21 cases (average 0.82 A). Furthermore, when a simple scoring function consisting of two statistical potentials was used to rank the structures, one or more structures with SSE-RMSD of at least 0.2 A lower than the starting value was found among the five best ranked structures in 11 out of the 21 cases. The average improvement in SSE-RMSD for the best models was 0.42 A. However, none of the scoring functions tested identified the structures with the lowest SSE-RMSD as the best models although all identified the native conformation as the one with lowest energy. This suggests that while the proposed protocol proved effective for the refinement of high-quality models of small proteins scoring functions remain one of the major limiting factors in structure refinement. This and other aspects by which the methodology could be further improved are discussed.
| S-EPMC2761145 | biostudies-literature
Project description:Motivation:Transcription factors bind regulatory DNA sequences in a combinatorial manner to modulate gene expression. Deep neural networks (DNNs) can learn the cis-regulatory grammars encoded in regulatory DNA sequences associated with transcription factor binding and chromatin accessibility. Several feature attribution methods have been developed for estimating the predictive importance of individual features (nucleotides or motifs) in any input DNA sequence to its associated output prediction from a DNN model. However, these methods do not reveal higher-order feature interactions encoded by the models. Results:We present a new method called Deep Feature Interaction Maps (DFIM) to efficiently estimate interactions between all pairs of features in any input DNA sequence. DFIM accurately identifies ground truth motif interactions embedded in simulated regulatory DNA sequences. DFIM identifies synergistic interactions between GATA1 and TAL1 motifs from in vivo TF binding models. DFIM reveals epistatic interactions involving nucleotides flanking the core motif of the Cbf1 TF in yeast from in vitro TF binding models. We also apply DFIM to regulatory sequence models of in vivo chromatin accessibility to reveal interactions between regulatory genetic variants and proximal motifs of target TFs as validated by TF binding quantitative trait loci. Our approach makes significant strides in improving the interpretability of deep learning models for genomics. Availability and implementation:Code is available at: https://github.com/kundajelab/dfim. Supplementary information:Supplementary data are available at Bioinformatics online.
| S-EPMC6129272 | biostudies-literature
Project description:Recent experimental and theoretical efforts have highlighted the fact that binding of transcription factors to DNA can be more accurately described by continuous measures of their binding affinities, rather than a discrete description in terms of binding sites. While the binding affinities can be predicted from a physical model, it is often desirable to know the distribution of binding affinities for specific sequence backgrounds. In this paper, we present a statistical approach to derive the exact distribution for sequence models with fixed GC content. We demonstrate that the affinity distribution of almost all known transcription factors can be effectively parametrized by a class of generalized extreme value distributions. Moreover, this parameterization also describes the affinity distribution for sequence backgrounds with variable GC content, such as human promoter sequences. Our approach is applicable to arbitrary sequences and all transcription factors with known binding preferences that can be described in terms of a motif matrix. The statistical treatment also provides a proper framework to directly compare transcription factors with very different affinity distributions. This is illustrated by our analysis of human promoters with known binding sites, for many of which we could identify the known regulators as those with the highest affinity. The combination of physical model and statistical normalization provides a quantitative measure which ranks transcription factors for a given sequence, and which can be compared directly with large-scale binding data. Its successful application to human promoter sequences serves as an encouraging example of how the method can be applied to other sequences.
| S-EPMC2266803 | biostudies-literature
Project description:Applications in structural biology and medicinal chemistry require protein-ligand scoring functions for two distinct tasks: (i) ranking different poses of a small molecule in a protein binding site and (ii) ranking different small molecules by their complementarity to a protein site. Using probability theory, we developed two atomic distance-dependent statistical scoring functions: PoseScore was optimized for recognizing native binding geometries of ligands from other poses and RankScore was optimized for distinguishing ligands from nonbinding molecules. Both scores are based on a set of 8,885 crystallographic structures of protein-ligand complexes but differ in the values of three key parameters. Factors influencing the accuracy of scoring were investigated, including the maximal atomic distance and non-native ligand geometries used for scoring, as well as the use of protein models instead of crystallographic structures for training and testing the scoring function. For the test set of 19 targets, RankScore improved the ligand enrichment (logAUC) and early enrichment (EF(1)) scores computed by DOCK 3.6 for 13 and 14 targets, respectively. In addition, RankScore performed better at rescoring than each of seven other scoring functions tested. Accepting both the crystal structure and decoy geometries with all-atom root-mean-square errors of up to 2 Å from the crystal structure as correct binding poses, PoseScore gave the best score to a correct binding pose among 100 decoys for 88% of all cases in a benchmark set containing 100 protein-ligand complexes. PoseScore accuracy is comparable to that of DrugScore(CSD) and ITScore/SE and superior to 12 other tested scoring functions. Therefore, RankScore can facilitate ligand discovery, by ranking complexes of the target with different small molecules; PoseScore can be used for protein-ligand complex structure prediction, by ranking different conformations of a given protein-ligand pair. The statistical potentials are available through the Integrative Modeling Platform (IMP) software package (http://salilab.org/imp) and the LigScore Web server (http://salilab.org/ligscore/).
| S-EPMC3246566 | biostudies-literature
Project description:This perspective outlines an approach to improve mechanistic understanding of macrophages in inflammation and tissue homeostasis, with a focus on human inflammatory bowel disease (IBD). The approach integrates wet-lab and in-silico experimentation, driven by mechanistic mathematical models of relevant biological processes. Although wet-lab experimentation with genetically modified mouse models and primary human cells and tissues have provided important insights, the role of macrophages in human IBD remains poorly understood. Key open questions include: (1) To what degree hyperinflammatory processes (e.g., gain of cytokine production) and immunodeficiency (e.g., loss of bacterial killing) intersect to drive IBD pathophysiology? and (2) What are the roles of macrophage heterogeneity in IBD onset and progression? Mathematical modeling offers a synergistic approach that can be used to address such questions. Mechanistic models are useful for informing wet-lab experimental designs and provide a knowledge constrained framework for quantitative analysis and interpretation of resulting experimental data. The majority of published mathematical models of macrophage function are based either on animal models, or immortalized human cell lines. These experimental models do not recapitulate important features of human gastrointestinal pathophysiology, and, therefore are limited in the extent to which they can fully inform understanding of human IBD. Thus, we envision a future where mechanistic mathematical models are based on features relevant to human disease and parametrized by richer human datasets, including biopsy tissues taken from IBD patients, human organ-on-a-chip systems and other high-throughput clinical data derived from experimental medicine studies and/or clinical trials on IBD patients.
| S-EPMC6563075 | biostudies-literature
Project description:We develop an R/Bioconductor package CancerInSilico to simulate bulk and single cell transcriptional data from a known ground truth obtained from mathematical models of cellular systems. This package contains a general R infrastructure for cell-based mathematical model, implemented for an off-lattice, cell-center Monte Carlo mathematical model. We also adapt this model to simulate the impact of growth suppression by targeted therapeutics in cancer and benchmark simulations against bulk in vitro experimental data. Sensitivity to parameters is evaluated and used to predict the relative impact of variation in cellular growth parameters and cell types on tumor heterogeneity in therapeutic response.
2020-07-09 | GSE114375 | GEO