A novel logical model of COVID-19 intracellular infection to support therapies development.

In this paper, a logical-based mathematical model of the cellular pathways involved in the COVID-19 infection has been developed to study various drug treatments (single or in combination), in different illness scenarios, providing insights into their mechanisms of action. Drug simulations suggest that the effects of single drugs are limited, or depending on the scenario counterproductive, whereas better results appear combining different treatments. Specifically, the combination of the anti-inflammatory Baricitinib and the anti-viral Remdesivir showed significant benefits while a stronger efficacy emerged from the triple combination of Baricitinib, Remdesivir, and the corticosteroid Dexamethasone. Together with a sensitivity analysis, we performed an analysis of the mechanisms of the drugs to reveal their impact on molecular pathways.

Computational Analysis of Cytokine Release Following Bispecific T-Cell Engager Therapy: Applications of a Logic-Based Model.

Bispecific T-cell engaging therapies harness the immune system to elicit an effective anticancer response. Modulating the immune activation avoiding potential adverse effects such as cytokine release syndrome (CRS) is a critical aspect to realizing the full potential of this therapy. The use of suitable exogenous intervention strategies to mitigate the CRS risk without compromising the antitumoral capability of bispecific antibody treatment is crucial. To this end, computational approaches can be instrumental to systematically exploring the effects of combining bispecific antibodies with CRS intervention strategies. Here, we employ a logical model to describe the action of bispecific antibodies and the complex interplay of various immune system components and use it to perform simulation experiments to improve the understanding of the factors affecting CRS. We performed a sensitivity analysis to identify the comedications that could ameliorate CRS without impairing tumor clearance. Our results agree with publicly available experimental data suggesting anti-TNF and anti-IL6 as possible co-treatments. Furthermore, we suggest anti-IFNγ as a suitable candidate for clinical studies.

Literature Mining and Mechanistic Graphical Modelling to Improve mRNA Vaccine Platforms.

RNA vaccines represent a milestone in the history of vaccinology. They provide several advantages over more traditional approaches to vaccine development, showing strong immunogenicity and an overall favorable safety profile. While preclinical testing has provided some key insights on how RNA vaccines interact with the innate immune system, their mechanism of action appears to be fragmented amid the literature, making it difficult to formulate new hypotheses to be tested in clinical settings and ultimately improve this technology platform. Here, we propose a systems biology approach, based on the combination of literature mining and mechanistic graphical modeling, to consolidate existing knowledge around mRNA vaccines mode of action and enhance the translatability of preclinical hypotheses into clinical evidence. A Natural Language Processing (NLP) pipeline for automated knowledge extraction retrieved key biological evidences that were joined into an interactive mechanistic graphical model representing the chain of immune events induced by mRNA vaccines administration. The achieved mechanistic graphical model will help the design of future experiments, foster the generation of new hypotheses and set the basis for the development of mathematical models capable of simulating and predicting the immune response to mRNA vaccines.

In Silico Logical Modelling to Uncover Cooperative Interactions in Cancer.

The multistep development of cancer involves the cooperation between multiple molecular lesions, as well as complex interactions between cancer cells and the surrounding tumour microenvironment. The search for these synergistic interactions using experimental models made tremendous contributions to our understanding of oncogenesis. Yet, these approaches remain labour-intensive and challenging. To tackle such a hurdle, an integrative, multidisciplinary effort is required. In this article, we highlight the use of logical computational models, combined with experimental validations, as an effective approach to identify cooperative mechanisms and therapeutic strategies in the context of cancer biology. In silico models overcome limitations of reductionist approaches by capturing tumour complexity and by generating powerful testable hypotheses. We review representative examples of logical models reported in the literature and their validation. We then provide further analyses of our logical model of Epithelium to Mesenchymal Transition (EMT), searching for additional cooperative interactions involving inputs from the tumour microenvironment and gain of function mutations in NOTCH.

A Novel Hybrid Logic-ODE Modeling Approach to Overcome Knowledge Gaps.

Mathematical modeling allows using different formalisms to describe, investigate, and understand biological processes. However, despite the advent of high-throughput experimental techniques, quantitative information is still a challenge when looking for data to calibrate model parameters. Furthermore, quantitative formalisms must cope with stiffness and tractability problems, more so if used to describe multicellular systems. On the other hand, qualitative models may lack the proper granularity to describe the underlying kinetic processes. We propose a hybrid modeling approach that integrates ordinary differential equations and logical formalism to describe distinct biological layers and their communication. We focused on a multicellular system as a case study by applying the hybrid formalism to the well-known Delta-Notch signaling pathway. We used a differential equation model to describe the intracellular pathways while the cell-cell interactions were defined by logic rules. The hybrid approach herein employed allows us to combine the pros of different modeling techniques by overcoming the lack of quantitative information with a qualitative description that discretizes activation and inhibition processes, thus avoiding complexity.

