QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration.

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration.

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

Structure-based prediction of Ras-effector binding affinities and design of "branchegetic" interface mutations.

Ras is a central cellular hub protein controlling multiple cell fates. How Ras interacts with a variety of potential effector proteins is relatively unexplored, with only some key effectors characterized in great detail. Here, we have used homology modeling based on X-ray and AlphaFold2 templates to build structural models for 54 Ras-effector complexes. These models were used to estimate binding affinities using a supervised learning regressor. Furthermore, we systematically introduced Ras "branch-pruning" (or branchegetic) mutations to identify 200 interface mutations that affect the binding energy with at least one of the model structures. The impacts of these branchegetic mutants were integrated into a mathematical model to assess the potential for rewiring interactions at the Ras hub on a systems level. These findings have provided a quantitative understanding of Ras-effector interfaces and their impact on systems properties of a key cellular hub.

Analysis of context-specific KRAS-effector (sub)complexes in Caco-2 cells.

Ras is a key switch controlling cell behavior. In the GTP-bound form, Ras interacts with numerous effectors in a mutually exclusive manner, where individual Ras-effectors are likely part of larger cellular (sub)complexes. The molecular details of these (sub)complexes and their alteration in specific contexts are not understood. Focusing on KRAS, we performed affinity purification (AP)-mass spectrometry (MS) experiments of exogenously expressed FLAG-KRAS WT and three oncogenic mutants ("genetic contexts") in the human Caco-2 cell line, each exposed to 11 different culture media ("culture contexts") that mimic conditions relevant in the colon and colorectal cancer. We identified four effectors present in complex with KRAS in all genetic and growth contexts ("context-general effectors"). Seven effectors are found in KRAS complexes in only some contexts ("context-specific effectors"). Analyzing all interactors in complex with KRAS per condition, we find that the culture contexts had a larger impact on interaction rewiring than genetic contexts. We investigated how changes in the interactome impact functional outcomes and created a Shiny app for interactive visualization. We validated some of the functional differences in metabolism and proliferation. Finally, we used networks to evaluate how KRAS-effectors are involved in the modulation of functions by random walk analyses of effector-mediated (sub)complexes. Altogether, our work shows the impact of environmental contexts on network rewiring, which provides insights into tissue-specific signaling mechanisms. This may also explain why KRAS oncogenic mutants may be causing cancer only in specific tissues despite KRAS being expressed in most cells and tissues.

Classification of GTP-dependent K-Ras4B active and inactive conformational states.

Classifying reliably active and inactive molecular conformations of wildtype (WT) and mutated oncogenic proteins is a key, ongoing challenge in molecular cancer studies. Here, we probe the GTP-bound K-Ras4B conformational dynamics using long-time atomistic molecular dynamics (MD) simulations. We extract and analyze the detailed underlying free energy landscape of WT K-Ras4B. We use two key reaction coordinates, labeled d1 and d2 (i.e., distances coordinating the Pß atom of the GTP ligand with two key residues, T35 and G60), shown to correlate closely with activities of WT and mutated K-Ras4B. However, our new K-Ras4B conformational kinetics study reveals a more complex network of equilibrium Markovian states. We show that a new reaction coordinate is required to account for the orientation of acidic K-Ras4B sidechains such as D38 with respect to the interface with binding effector RAF1 and rationalize the activation/inactivation propensities and the corresponding molecular binding mechanisms. We use this understanding to unveil how a relatively conservative mutation (i.e., D33E, in the switch I region) can lead to significantly different activation propensities compared with WT K-Ras4B. Our study sheds new light on the ability of residues near the K-Ras4B-RAF1 interface to modulate the network of salt bridges at the binding interface with the RAF1 downstream effector and, thus, to influence the underlying GTP-dependent activation/inactivation mechanism. Altogether, our hybrid MD-docking modeling approach enables the development of new in silico methods for quantitative assessment of activation propensity changes (e.g., due to mutations or local binding environment). It also unveils the underlying molecular mechanisms and facilitates the rational design of new cancer drugs.

