Search

The angiotensin-converting enzyme (ACE) gene (ACE) insertion/deletion (I/D) polymorphism raises the possibility of personalising ACE inhibitor therapy to optimise its efficiency and reduce side effects in genetically distinct subgroups. However, the extent of its influence among these subgroups is unknown. Therefore, we extended our computational model of blood pressure regulation to investigate the effect of the ACE I/D polymorphism on haemodynamic parameters in humans undergoing antihypertensive therapy. The model showed that the dependence of blood pressure on serum ACE activity is a function of saturation and therefore, the lack of association between ACE I/D and blood pressure levels may be due to high ACE activity in specific populations. Additionally, in an extended model simulating the effects of different classes of antihypertensive drugs, we explored the relationship between ACE I/D and the efficacy of inhibitors of the renin–angiotensin–aldosterone system. The model predicted that the response of cardiovascular and renal parameters to treatment directly depends on ACE activity. However, significant differences in parameter changes were observed only between groups with high and low ACE levels, while different ACE I/D genotypes within the same group had similar changes in absolute values. We conclude that a single genetic variant is responsible for only a small fraction of heredity in treatment success and its predictive value is limited.

Hypertension is a multifactorial disease arising from complex pathophysiological pathways. Individual characteristics of patients result in different responses to various classes of antihypertensive medications. Therefore, evaluating the efficacy of therapy based on in silico predictions is an important task. This study is a continuation of research on the modular agent-based model of the cardiovascular and renal systems (presented in the previously published article). In the current work, we included in the model equations simulating the response to antihypertensive therapies with different mechanisms of action. For this, we used the pharmacodynamic effects of the angiotensin II receptor blocker losartan, the calcium channel blocker amlodipine, the angiotensin-converting enzyme inhibitor enalapril, the direct renin inhibitor aliskiren, the thiazide diuretic hydrochlorothiazide, and the β-blocker bisoprolol. We fitted therapy parameters based on known clinical trials for all considered medications, and then tested the model’s ability to show reasonable dynamics (expected by clinical observations) after treatment with individual drugs and their dual combinations in a group of virtual patients with hypertension. The extended model paves the way for the next step in personalized medicine that is adapting the model parameters to a real patient and predicting his response to antihypertensive therapy. The model is implemented in the BioUML software and is available at https://gitlab.sirius-web.org/virtual-patient/antihypertensive-treatment-modeling.

Single-nucleotide variants in enhancers or promoters may affect gene transcription by altering transcription factor binding sites. Here the authors present a meta-analysis empowered by a new statistical method covering thousands of ChIP-Seq experiments resulting in the identification of more than 500 thousand allele-specific binding (ASB) events in the human genome.

Abstract Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction‐based models and packages that extend the core with features suited to other model types including constraint‐based models, reaction‐diffusion models, logical network models, and rule‐based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single‐cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

Here we present a modular agent-based mathematical model of the human cardiovascular and renal systems. It integrates the previous models primarily developed by A. C. Guyton, F. Karaaslan, K. M. Hallow, and Y. V. Solodyannikov. We performed the model calibration to find an equilibrium state within the normal vital sign ranges for a healthy adult. We verified the model’s abilities to reproduce equilibrium states with abnormal physiological values related to different combinations of cardiovascular diseases (such as systemic hypertension, chronic heart failure, pulmonary hypertension, etc.). For the model creation and validation, we involved over 200 scientific studies covering known models of the human cardiovascular and renal functions, biosimulation platforms, and clinical measurements of physiological quantities in normal and pathological conditions. We compiled detailed documentation describing all equations, parameters and variables of the model with justification of all formulas and values. The model is implemented in BioUML and available in the web-version of the software.