Toward in vivo-relevant hERG safety assessment and mitigation strategies based on relationships between non-equilibrium blocker binding, three-dimensional channel-blocker interactions, dynamic occupancy, dynamic exposure, and cellular arrhythmia.

The human ether-a-go-go-related voltage-gated cardiac ion channel (commonly known as hERG) conducts the rapid outward repolarizing potassium current in cardiomyocytes (IKr). Inadvertent blockade of this channel by drug-like molecules represents a key challenge in pharmaceutical R&D due to frequent overlap between the structure-activity relationships of hERG and many primary targets. Building on our previous work, together with recent cryo-EM structures of hERG, we set about to better understand the energetic and structural basis of promiscuous blocker-hERG binding in the context of Biodynamics theory. We propose a two-step blocker binding process consisting of: The initial capture step: diffusion of a single fully solvated blocker copy into a large cavity lined by the intra-cellular cyclic nucleotide binding homology domain (CNBHD). Occupation of this cavity is a necessary but insufficient condition for ion current disruption.The IKr disruption step: translocation of the captured blocker along the channel axis, such that: The head group, consisting of a quasi-rod-shaped moiety, projects into the open pore, accompanied by partial de-solvation of the binding interface.One tail moiety packs along a kink between the S6 helix and proximal C-linker helix adjacent to the intra-cellular entrance of the pore, likewise accompanied by mutual de-solvation of the binding interface (noting that the association barrier is comprised largely of the total head + tail group de-solvation cost).Blockers containing a highly planar moiety that projects into a putative constriction zone within the closed channel become trapped upon closing, as do blockers terminating prior to this region.A single captured blocker copy may conceivably associate and dissociate to/from the pore many times before exiting the CNBHD cavity. Lastly, we highlight possible flaws in the current hERG safety index (SI), and propose an alternate in vivo-relevant strategy factoring in: Benefit/risk.The predicted arrhythmogenic fractional hERG occupancy (based on action potential (AP) simulations of the undiseased human ventricular cardiomyocyte).Alteration of the safety threshold due to underlying disease.Risk of exposure escalation toward the predicted arrhythmic limit due to patient-to-patient pharmacokinetic (PK) variability, drug-drug interactions, overdose, and use for off-label indications in which the hERG safety parameters may differ from their on-label counterparts.

Hybrid Epithelial-Mesenchymal Phenotypes Are Controlled by Microenvironmental Factors.

Epithelial-to-mesenchymal transition (EMT) has been associated with cancer cell heterogeneity, plasticity, and metastasis. However, the extrinsic signals supervising these phenotypic transitions remain elusive. To assess how selected microenvironmental signals control cancer-associated phenotypes along the EMT continuum, we defined a logical model of the EMT cellular network that yields qualitative degrees of cell adhesions by adherens junctions and focal adhesions, two features affected during EMT. The model attractors recovered epithelial, mesenchymal, and hybrid phenotypes. Simulations showed that hybrid phenotypes may arise through independent molecular paths involving stringent extrinsic signals. Of particular interest, model predictions and their experimental validations indicated that: (i) stiffening of the extracellular matrix was a prerequisite for cells overactivating FAK_SRC to upregulate SNAIL and acquire a mesenchymal phenotype and (ii) FAK_SRC inhibition of cell-cell contacts through the receptor-type tyrosine-protein phosphatases kappa led to acquisition of a full mesenchymal, rather than a hybrid, phenotype. Altogether, these computational and experimental approaches allow assessment of critical microenvironmental signals controlling hybrid EMT phenotypes and indicate that EMT involves multiple molecular programs. SIGNIFICANCE: A multidisciplinary study sheds light on microenvironmental signals controlling cancer cell plasticity along EMT and suggests that hybrid and mesenchymal phenotypes arise through independent molecular paths.

Biodynamics: A novel quasi-first principles theory on the fundamental mechanisms of cellular function/dysfunction and the pharmacological modulation thereof.

