Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.

Introduction: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Methods: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors. Results: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19. Discussion: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

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

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.

Human intestinal lipid storage through sequential meals reveals faster dinner appearance is associated with hyperlipidemia.

BackgroundIt is increasingly recognized that intestinal cells can store lipids after a meal, yet the effect of this phenomenon on lipid absorption patterns in insulin resistance remains unknown.MethodsThe kinetics of meal fat appearance were measured in insulin-sensitive (IS, n = 8) and insulin-resistant (IR, n = 8) subjects after sequential, isotopically labeled lunch and dinner meals. Plasma dynamics on triacylglycerol-rich (TAG-rich) lipoproteins and plasma hormones were analyzed using a nonlinear, non-steady state kinetic model.ResultsAt the onset of dinner, IS subjects showed an abrupt plasma appearance of lunch lipid consistent with the "second-meal effect," followed by slower appearance of dinner fat in plasma, resulting in reduced accumulation of dinner TAG of 48% compared with lunch. By contrast, IR subjects exhibited faster meal TAG appearance rates after both lunch and dinner. This effect of lower enterocyte storage between meals was associated with greater nocturnal and next-morning hyperlipidemia. The biochemical data and the kinetic analysis of second-meal effect dynamics are consistent with rapid secretion of stored TAG bypassing lipolysis and resynthesis. In addition, the data are consistent with a role for the diurnal pattern of plasma leptin in regulating the processing of dietary lipid.ConclusionThese data support the concept that intestinal lipid storage may be physiologically beneficial in IS subjects.Trial registrationClinicalTrials.gov NCT02020343.FundingThis study was supported by a grant from the American Diabetes Association (grant 1-13-TS-12).

The IDO Metabolic Trap Hypothesis for the Etiology of ME/CFS.

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a debilitating noncommunicable disease brandishing an enormous worldwide disease burden with some evidence of inherited genetic risk. Absence of measurable changes in patients' standard blood work has necessitated ad hoc symptom-driven therapies and a dearth of mechanistic hypotheses regarding its etiology and possible cure. A new hypothesis, the indolamine-2,3-dioxygenase (IDO) metabolic trap, was developed and formulated as a mathematical model. The historical occurrence of ME/CFS outbreaks is a singular feature of the disease and implies that any predisposing genetic mutation must be common. A database search for common damaging mutations in human enzymes produces 208 hits, including IDO2 with four such mutations. Non-functional IDO2, combined with well-established substrate inhibition of IDO1 and kinetic asymmetry of the large neutral amino acid transporter, LAT1, yielded a mathematical model of tryptophan metabolism that displays both physiological and pathological steady-states. Escape from the pathological one requires an exogenous perturbation. This model also identifies a critical point in cytosolic tryptophan abundance beyond which descent into the pathological steady-state is inevitable. If, however, means can be discovered to return cytosolic tryptophan below the critical point, return to the normal physiological steady-state is assured. Testing this hypothesis for any cell type requires only labelled tryptophan, a means to measure cytosolic tryptophan and kynurenine, and the standard tools of tracer kinetics.

The Development and Enhancement of FRAP as a Key Tool for Investigating Protein Dynamics.

The saga of fluorescence recovery after photobleaching (FRAP) illustrates how disparate technical developments impact science. Starting with the classic 1976 Axelrod et al. work in Biophysical Journal, FRAP (originally fluorescence photobleaching recovery) opened the door to extraction of quantitative information from photobleaching experiments, laying the experimental and theoretical groundwork for quantifying both the mobility and the mobile fraction of a labeled population of proteins. Over the ensuing years, FRAP's reach dramatically expanded, with new developments in GFP technology and turn-key confocal microscopy, which enabled measurement of protein diffusion and binding/dissociation rates in virtually every compartment within the cell. The FRAP technique and data catalyzed an exchange of ideas between biophysicists studying membrane dynamics, cell biologists focused on intracellular dynamics, and systems biologists modeling the dynamics of cell activity. The outcome transformed the field of cellular biology, leading to a fundamental rethinking of long-held theories of cellular dynamism. Here, we review the pivotal FRAP studies that made these developments and conceptual changes possible, which gave rise to current models of complex cell dynamics.

Differential equation methods for simulation of GFP kinetics in non-steady state experiments.