Whole-cell energy modeling reveals quantitative changes of predicted energy flows in RAS mutant cancer cell lines.

Cellular utilization of available energy flows to drive a multitude of forms of cellular "work" is a major biological constraint. Cells steer metabolism to address changing phenotypic states but little is known as to how bioenergetics couples to the richness of processes in a cell as a whole. Here, we outline a whole-cell energy framework that is informed by proteomic analysis and an energetics-based gene ontology. We separate analysis of metabolic supply and the capacity to generate high-energy phosphates from a representation of demand that is built on the relative abundance of ATPases and GTPases that deliver cellular work. We employed mouse embryonic fibroblast cell lines that express wild-type KRAS or oncogenic mutations and with distinct phenotypes. We observe shifts between energy-requiring processes. Calibrating against Seahorse analysis, we have created a whole-cell energy budget with apparent predictive power, for instance in relation to protein synthesis.

Vitronectin and Its Interaction with PAI-1 Suggests a Functional Link to Vascular Changes in AMD Pathobiology.

The pathogenesis of age-related macular degeneration (AMD), a frequent disorder of the central retina, is incompletely understood. Genome-wide association studies (GWAS) suggest a strong contribution of genomic variation in AMD susceptibility. Nevertheless, little is known about biological mechanisms of the disease. We reported previously that the AMD-associated polymorphism rs704C > T in the vitronectin (VTN) gene influences protein expression and functional aspects of encoded vitronectin, a human blood and extracellular matrix (ECM) protein. Here, we refined the association of rs704 with AMD in 16,144 cases and 17,832 controls and noted that rs704 is carried exclusively by the neovascular AMD subtype. Interaction studies demonstrate that rs704 affects the ability of vitronectin to bind the angiogenic regulator plasminogen activator inhibitor 1 (PAI-1) but has no influence on stabilizing its active state. Western blot analysis and confocal imaging reveal a strong enrichment of PAI-1 in the ECM of cultured endothelial cells and RPE cell line ARPE-19 exposed to vitronectin. Large-scale gene expression of VTN and PAI-1 showed positive correlations and a statistically significant increase in human retinal and blood tissues aged 60 years and older. Our results suggest a mechanism by which the AMD-associated rs704 variant in combination with ageing may contribute to the vascular complications in AMD.

Reconstruction and analysis of a large-scale binary Ras-effector signaling network.

BACKGROUND: Ras is a key cellular signaling hub that controls numerous cell fates via multiple downstream effector pathways. While pathways downstream of effectors such as Raf, PI3K and RalGDS are extensively described in the literature, how other effectors signal downstream of Ras is often still enigmatic. METHODS: A comprehensive and unbiased Ras-effector network was reconstructed downstream of 43 effector proteins (converging onto 12 effector classes) using public pathway and protein-protein interaction (PPI) databases. The output is an oriented graph of pairwise interactions defining a 3-layer signaling network downstream of Ras. The 2290 proteins comprising the network were studied for their implication in signaling crosstalk and feedbacks, their subcellular localizations, and their cellular functions. RESULTS: The final Ras-effector network consists of 2290 proteins that are connected via 19,080 binary PPIs, increasingly distributed across the downstream layers, with 441 PPIs in layer 1, 1660 in layer 2, and 16,979 in layer 3. We identified a high level of crosstalk among proteins of the 12 effector classes. A class-specific Ras sub-network was generated in CellDesigner (.xml file) and a functional enrichment analysis thereof shows that 58% of the processes have previously been associated to a respective effector pathway, with the remaining providing insights into novel and unexplored functions of specific effector pathways. CONCLUSIONS: Our large-scale and cell general Ras-effector network is a crucial steppingstone towards defining the network boundaries. It constitutes a 'reference interactome' and can be contextualized for specific conditions, e.g. different cell types or biopsy material obtained from cancer patients. Further, it can serve as a basis for elucidating systems properties, such as input-output relationships, crosstalk, and pathway redundancy. Video Abstract.

