BioUML-towards a universal research platform.

BioUML (https://www.biouml.org)-is a web-based integrated platform for systems biology and data analysis. It supports visual modelling and construction of hierarchical biological models that allow us to construct the most complex modular models of blood pressure regulation, skeletal muscle metabolism, COVID-19 epidemiology. BioUML has been integrated with git repositories where users can store their models and other data. We have also expanded the capabilities of BioUML for data analysis and visualization of biomedical data: (i) any programs and Jupyter kernels can be plugged into the BioUML platform using Docker technology; (ii) BioUML is integrated with the Galaxy and Galaxy Tool Shed; (iii) BioUML provides two-way integration with R and Python (Jupyter notebooks): scripts can be executed on the BioUML web pages, and BioUML functions can be called from scripts; (iv) using plug-in architecture, specialized viewers and editors can be added. For example, powerful genome browsers as well as viewers for molecular 3D structure are integrated in this way; (v) BioUML supports data analyses using workflows (own format, Galaxy, CWL, BPMN, nextFlow). Using these capabilities, we have initiated a new branch of the BioUML development-u-science-a universal scientific platform that can be configured for specific research requirements.

Physiologically Based Pharmacokinetic Modeling of Nanoparticle Biodistribution: A Review of Existing Models, Simulation Software, and Data Analysis Tools.

Cancer treatment and pharmaceutical development require targeted treatment and less toxic therapeutic intervention to achieve real progress against this disease. In this scenario, nanomedicine emerged as a reliable tool to improve drug pharmacokinetics and to translate to the clinical biologics based on large molecules. However, the ability of our body to recognize foreign objects together with carrier transport heterogeneity derived from the combination of particle physical and chemical properties, payload and surface modification, make the designing of effective carriers very difficult. In this scenario, physiologically based pharmacokinetic modeling can help to design the particles and eventually predict their ability to reach the target and treat the tumor. This effort is performed by scientists with specific expertise and skills and familiarity with artificial intelligence tools such as advanced software that are not usually in the "cords" of traditional medical or material researchers. The goal of this review was to highlight the advantages that computational modeling could provide to nanomedicine and bring together scientists with different background by portraying in the most simple way the work of computational developers through the description of the tools that they use to predict nanoparticle transport and tumor targeting in our body.

A Modular Mathematical Model of Exercise-Induced Changes in Metabolism, Signaling, and Gene Expression in Human Skeletal Muscle.

Skeletal muscle is the principal contributor to exercise-induced changes in human metabolism. Strikingly, although it has been demonstrated that a lot of metabolites accumulating in blood and human skeletal muscle during an exercise activate different signaling pathways and induce the expression of many genes in working muscle fibres, the systematic understanding of signaling-metabolic pathway interrelations with downstream genetic regulation in the skeletal muscle is still elusive. Herein, a physiologically based computational model of skeletal muscle comprising energy metabolism, Ca2+, and AMPK (AMP-dependent protein kinase) signaling pathways and the expression regulation of genes with early and delayed responses was developed based on a modular modeling approach and included 171 differential equations and more than 640 parameters. The integrated modular model validated on diverse including original experimental data and different exercise modes provides a comprehensive in silico platform in order to decipher and track cause-effect relationships between metabolic, signaling, and gene expression levels in skeletal muscle.

Multisensing Wearables for Real-Time Monitoring of Sweat Electrolyte Biomarkers During Exercise and Analysis on Their Correlation With Core Body Temperature.

Sweat secreted by the human eccrine sweat glands can provide valuable biomarker information during exercise. Real-time non-invasive biomarker recordings are therefore useful for evaluating the physiological conditions of an athlete such as their hydration status during endurance exercise. This work describes a wearable sweat biomonitoring patch incorporating printed electrochemical sensors into a plastic microfluidic sweat collector and data analysis that shows the real-time recorded sweat biomarkers can be used to predict a physiological biomarker. The system was placed on subjects carrying out an hour-long exercise session and results were compared to a wearable system using potentiometric robust silicon-based sensors and to commercially available HORIBA-LAQUAtwin devices. Both prototypes were applied to the real-time monitoring of sweat during cycling sessions and showed stable readings for around an hour. Analysis of the sweat biomarkers collected from the printed patch prototype shows that their real-time measurements correlate well (correlation coefficient ≥ 0.65) with other physiological biomarkers such as heart rate and regional sweat rate collected in the same session. We show for the first time, that the real-time sweat sodium and potassium concentration biomarker measurements from the printed sensors can be used to predict the core body temperature with root mean square error (RMSE) of 0.02 °C which is 71% lower compared to the use of only the physiological biomarkers. These results show that these wearable patch technologies are promising for real-time portable sweat monitoring analytical platforms, especially for athletes performing endurance exercise.

Mathematical modeling of antihypertensive therapy.

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.

Thoroughly Calibrated Modular Agent-Based Model of the Human Cardiovascular and Renal Systems for Blood Pressure Regulation in Health and Disease.

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.

Mathematical Modeling and Analyses of Interspike-Intervals of Spontaneous Activity in Afferent Neurons of the Zebrafish Lateral Line.

Without stimuli, hair cells spontaneously release neurotransmitter leading to spontaneous generation of action potentials (spikes) in innervating afferent neurons. We analyzed spontaneous spike patterns recorded from the lateral line of zebrafish and found that distributions of interspike intervals (ISIs) either have an exponential shape or an "L" shape that is characterized by a sharp decay but wide tail. ISI data were fitted to renewal-process models that accounted for the neuron refractory periods and hair-cell synaptic release. Modeling the timing of synaptic release using a mixture of two exponential distributions yielded the best fit for our ISI data. Additionally, lateral line ISIs displayed positive serial correlation and appeared to exhibit switching between faster and slower modes of spike generation. This pattern contrasts with previous findings from the auditory system where ISIs tended to have negative serial correlation due to synaptic depletion. We propose that afferent neuron innervation with multiple and heterogenous hair-cells synapses, each influenced by changes in calcium domains, can serve as a mechanism for the random switching behavior. Overall, our analyses provide evidence of how physiological similarities and differences between synapses and innervation patterns in the auditory, vestibular, and lateral line systems can lead to variations in spontaneous activity.

Computational master-regulator search reveals mTOR and PI3K pathways responsible for low sensitivity of NCI-H292 and A427 lung cancer cell lines to cytotoxic action of p53 activator Nutlin-3.

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.

Toward Community Standards and Software for Whole-Cell Modeling.

OBJECTIVE: Whole-cell (WC) modeling is a promising tool for biological research, bioengineering, and medicine. However, substantial work remains to create accurate comprehensive models of complex cells. METHODS: We organized the 2015 Whole-Cell Modeling Summer School to teach WC modeling and evaluate the need for new WC modeling standards and software by recoding a recently published WC model in the Systems Biology Markup Language. RESULTS: Our analysis revealed several challenges to representing WC models using the current standards. CONCLUSION: We, therefore, propose several new WC modeling standards, software, and databases. SIGNIFICANCE: We anticipate that these new standards and software will enable more comprehensive models.