ABSTRACT
The currently accepted intestinal epithelial cell organization model proposes that Lgr5+ crypt-base columnar (CBC) cells represent the sole intestinal stem cell (ISC) compartment. However, previous studies have indicated that Lgr5+ cells are dispensable for intestinal regeneration, leading to two major hypotheses: one favoring the presence of a quiescent reserve ISC and the other calling for differentiated cell plasticity. To investigate these possibilities, we studied crypt epithelial cells in an unbiased fashion via high-resolution single-cell profiling. These studies, combined with in vivo lineage tracing, show that Lgr5 is not a specific ISC marker and that stemness potential exists beyond the crypt base and resides in the isthmus region, where undifferentiated cells participate in intestinal homeostasis and regeneration following irradiation (IR) injury. Our results provide an alternative model of intestinal epithelial cell organization, suggesting that stemness potential is not restricted to CBC cells, and neither de-differentiation nor reserve ISC are drivers of intestinal regeneration.
Subject(s)
Homeostasis , Intestinal Mucosa , Receptors, G-Protein-Coupled , Regeneration , Stem Cells , Animals , Stem Cells/metabolism , Stem Cells/cytology , Mice , Intestinal Mucosa/metabolism , Receptors, G-Protein-Coupled/metabolism , Intestines/cytology , Cell Differentiation , Mice, Inbred C57BL , Epithelial Cells/metabolism , Single-Cell Analysis , MaleABSTRACT
Despite considerable efforts, the mechanisms linking genomic alterations to the transcriptional identity of cancer cells remain elusive. Integrative genomic analysis, using a network-based approach, identified 407 master regulator (MR) proteins responsible for canalizing the genetics of individual samples from 20 cohorts in The Cancer Genome Atlas (TCGA) into 112 transcriptionally distinct tumor subtypes. MR proteins could be further organized into 24 pan-cancer, master regulator block modules (MRBs), each regulating key cancer hallmarks and predictive of patient outcome in multiple cohorts. Of all somatic alterations detected in each individual sample, >50% were predicted to induce aberrant MR activity, yielding insight into mechanisms linking tumor genetics and transcriptional identity and establishing non-oncogene dependencies. Genetic and pharmacological validation assays confirmed the predicted effect of upstream mutations and MR activity on downstream cellular identity and phenotype. Thus, co-analysis of mutational and gene expression profiles identified elusive subtypes and provided testable hypothesis for mechanisms mediating the effect of genetic alterations.
Subject(s)
Neoplasms/genetics , Transcription, Genetic , Adenocarcinoma/genetics , Animals , Cell Line, Tumor , Colonic Neoplasms/genetics , Gene Expression Regulation, Neoplastic , Gene Regulatory Networks , Genome, Human , HEK293 Cells , Humans , Mice, Nude , Mutation/genetics , Reproducibility of ResultsABSTRACT
Understanding the organizational logic of neural circuits requires deciphering the biological basis of neuronal diversity and identity, but there is no consensus on how neuron types should be defined. We analyzed single-cell transcriptomes of a set of anatomically and physiologically characterized cortical GABAergic neurons and conducted a computational genomic screen for transcriptional profiles that distinguish them from one another. We discovered that cardinal GABAergic neuron types are delineated by a transcriptional architecture that encodes their synaptic communication patterns. This architecture comprises 6 categories of â¼40 gene families, including cell-adhesion molecules, transmitter-modulator receptors, ion channels, signaling proteins, neuropeptides and vesicular release components, and transcription factors. Combinatorial expression of select members across families shapes a multi-layered molecular scaffold along the cell membrane that may customize synaptic connectivity patterns and input-output signaling properties. This molecular genetic framework of neuronal identity integrates cell phenotypes along multiple axes and provides a foundation for discovering and classifying neuron types.
Subject(s)
GABAergic Neurons/cytology , Gene Expression Profiling , Single-Cell Analysis , Animals , Cell Adhesion Molecules, Neuronal/metabolism , Extracellular Matrix/metabolism , GABAergic Neurons/metabolism , Mice , Receptors, GABA/metabolism , Receptors, Ionotropic Glutamate/metabolism , Signal Transduction , Synapses , Transcription, Genetic , Zinc/metabolism , gamma-Aminobutyric Acid/metabolismABSTRACT
Regulation of gene expression is central to many biological processes. Although reconstruction of regulatory circuits from genomic data alone is therefore desirable, this remains a major computational challenge. Comparative approaches that examine the conservation and divergence of circuits and their components across strains and species can help reconstruct circuits as well as provide insights into the evolution of gene regulatory processes and their adaptive contribution. In recent years, advances in genomic and computational tools have led to a wealth of methods for such analysis at the sequence, expression, pathway, module, and entire network level. Here, we review computational methods developed to study transcriptional regulatory networks using comparative genomics, from sequence to functional data. We highlight how these methods use evolutionary conservation and divergence to reliably detect regulatory components as well as estimate the extent and rate of divergence. Finally, we discuss the promise and open challenges in linking regulatory divergence to phenotypic divergence and adaptation.
Subject(s)
Gene Expression Regulation/genetics , Gene Regulatory Networks/genetics , Adaptation, Physiological/genetics , Animals , Computational Biology/methods , Evolution, Molecular , Genome/genetics , Genomics/methods , HumansABSTRACT
We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-ß dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.
Subject(s)
Genomics/methods , Neoplasms , Adolescent , Adult , Aged , Aged, 80 and over , Child , Female , Humans , Interferon-gamma/genetics , Interferon-gamma/immunology , Macrophages/immunology , Male , Middle Aged , Neoplasms/classification , Neoplasms/genetics , Neoplasms/immunology , Prognosis , Th1-Th2 Balance/physiology , Transforming Growth Factor beta/genetics , Transforming Growth Factor beta/immunology , Wound Healing/genetics , Wound Healing/immunology , Young AdultABSTRACT
Global migrations of diverse animal species often converge along the same routes, bringing together seasonal assemblages of animals that may compete, prey on each other, and share information or pathogens. These interspecific interactions, when energetic demands are high and the time to complete journeys is short, may influence survival, migratory success, stopover ecology, and migratory routes. Numerous accounts suggest that interspecific co-migrations are globally distributed in aerial, aquatic, and terrestrial systems, although the study of migration to date has rarely investigated species interactions among migrating animals. Here, we test the hypothesis that migrating animals are communities engaged in networks of ecological interactions. We leverage over half a million records of 50 bird species from five bird banding sites collected over 8 to 23 y to test for species associations using social network analyses. We find strong support for persistent species relationships across sites and between spring and fall migration. These relationships may be ecologically meaningful: They are often stronger among phylogenetically related species with similar foraging behaviors and nonbreeding ranges even after accounting for the nonsocial contributions to associations, including overlap in migration timing and habitat use. While interspecific interactions could result in costly competition or beneficial information exchange, we find that relationships are largely positive, suggesting limited competitive exclusion at the scale of a banding station during migratory stopovers. Our findings support an understanding of animal migrations that consist of networked communities rather than random assemblages of independently migrating species, encouraging future studies of the nature and consequences of co-migrant interactions.
Subject(s)
Animal Migration , Birds , Ecosystem , Seasons , Animals , Animal Migration/physiology , Birds/physiologyABSTRACT
Untargeted metabolomics based on liquid chromatography-mass spectrometry technology is quickly gaining widespread application, given its ability to depict the global metabolic pattern in biological samples. However, the data are noisy and plagued by the lack of clear identity of data features measured from samples. Multiple potential matchings exist between data features and known metabolites, while the truth can only be one-to-one matches. Some existing methods attempt to reduce the matching uncertainty, but are far from being able to remove the uncertainty for most features. The existence of the uncertainty causes major difficulty in downstream functional analysis. To address these issues, we develop a novel approach for Bayesian Analysis of Untargeted Metabolomics data (BAUM) to integrate previously separate tasks into a single framework, including matching uncertainty inference, metabolite selection and functional analysis. By incorporating the knowledge graph between variables and using relatively simple assumptions, BAUM can analyze datasets with small sample sizes. By allowing different confidence levels of feature-metabolite matching, the method is applicable to datasets in which feature identities are partially known. Simulation studies demonstrate that, compared with other existing methods, BAUM achieves better accuracy in selecting important metabolites that tend to be functionally consistent and assigning confidence scores to feature-metabolite matches. We analyze a COVID-19 metabolomics dataset and a mouse brain metabolomics dataset using BAUM. Even with a very small sample size of 16 mice per group, BAUM is robust and stable. It finds pathways that conform to existing knowledge, as well as novel pathways that are biologically plausible.
Subject(s)
Metabolomics , Mice , Animals , Bayes Theorem , Sample Size , Uncertainty , Metabolomics/methods , Computer SimulationABSTRACT
Veterinary systems biology is an innovative approach that integrates biological data at the molecular and cellular levels, allowing for a more extensive understanding of the interactions and functions of complex biological systems in livestock and veterinary science. It has tremendous potential to integrate multi-omics data with the support of vetinformatics resources for bridging the phenotype-genotype gap via computational modeling. To understand the dynamic behaviors of complex systems, computational models are frequently used. It facilitates a comprehensive understanding of how a host system defends itself against a pathogen attack or operates when the pathogen compromises the host's immune system. In this context, various approaches, such as systems immunology, network pharmacology, vaccinology and immunoinformatics, can be employed to effectively investigate vaccines and drugs. By utilizing this approach, we can ensure the health of livestock. This is beneficial not only for animal welfare but also for human health and environmental well-being. Therefore, the current review offers a detailed summary of systems biology advancements utilized in veterinary sciences, demonstrating the potential of the holistic approach in disease epidemiology, animal welfare and productivity.
Subject(s)
Animal Welfare , Systems Biology , Animals , Computational Biology , Computer Simulation , Genotype , PhenotypeABSTRACT
Micro ribonucleic acids (miRNAs) play a pivotal role in governing the human transcriptome in various biological phenomena. Hence, the accumulation of miRNA expression dysregulation frequently assumes a noteworthy role in the initiation and progression of complex diseases. However, accurate identification of dysregulated miRNAs still faces challenges at the current stage. Several bioinformatics tools have recently emerged for forecasting the associations between miRNAs and diseases. Nonetheless, the existing reference tools mainly identify the miRNA-disease associations in a general state and fall short of pinpointing dysregulated miRNAs within a specific disease state. Additionally, no studies adequately consider miRNA-miRNA interactions (MMIs) when analyzing the miRNA-disease associations. Here, we introduced a systematic approach, called IDMIR, which enabled the identification of expression dysregulated miRNAs through an MMI network under the gene expression context, where the network's architecture was designed to implicitly connect miRNAs based on their shared biological functions within a particular disease context. The advantage of IDMIR is that it uses gene expression data for the identification of dysregulated miRNAs by analyzing variations in MMIs. We illustrated the excellent predictive power for dysregulated miRNAs of the IDMIR approach through data analysis on breast cancer and bladder urothelial cancer. IDMIR could surpass several existing miRNA-disease association prediction approaches through comparison. We believe the approach complements the deficiencies in predicting miRNA-disease association and may provide new insights and possibilities for diagnosing and treating diseases. The IDMIR approach is now available as a free R package on CRAN (https://CRAN.R-project.org/package=IDMIR).
Subject(s)
Computational Biology , Gene Regulatory Networks , MicroRNAs , Urinary Bladder Neoplasms , Humans , MicroRNAs/genetics , MicroRNAs/metabolism , Computational Biology/methods , Urinary Bladder Neoplasms/genetics , Urinary Bladder Neoplasms/metabolism , Breast Neoplasms/genetics , Breast Neoplasms/metabolism , Gene Expression Profiling , Female , Gene Expression Regulation, NeoplasticABSTRACT
Distinct molecular pathways govern the differentiation of CD8+ effector T cells into memory or exhausted T cells during acute and chronic viral infection, but these are not well studied in humans. Here, we employed an integrative systems immunology approach to identify transcriptional commonalities and differences between virus-specific CD8+ T cells from patients with persistent and spontaneously resolving hepatitis C virus (HCV) infection during the acute phase. We observed dysregulation of metabolic processes during early persistent infection that was linked to changes in expression of genes related to nucleosomal regulation of transcription, T cell differentiation, and the inflammatory response and correlated with subject age, sex, and the presence of HCV-specific CD4+ T cell populations. These early changes in HCV-specific CD8+ T cell transcription preceded the overt establishment of T cell exhaustion, making this signature a prime target in the search for the regulatory origins of T cell dysfunction in chronic viral infection.
Subject(s)
CD8-Positive T-Lymphocytes/immunology , Hepacivirus/immunology , Hepatitis C, Chronic/immunology , Transcription, Genetic/immunology , Acute Disease , Adaptive Immunity/genetics , Adaptive Immunity/immunology , Adult , Aged , CD4-Positive T-Lymphocytes/immunology , CD4-Positive T-Lymphocytes/metabolism , CD4-Positive T-Lymphocytes/virology , CD8-Positive T-Lymphocytes/metabolism , CD8-Positive T-Lymphocytes/virology , Cluster Analysis , Female , Gene Expression Profiling/methods , Gene Regulatory Networks/immunology , Genetic Variation/immunology , Hepacivirus/physiology , Hepatitis C, Chronic/genetics , Hepatitis C, Chronic/virology , Humans , Lymphocyte Activation/genetics , Lymphocyte Activation/immunology , Male , Middle Aged , Multivariate Analysis , Time Factors , Young AdultABSTRACT
Notch signaling plays a critical role in cell fate decisions in all cell types. Furthermore, gain-of-function mutations in NOTCH1 have been uncovered in many human cancers. Disruption of Notch signaling has recently emerged as an attractive disease treatment strategy. However, the nuclear interaction landscape of the oncoprotein NOTCH1 remains largely unexplored. We therefore employed here a proximity-dependent biotin identification approach to identify in vivo protein associations with the nuclear Notch1 intracellular domain in live cells. We identified a large set of previously reported and unreported proteins that associate with NOTCH1, including general transcription and elongation factors, DNA repair and replication factors, coactivators, corepressors, and components of the NuRD and SWI/SNF chromatin remodeling complexes. We also found that Notch1 intracellular domain associates with protein modifiers and components of other signaling pathways that may influence Notch signal transduction and protein stability such as USP7. We further validated the interaction of NOTCH1 with histone deacetylase 1 or GATAD2B using protein network analysis, proximity-based ligation, in vivo cross-linking and coimmunoprecipitation assays in several Notch-addicted cancer cell lines. Through data mining, we also revealed potential drug targets for the inhibition of Notch signaling. Collectively, these results provide a valuable resource to uncover the mechanisms that fine-tune Notch signaling in tumorigenesis and inform therapeutic targets for Notch-addicted tumors.
Subject(s)
Carcinogenesis , Neoplasms , Oncogene Proteins , Receptor, Notch1 , Humans , Cell Differentiation , Cell Line , Oncogene Proteins/genetics , Oncogene Proteins/metabolism , Receptor, Notch1/genetics , Receptor, Notch1/metabolism , Receptors, Notch/metabolism , Signal Transduction , Ubiquitin-Specific Peptidase 7/metabolism , Carcinogenesis/genetics , Carcinogenesis/metabolism , Neoplasms/genetics , Neoplasms/metabolismABSTRACT
Frontotemporal dementia (FTD) is a primary cause of dementia encompassing a broad range of clinical phenotypes and cellular pathologies. Genetic discoveries in FTD have largely been driven by linkage studies in well-documented extended families, explaining most of the patients with a known pathogenic mutation. In the context of complex diseases, it is hypothesized that mutations with reduced penetrance or a combination of low-effect size variants with environmental factors drive disease. Furthermore, these genes are likely to be part of the interaction networks of known FTD genes, contributing to converging cellular processes. In this review, we examine gene discovery approaches in FTD and introduce network biology concepts as tools to assist gene identification studies in genetically complex disease.
Subject(s)
Frontotemporal Dementia , Frontotemporal Dementia/genetics , Frontotemporal Dementia/pathology , Genetic Linkage , Humans , Mutation , PhenotypeABSTRACT
Understanding drug selectivity mechanism is a long-standing issue for helping design drugs with high specificity. Designing drugs targeting cyclin-dependent kinases (CDKs) with high selectivity is challenging because of their highly conserved binding pockets. To reveal the underlying general selectivity mechanism, we carried out comprehensive analyses from both the thermodynamics and kinetics points of view on a representative CDK12 inhibitor. To fully capture the binding features of the drug-target recognition process, we proposed to use kinetic residue energy analysis (KREA) in conjunction with the community network analysis (CNA) to reveal the underlying cooperation effect between individual residues/protein motifs to the binding/dissociating process of the ligand. The general mechanism of drug selectivity in CDKs can be summarized as that the difference of structural cooperation between the ligand and the protein motifs leads to the difference of the energetic contribution of the key residues to the ligand. The proposed mechanisms may be prevalent in drug selectivity issues, and the insights may help design new strategies to overcome/attenuate the drug selectivity associated problems.
Subject(s)
Cyclin-Dependent Kinases , Molecular Dynamics Simulation , Cyclin-Dependent Kinases/metabolism , Ligands , Protein Binding , ThermodynamicsABSTRACT
BACKGROUND: It is valuable to analyze the genome-wide association studies (GWAS) data for a complex disease phenotype in the context of the protein-protein interaction (PPI) network, as the related pathophysiology results from the function of interacting polyprotein pathways. The analysis may include the design and curation of a phenotype-specific GWAS meta-database incorporating genotypic and eQTL data linking to PPI and other biological datasets, and the development of systematic workflows for PPI network-based data integration toward protein and pathway prioritization. Here, we pursued this analysis for blood pressure (BP) regulation. METHODS: The relational scheme of the implemented in Microsoft SQL Server BP-GWAS meta-database enabled the combined storage of: GWAS data and attributes mined from GWAS Catalog and the literature, Ensembl-defined SNP-transcript associations, and GTEx eQTL data. The BP-protein interactome was reconstructed from the PICKLE PPI meta-database, extending the GWAS-deduced network with the shortest paths connecting all GWAS-proteins into one component. The shortest-path intermediates were considered as BP-related. For protein prioritization, we combined a new integrated GWAS-based scoring scheme with two network-based criteria: one considering the protein role in the reconstructed by shortest-path (RbSP) interactome and one novel promoting the common neighbors of GWAS-prioritized proteins. Prioritized proteins were ranked by the number of satisfied criteria. RESULTS: The meta-database includes 6687 variants linked with 1167 BP-associated protein-coding genes. The GWAS-deduced PPI network includes 1065 proteins, with 672 forming a connected component. The RbSP interactome contains 1443 additional, network-deduced proteins and indicated that essentially all BP-GWAS proteins are at most second neighbors. The prioritized BP-protein set was derived from the union of the most BP-significant by any of the GWAS-based or the network-based criteria. It included 335 proteins, with ~ 2/3 deduced from the BP PPI network extension and 126 prioritized by at least two criteria. ESR1 was the only protein satisfying all three criteria, followed in the top-10 by INSR, PTN11, CDK6, CSK, NOS3, SH2B3, ATP2B1, FES and FINC, satisfying two. Pathway analysis of the RbSP interactome revealed numerous bioprocesses, which are indeed functionally supported as BP-associated, extending our understanding about BP regulation. CONCLUSIONS: The implemented workflow could be used for other multifactorial diseases.
Subject(s)
Genome-Wide Association Study , Protein Interaction Maps , Humans , Protein Interaction Maps/genetics , Genome-Wide Association Study/methods , Blood Pressure/genetics , Genotype , Databases, Factual , Plasma Membrane Calcium-Transporting ATPasesABSTRACT
Monogenic diseases are well-suited paradigms for the causal analysis of disease-driving molecular patterns. Spinal Muscular Atrophy (SMA) is one such monogenic model caused by mutation or deletion of the Survival of motor neuron 1 (SMN1) gene. Although several functions of the SMN protein have been studied, single functions and pathways alone do not allow to identify critical disease-driving molecules. Here, we analyzed the systemic characteristics of SMA employing proteomics, phosphoproteomics, translatomics and interactomics from two mouse models with different disease-severities and genetics. This systems approach revealed sub-networks and proteins characterizing commonalities and differences of both models. To link the identified molecular networks with the disease-causing SMN protein, we combined SMN-interactome data with both proteomes creating a comprehensive representation of SMA. By this approach, disease hubs and bottlenecks between SMN and downstream pathways could be identified. Linking a disease-causing molecule with widespread molecular dysregulations via multiomics is a concept for analyses of monogenic diseases.
ABSTRACT
Working memory, which is foundational to higher cognitive function, is the "sketchpad of volitional control." Successful working memory is the inevitable outcome of the individual's active control and manipulation of thoughts and turning them into internal goals during which the causal brain processes information in real time. However, little is known about the dynamic causality among distributed brain regions behind thought control that underpins successful working memory. In our present study, given that correct responses and incorrect ones did not differ in either contralateral delay activity or alpha suppression, further rooting on the high-temporal-resolution EEG time-varying directed network analysis, we revealed that successful working memory depended on both much stronger top-down connections from the frontal to the temporal lobe and bottom-up linkages from the occipital to the temporal lobe, during the early maintenance period, as well as top-down flows from the frontal lobe to the central areas as the delay behavior approached. Additionally, the correlation between behavioral performance and casual interactions increased over time, especially as memory-guided delayed behavior approached. Notably, when using the network metrics as features, time-resolved multiple linear regression of overall behavioral accuracy was exactly achieved as delayed behavior approached. These results indicate that accurate memory depends on dynamic switching of causal network connections and shifting to more task-related patterns during which the appropriate intervention may help enhance memory.
Subject(s)
Brain , Memory, Short-Term , Memory, Short-Term/physiology , Brain/physiology , Temporal Lobe/physiology , Frontal Lobe/physiology , Brain MappingABSTRACT
Identifying novel breast cancer biomarkers will improve patient stratification, enhance therapeutic outcomes, and help develop non-invasive diagnostics. We compared the proteomic profiles of whole-cell and exosomal samples of representative breast cancer cell subtypes to evaluate the potential of extracellular vesicles as non-invasive disease biomarkers in liquid biopsies. Overall, differentially-expressed proteins in whole-cell and exosome samples (which included markers for invasion, metastasis, angiogenesis, and drug resistance) effectively discriminated subtypes; furthermore, our results confirmed that the proteomic profile of exosomes reflects breast cancer cell-of-origin, which underscores their potential as disease biomarkers. Our study will contribute to identifying biomarkers that support breast cancer patient stratification and developing novel therapeutic strategies. We include an open, interactive web tool to explore the data as a molecular resource that can explain the role of these protein signatures in breast cancer classification.
Subject(s)
Biomarkers, Tumor , Breast Neoplasms , Exosomes , Proteomics , Humans , Exosomes/metabolism , Female , Breast Neoplasms/pathology , Breast Neoplasms/metabolism , Breast Neoplasms/classification , Breast Neoplasms/genetics , Biomarkers, Tumor/metabolism , Proteomics/methods , Cell Line, Tumor , Proteome/metabolismABSTRACT
Biological networks have been widely used in many different diseases to identify potential biomarkers and design drug targets. In the present review, we describe the main computational techniques for reconstructing and analyzing different types of protein networks and summarize the previous applications of such techniques in cardiovascular diseases. Existing tools are critically compared, discussing when each method is preferred such as the use of co-expression networks for functional annotation of protein clusters and the use of directed networks for inferring regulatory associations. Finally, we are presenting examples of reconstructing protein networks of different types (regulatory, co-expression, and protein-protein interaction networks). We demonstrate the necessity to reconstruct networks separately for each cardiovascular tissue type and disease entity and provide illustrative examples of the importance of taking into consideration relevant post-translational modifications. Finally, we demonstrate and discuss how the findings of protein networks could be interpreted using single-cell RNA-sequencing data.
Subject(s)
Gene Regulatory Networks , Proteomics , Protein Interaction Maps , Proteins , Computational Biology/methodsABSTRACT
Topological analysis of protein residue networks (PRNs) is a common method that can help to understand the roles of individual residues. Here, we used protein kinase A as a study object and asked what already known functionally important residues can be detected by network analysis. Along several traditional approaches to weight edges in PRNs we used local spatial pattern (LSP) alignment that assigns high weights to edges only if CαCß vectors for the corresponding residues retain their mutual positions and orientation. Our results show that even short molecular dynamic simulations of 10 to 20 ns can give convergent values for betweenness and degree centralities calculated from the LSP-based PRNs. Using these centralities, we were able to clearly distinguish a group of residues that are highly conserved in protein kinases and play important functional and regulatory roles. In comparison, traditional methods based on cross-correlation and linear mutual information were much less efficient for this particular task. These results call for reevaluation of the current methods to generate PRNs.
Subject(s)
Cyclic AMP-Dependent Protein Kinases , Molecular Dynamics SimulationABSTRACT
Networks are vital tools for understanding and modeling interactions in complex systems in science and engineering, and direct and indirect interactions are pervasive in all types of networks. However, quantitatively disentangling direct and indirect relationships in networks remains a formidable task. Here, we present a framework, called iDIRECT (Inference of Direct and Indirect Relationships with Effective Copula-based Transitivity), for quantitatively inferring direct dependencies in association networks. Using copula-based transitivity, iDIRECT eliminates/ameliorates several challenging mathematical problems, including ill-conditioning, self-looping, and interaction strength overflow. With simulation data as benchmark examples, iDIRECT showed high prediction accuracies. Application of iDIRECT to reconstruct gene regulatory networks in Escherichia coli also revealed considerably higher prediction power than the best-performing approaches in the DREAM5 (Dialogue on Reverse Engineering Assessment and Methods project, #5) Network Inference Challenge. In addition, applying iDIRECT to highly diverse grassland soil microbial communities in response to climate warming showed that the iDIRECT-processed networks were significantly different from the original networks, with considerably fewer nodes, links, and connectivity, but higher relative modularity. Further analysis revealed that the iDIRECT-processed network was more complex under warming than the control and more robust to both random and target species removal (P < 0.001). As a general approach, iDIRECT has great advantages for network inference, and it should be widely applicable to infer direct relationships in association networks across diverse disciplines in science and engineering.