Abstract Background The search for molecular biomarkers of early-onset colorectal cancer (CRC) is an important but still quite challenging and unsolved task. Detection of CpG methylation in human DNA obtained from blood or stool has been proposed as a promising approach to a noninvasive early diagnosis of CRC. Thousands of abnormally methylated CpG positions in CRC genomes are often located in non-coding parts of genes. Novel bioinformatic methods are thus urgently needed for multi-omics data analysis to reveal causative biomarkers with a potential driver role in early stages of cancer. Methods We have developed a method for finding potential causal relationships between epigenetic changes (DNA methylations) in gene regulatory regions that affect transcription factor binding sites (TFBS) and gene expression changes. This method also considers the topology of the involved signal transduction pathways and searches for positive feedback loops that may cause the carcinogenic aberrations in gene expression. We call this method “Walking pathways”, since it searches for potential rewiring mechanisms in cancer pathways due to dynamic changes in the DNA methylation status of important gene regulatory regions (“epigenomic walking”). Results In this paper, we analysed an extensive collection of full genome gene-expression data (RNA-seq) and DNA methylation data of genomic CpG islands (using Illumina methylation arrays) generated from a sample of tumor and normal gut epithelial tissues of 300 patients with colorectal cancer (at different stages of the disease) (data generated in the EU-supported SysCol project). Identification of potential epigenetic biomarkers of DNA methylation was performed using the fully automatic multi-omics analysis web service “My Genome Enhancer” (MGE) (my-genome-enhancer.com). MGE uses the database on gene regulation TRANSFAC®, the signal transduction pathways database TRANSPATH®, and software that employs AI (artificial intelligence) methods for the analysis of cancer-specific enhancers. Conclusions The identified biomarkers underwent experimental testing on an independent set of blood samples from patients with colorectal cancer. As a result, using advanced methods of statistics and machine learning, a minimum set of 6 biomarkers was selected, which together achieve the best cancer detection potential. The markers include hypermethylated positions in regulatory regions of the following genes: CALCA, ENO1, MYC, PDX1, TCF7, ZNF43.

Here we describe a toolkit for presenting physiologically based pharmacokinetic (PBPK) models in a modular graphical view in the BioUML platform. Firstly, we demonstrate the BioUML capabilities for PBPK modeling tested on an existing model of nanoparticles delivery to solid tumors in mice. Secondly, we provide guidance on the conversion of the PBPK model code from a text modeling language like Berkeley Madonna to a visual modular diagram in the BioUML. We give step-by-step explanations of the model transformation and demonstrate that simulation results from the original model are exactly the same as numerical results obtained for the transformed model. The main advantage of the proposed approach is its clarity and ease of perception. Additionally, the modular representation serves as a simplified and convenient base for in silico investigation of the model and reduces the risk of technical errors during its reuse and extension by concomitant biochemical processes. In summary, this article demonstrates that BioUML can be used as an alternative and robust tool for PBPK modeling.

Abstract Background Small molecule Nutlin-3 reactivates p53 in cancer cells by interacting with the complex between p53 and its repressor Mdm-2 and causing an increase in cancer cell apoptosis. Therefore, Nutlin-3 has potent anticancer properties. Clinical and experimental studies of Nutlin-3 showed that some cancer cells may lose sensitivity to this compound. Here we analyze possible mechanisms for insensitivity of cancer cells to Nutlin-3. Methods We applied upstream analysis approach implemented in geneXplain platform (genexplain.com) using TRANSFAC® database of transcription factors and their binding sites in genome and using TRANSPATH® database of signal transduction network with associated software such as Match™ and Composite Module Analyst (CMA). Results Using genome-wide gene expression profiling we compared several lung cancer cell lines and showed that expression programs executed in Nutlin-3 insensitive cell lines significantly differ from that of Nutlin-3 sensitive cell lines. Using artificial intelligence approach embed in CMA software, we identified a set of transcription factors cooperatively binding to the promoters of genes up-regulated in the Nutlin-3 insensitive cell lines. Graph analysis of signal transduction network upstream of these transcription factors allowed us to identify potential master-regulators responsible for maintaining such low sensitivity to Nutlin-3 with the most promising candidate mTOR, which acts in the context of activated PI3K pathway. These finding were validated experimentally using an array of chemical inhibitors. Conclusions We showed that the Nutlin-3 insensitive cell lines are actually highly sensitive to the dual PI3K/mTOR inhibitor NVP-BEZ235, while no responding to either PI3K –specific LY294002 nor Bcl-XL specific 2,3-DCPE compounds.