Guido Santos

New Targets for Drug Discovery against Malaria

A mathematical model which predicts the intraerythrocytic stages of Plasmodium falciparum infection was developed using data from malaria-infected mice. Variables selected accounted for levels of healthy red blood cells, merozoite (Plasmodium asexual phase) infected red blood cells, gametocyte (Plasmodium sexual phase) infected red blood cells and a phenomenological variable which accounts for the mean activity of the immune system of the host. The model built was able to reproduce the behavior of three different scenarios of malaria. It predicts the later dynamics of malaria-infected humans well after the first peak of parasitemia, the qualitative response of malaria-infected monkeys to vaccination and the changes observed in malaria-infected mice when they are treated with antimalarial drugs. The mathematical model was used to identify new targets to be focused on drug design. Optimization methodologies were applied to identify five targets for minimizing the parasite load; four of the targets thus identified have never before been taken into account in drug design. The potential targets include: 1) increasing the death rate of the gametocytes, 2) decreasing the invasion rate of the red blood cells by the merozoites, 3) increasing the transformation of merozoites into gametocytes, 4) decreasing the activation of the immune system by the gametocytes, and finally 5) a combination of the previous target with decreasing the recycling rate of the red blood cells. The first target is already used in current therapies, whereas the remainders are proposals for potential new targets. Furthermore, the combined target (the simultaneous decrease of the activation of IS by gRBC and the decrease of the influence of IS on the recycling of hRBC) is interesting, since this combination does not affect the parasite directly. Thus, it is not expected to generate selective pressure on the parasites, which means that it would not produce resistance in Plasmodium.

Model-based genotype-phenotype mapping used to investigate gene signatures of immune sensitivity and resistance in melanoma micrometastasis

In this paper, we combine kinetic modelling and patient gene expression data analysis to elucidate biological mechanisms by which melanoma becomes resistant to the immune system and to immunotherapy. To this end, we systematically perturbed the parameters in a kinetic model and performed a mathematical analysis of their impact, thereby obtaining signatures associated with the emergence of phenotypes of melanoma immune sensitivity and resistance. Our phenotypic signatures were compared with published clinical data on pretreatment tumor gene expression in patients subjected to immunotherapy against metastatic melanoma. To this end, the differentially expressed genes were annotated with standard gene ontology terms and aggregated into metagenes. Our method sheds light on putative mechanisms by which melanoma may develop immunoresistance. Precisely, our results and the clinical data point to the existence of a signature of intermediate expression levels for genes related to antigen presentation that constitutes an intriguing resistance mechanism, whereby micrometastases are able to minimize the combined anti-tumor activity of complementary responses mediated by cytotoxic T cells and natural killer cells, respectively. Finally, we computationally explored the efficacy of cytokines used as low-dose co-adjuvants for the therapeutic anticancer vaccine to overcome tumor immunoresistance.

The (Mathematical) Modeling Process in Biosciences

In this communication, we introduce a general framework and discussion on the role of models and the modeling process in the field of biosciences. The objective is to sum up the common procedures during the formalization and analysis of a biological problem from the perspective of Systems Biology, which approaches the study of biological systems as a whole. We begin by presenting the definitions of (biological) system and model. Particular attention is given to the meaning of mathematical model within the context of biology. Then, we present the process of modeling and analysis of biological systems. Three stages are described in detail: conceptualization of the biological system into a model, mathematical formalization of the previous conceptual model and optimization and system management derived from the analysis of the mathematical model. All along this work the main features and shortcomings of the process are analyzed and a set of rules that could help in the task of modeling any biological system are presented. Special regard is given to the formative requirements and the interdisciplinary nature of this approach. We conclude with some general considerations on the challenges that modeling is posing to current biology.

Quantitative analysis of the processes and signaling events involved in early HIV-1 infection of T cells.

