Deep models of integrated multiscale molecular data decipher the endothelial cell response to ionizing radiation.

The vascular endothelium is a hot spot in the response to radiation therapy for both tumors and normal tissues. To improve patient outcomes, interpretable systemic hypotheses are needed to help radiobiologists and radiation oncologists propose endothelial targets that could protect normal tissues from the adverse effects of radiation therapy and/or enhance its antitumor potential. To this end, we captured the kinetics of multi-omics layers-i.e. miRNome, targeted transcriptome, proteome, and metabolome-in irradiated primary human endothelial cells cultured in vitro. We then designed a strategy of deep learning as in convolutional graph networks that facilitates unsupervised high-level feature extraction of important omics data to learn how ionizing radiation-induced endothelial dysfunction may evolve over time. Last, we present experimental data showing that some of the features identified using our approach are involved in the alteration of angiogenesis by ionizing radiation.

Deciphering the Dynamic Molecular Program of Radiation-Induced Endothelial Senescence.

PURPOSE: Radiation-induced cellular senescence is a double-edged sword, acting as both a tumor suppression process limiting tumor proliferation, and a crucial process contributing to normal tissue injury. Endothelial cells play a role in normal tissue injury after radiation therapy. Recently, a study observed an accumulation of senescent endothelial cells (ECs) around radiation-induced lung focal lesions following stereotactic radiation injury in mice. However, the effect of radiation on EC senescence remains unclear because it depends on dose and fractionation, and because the senescent phenotype is heterogeneous and dynamic. METHODS AND MATERIALS: Using a systems biology approach in vitro, we deciphered the dynamic senescence-associated transcriptional program induced by irradiation. RESULTS: Flow cytometry and single-cell RNA sequencing experiments revealed the heterogeneous senescent status of irradiated ECs and allowed to deciphered the molecular program involved in this status. We identified the Interleukin-1 signaling pathway as a key player in the radiation-induced premature senescence of ECs, as well as the endothelial-to-mesenchymal transition process, which shares strong hallmarks of senescence. CONCLUSIONS: Our work provides crucial information on the dynamics of the radiation-induced premature senescence process, the effect of the radiation dose, as well as the molecular program involved in the heterogeneous senescent status of ECs.

Singular manifolds of proteomic drivers to model the evolution of inflammatory bowel disease status.

The conditions used to describe the presence of an immune disease are often represented by interaction graphs. These informative, but intricate structures are susceptible to perturbations at different levels. The mode in which that perturbation occurs is still of utmost importance in areas such as cell reprogramming and therapeutics models. In this sense, module identification can be useful to well characterise the global graph architecture. To help us with this identification, we perform topological overlap-related measures. Thanks to these measures, the location of highly disease-specific module regulators is possible. Such regulators can perturb other nodes, potentially causing the entire system to change behaviour or collapse. We provide a geometric framework explaining such situations in the context of inflammatory bowel diseases (IBD). IBD are severe chronic disorders of the gastrointestinal tract whose incidence is dramatically increasing worldwide. Our approach models different IBD status as Riemannian manifolds defined by the graph Laplacian of two high throughput proteome screenings. It also identifies module regulators as singularities within the manifolds (the so-called singular manifolds). Furthermore, it reinterprets the characteristic nonlinear dynamics of IBD as compensatory responses to perturbations on those singularities. Then, particular reconfigurations of the immune system could make the disease status move towards an innocuous target state.

Mathematical deconvolution uncovers the genetic regulatory signal of cancer cellular heterogeneity on resistance to paclitaxel.

Drug resistance remains a major problem in combating malignancies, resulting critical the resistance to paclitaxel used in the treatment of many different cancers. Elucidating the cellular heterogeneity composition of tumours may be relevant to designing more effective treatment strategies on drug resistance. In particular, such heterogeneity correlates with the measurement of gene expression below the population level. However, experimental assays capturing differential response are limited and cannot discern the variation in gene expression specific to different cellular types in tumour populations. These limitations led us to consider a mathematical modelling approach, in which the gene expression of cellular subpopulations is recovered by deconvolution. Mathematically, the deconvolution is a multi-linear regression-based problem. We combined herein data on cellular subpopulation frequency composition with gene expression values from 16 tumour lines (8 resistant and 8 sensitive to paclitaxel treatment) to find genes that are differentially expressed between paclitaxel resistant and paclitaxel sensitive tumour lines in different cellular subpopulations. The results indicate that many genes differentially expressed between paclitaxel resistant and sensitive cancer lines are only detected when considering their heterogeneous cellular composition. Overall, our methodology is thought to keep in mind phenotypic heterogeneity improving our resolution in the identification of biomarkers on resistance to chemo-therapeutic agents.