Genetically encoded fluorescent proteins, combined with fluorescence microscopy, are widely used in cell biology to collect kinetic data on intracellular trafficking. Methods for extraction of quantitative information from these data are based on the mathematics of diffusion and tracer kinetics. Current methods, although useful and powerful, depend on the assumption that the cellular system being studied is in a steady state, that is, the assumption that all the molecular concentrations and fluxes are constant for the duration of the experiment. Here, we derive new tracer kinetic analytical methods for non-steady state biological systems by constructing mechanistic nonlinear differential equation models of the underlying cell biological processes and linking them to a separate set of differential equations governing the kinetics of the fluorescent tracer. Linking the two sets of equations is based on a new application of the fundamental tracer principle of indistinguishability and, unlike current methods, supports correct dependence of tracer kinetics on cellular dynamics. This approach thus provides a general mathematical framework for applications of GFP fluorescence microscopy (including photobleaching [FRAP, FLIP] and photoactivation to frequently encountered experimental protocols involving physiological or pharmacological perturbations (e.g., growth factors, neurotransmitters, acute knockouts, inhibitors, hormones, cytokines, and metabolites) that initiate mechanistically informative intracellular transients. When a new steady state is achieved, these methods automatically reduce to classical steady state tracer kinetic analysis.

Myc Regulates Chromatin Decompaction and Nuclear Architecture during B Cell Activation.

50 years ago, Vincent Allfrey and colleagues discovered that lymphocyte activation triggers massive acetylation of chromatin. However, the molecular mechanisms driving epigenetic accessibility are still unknown. We here show that stimulated lymphocytes decondense chromatin by three differentially regulated steps. First, chromatin is repositioned away from the nuclear periphery in response to global acetylation. Second, histone nanodomain clusters decompact into mononucleosome fibers through a mechanism that requires Myc and continual energy input. Single-molecule imaging shows that this step lowers transcription factor residence time and non-specific collisions during sampling for DNA targets. Third, chromatin interactions shift from long range to predominantly short range, and CTCF-mediated loops and contact domains double in numbers. This architectural change facilitates cognate promoter-enhancer contacts and also requires Myc and continual ATP production. Our results thus define the nature and transcriptional impact of chromatin decondensation and reveal an unexpected role for Myc in the establishment of nuclear topology in mammalian cells.

Evaluation of CETP activity in vivo under non-steady-state conditions: influence of anacetrapib on HDL-TG flux.

Studies in lipoprotein kinetics almost exclusively rely on steady-state approaches to modeling. Herein, we have used a non-steady-state experimental design to examine the role of cholesteryl ester transfer protein (CETP) in mediating HDL-TG flux in vivo in rhesus macaques, and therefore, we developed an alternative strategy to model the data. Two isotopomers ([(2)H11] and [(13)C18]) of oleic acid were administered (orally and intravenously, respectively) to serve as precursors for labeling TGs in apoB-containing lipoproteins. The flux of a specific TG (52:2) from these donor lipoproteins to HDL was used as the measure of CETP activity; calculations are also presented to estimate total HDL-TG flux. Based on our data, we estimate that the peak total postprandial TG flux to HDL via CETP is ∼ 13 mg · h(-1) · kg(-1) and show that this transfer was inhibited by 97% following anacetrapib treatment. Collectively, these data demonstrate that HDL TG flux can be used as a measure of CETP activity in vivo. The fact that the donor lipoproteins can be labeled in situ using well-established stable isotope tracer techniques suggests ways to measure this activity for native lipoproteins in free-living subjects under any physiological conditions.

Mechanistic modeling confronts the complexity of molecular cell biology.

Mechanistic modeling has the potential to transform how cell biologists contend with the inescapable complexity of modern biology. I am a physiologist-electrical engineer-systems biologist who has been working at the level of cell biology for the past 24 years. This perspective aims 1) to convey why we build models, 2) to enumerate the major approaches to modeling and their philosophical differences, 3) to address some recurrent concerns raised by experimentalists, and then 4) to imagine a future in which teams of experimentalists and modelers build-and subject to exhaustive experimental tests-models covering the entire spectrum from molecular cell biology to human pathophysiology. There is, in my view, no technical obstacle to this future, but it will require some plasticity in the biological research mind-set.

A genome-wide map of CTCF multivalency redefines the CTCF code.

The "CTCF code" hypothesis posits that CTCF pleiotropic functions are driven by recognition of diverse sequences through combinatorial use of its 11 zinc fingers (ZFs). This model, however, is supported by in vitro binding studies of a limited number of sequences. To study CTCF multivalency in vivo, we define ZF binding requirements at ∼50,000 genomic sites in primary lymphocytes. We find that CTCF reads sequence diversity through ZF clustering. ZFs 4-7 anchor CTCF to ∼80% of targets containing the core motif. Nonconserved flanking sequences are recognized by ZFs 1-2 and ZFs 8-11 clusters, which also stabilize CTCF broadly. Alternatively, ZFs 9-11 associate with a second phylogenetically conserved upstream motif at ∼15% of its sites. Individually, ZFs increase overall binding and chromatin residence time. Unexpectedly, we also uncovered a conserved downstream DNA motif that destabilizes CTCF occupancy. Thus, CTCF associates with a wide array of DNA modules via combinatorial clustering of its 11 ZFs.