Lymphocyte invasion by HIV-1 is a complex, highly regulated process involving many different types of molecules that is prompted by the viruss association with viral receptors located at the cell-surface membrane that culminates in the formation of a fusion pore through which the virus enters the cell. A great deal of work has been done to identify the key actors in the process and determine the regulatory interactions; however, there have been no reports to date of attempts being made to fully understand the system dynamics through a systemic, quantitative modeling approach. In this paper, we introduce a dynamic mathematical model that integrates the available information on the molecular events involved in lymphocyte invasion. Our model shows that moesin activation is induced by virus signaling, while filamin-A is mobilized by the receptor capping. Actin disaggregation from the cap is facilitated by cofilin. Cofilin is inactivated by HIV-1 signaling in activated lymphocytes, while in resting lymphocytes another signal is required to activate cofilin in the later stages in order to accelerate the decay of the aggregated actin as a restriction factor for the viral entry. Furthermore, stopping the activation signaling of moesin is sufficient to liberate the actin filaments from the cap. The model also shows the positive effect of gelsolin on actin capping by means of the nucleation effect. These findings allow us to propose novel approaches in the search for new therapeutic strategies. In particular, gelsolin inhibition is seen as a promising target for preventing HIV-1 entry into lymphocytes, due to its role in facilitating the capping needed for the invasion. Also it is shown that HIV-1 should overcome the cortical actin barrier during early infection and predicts the different susceptibility of CD4+ T cells to be infected in terms of actin cytoskeleton dynamics driven by associated cellular factors.

Lipid Raft Size and Lipid Mobility in Non-raft Domains Increase during Aging and Are Exacerbated in APP/PS1 Mice Model of Alzheimer's Disease. Predictions from an Agent-Based Mathematical Model

A connection between lipid rafts and Alzheimers disease has been studied during the last decades. Mathematical modeling approaches have recently been used to correlate the effects of lipid composition changes in the physicochemical properties of raft-like membranes. Here we propose an agent based model to assess the effect of lipid changes in lipid rafts on the evolution and progression of Alzheimers disease using lipid profile data obtained in an established model of familial Alzheimers disease. We have observed that lipid raft size and lipid mobility in non-raft domains are two main factors that increase during age and are accelerated in the transgenic Alzheimers disease mouse model. The consequences of these changes are discussed in the context of neurotoxic amyloid β production. Our agent based model predicts that increasing sterols (mainly cholesterol) and long-chain polyunsaturated fatty acids (LCPUFA) (mainly DHA, docosahexaenoic acid) proportions in the membrane composition might delay the onset and progression of the disease.

Bacterial Adherence and Dwelling Probability: Two Drivers of Early Alveolar Infection by Streptococcus pneumoniae Identified in Multi-Level Mathematical Modeling

Pneumococcal infection is the most frequent cause of pneumonia, and one of the most prevalent diseases worldwide. The population groups at high risk of death from bacterial pneumonia are infants, elderly and immunosuppressed people. These groups are more vulnerable because they have immature or impaired immune systems, the efficacy of their response to vaccines is lower, and antibiotic treatment often does not take place until the inflammatory response triggered is already overwhelming. The immune response to bacterial lung infections involves dynamic interactions between several types of cells whose activation is driven by intracellular molecular networks. A feasible approach to the integration of knowledge and data linking tissue, cellular and intracellular events and the construction of hypotheses in this area is the use of mathematical modeling. For this paper, we used a multi-level computational model to analyse the role of cellular and molecular interactions during the first 10 h after alveolar invasion of Streptococcus pneumoniae bacteria. By “multi-level” we mean that we simulated the interplay between different temporal and spatial scales in a single computational model. In this instance, we included the intracellular scale of processes driving lung epithelial cell activation together with the scale of cell-to-cell interactions at the alveolar tissue. In our analysis, we combined systematic model simulations with logistic regression analysis and decision trees to find genotypic-phenotypic signatures that explain differences in bacteria strain infectivity. According to our simulations, pneumococci benefit from a high dwelling probability and a high proliferation rate during the first stages of infection. In addition to this, the model predicts that during the very early phases of infection the bacterial capsule could be an impediment to the establishment of the alveolar infection because it impairs bacterial colonization.

Model-Based Phenotypic Signatures Governing the Dynamics of the Stem and Semi-differentiated Cell Populations in Dysplastic Colonic Crypts