Integrative Network-based Analysis of Colonic Detoxification Gene Expression in Ulcerative Colitis According to Smoking Status.

BACKGROUNDS AND AIMS: The effect of cigarette smoking [CS] is ambivalent since smoking improves ulcerative colitis [UC] while it worsens Crohn's disease [CD]. Although this clinical relationship between inflammatory bowel disease [IBD] and tobacco is well established, only a few experimental works have investigated the effect of smoking on the colonic barrier homeostasis focusing on xenobiotic detoxification genes. METHODS: A comprehensive and integrated comparative analysis of the global xenobiotic detoxification capacity of the normal colonic mucosa of healthy smokers [n = 8] and non-smokers [n = 9] versus the non-affected colonic mucosa of UC patients [n = 19] was performed by quantitative real-time polymerase chain reaction [qRT PCR]. The detoxification gene expression profile was analysed in CD patients [n = 18], in smoking UC patients [n = 5], and in biopsies from non-smoking UC patients cultured or not with cigarette smoke extract [n = 8]. RESULTS: Of the 244 detoxification genes investigated, 65 were dysregulated in UC patients in comparison with healthy controls or CD patients. The expression of ≥ 45/65 genes was inversed by CS in biopsies of smoking UC patients in remission and in colonic explants of UC patients exposed to cigarette smoke extract. We devised a network-based data analysis approach for differentially assessing changes in genetic interactions, allowing identification of unexpected regulatory detoxification genes that may play a major role in the beneficial effect of smoking on UC. CONCLUSIONS: Non-inflamed colonic mucosa in UC is characterised by a specifically altered detoxification gene network, which is partially restored by tobacco. These mucosal signatures could be useful for developing new therapeutic strategies and biomarkers of drug response in UC.

Exploring the interactions of the RAS family in the human protein network and their potential implications in RAS-directed therapies.

RAS proteins are the founding members of the RAS superfamily of GTPases. They are involved in key signaling pathways regulating essential cellular functions such as cell growth and differentiation. As a result, their deregulation by inactivating mutations often results in aberrant cell proliferation and cancer. With the exception of the relatively well-known KRAS, HRAS and NRAS proteins, little is known about how the interactions of the other RAS human paralogs affect cancer evolution and response to treatment. In this study we performed a comprehensive analysis of the relationship between the phylogeny of RAS proteins and their location in the protein interaction network. This analysis was integrated with the structural analysis of conserved positions in available 3D structures of RAS complexes. Our results show that many RAS proteins with divergent sequences are found close together in the human interactome. We found specific conserved amino acid positions in this group that map to the binding sites of RAS with many of their signaling effectors, suggesting that these pairs could share interacting partners. These results underscore the potential relevance of cross-talking in the RAS signaling network, which should be taken into account when considering the inhibitory activity of drugs targeting specific RAS oncoproteins. This study broadens our understanding of the human RAS signaling network and stresses the importance of considering its potential cross-talk in future therapies.

Protein-protein interfaces from cytochrome c oxidase I evolve faster than nonbinding surfaces, yet negative selection is the driving force.

Respiratory complexes are encoded by two genomes (mitochondrial DNA [mtDNA] and nuclear DNA [nDNA]). Although the importance of intergenomic coadaptation is acknowledged, the forces and constraints shaping such coevolution are largely unknown. Previous works using cytochrome c oxidase (COX) as a model enzyme have led to the so-called "optimizing interaction" hypothesis. According to this view, mtDNA-encoded residues close to nDNA-encoded residues evolve faster than the rest of positions, favoring the optimization of protein-protein interfaces. Herein, using evolutionary data in combination with structural information of COX, we show that failing to discern the effects of interaction from other structural and functional effects can lead to deceptive conclusions such as the "optimizing hypothesis." Once spurious factors have been accounted for, data analysis shows that mtDNA-encoded residues engaged in contacts are, in general, more constrained than their noncontact counterparts. Nevertheless, noncontact residues from the surface of COX I subunit are a remarkable exception, being subjected to an exceptionally high purifying selection that may be related to the maintenance of a suitable heme environment. We also report that mtDNA-encoded residues involved in contacts with other mtDNA-encoded subunits are more constrained than mtDNA-encoded residues interacting with nDNA-encoded polypeptides. This differential behavior cannot be explained on the basis of predicted thermodynamic stability, as interactions between mtDNA-encoded subunits contribute more weakly to the complex stability than those interactions between subunits encoded by different genomes. Therefore, the higher conservation observed among mtDNA-encoded residues involved in intragenome interactions is likely due to factors other than structural stability.