Insulin activation of plasma nonesterified fatty acid uptake in metabolic syndrome.

OBJECTIVE: Insulin control of fatty acid metabolism has long been deemed dominated by suppression of adipose lipolysis. The goal of the present study was to test the hypothesis that this single role of insulin is insufficient to explain observed fatty acid dynamics. METHODS AND RESULTS: Fatty acid kinetics were measured during a meal tolerance test and insulin sensitivity assessed by intravenous glucose tolerance test in overweight human subjects (n=15; body mass index, 35.8 ± 7.1 kg/m(2)). Non-steady state tracer kinetic models were formulated and tested using ProcessDB software. Suppression of adipose fatty acid release, by itself, could not account for postprandial nonesterified fatty acid concentration changes, but adipose suppression combined with insulin activation of fatty acid uptake was consistent with the measured data. The observed insulin K(m) for nonesterified fatty acid uptake was inversely correlated with both insulin sensitivity of glucose uptake (intravenous glucose tolerance test insulin sensitivity; r=-0.626; P=0.01) and whole body fat oxidation after the meal (r=-0.538; P=0.05). CONCLUSIONS: These results support insulin regulation of fatty acid turnover by both release and uptake mechanisms. Activation of fatty acid uptake is consistent with the human data, has mechanistic precedent in cell culture, and highlights a new potential target for therapies aimed at improving the control of fatty acid metabolism in insulin-resistant disease states.

Why and how to expand the role of systems biology in pharmaceutical research and development.

Seen from the perspective of funding organizations, investors, and the general public, the productivity of our world-wide biomedical research enterprise is declining despite increased investment. This opinion piece suggests a cause and a solution. The cause is the enormous complexity of human biology and pathophysiology. The unsolved human diseases involve so many interacting variables that single research laboratories headed by skilled principal investigators doing innovative experimental work cannot be expected to assemble the reductionist pieces into an integrated working model. Systems biology offers a solution, but it will require teamwork. Co-equal teams of experimental and computational biologists can construct multiscale differential equation models and test them against experimental data. A successful model provides actionable evidence-based guidance to the entire research and development team. These integrative biology teams may, for historical and cultural reasons, be unsustainable in academia, but they seem naturally suited to modern pharmaceutical research and development. One way to organize such teams and their workflow is described in detail.

Lipids and cholesterol as regulators of traffic in the endomembrane system.

The endomembrane system of eukaryotic cells uses membrane-enclosed carriers to move diverse macromolecules among different membrane-bound compartments, a requirement for cells to secrete and take up molecules from their environment. Two recycling pathways-biosynthetic and endocytic, each with specific lipid components-make up this system, with the Golgi apparatus mediating transport between the two. Here, we integrate lipid-based mechanisms into the description of this system. A partitioning model of the Golgi apparatus is discussed as a working hypothesis to explain how membrane lipids and proteins that are segregated based on lateral lipid partitioning support the unique composition of the biosynthetic and endocytic recycling pathways in the face of constant trafficking of molecular constituents. We further discuss how computational modeling can allow for interpretation of experimental findings and provide mechanistic insight into these important cellular pathways.

Salivary gland branching morphogenesis: a quantitative systems analysis of the Eda/Edar/NFkappaB paradigm.

BACKGROUND: Ectodysplasin-A appears to be a critical component of branching morphogenesis. Mutations in mouse Eda or human EDA are associated with absent or hypoplastic sweat glands, sebaceous glands, lacrimal glands, salivary glands (SMGs), mammary glands and/or nipples, and mucous glands of the bronchial, esophageal and colonic mucosa. In this study, we utilized EdaTa (Tabby) mutant mice to investigate how a marked reduction in functional Eda propagates with time through a defined genetic subcircuit and to test the proposition that canonical NFkappaB signaling is sufficient to account for the differential expression of developmentally regulated genes in the context of Eda polymorphism. RESULTS: The quantitative systems analyses do not support the stated hypothesis. For most NFkappaB-regulated genes, the observed time course of gene expression is nearly unchanged in Tabby (EdaTa) as compared to wildtype mice, as is NFkappaB itself. Importantly, a subset of genes is dramatically differentially expressed in Tabby (Edar, Fgf8, Shh, Egf, Tgfa, Egfr), strongly suggesting the existence of an alternative Eda-mediated transcriptional pathway pivotal for SMG ontogeny. Experimental and in silico investigations have identified C/EBPalpha as a promising candidate. CONCLUSION: In Tabby SMGs, upregulation of the Egf/Tgfalpha/Egfr pathway appears to mitigate the potentially severe abnormal phenotype predicted by the downregulation of Fgf8 and Shh. Others have suggested that the buffering of the phenotypic outcome that is coincident with variant Eda signaling could be a common mechanism that permits viable and diverse phenotypes, normal and abnormal. Our results support this proposition. Further, if branching epithelia use variations of a canonical developmental program, our results are likely applicable to understanding the phenotypes of other branching organs affected by Eda (EDA) mutation.