Mathematical modeling of cell differentiated in colonic crypts can contribute to a better understanding of basic mechanisms underlying colonic tissue organization, but also its deregulation during carcinogenesis and tumor progression. Here, we combined bifurcation analysis to assess the effect that time delay has in the complex interplay of stem cells and semi-differentiated cells at the niche of colonic crypts, and systematic model perturbation and simulation to find model-based phenotypes linked to cancer progression. The models suggest that stem cell and semi-differentiated cell population dynamics in colonic crypts can display chaotic behavior. In addition, we found that clinical profiling of colorectal cancer correlates with the in silico phenotypes proposed by the mathematical model. Further, potential therapeutic targets for chemotherapy resistant phenotypes are proposed, which in any case will require experimental validation.

Modelling and Analysis of Central Metabolism Operating Regulatory Interactions in Salt Stress Conditions in a L-Carnitine Overproducing E. coli Strain

Based on experimental data from E. coli cultures, we have devised a mathematical model in the GMA-power law formalism that describes the central and L-carnitine metabolism in and between two steady states, non-osmotic and hyperosmotic (0.3 M NaCl). A key feature of this model is the introduction of type of kinetic order, the osmotic stress kinetic orders (gOSn), derived from the power law general formalism, which represent the effect of osmotic stress in each metabolic process of the model. By considering the values of the gOSn linked to each metabolic process we found that osmotic stress has a positive and determinant influence on the increase in flux in energetic metabolism (glycolysis); L-carnitine biosynthesis production; the transformation/excretion of Acetyl-CoA into acetate and ethanol; the input flux of peptone into the cell; the anabolic use of pyruvate and biomass decomposition. In contrast, we found that although the osmotic stress has an inhibitory effect on the transformation flux from the glycolytic end products (pyruvate) to Acetyl-CoA, this inhibition is counteracted by other effects (the increase in pyruvate concentration) to the extent that the whole flux increases. In the same vein, the down regulation exerted by osmotic stress on fumarate uptake and its oxidation and the production and export of lactate and pyruvate are reversed by other factors up to the point that the first increased and the second remained unchanged. The model analysis shows that in osmotic conditions the energy and fermentation pathways undergo substantial rearrangement. This is illustrated by the observation that the increase in the fermentation fluxes is not connected with fluxes towards the tricaboxylic acid intermediates and the synthesis of biomass. The osmotic stress associated with these fluxes reflects these changes. All these observations support that the responses to salt stress observed in E. coli might be conserved in halophiles. Flux evolution during osmotic adaptations showed a hyperbolic (increasing or decreasing) pattern except in the case of peptone and fumarate uptake by the cell, which initially decreased. Finally, the model also throws light on the role of L-carnitine as osmoprotectant.

FPtool a software tool to obtain in silico genotype-phenotype signatures and fingerprints based on massive model simulations

Fptool is an intuitive tool that provides to the user a preliminary fingerprint of the behaviour simulated by a mathematical model of a biochemical network when comparing two biological scenarios defined by the user. Here we present the tool and we applied to an already published mathematical model of lung legionella infection. The fingerprint obtained correlates with the results obtained in the original article. This tool is optimal for the users that would like to obtain a fast and preliminary view of the qualitative behaviour of a mathematical model before deciding for more elaborate analyses.

A Mathematical model to predict Trichodesmiun spp. bloom during a sewage burst

A modified version of the Monod equation is proposed to understand the bloom emergence of Trichodesmiun spp. in presence of sewage outburst. The mathematical model predicts that the bacterial population does not affect the bloom growth. Also, the quotient between the bacteria during the sewage outburst and before is proposed to be equal to the quotient between the influx of limiting resources during burst and the influx before the burst.Abstract A modified version of the Monod equation is proposed to understand the bloom emergence of Trichodesmiun spp. in presence of sewage outburst. The mathematical model predicts that the bacterial population does not affect the bloom growth. Also, the quotient between the bacteria during the sewage outburst and before is proposed to be equal to the quotient between the influx of limiting resources during burst and the influx before the burst.A modified version of the Monod equation is proposed to understand the bloom emergence of Trichodesmiun spp. in presence of sewage outburst. The mathematical model predicts that the bacterial population does not affect the bloom growth. Also, the quotient between the bacteria during the sewage outburst and before is proposed to be equal to the quotient between the influx of limiting resources during burst and the influx before the burst.