Seed sequence polymorphism rs2168518 and allele-specific target gene regulation of hsa-miR-4513.

MicroRNAs (miRNAs) are small post-transcriptional regulators that offer promising targets for treating complex diseases. To this end, hsa-miR-4513 is an excellent candidate as this gene harbors within its conserved heptametrical seed sequence a frequent polymorphism (rs2168518), which has previously been associated with several complex phenotypes. So far, little is known about the biological mechanism(s) underlying these associations. In an initial step, we now aimed to identify allele-specific target genes of hsa-miR-4513. We performed RNA sequencing in a miRNA overexpression model in human umbilical vein endothelial cells transfected with separated hsa-miR-4513 alleles at rs2168518, namely hsa-miR-4513-G and hsa-miR-4513-A. Genes specifically regulated by the rs2168518 alleles were independently verified by quantitative reverse transcription PCR (qRT-PCR), western blot analysis and allele-specific miRNA binding via a luciferase reporter assay. By a text-based search publicly available databases such as Online Mendelian Inheritance in Man and Mouse Genome Informatics were utilized to link target genes of hsa-miR-4513 to previously described phenotypes. Overall, we identified 23 allele-specific hsa-miR-4513 target genes and replicated 19 of those independently via qRT-PCR. Western blot analysis and luciferase reporter assays conducted for an exemplary subsample further confirmed the allele-specific regulation of these genes by hsa-miR-4513. Remarkably, multiple allele-specific target genes identified are linked via text retrieval to several phenotypes previously reported to be associated with hsa-miR-4513. These genes offer promising candidates for ongoing research on the functional pathobiological impact of hsa-miR-4513 and its seed polymorphism rs2168518. This could give rise to therapeutic applications targeting this miRNA.

Network pharmacology: curing causal mechanisms instead of treating symptoms.

For complex diseases, most drugs are highly ineffective, and the success rate of drug discovery is in constant decline. While low quality, reproducibility issues, and translational irrelevance of most basic and preclinical research have contributed to this, the current organ-centricity of medicine and the 'one disease-one target-one drug' dogma obstruct innovation in the most profound manner. Systems and network medicine and their therapeutic arm, network pharmacology, revolutionize how we define, diagnose, treat, and, ideally, cure diseases. Descriptive disease phenotypes are replaced by endotypes defined by causal, multitarget signaling modules that also explain respective comorbidities. Precise and effective therapeutic intervention is achieved by synergistic multicompound network pharmacology and drug repurposing, obviating the need for drug discovery and speeding up clinical translation.

HOMELETTE: a unified interface to homology modelling software.