Transport through the Golgi apparatus by rapid partitioning within a two-phase membrane system.

The prevailing view of intra-Golgi transport is cisternal progression, which has a key prediction--that newly arrived cargo exhibits a lag or transit time before exiting the Golgi. Instead, we find that cargo molecules exit at an exponential rate proportional to their total Golgi abundance with no lag. Incoming cargo molecules rapidly mix with those already in the system and exit from partitioned domains with no cargo privileged for export based on its time of entry into the system. Given these results, we constructed a new model of intra-Golgi transport that involves rapid partitioning of enzymes and transmembrane cargo between two lipid phases combined with relatively rapid exchange among cisternae. Simulation and experimental testing of this rapid partitioning model reproduced all the key characteristics of the Golgi apparatus, including polarized lipid and protein gradients, exponential cargo export kinetics, and cargo waves.

In vivo dynamics of RNA polymerase II transcription.

We imaged transcription in living cells using a locus-specific reporter system, which allowed precise, single-cell kinetic measurements of promoter binding, initiation and elongation. Photobleaching of fluorescent RNA polymerase II revealed several kinetically distinct populations of the enzyme interacting with a specific gene. Photobleaching and photoactivation of fluorescent MS2 proteins used to label nascent messenger RNAs provided sensitive elongation measurements. A mechanistic kinetic model that fits our data was validated using specific inhibitors. Polymerases elongated at 4.3 kilobases min(-1), much faster than previously documented, and entered a paused state for unexpectedly long times. Transcription onset was inefficient, with only 1% of polymerase-gene interactions leading to completion of an mRNA. Our systems approach, quantifying both polymerase and mRNA kinetics on a defined DNA template in vivo with high temporal resolution, opens new avenues for studying regulation of transcriptional processes in vivo.

The ATM repair pathway inhibits RNA polymerase I transcription in response to chromosome breaks.

DNA lesions interfere with DNA and RNA polymerase activity. Cyclobutane pyrimidine dimers and photoproducts generated by ultraviolet irradiation cause stalling of RNA polymerase II, activation of transcription-coupled repair enzymes, and inhibition of RNA synthesis. During the S phase of the cell cycle, collision of replication forks with damaged DNA blocks ongoing DNA replication while also triggering a biochemical signal that suppresses the firing of distant origins of replication. Whether the transcription machinery is affected by the presence of DNA double-strand breaks remains a long-standing question. Here we monitor RNA polymerase I (Pol I) activity in mouse cells exposed to genotoxic stress and show that induction of DNA breaks leads to a transient repression in Pol I transcription. Surprisingly, we find Pol I inhibition is not itself the direct result of DNA damage but is mediated by ATM kinase activity and the repair factor proteins NBS1 (also known as NLRP2) and MDC1. Using live-cell imaging, laser micro-irradiation, and photobleaching technology we demonstrate that DNA lesions interfere with Pol I initiation complex assembly and lead to a premature displacement of elongating holoenzymes from ribosomal DNA. Our data reveal a novel ATM/NBS1/MDC1-dependent pathway that shuts down ribosomal gene transcription in response to chromosome breaks.

ArfGAP1 dynamics and its role in COPI coat assembly on Golgi membranes of living cells.

Secretory protein trafficking relies on the COPI coat, which by assembling into a lattice on Golgi membranes concentrates cargo at specific sites and deforms the membranes at these sites into coated buds and carriers. The GTPase-activating protein (GAP) responsible for catalyzing Arf1 GTP hydrolysis is an important part of this system, but the mechanism whereby ArfGAP is recruited to the coat, its stability within the coat, and its role in maintenance of the coat are unclear. Here, we use FRAP to monitor the membrane turnover of GFP-tagged versions of ArfGAP1, Arf1, and coatomer in living cells. ArfGAP1 underwent fast cytosol/Golgi exchange with approximately 40% of the exchange dependent on engagement of ArfGAP1 with coatomer and Arf1, and affected by secretory cargo load. Permanent activation of Arf1 resulted in ArfGAP1 being trapped on the Golgi in a coatomer-dependent manner. These data suggest that ArfGAP1, coatomer and Arf1 play interdependent roles in the assembly-disassembly cycle of the COPI coat in vivo.