Cellular function depends on heterogeneous dynamic intra-, inter-, and supramolecular structure-function relationships. However, the specific mechanisms by which cellular function is transduced from molecular systems, and by which cellular dysfunction arises from molecular dysfunction are poorly understood. We proposed previously that cellular function manifests as a molecular form of analog computing, in which specific time-dependent state transition fluxes within sets of molecular species ("molecular differential equations" (MDEs)) are sped and slowed in response to specific perturbations (inputs). In this work, we offer a theoretical treatment of the molecular mechanisms underlying cellular analog computing (which we refer to as "biodynamics"), focusing primarily on non-equilibrium (dynamic) intermolecular state transitions that serve as the principal means by which MDE systems are solved (the molecular equivalent of mathematical "integration"). Under these conditions, bound state occupancy is governed by kon and koff, together with the rates of binding partner buildup and decay. Achieving constant fractional occupancy over time depends on: 1) equivalence between kon and the rate of binding site buildup); 2) equivalence between koff and the rate of binding site decay; and 3) free ligand concentration relative to koff/kon (n · Kd, where n is the fold increase in binding partner concentration needed to achieve a given fractional occupancy). Failure to satisfy these conditions results in fractional occupancy well below that corresponding to n · Kd. The implications of biodynamics for cellular function/dysfunction and drug discovery are discussed.

Mapping the phenotypic repertoire of the cytoplasmic 2-Cys peroxiredoxin - Thioredoxin system. 1. Understanding commonalities and differences among cell types.

The system (PTTRS) formed by typical 2-Cys peroxiredoxins (Prx), thioredoxin (Trx), Trx reductase (TrxR), and sulfiredoxin (Srx) is central in antioxidant protection and redox signaling in the cytoplasm of eukaryotic cells. Understanding how the PTTRS integrates these functions requires tracing phenotypes to molecular properties, which is non-trivial. Here we analyze this problem based on a model that captures the PTTRS' conserved features. We have mapped the conditions that generate each distinct response to H2O2 supply rates (vsup), and estimated the parameters for thirteen human cell types and for Saccharomyces cerevisiae. The resulting composition-to-phenotype map yielded the following experimentally testable predictions. The PTTRS permits many distinct responses including ultra-sensitivity and hysteresis. However, nearly all tumor cell lines showed a similar response characterized by limited Trx-S- depletion and a substantial but self-limited gradual accumulation of hyperoxidized Prx at high vsup. This similarity ensues from strong correlations between the TrxR, Srx and Prx activities over cell lines, which contribute to maintain the Prx-SS reduction capacity in slight excess over the maximal steady state Prx-SS production. In turn, in erythrocytes, hepatocytes and HepG2 cells high vsup depletes Trx-S- and oxidizes Prx mainly to Prx-SS. In all nucleated human cells the Prx-SS reduction capacity defined a threshold separating two different regimes. At sub-threshold vsup the cytoplasmic H2O2 concentration is determined by Prx, nM-range and spatially localized, whereas at supra-threshold vsup it is determined by much less active alternative sinks and µM-range throughout the cytoplasm. The yeast shows a distinct response where the Prx Tsa1 accumulates in sulfenate form at high vsup. This is mainly due to an exceptional stability of Tsa1's sulfenate. The implications of these findings for thiol redox regulation and cell physiology are discussed. All estimates were thoroughly documented and provided, together with analytical approximations for system properties, as a resource for quantitative redox biology.

Synthetic mixed-signal computation in living cells.

Living cells implement complex computations on the continuous environmental signals that they encounter. These computations involve both analogue- and digital-like processing of signals to give rise to complex developmental programs, context-dependent behaviours and homeostatic activities. In contrast to natural biological systems, synthetic biological systems have largely focused on either digital or analogue computation separately. Here we integrate analogue and digital computation to implement complex hybrid synthetic genetic programs in living cells. We present a framework for building comparator gene circuits to digitize analogue inputs based on different thresholds. We then demonstrate that comparators can be predictably composed together to build band-pass filters, ternary logic systems and multi-level analogue-to-digital converters. In addition, we interface these analogue-to-digital circuits with other digital gene circuits to enable concentration-dependent logic. We expect that this hybrid computational paradigm will enable new industrial, diagnostic and therapeutic applications with engineered cells.

Hydrogen peroxide metabolism and sensing in human erythrocytes: a validated kinetic model and reappraisal of the role of peroxiredoxin II.

Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10(7)-10(8) M(-1) s(-1)) for H2O2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H2O2. However, this is inconsistent with extensive evidence that Prx2's contribution to H2O2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10(5) M(-1) s(-1) agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H2O2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Is Peroxiredoxin II's peroxidase activity strongly inhibited in human erythrocytes?