SUMMARY: Homology modelling, the technique of generating models of 3D protein structures based on experimental structures from related proteins, has become increasingly popular over the years. An abundance of different tools for model generation and model evaluation is available from various research groups. We present HOMELETTE, an interface which implements a unified programmatic access to these tools. This allows for the assemble of custom pipelines from pre- or self-implemented building blocks. AVAILABILITY AND IMPLEMENTATION: HOMELETTE is implemented in Python, compatible with version 3.6 and newer. It is distributed under the MIT license. Documentation and tutorials are available at Read the Docs (https://homelette.readthedocs.io/). The latest version of HOMELETTE is available on PyPI (https://pypi.org/project/homelette/) and GitHub (https://github.com/PhilippJunk/homelette). A full installation of the latest version of HOMELETTE with all dependencies is also available as a Docker container (https://hub.docker.com/r/philippjunk/homelette_template). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Disease-Gene Networks of Skin Pigmentation Disorders and Reconstruction of Protein-Protein Interaction Networks.

Melanin, a light and free radical absorbing pigment, is produced in melanocyte cells that are found in skin, but also in hair follicles, eyes, the inner ear, heart, brain and other organs. Melanin synthesis is the result of a complex network of signaling and metabolic reactions. It therefore comes as no surprise that mutations in many of the genes involved are associated with various types of pigmentation diseases and phenotypes ('pigmentation genes'). Here, we used bioinformatics tools to first reconstruct gene-disease/phenotype associations for all pigmentation genes. Next, we reconstructed protein-protein interaction (PPI) networks centered around pigmentation gene products ('pigmentation proteins') and supplemented the PPI networks with protein expression information obtained by mass spectrometry in a panel of melanoma cell lines (both pigment producing and non-pigment producing cells). The analysis provides a systems network representation of all genes/ proteins centered around pigmentation and melanin biosynthesis pathways ('pigmentation network map'). Our work will enable the pigmentation research community to experimentally test new hypothesis arising from the pigmentation network map and to identify new targets for drug discovery.

The metabolic environment of the developing embryo: A multidisciplinary approach on oilseed rapeseed.

Brassicaceae seeds consist of three genetically distinct structures: the embryo, endosperm and seed coat, all of which are involved in assimilate allocation during seed development. The complexity of their metabolic interrelations remains unresolved to date. In the present study, we apply state-of-the-art imaging and analytical approaches to assess the metabolic environment of the Brassica napus embryo. Nuclear magnetic resonance imaging (MRI) provided volumetric data on the living embryo and endosperm, revealing how the endosperm envelops the embryo, determining endosperm's priority in assimilate uptake from the seed coat during early development. MRI analysis showed higher levels of sugars in the peripheral endosperm facing the seed coat, but a lower sugar content within the central vacuole and the region surrounding the embryo. Feeding intact siliques with 13C-labeled sucrose allowed tracing of the post-phloem route of sucrose transfer within the seed at the heart stage of embryogenesis, by means of mass spectrometry imaging. Quantification of over 70 organic and inorganic compounds in the endosperm revealed shifts in their abundance over different stages of development, while sugars and potassium were the main determinants of osmolality throughout these stages. Our multidisciplinary approach allows access to the hidden aspects of endosperm metabolism, a task which remains unattainable for the small-seeded model plant Arabidopsis thaliana.

Assigning Co-Regulated Human Genes and Regulatory Gene Clusters.

Elucidating the role of genetic variation in the regulation of gene expression is key to understanding the pathobiology of complex diseases which, in consequence, is crucial in devising targeted treatment options. Expression quantitative trait locus (eQTL) analysis correlates a genetic variant with the strength of gene expression, thus defining thousands of regulated genes in a multitude of human cell types and tissues. Some eQTL may not act independently of each other but instead may be regulated in a coordinated fashion by seemingly independent genetic variants. To address this issue, we combined the approaches of eQTL analysis and colocalization studies. Gene expression was determined in datasets comprising 49 tissues from the Genotype-Tissue Expression (GTEx) project. From about 33,000 regulated genes, over 14,000 were found to be co-regulated in pairs and were assembled across all tissues to almost 15,000 unique clusters containing up to nine regulated genes affected by the same eQTL signal. The distance of co-regulated eGenes was, on average, 112 kilobase pairs. Of 713 genes known to express clinical symptoms upon haploinsufficiency, 231 (32.4%) are part of at least one of the identified clusters. This calls for caution should treatment approaches aim at an upregulation of a haploinsufficient gene. In conclusion, we present an unbiased approach to identifying co-regulated genes in and across multiple tissues. Knowledge of such common effects is crucial to appreciate implications on biological pathways involved, specifically when a treatment option targets a co-regulated disease gene.

Cell Adhesion Molecules in Normal Skin and Melanoma.

Cell adhesion molecules (CAMs) of the cadherin, integrin, immunoglobulin, and selectin protein families are indispensable for the formation and maintenance of multicellular tissues, especially epithelia. In the epidermis, they are involved in cell-cell contacts and in cellular interactions with the extracellular matrix (ECM), thereby contributing to the structural integrity and barrier formation of the skin. Bulk and single cell RNA sequencing data show that >170 CAMs are expressed in the healthy human skin, with high expression levels in melanocytes, keratinocytes, endothelial, and smooth muscle cells. Alterations in expression levels of CAMs are involved in melanoma propagation, interaction with the microenvironment, and metastasis. Recent mechanistic analyses together with protein and gene expression data provide a better picture of the role of CAMs in the context of skin physiology and melanoma. Here, we review progress in the field and discuss molecular mechanisms in light of gene expression profiles, including recent single cell RNA expression information. We highlight key adhesion molecules in melanoma, which can guide the identification of pathways and strategies for novel anti-melanoma therapies.

Engineering of Biological Pathways: Complex Formation and Signal Transduction.

The rational in silico design of interface mutations within protein complexes is a synthetic biology tool that enables-when introduced into biological systems-the artificial rewiring of biological pathways. Here we describe the three-dimensional structure-based design of "rewiring" mutations using the FoldX force field. Specifically, we provide the protocol for the design and selection of interface mutations in three Ras-effector complex structures (PDB entries 3KUD, 4K81, and 6AMB). Ras mutations that impair binding to some but not all interacting partners are selected.

Epistatic interactions of genetic loci associated with age-related macular degeneration.

The currently largest genome-wide association study (GWAS) for age-related macular degeneration (AMD) defines disease association with genome-wide significance for 52 independent common and rare genetic variants across 34 chromosomal loci. Overall, these loci contain over 7200 variants and are enriched for genes with functions indicating several shared cellular processes. Still, the precise mechanisms leading to AMD pathology are largely unknown. Here, we exploit the phenomenon of epistatic interaction to identify seemingly independent AMD-associated variants that reveal joint effects on gene expression. We focus on genetic variants associated with lipid metabolism, organization of extracellular structures, and innate immunity, specifically the complement cascade. Multiple combinations of independent variants were used to generate genetic risk scores allowing gene expression in liver to be compared between low and high-risk AMD. We identified genetic variant combinations correlating significantly with expression of 26 genes, of which 19 have not been associated with AMD before. This study defines novel targets and allows prioritizing further functional work into AMD pathobiology.

Signaling pathways in intestinal homeostasis and colorectal cancer: KRAS at centre stage.

The intestinal epithelium acts as a physical barrier that separates the intestinal microbiota from the host and is critical for preserving intestinal homeostasis. The barrier is formed by tightly linked intestinal epithelial cells (IECs) (i.e. enterocytes, goblet cells, neuroendocrine cells, tuft cells, Paneth cells, and M cells), which constantly self-renew and shed. IECs also communicate with microbiota, coordinate innate and adaptive effector cell functions. In this review, we summarize the signaling pathways contributing to intestinal cell fates and homeostasis functions. We focus especially on intestinal stem cell proliferation, cell junction formation, remodelling, hypoxia, the impact of intestinal microbiota, the immune system, inflammation, and metabolism. Recognizing the critical role of KRAS mutants in colorectal cancer, we highlight the connections of KRAS signaling pathways in coordinating these functions. Furthermore, we review the impact of KRAS colorectal cancer mutants on pathway rewiring associated with disruption and dysfunction of the normal intestinal homeostasis. Given that KRAS is still considered undruggable and the development of treatments that directly target KRAS are unlikely, we discuss the suitability of targeting pathways downstream of KRAS as well as alterations of cell extrinsic/microenvironmental factors as possible targets for modulating signaling pathways in colorectal cancer. Video Abstract.

The Ins and Outs of RAS Effector Complexes.

RAS oncogenes are among the most commonly mutated proteins in human cancers. They regulate a wide range of effector pathways that control cell proliferation, survival, differentiation, migration and metabolic status. Including aberrations in these pathways, RAS-dependent signaling is altered in more than half of human cancers. Targeting mutant RAS proteins and their downstream oncogenic signaling pathways has been elusive. However, recent results comprising detailed molecular studies, large scale omics studies and computational modeling have painted a new and more comprehensive portrait of RAS signaling that helps us to understand the intricacies of RAS, how its physiological and pathophysiological functions are regulated, and how we can target them. Here, we review these efforts particularly trying to relate the detailed mechanistic studies with global functional studies. We highlight the importance of computational modeling and data integration to derive an actionable understanding of RAS signaling that will allow us to design new mechanism-based therapies for RAS mutated cancers.