H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione peroxidase (GPx1) and the more recently discovered peroxiredoxin 2 (Prx2). However, the contribution of Prx2 to H2O2 consumption is still unclear. Prx2's high reactivity with H2O2 (kPrx2=10×10(7) M(-1)s(-1), kCat =7×10(7) M(-1)s(-1), kGPx1 =4×10(7) M(-1)s(-1)) and high abundance ([Prx2]= 570µM, [Cat]= 32µM, [GPx1]= 1µM) suggest that under low H2O2 supply rates it should consume >99% of the H2O2. However, extensive evidence indicates that in intact erythrocytes Prx2 contributes no more than Cat to H2O2 consumption. In order for this to be attained, Prx2's effective rate constant with H2O2would have to be just ~10(5) M(-1)s(-1), much lower than that determined in multiple experiments with the purified proteins. Nevertheless, nearly all Prx2 is oxidized within 1min of exposing erythrocytes to a H2O2 bolus, which is inconsistent with an irreversible inhibition. A mathematical model of the H2O2 metabolism in human erythrocytes [Benfeitas et al. (2014) Free Radic. Biol. Med.] where Prx2 either has a low kPrx2 or is subject to a strong (>99%) but readily reversible inhibition achieves quantitative agreement with detailed experimental observations of the responses of the redox status of Prx2 in human erythrocytes and suggests functional advantages of this design (see companion abstract). By contrast, a variant where Prx2 is fully active with kPrx2=10(8) M(-1)s(-1) shows important qualitative discrepancies. Altogether, these results suggest that Prx2's peroxidase activity is strongly inhibited in human erythrocytes. We acknowledge fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) co-financed by FEDER through the COMPETE program and by FCT.

Is the Peroxiredoxin 2/Thioredoxin/Thioredoxin Reductase system in human erythrocytes designed for redox signaling?

In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredoxin 2 (Prx2). Our previous analyses indicate that Prx2's peroxidase activity is subjected to a strong but quickly reversible inhibition (see companion abstract). If this activity is inhibited then the main role of Prx2 cannot be to eliminate H2O2. What functional advantages could then such an inhibition confer?We set up and validated a kinetic model of H2O2 metabolism human erythrocytes that shows quantitative agreement with extensive experimental observations. We then applied it to analyze the behavior of Prx2 and Trx under the H2O2 exposure dynamics that erythrocytes face in circulation. The significance of Prx2 inhibition was assessed by comparing the behavior of this model with that of an otherwise identical model lacking inhibition.Our analysis shows that Prx2 inhibition leads to 25-40% lower NADPH consumption under low to moderately high H2O2 supply (<0.8µM H2O2/s). Further, the inhibition extends the range where the concentrations of potential redox signaling readouts - H2O2, Prx2 sulfenic acid, Prx2 disulfide and Trx disulfide- show a proportional response to changes in H2O2 supply, covering practically the whole physiological range of the latter. This is desirable for analogic signal transduction and allows the Prx2/Trx/TrxR system to reliably transduce changes in H2O2 supply as changes in thiol oxidation. Finally, the inhibition allows other less abundant peroxiredoxins in the erythrocyte to be oxidized by H2O2 at physiological H2O2 supplies.Altogether, these results suggest that the postulated reversible inhibition of Prx2's peroxidase facilitates signal transduction by the peroxiredoxins and spares NADPH.We acknowledge: fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 co-financed by FEDER through the COMPETE program and by FCT (project PTDC/QUI-BIQ/119657/2010).

Integration between anticipatory blocking and redox signaling by the peroxiredoxin/thioredoxin/thioredoxin-reductase system.

Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other electrophylic compounds. Both H2O2 and these compounds can irreversibly modify protein thiols, with deleterious consequences. Induction of enzymatic defenses against those agents is too slow to avoid significant damage. Cells may solve this conundrum by reversibly "blocking" the thiols once H2O2 concentrations begin to increase. We term this mechanism "anticipatory blocking" because it acts in anticipation of irreversible damage upon detection of early signs of stress. Here we examine the design requirements for the Peroxiredoxin/Thioredoxin/Thioredoxin-Reductase/Protein-Dithiol System (PTTRDS) to effectively integrate H2O2 signaling and anticipatory blocking of protein dithiols as disulfides, and we compared them to the designs found in cells. To that effect, we developed a minimal model of the PTTRDS, and we defined a set of quantitative performance criteria that embody the requirements for (a) efficient scavenging capacity, (b) low NADPH consumption, (c) effective signal propagation, and (d) effective anticipatory blocking. We then sought the design principles (relationships among rate constants and species concentrations) that warrant fulfillment of all these criteria. Experimental data indicates that the design of the PTTRDS in human erythrocytes fulfills these principles and thus accomplishes effective integration between anticipatory blocking, antioxidant protection and redox signaling. A more general analysis suggests that the same principles hold in a wide variety of cell types and organisms. We acknowledge grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) financed by FEDER through the "Programa Operacional Factores de Competitividade, COMPETE" and by national funds through "FCT, Fundação para a Ciência e a Tecnologia".