Computational modeling of DLBCL predicts response to BH3-mimetics

In healthy cells, pro- and anti-apoptotic BCL2 family and BH3-only proteins are expressed in a delicate equilibrium. In contrast, this homeostasis is frequently perturbed in cancer cells due to the overexpression of anti-apoptotic BCL2 family proteins. Variability in the expression and sequestration of these proteins in Diffuse Large B cell Lymphoma (DLBCL) likely contributes to variability in response to BH3-mimetics. Successful deployment of BH3-mimetics in DLBCL requires reliable predictions of which lymphoma cells will respond. Here we show that a computational systems biology approach enables accurate prediction of the sensitivity of DLBCL cells to BH3-mimetics. We found that fractional killing of DLBCL, can be explained by cell-to-cell variability in the molecular abundances of signaling proteins. Importantly, by combining protein interaction data with a knowledge of genetic lesions in DLBCL cells, our in silico models accurately predict in vitro response to BH3-mimetics. Furthermore, through virtual DLBCL cells we predict synergistic combinations of BH3-mimetics, which we then experimentally validated. These results show that computational systems biology models of apoptotic signaling, when constrained by experimental data, can facilitate the rational assignment of efficacious targeted inhibitors in B cell malignancies, paving the way for development of more personalized approaches to treatment.

Non-Antigenic Modulation of Antigen Receptor (TCR) Cβ-FG Loop Modulates Signalling: Implications of External Factors Influencing T-Cell Responses

T-cell recognition of antigens is complex, leading to biochemical and cellular events that impart both specific and targeted immune responses. The end result is an array of cytokines that facilitate the direction and intensity of the immune reaction-such as T-cell proliferation, differentiation, macrophage activation, and B-cell isotype switching-all of which may be necessary and appropriate to eliminate the antigen and induce adaptive immunity. Using in silico docking to identify small molecules that putatively bind to the T-cell Cβ-FG loop, we have shown in vitro using an antigen presentation assay that T-cell signalling is altered. The idea of modulating T-cell signalling independently of antigens by directly targeting the FG loop is novel and warrants further study.

Mapping the use of computational modelling and simulation in clinics: A survey

In silico medicine describes the application of computational modelling and simulation (CM&S) to the study, diagnosis, treatment or prevention of a disease. Tremendous research advances have been achieved to facilitate the use of CM&S in clinical applications. Nevertheless, the uptake of CM&S in clinical practice is not always timely and accurately reflected in the literature. A clear view on the current awareness, actual usage and opinions from the clinicians is needed to identify barriers and opportunities for the future of in silico medicine. The aim of this study was capturing the state of CM&S in clinics by means of a survey toward the clinical community. Responses were collected online using the Virtual Physiological Human institute communication channels, engagement with clinical societies, hospitals and individual contacts, between 2020 and 2021. Statistical analyses were done with R. Participants (n = 163) responded from all over the world. Clinicians were mostly aged between 35 and 64 years-old, with heterogeneous levels of experience and areas of expertise (i.e., 48% cardiology, 13% musculoskeletal, 8% general surgery, 5% paediatrics). The CM&S terms "Personalised medicine" and "Patient-specific modelling" were the most well-known within the respondents. "In silico clinical trials" and "Digital Twin" were the least known. The familiarity with different methods depended on the medical specialty. CM&S was used in clinics mostly to plan interventions. To date, the usage frequency is still scarce. A well-recognized benefit associated to CM&S is the increased trust in planning procedures. Overall, the recorded level of trust for CM&S is high and not proportional to awareness level. The main barriers appear to be access to computing resources, perception that CM&S is slow. Importantly, clinicians see a role for CM&S expertise in their team in the future. This survey offers a snapshot of the current situation of CM&S in clinics. Although the sample size and representativity could be increased, the results provide the community with actionable data to build a responsible strategy for accelerating a positive uptake of in silico medicine. New iterations and follow-up activities will track the evolution of responses over time and contribute to strengthen the engagement with the medical community.

Understanding repertoire sequencing data through a multiscale computational model of the germinal center

Sequencing of B-cell and T-cell immune receptor repertoires helps us to understand the adaptive immune response, although it only provides information about the clonotypes (lineages) and their frequencies and not about, for example, their affinity or antigen (Ag) specificity. To further characterize the identified clones, usually with special attention to the particularly abundant ones (dominant), additional time-consuming or expensive experiments are generally required. Here, we present an extension of a multiscale model of the germinal center (GC) that we previously developed to gain more insight in B-cell repertoires. We compare the extent that these simulated repertoires deviate from experimental repertoires established from single GCs, blood, or tissue. Our simulations show that there is a limited correlation between clonal abundance and affinity and that there is large affinity variability among same-ancestor (same-clone) subclones. Our simulations suggest that low-abundance clones and subclones, might also be of interest since they may have high affinity for the Ag. We show that the fraction of plasma cells (PCs) with high B-cell receptor (BcR) mRNA content in the GC does not significantly affect the number of dominant clones derived from single GCs by sequencing BcR mRNAs. Results from these simulations guide data interpretation and the design of follow-up experiments.

Optimization of nutrient utilization efficiency and productivity for algal cultures under light and dark cycles using genome-scale model process control

Algal cultivations are strongly influenced by light and dark cycles. In this study, genome-scale metabolic models were applied to optimize nutrient supply during alternating light and dark cycles of Chlorella vulgaris. This approach lowered the glucose requirement by 75% and nitrate requirement by 23%, respectively, while maintaining high final biomass densities that were more than 80% of glucose-fed heterotrophic culture. Furthermore, by strictly controlling glucose feeding during the alternating cycles based on model-input, yields of biomass, lutein, and fatty acids per gram of glucose were more than threefold higher with cycling compared to heterotrophic cultivation. Next, the model was incorporated into open-loop and closed-loop control systems and compared with traditional fed-batch systems. Closed-loop systems which incorporated a feed-optimizing algorithm increased biomass yield on glucose more than twofold compared to standard fed-batch cultures for cycling cultures. Finally, the performance was compared to conventional proportional-integral-derivative (PID) controllers. Both simulation and experimental results exhibited superior performance for genome-scale model process control (GMPC) compared to traditional PID systems, reducing the overall measured value and setpoint error by 80% over 8 h. Overall, this approach provides researchers with the capability to enhance nutrient utilization and productivity of cell factories systematically by combining genome-scale models and controllers into an integrated platform with superior performance to conventional fed-batch and PID methodologies.

An agent-based approach for modelling and simulation of glycoprotein VI receptor diffusion, localisation and dimerisation in platelet lipid rafts

Receptor diffusion plays an essential role in cellular signalling via the plasma membrane microenvironment and receptor interactions, but the regulation is not well understood. To aid in understanding of the key determinants of receptor diffusion and signalling, we developed agent-based models (ABMs) to explore the extent of dimerisation of the platelet- and megakaryocyte-specific receptor for collagen glycoprotein VI (GPVI). This approach assessed the importance of glycolipid enriched raft-like domains within the plasma membrane that lower receptor diffusivity. Our model simulations demonstrated that GPVI dimers preferentially concentrate in confined domains and, if diffusivity within domains is decreased relative to outside of domains, dimerisation rates are increased. While an increased amount of confined domains resulted in further dimerisation, merging of domains, which may occur upon membrane rearrangements, was without effect. Modelling of the proportion of the cell membrane which constitutes lipid rafts indicated that dimerisation levels could not be explained by these alone. Crowding of receptors by other membrane proteins was also an important determinant of GPVI dimerisation. Together, these results demonstrate the value of ABM approaches in exploring the interactions on a cell surface, guiding the experimentation for new therapeutic avenues.

Modulation of transcription factor dynamics allows versatile information transmission

Cells detect changes in their environment and generate responses, often involving changes in gene expression. In this paper we use information theory and a simple transcription model to analyze whether the resulting gene expression serves to identify extracellular stimuli and assess their intensity when they are encoded in the amplitude, duration or frequency of pulses of a transcription factor's nuclear concentration (or activation state). We find, for all cases, that about three ranges of input strengths can be distinguished and that maximum information transmission occurs for fast and high activation threshold promoters. The three input modulation modes differ in the sensitivity to changes in the promoters parameters. Frequency modulation is the most sensitive and duration modulation, the least. This is key for signal identification: there are promoter parameters that yield a relatively high information transmission for duration or amplitude modulation and a much smaller value for frequency modulation. The reverse situation cannot be found with a single promoter transcription model. Thus, pulses of transcription factors can selectively activate the "frequency-tuned" promoter while prolonged nuclear accumulation would activate promoters of all three modes simultaneously. Frequency modulation is therefore highly selective and better suited than the other encoding modes for signal identification without requiring other mediators of the transduction process.

Genome-scale model of Pseudomonas aeruginosa metabolism unveils virulence and drug potentiation

Pseudomonas aeruginosa is one of the leading causes of hospital-acquired infections. To decipher the metabolic mechanisms associated with virulence and antibiotic resistance, we have developed an updated genome-scale model (GEM) of P. aeruginosa. The model (iSD1509) is an extensively curated, three-compartment, and mass-and-charge balanced BiGG model containing 1509 genes, the largest gene content for any P. aeruginosa GEM to date. It is the most accurate with prediction accuracies as high as 92.4% (gene essentiality) and 93.5% (substrate utilization). In iSD1509, we newly added a recently discovered pathway for ubiquinone-9 biosynthesis which is required for anaerobic growth. We used a modified iSD1509 to demonstrate the role of virulence factor (phenazines) in the pathogen survival within biofilm/oxygen-limited condition. Further, the model can mechanistically explain the overproduction of a drug susceptibility biomarker in the P. aeruginosa mutants. Finally, we use iSD1509 to demonstrate the drug potentiation by metabolite supplementation, and elucidate the mechanisms behind the phenotype, which agree with experimental results.

Software JimenaE allows efficient dynamic simulations of Boolean networks, centrality and system state analysis

The signal modelling framework JimenaE simulates dynamically Boolean networks. In contrast to SQUAD, there is systematic and not just heuristic calculation of all system states. These specific features are not present in CellNetAnalyzer and BoolNet. JimenaE is an expert extension of Jimena, with new optimized code, network conversion into different formats, rapid convergence both for system state calculation as well as for all three network centralities. It allows higher accuracy in determining network states and allows to dissect networks and identification of network control type and amount for each protein with high accuracy. Biological examples demonstrate this: (i) High plasticity of mesenchymal stromal cells for differentiation into chondrocytes, osteoblasts and adipocytes and differentiation-specific network control focusses on wnt-, TGF-beta and PPAR-gamma signaling. JimenaE allows to study individual proteins, removal or adding interactions (or autocrine loops) and accurately quantifies effects as well as number of system states. (ii) Dynamical modelling of cell-cell interactions of plant Arapidopsis thaliana against Pseudomonas syringae DC3000: We analyze for the first time the pathogen perspective and its interaction with the host. We next provide a detailed analysis on how plant hormonal regulation stimulates specific proteins and who and which protein has which type and amount of network control including a detailed heatmap of the A.thaliana response distinguishing between two states of the immune response. (iii) In an immune response network of dendritic cells confronted with Aspergillus fumigatus, JimenaE calculates now accurately the specific values for centralities and protein-specific network control including chemokine and pattern recognition receptors.

Hunting behavior of a solitary sailfish Istiophorus platypterus and estimated energy gain after prey capture

Foraging behavior and interaction with prey is an integral component of the ecological niche of predators but is inherently difficult to observe for highly mobile animals in the marine environment. Billfishes have been described as energy speculators, expending a large amount of energy foraging, expecting to offset high costs with periodic high energetic gain. Surface-based group feeding of sailfish, Istiophorus platypterus, is commonly observed, yet sailfish are believed to be largely solitary roaming predators with high metabolic requirements, suggesting that individual foraging also represents a major component of predator-prey interactions. Here, we use biologging data and video to examine daily activity levels and foraging behavior, estimate metabolic costs, and document a solitary predation event for a 40 kg sailfish. We estimate a median active metabolic rate of 218.9 ± 70.5 mgO₂ kg^-1 h^-1 which increased to 518.8 ± 586.3 mgO₂ kg^-1 h^-1 during prey pursuit. Assuming a successful predation, we estimate a daily net energy gain of 2.4 MJ (5.1 MJ acquired, 2.7 MJ expended), supporting the energy speculator model. While group hunting may be a common activity used by sailfish to acquire energy, our calculations indicate that opportunistic individual foraging events offer a net energy return that contributes to the fitness of these highly mobile predators.

Metabolic modelling of the human gut microbiome in type 2 diabetes patients in response to metformin treatment

The human gut microbiome has been associated with several metabolic disorders including type 2 diabetes mellitus. Understanding metabolic changes in the gut microbiome is important to elucidate the role of gut bacteria in regulating host metabolism. Here, we used available metagenomics data from a metformin study, together with genome-scale metabolic modelling of the key bacteria in individual and community-level to investigate the mechanistic role of the gut microbiome in response to metformin. Individual modelling predicted that species that are increased after metformin treatment have higher growth rates in comparison to species that are decreased after metformin treatment. Gut microbial enrichment analysis showed prior to metformin treatment pathways related to the hypoglycemic effect were enriched. Our observations highlight how the key bacterial species after metformin treatment have commensal and competing behavior, and how their cellular metabolism changes due to different nutritional environment. Integrating different diets showed there were specific microbial alterations between different diets. These results show the importance of the nutritional environment and how dietary guidelines may improve drug efficiency through the gut microbiota.

Jag1-Notch cis-interaction determines cell fate segregation in pancreatic development

The Notch ligands Jag1 and Dll1 guide differentiation of multipotent pancreatic progenitor cells (MPCs) into unipotent pro-acinar cells (PACs) and bipotent duct/endocrine progenitors (BPs). Ligand-mediated trans-activation of Notch receptors induces oscillating expression of the transcription factor Hes1, while ligand-receptor cis-interaction indirectly represses Hes1 activation. Despite Dll1 and Jag1 both displaying cis- and trans-interactions, the two mutants have different phenotypes for reasons not fully understood. Here, we present a mathematical model that recapitulates the spatiotemporal differentiation of MPCs into PACs and BPs. The model correctly captures cell fate changes in Notch pathway knockout mice and small molecule inhibitor studies, and a requirement for oscillatory Hes1 expression to maintain the multipotent state. Crucially, the model entails cell-autonomous attenuation of Notch signaling by Jag1-mediated cis-inhibition in MPC differentiation. The model sheds light on the underlying mechanisms, suggesting that cis-interaction is crucial for exiting the multipotent state, while trans-interaction is required for adopting the bipotent fate.

Position-dependent effects of RNA-binding proteins in the context of co-transcriptional splicing

Alternative splicing is an important step in eukaryotic mRNA pre-processing which increases the complexity of gene expression programs, but is frequently altered in disease. Previous work on the regulation of alternative splicing has demonstrated that splicing is controlled by RNA-binding proteins (RBPs) and by epigenetic DNA/histone modifications which affect splicing by changing the speed of polymerase-mediated pre-mRNA transcription. The interplay of these different layers of splicing regulation is poorly understood. In this paper, we derived mathematical models describing how splicing decisions in a three-exon gene are made by combinatorial spliceosome binding to splice sites during ongoing transcription. We additionally take into account the effect of a regulatory RBP and find that the RBP binding position within the sequence is a key determinant of how RNA polymerase velocity affects splicing. Based on these results, we explain paradoxical observations in the experimental literature and further derive rules explaining why the same RBP can act as inhibitor or activator of cassette exon inclusion depending on its binding position. Finally, we derive a stochastic description of co-transcriptional splicing regulation at the single-cell level and show that splicing outcomes show little noise and follow a binomial distribution despite complex regulation by a multitude of factors. Taken together, our simulations demonstrate the robustness of splicing outcomes and reveal that quantitative insights into kinetic competition of co-transcriptional events are required to fully understand this important mechanism of gene expression diversity.

Genetic complementation fosters evolvability in complex fitness landscapes

The ability of natural selection to optimize traits depends on the topology of the genotype-fitness map (fitness landscape). Epistatic interactions produce rugged fitness landscapes, where adaptation is constrained by the presence of low-fitness intermediates. Here, we used simulations to explore how evolvability in rugged fitness landscapes is influenced by genetic complementation, a process whereby different sequence variants mutually compensate for their deleterious mutations. We designed our model inspired by viral populations, in which genetic variants are known to interact frequently through coinfection. Our simulations indicate that genetic complementation enables a more efficient exploration of rugged fitness landscapes. Although this benefit may be undermined by genetic parasites, its overall effect on evolvability remains positive in populations that exhibit strong relatedness between interacting sequences. Similar processes could operate in contexts other than viral coinfection, such as in the evolution of ploidy.

Evidence that endosperm turgor pressure both promotes and restricts seed growth and size

In plants, as in animals, organ growth depends on mechanical interactions between cells and tissues, and is controlled by both biochemical and mechanical cues. Here, we investigate the control of seed size, a key agronomic trait, by mechanical interactions between two compartments: the endosperm and the testa. By combining experiments with computational modelling, we present evidence that endosperm pressure plays two antagonistic roles: directly driving seed growth, but also indirectly inhibiting it through tension it generates in the surrounding testa, which promotes wall stiffening. We show that our model can recapitulate wild type growth patterns, and is consistent with the small seed phenotype of the haiku2 mutant, and the results of osmotic treatments. Our work suggests that a developmental regulation of endosperm pressure is required to prevent a precocious reduction of seed growth rate induced by force-dependent seed coat stiffening.

Three-dimensional unsupervised probabilistic pose reconstruction (3D-UPPER) for freely moving animals

A key step in understanding animal behaviour relies in the ability to quantify poses and movements. Methods to track body landmarks in 2D have made great progress over the last few years but accurate 3D reconstruction of freely moving animals still represents a challenge. To address this challenge here we develop the 3D-UPPER algorithm, which is fully automated, requires no a priori knowledge of the properties of the body and can also be applied to 2D data. We find that 3D-UPPER reduces by [Formula: see text] fold the error in 3D reconstruction of mouse body during freely moving behaviour compared with the traditional triangulation of 2D data. To achieve that, 3D-UPPER performs an unsupervised estimation of a Statistical Shape Model (SSM) and uses this model to constrain the viable 3D coordinates. We show, by using simulated data, that our SSM estimator is robust even in datasets containing up to 50% of poses with outliers and/or missing data. In simulated and real data SSM estimation converges rapidly, capturing behaviourally relevant changes in body shape associated with exploratory behaviours (e.g. with rearing and changes in body orientation). Altogether 3D-UPPER represents a simple tool to minimise errors in 3D reconstruction while capturing meaningful behavioural parameters.

A new method of ventilation inhomogeneity assessment based on a simulation study using clinical data on congenital diaphragmatic hernia cases

Congenital Diaphragmatic Hernia (CDH) is a diaphragm defect associated with lung hypoplasia and ventilation inhomogeneity (VI). The affected neonates are usually born with respiratory failure and require mechanical ventilation after birth. However, significant interindividual VI differences make ventilation difficult. So far, there are no clinical methods of VI assessment that could be applied to optimize ventilation at the bedside. A new VI index is a ratio of time constants T₁/T₂ of gas flows in both lungs. Pressure-controlled ventilation simulations were conducted using an infant hybrid (numerical-physical) respiratory simulator connected to a ventilator. The parameters of the respiratory system model and ventilator settings were based on retrospective clinical data taken from three neonates (2, 2.6, 3.6 kg) treated in the Paediatric Teaching Clinical Hospital of the Medical University of Warsaw. We searched for relationships between respiratory system impedance (Z) and ventilation parameters: work of breathing (WOB), peak inspiratory pressure (PIP), and mean airway pressure (MAP). The study showed the increased VI described by the T₁/T₂ index value highly correlated with elevated Z, WOB, PIP and MAP (0.8-0.9, the Spearman correlation coefficients were significant at P < 0.001). It indicates that the T₁/T₂ index may help to improve the ventilation therapy of CDH neonates.

Mechanism based therapies enable personalised treatment of hypertrophic cardiomyopathy

Cardiomyopathies have unresolved genotype-phenotype relationships and lack disease-specific treatments. Here we provide a framework to identify genotype-specific pathomechanisms and therapeutic targets to accelerate the development of precision medicine. We use human cardiac electromechanical in-silico modelling and simulation which we validate with experimental hiPSC-CM data and modelling in combination with clinical biomarkers. We select hypertrophic cardiomyopathy as a challenge for this approach and study genetic variations that mutate proteins of the thick (MYH7^R403Q/+) and thin filaments (TNNT2^R92Q/+, TNNI3^R21C/+) of the cardiac sarcomere. Using in-silico techniques we show that the destabilisation of myosin super relaxation observed in hiPSC-CMs drives disease in virtual cells and ventricles carrying the MYH7^R403Q/+ variant, and that secondary effects on thin filament activation are necessary to precipitate slowed relaxation of the cell and diastolic insufficiency in the chamber. In-silico modelling shows that Mavacamten corrects the MYH7^R403Q/+ phenotype in agreement with hiPSC-CM experiments. Our in-silico model predicts that the thin filament variants TNNT2^R92Q/+ and TNNI3^R21C/+ display altered calcium regulation as central pathomechanism, for which Mavacamten provides incomplete salvage, which we have corroborated in TNNT2^R92Q/+ and TNNI3^R21C/+ hiPSC-CMs. We define the ideal characteristics of a novel thin filament-targeting compound and show its efficacy in-silico. We demonstrate that hybrid human-based hiPSC-CM and in-silico studies accelerate pathomechanism discovery and classification testing, improving clinical interpretation of genetic variants, and directing rational therapeutic targeting and design.

Predicting stress response and improved protein overproduction in Bacillus subtilis

Bacillus subtilis is a well-characterized microorganism and a model for the study of Gram-positive bacteria. The bacterium can produce proteins at high densities and yields, which has made it valuable for industrial bioproduction. Like other cell factories, metabolic modeling of B. subtilis has discovered ways to optimize its metabolism toward various applications. The first genome-scale metabolic model (M-model) of B. subtilis was published more than a decade ago and has been applied extensively to understand metabolism, to predict growth phenotypes, and served as a template to reconstruct models for other Gram-positive bacteria. However, M-models are ill-suited to simulate the production and secretion of proteins as well as their proteomic response to stress. Thus, a new generation of metabolic models, known as metabolism and gene expression models (ME-models), has been initiated. Here, we describe the reconstruction and validation of a ME model of B. subtilis, iJT964-ME. This model achieved higher performance scores on the prediction of gene essentiality as compared to the M-model. We successfully validated the model by integrating physiological and omics data associated with gene expression responses to ethanol and salt stress. The model further identified the mechanism by which tryptophan synthesis is upregulated under ethanol stress. Further, we employed iJT964-ME to predict amylase production rates under two different growth conditions. We analyzed these flux distributions and identified key metabolic pathways that permitted the increase in amylase production. Models like iJT964-ME enable the study of proteomic response to stress and the illustrate the potential for optimizing protein production in bacteria.

In Silico Modelling to Assess the Electrical and Thermal Disturbance Provoked by a Metal Intracoronary Stent during Epicardial Pulsed Electric Field Ablation

Background: Pulsed Electric Field (PEF) ablation has been recently proposed to ablate cardiac ganglionic plexi (GP) aimed to treat atrial fibrillation. The effect of metal intracoronary stents in the vicinity of the ablation electrode has not been yet assessed. Methods: A 2D numerical model was developed accounting for the different tissues involved in PEF ablation with an irrigated ablation device. A coronary artery (with and without a metal intracoronary stent) was considered near the ablation source (0.25 and 1 mm separation). The 1000 V/cm threshold was used to estimate the ‘PEF-zone’. Results: The presence of the coronary artery (with or without stent) distorts the E-field distribution, creating hot spots (higher E-field values) in the front and rear of the artery, and cold spots (lower E-field values) on the sides of the artery. The value of the E-field inside the coronary artery is very low (~200 V/cm), and almost zero with a metal stent. Despite this distortion, the PEF-zone contour is almost identical with and without artery/stent, remaining almost completely confined within the fat layer in any case. The mentioned hot spots of E-field translate into a moderate temperature increase (<48 °C) in the area between the artery and electrode. These thermal side effects are similar for pulse intervals of 10 and 100 μs. Conclusions: The presence of a metal intracoronary stent near the ablation device during PEF ablation simply ‘amplifies’ the E-field distortion already caused by the presence of the vessel. This distortion may involve moderate heating (<48 °C) in the tissue between the artery and ablation electrode without associated thermal damage.

Molecular models of multiple sclerosis severity identify heterogeneity of pathogenic mechanisms

While autopsy studies identify many abnormalities in the central nervous system (CNS) of subjects dying with neurological diseases, without their quantification in living subjects across the lifespan, pathogenic processes cannot be differentiated from epiphenomena. Using machine learning (ML), we searched for likely pathogenic mechanisms of multiple sclerosis (MS). We aggregated cerebrospinal fluid (CSF) biomarkers from 1305 proteins, measured blindly in the training dataset of untreated MS patients (N = 129), into models that predict past and future speed of disability accumulation across all MS phenotypes. Healthy volunteers (N = 24) data differentiated natural aging and sex effects from MS-related mechanisms. Resulting models, validated (Rho 0.40-0.51, p < 0.0001) in an independent longitudinal cohort (N = 98), uncovered intra-individual molecular heterogeneity. While candidate pathogenic processes must be validated in successful clinical trials, measuring them in living people will enable screening drugs for desired pharmacodynamic effects. This will facilitate drug development making, it hopefully more efficient and successful.

[Review on the speciation analysis methods of platinum group metals in environmental media]

Platinum group metals (PGMs) present a variety of forms in the environment, and analysis of speciation is essential for identifying their ecological risk. Here, we reviewed the methods for the morphological analysis of three major PGMs (platinum, palladium and rhodium) in the environment, including chemical sequential extraction, hyphenated techniques for instruments, computer simulations. We outlined the types, characteristics and applications of these methods, elaborated the weaknesses, and provided prospects for future development. Among them, chemical sequential extraction is universally applied in the morphological analysis of solid-phase samples, with diverse extraction conditions and procedures proposed in the current study. However, it has not been well standardized. The hyphenated techniques for instruments have significant advantages for the determination of elemental forms in solution, of which capillary electrophoresis system can separate similar substances with the same electrophoresis ability. Liquid chromatography systems have better performance in terms of separation capacity and detection limit. The computer simulations further expand the access to morphological analysis, enabling complex morphological calculations. It was proposed to combine multiple methods in the future to continuously improve the accuracy of analytical techniques by complementing and optimizing each other.

Pulsed Electric Field Ablation of Epicardial Autonomic Ganglia: Computer Analysis of Monopolar Electric Field across the Tissues Involved

BACKGROUND AND OBJECTIVES

Pulsed Electric Field (PEF) ablation has been proposed as a non-thermal energy to treat atrial fibrillation (AF) by epicardial ablation of ganglionated plexi (GP), which are embedded within epicardial fat. Our objective was to study the distribution of the electric field through the involved tissues (fat, GPs, myocardium and blood) during epicardial PEF ablation.

METHODS

A two-dimensional model was built considering different tissue layers below the ablation device which consists of an irrigated electrode. The 1000 V/cm threshold was used to estimate the 'PEF-zone'.

RESULTS

The PEF-zone was almost 100% circumscribed in the epicardial fat layer, with very little incidence in the myocardium. The presence of the saline on the epicardial fat causes the PEF-zone to spread laterally around the electrode from ~5 mm to ~15 mm, relatively independently of how embedded the electrode is in the saline layer. For a saline layer well spread over the tissue surface and an electrode fully embedded in the saline layer, the PEF-zone width decreases as the fat layer thickens: from ~15 mm for fat thickness of 1 and 2 mm, down to ~10 mm for fat thickness of 5 mm. The presence of a GP in the center of the fat layer hardly affects the size of the PEF-zone, but significantly alters the distribution of the electric field around the GP, resulting in progressively lower values than in the surrounding adipose tissue as the fat layer thickness increased.

CONCLUSIONS

Our results suggest how some procedural (irrigation) and anatomical parameters (fat thicknesses and presence of GPs) could be relevant in terms of the size of the tissue area affected by pulsed field ablation.

Model-based exploration of hypokalemia in dairy cows

Hypokalemia in dairy cows, which is characterized by too low serum potassium levels, is a severe mineral disorder that can be life threatening. In this paper, we explore different originating conditions of hypokalemia-reduced potassium intake, increased excretion, acid-base disturbances, and increased insulin-by using a dynamic mathematical model for potassium balance in non-lactating and lactating cows. The simulations confirm observations described in literature. They illustrate, for example, that changes in dietary intake or excretion highly effect intracellular potassium levels, whereas extracellular levels vary only slightly. Simulations also show that the higher the potassium content in the diet, the more potassium is excreted with urine. Application of the mathematical model assists in experimental planning and therefore contributes to the 3R strategy: reduction, refinement and replacement of animal experiments.

Integrase-mediated differentiation circuits improve evolutionary stability of burdensome and toxic functions in E. coli

Advances in synthetic biology, bioengineering, and computation allow us to rapidly and reliably program cells with increasingly complex and useful functions. However, because the functions we engineer cells to perform are typically burdensome to cell growth, they can be rapidly lost due to the processes of mutation and natural selection. Here, we show that a strategy of terminal differentiation improves the evolutionary stability of burdensome functions in a general manner by realizing a reproductive and metabolic division of labor. To implement this strategy, we develop a genetic differentiation circuit in Escherichia coli using unidirectional integrase-recombination. With terminal differentiation, differentiated cells uniquely express burdensome functions driven by the orthogonal T7 RNA polymerase, but their capacity to proliferate is limited to prevent the propagation of advantageous loss-of-function mutations that inevitably occur. We demonstrate computationally and experimentally that terminal differentiation increases duration and yield of high-burden expression and that its evolutionary stability can be improved with strategic redundancy. Further, we show this strategy can even be applied to toxic functions. Overall, this study provides an effective, generalizable approach for protecting burdensome engineered functions from evolutionary degradation.

Nonlinear model of infection wavy oscillation of COVID-19 in Japan based on diffusion kinetics

The infectious propagation of SARS-CoV-2 is continuing worldwide, and specifically, Japan is facing severe circumstances. Medical resource maintenance and action limitations remain the central measures. An analysis of long-term follow-up reports in Japan shows that the infection number follows a unique wavy oscillation, increasing and decreasing over time. However, only a few studies explain the infection wavy oscillation. This study introduces a novel nonlinear mathematical model of the new infection wavy oscillation by applying the macromolecule diffusion theory. In this model, the diffusion coefficient that depends on population density gives nonlinearity in infection propagation. As a result, our model accurately simulated infection wavy oscillations, and the infection wavy oscillation frequency and amplitude were closely linked with the recovery rate of infected individuals. In conclusion, our model provides a novel nonlinear contact infection analysis framework.

Network analysis reveals that the tumor suppressor lncRNA GAS5 acts as a double-edged sword in response to DNA damage in gastric cancer

The lncRNA GAS5 acts as a tumor suppressor and is downregulated in gastric cancer (GC). In contrast, E2F1, an important transcription factor and tumor promoter, directly inhibits miR-34c expression in GC cell lines. Furthermore, in the corresponding GC cell lines, lncRNA GAS5 directly targets E2F1. However, lncRNA GAS5 and miR-34c remain to be studied in conjunction with GC. Here, we present a dynamic Boolean network to classify gene regulation between these two non-coding RNAs (ncRNAs) in GC. This is the first study to show that lncRNA GAS5 can positively regulate miR-34c in GC through a previously unknown molecular pathway coupling lncRNA/miRNA. We compared our network to several in-vivo/in-vitro experiments and obtained an excellent agreement. We revealed that lncRNA GAS5 regulates miR-34c by targeting E2F1. Additionally, we found that lncRNA GAS5, independently of p53, inhibits GC proliferation through the ATM/p38 MAPK signaling pathway. Accordingly, our results support that E2F1 is an engaging target of drug development in tumor growth and aggressive proliferation of GC, and favorable results can be achieved through tumor suppressor lncRNA GAS5/miR-34c axis in GC. Thus, our findings unlock a new avenue for GC treatment in response to DNA damage by these ncRNAs.

Drug transporters OAT1 and OAT3 have specific effects on multiple organs and gut microbiome as revealed by contextualized metabolic network reconstructions

In vitro and in vivo studies have established the organic anion transporters OAT1 (SLC22A6, NKT) and OAT3 (SLC22A8) among the main multi-specific "drug" transporters. They also transport numerous endogenous metabolites, raising the possibility of drug-metabolite interactions (DMI). To help understand the role of these drug transporters on metabolism across scales ranging from organ systems to organelles, a formal multi-scale analysis was performed. Metabolic network reconstructions of the omics-alterations resulting from Oat1 and Oat3 gene knockouts revealed links between the microbiome and human metabolism including reactions involving small organic molecules such as dihydroxyacetone, alanine, xanthine, and p-cresol-key metabolites in independent pathways. Interestingly, pairwise organ-organ interactions were also disrupted in the two Oat knockouts, with altered liver, intestine, microbiome, and skin-related metabolism. Compared to older models focused on the "one transporter-one organ" concept, these more sophisticated reconstructions, combined with integration of a multi-microbial model and more comprehensive metabolomics data for the two transporters, provide a considerably more complex picture of how renal "drug" transporters regulate metabolism across the organelle (e.g. endoplasmic reticulum, Golgi, peroxisome), cellular, organ, inter-organ, and inter-organismal scales. The results suggest that drugs interacting with OAT1 and OAT3 can have far reaching consequences on metabolism in organs (e.g. skin) beyond the kidney. Consistent with the Remote Sensing and Signaling Theory (RSST), the analysis demonstrates how transporter-dependent metabolic signals mediate organ crosstalk (e.g., gut-liver-kidney) and inter-organismal communication (e.g., gut microbiome-host).

Integrated network-based multiple computational analyses for identification of co-expressed candidate genes associated with neurological manifestations of COVID-19

'Tripartite network' (TN) and 'combined gene network' (CGN) were constructed and their hub-bottleneck and driver nodes (44 genes) were evaluated as 'target genes' (TG) to identify 21 'candidate genes' (CG) and their relationship with neurological manifestations of COVID-19. TN was developed using neurological symptoms of COVID-19 found in literature. Under query genes (TG of TN), co-expressed genes were identified using pair-wise mutual information to genes available in RNA-Seq autopsy data of frontal cortex of COVID-19 victims. CGN was constructed with genes selected from TN and co-expressed in COVID-19. TG and their connecting genes of respective networks underwent functional analyses through findings of their enrichment terms and pair-wise 'semantic similarity scores' (SSS). A new integrated 'weighted harmonic mean score' was formulated assimilating values of SSS and STRING-based 'combined score' of the selected TG-pairs, which provided CG-pairs with properties of CGs as co-expressed and 'indispensable nodes' in CGN. Finally, six pairs sharing seven 'prevalent CGs' (ADAM10, ADAM17, AKT1, CTNNB1, ESR1, PIK3CA, FGFR1) showed linkages with the phenotypes (a) directly under neurodegeneration, neurodevelopmental diseases, tumour/cancer and cellular signalling, and (b) indirectly through other CGs under behavioural/cognitive and motor dysfunctions. The pathophysiology of 'prevalent CGs' has been discussed to interpret neurological phenotypes of COVID-19.

Cross clinical-experimental-computational qualification of in silico drug trials on human cardiac purkinje cells for proarrhythmia risk prediction

The preclinical identification of drug-induced cardiotoxicity and its translation into human risk are still major challenges in pharmaceutical drug discovery. The ICH S7B Guideline and Q&A on Clinical and Nonclinical Evaluation of QT/QTc Interval Prolongation and Proarrhythmic Potential promotes human in silico drug trials as a novel tool for proarrhythmia risk assessment. To facilitate the use of in silico data in regulatory submissions, explanatory control compounds should be tested and documented to demonstrate consistency between predictions and the historic validation data. This study aims to quantify drug-induced electrophysiological effects on in silico cardiac human Purkinje cells, to compare them with existing in vitro rabbit data, and to assess their accuracy for clinical pro-arrhythmic risk predictions. The effects of 14 reference compounds were quantified in simulations with a population of in silico human cardiac Purkinje models. For each drug dose, five electrophysiological biomarkers were quantified at three pacing frequencies, and results compared with available in vitro experiments and clinical proarrhythmia reports. Three key results were obtained: 1) In silico, repolarization abnormalities in human Purkinje simulations predicted drug-induced arrhythmia for all risky compounds, showing higher predicted accuracy than rabbit experiments; 2) Drug-induced electrophysiological changes observed in human-based simulations showed a high degree of consistency with in vitro rabbit recordings at all pacing frequencies, and depolarization velocity and action potential duration were the most consistent biomarkers; 3) discrepancies observed for dofetilide, sotalol and terfenadine are mainly caused by species differences between humans and rabbit. Taken together, this study demonstrates higher accuracy of in silico methods compared to in vitro animal models for pro-arrhythmic risk prediction, as well as a high degree of consistency with in vitro experiments commonly used in safety pharmacology, supporting the potential for industrial and regulatory adoption of in silico trials for proarrhythmia prediction.

HOXA9 has the hallmarks of a biological switch with implications in blood cancers

Blood malignancies arise from the dysregulation of haematopoiesis. The type of blood cell and the specific order of oncogenic events initiating abnormal growth ultimately determine the cancer subtype and subsequent clinical outcome. HOXA9 plays an important role in acute myeloid leukaemia (AML) prognosis by promoting blood cell expansion and altering differentiation; however, the function of HOXA9 in other blood malignancies is still unclear. Here, we highlight the biological switch and prognosis marker properties of HOXA9 in AML and chronic myeloproliferative neoplasms (MPN). First, we establish the ability of HOXA9 to stratify AML patients with distinct cellular and clinical outcomes. Then, through the use of a computational network model of MPN, we show that the self-activation of HOXA9 and its relationship to JAK2 and TET2 can explain the branching progression of JAK2/TET2 mutant MPN patients towards divergent clinical characteristics. Finally, we predict a connection between the RUNX1 and MYB genes and a suppressive role for the NOTCH pathway in MPN diseases.

Molecular pathways enhance drug response prediction using transfer learning from cell lines to tumors and patient-derived xenografts

Computational models have been successful in predicting drug sensitivity in cancer cell line data, creating an opportunity to guide precision medicine. However, translating these models to tumors remains challenging. We propose a new transfer learning workflow that transfers drug sensitivity predicting models from large-scale cancer cell lines to both tumors and patient derived xenografts based on molecular pathways derived from genomic features. We further compute feature importance to identify pathways most important to drug response prediction. We obtained good performance on tumors (AUROC = 0.77) and patient derived xenografts from triple negative breast cancers (RMSE = 0.11). Using feature importance, we highlight the association between ER-Golgi trafficking pathway in everolimus sensitivity within breast cancer patients and the role of class II histone deacetylases and interlukine-12 in response to drugs for triple-negative breast cancer. Pathway information support transfer of drug response prediction models from cell lines to tumors and can provide biological interpretation underlying the predictions, serving as a steppingstone towards usage in clinical setting.

An expanded whole-cell model of E. coli links cellular physiology with mechanisms of growth rate control

Growth and environmental responses are essential for living organisms to survive and adapt to constantly changing environments. In order to simulate new conditions and capture dynamic responses to environmental shifts in a developing whole-cell model of E. coli, we incorporated additional regulation, including dynamics of the global regulator guanosine tetraphosphate (ppGpp), along with dynamics of amino acid biosynthesis and translation. With the model, we show that under perturbed ppGpp conditions, small molecule feedback inhibition pathways, in addition to regulation of expression, play a role in ppGpp regulation of growth. We also found that simulations with dysregulated amino acid synthesis pathways provide average amino acid concentration predictions that are comparable to experimental results but on the single-cell level, concentrations unexpectedly show regular fluctuations. Additionally, during both an upshift and downshift in nutrient availability, the simulated cell responds similarly with a transient increase in the mRNA:rRNA ratio. This additional simulation functionality should support a variety of new applications and expansions of the E. coli Whole-Cell Modeling Project.

Dynamic cybergenetic control of bacterial co-culture composition via optogenetic feedback

Communities of microbes play important roles in natural environments and hold great potential for deploying division-of-labor strategies in synthetic biology and bioproduction. However, the difficulty of controlling the composition of microbial consortia over time hinders their optimal use in many applications. Here, we present a fully automated, high-throughput platform that combines real-time measurements and computer-controlled optogenetic modulation of bacterial growth to implement precise and robust compositional control of a two-strain E. coli community. In addition, we develop a general framework for dynamic modeling of synthetic genetic circuits in the physiological context of E. coli and use a host-aware model to determine the optimal control parameters of our closed-loop compositional control system. Our platform succeeds in stabilizing the strain ratio of multiple parallel co-cultures at arbitrary levels and in changing these targets over time, opening the door for the implementation of dynamic compositional programs in synthetic bacterial communities.

Reconstruction of a catalogue of genome-scale metabolic models with enzymatic constraints using GECKO 2.0

Genome-scale metabolic models (GEMs) have been widely used for quantitative exploration of the relation between genotype and phenotype. Streamlined integration of enzyme constraints and proteomics data into such models was first enabled by the GECKO toolbox, allowing the study of phenotypes constrained by protein limitations. Here, we upgrade the toolbox in order to enhance models with enzyme and proteomics constraints for any organism with a compatible GEM reconstruction. With this, enzyme-constrained models for the budding yeasts Saccharomyces cerevisiae, Yarrowia lipolytica and Kluyveromyces marxianus are generated to study their long-term adaptation to several stress factors by incorporation of proteomics data. Predictions reveal that upregulation and high saturation of enzymes in amino acid metabolism are common across organisms and conditions, suggesting the relevance of metabolic robustness in contrast to optimal protein utilization as a cellular objective for microbial growth under stress and nutrient-limited conditions. The functionality of GECKO is expanded with an automated framework for continuous and version-controlled update of enzyme-constrained GEMs, also producing such models for Escherichia coli and Homo sapiens. In this work, we facilitate the utilization of enzyme-constrained GEMs in basic science, metabolic engineering and synthetic biology purposes.

Artificial neural networks enable genome-scale simulations of intracellular signaling

Mammalian cells adapt their functional state in response to external signals in form of ligands that bind receptors on the cell-surface. Mechanistically, this involves signal-processing through a complex network of molecular interactions that govern transcription factor activity patterns. Computer simulations of the information flow through this network could help predict cellular responses in health and disease. Here we develop a recurrent neural network framework constrained by prior knowledge of the signaling network with ligand-concentrations as input and transcription factor-activity as output. Applied to synthetic data, it predicts unseen test-data (Pearson correlation r = 0.98) and the effects of gene knockouts (r = 0.8). We stimulate macrophages with 59 different ligands, with and without the addition of lipopolysaccharide, and collect transcriptomics data. The framework predicts this data under cross-validation (r = 0.8) and knockout simulations suggest a role for RIPK1 in modulating the lipopolysaccharide response. This work demonstrates the feasibility of genome-scale simulations of intracellular signaling.

Dynamic modeling of Nrf2 pathway activation in liver cells after toxicant exposure

Cells are exposed to oxidative stress and reactive metabolites every day. The Nrf2 signaling pathway responds to oxidative stress by upregulation of antioxidants like glutathione (GSH) to compensate the stress insult and re-establish homeostasis. Although mechanisms describing the interaction between the key pathway constituents Nrf2, Keap1 and p62 are widely reviewed and discussed in literature, quantitative dynamic models bringing together these mechanisms with time-resolved data are limited. Here, we present an ordinary differential equation (ODE) based dynamic model to describe the dynamic response of Nrf2, Keap1, Srxn1 and GSH to oxidative stress caused by the soft-electrophile diethyl maleate (DEM). The time-resolved data obtained by single-cell confocal microscopy of green fluorescent protein (GFP) reporters and qPCR of the Nrf2 pathway components complemented with siRNA knock down experiments, is accurately described by the calibrated mathematical model. We show that the quantitative model can describe the activation of the Nrf2 pathway by compounds with a different mechanism of activation, including drugs which are known for their ability to cause drug induced liver-injury (DILI) i.e., diclofenac (DCF) and omeprazole (OMZ). Finally, we show that our model can reveal differences in the processes leading to altered activation dynamics amongst DILI inducing drugs.

Mini-batch optimization enables training of ODE models on large-scale datasets

Quantitative dynamic models are widely used to study cellular signal processing. A critical step in modelling is the estimation of unknown model parameters from experimental data. As model sizes and datasets are steadily growing, established parameter optimization approaches for mechanistic models become computationally extremely challenging. Mini-batch optimization methods, as employed in deep learning, have better scaling properties. In this work, we adapt, apply, and benchmark mini-batch optimization for ordinary differential equation (ODE) models, thereby establishing a direct link between dynamic modelling and machine learning. On our main application example, a large-scale model of cancer signaling, we benchmark mini-batch optimization against established methods, achieving better optimization results and reducing computation by more than an order of magnitude. We expect that our work will serve as a first step towards mini-batch optimization tailored to ODE models and enable modelling of even larger and more complex systems than what is currently possible.

Boolean modelling in plant biology

Signalling and genetic networks underlie most biological processes and are often complex, containing many highly connected components. Modelling these networks can provide insight into mechanisms but is challenging given that rate parameters are often not well defined. Boolean modelling, in which components can only take on a binary value with connections encoded by logic equations, is able to circumvent some of these challenges, and has emerged as a viable tool to probe these complex networks. In this review, we will give an overview of Boolean modelling, with a specific emphasis on its use in plant biology. We review how Boolean modelling can be used to describe biological networks and then discuss examples of its applications in plant genetics and plant signalling.

Complementary resource preferences spontaneously emerge in diauxic microbial communities

Many microbes grow diauxically, utilizing the available resources one at a time rather than simultaneously. The properties of communities of microbes growing diauxically remain poorly understood, largely due to a lack of theory and models of such communities. Here, we develop and study a minimal model of diauxic microbial communities assembling in a serially diluted culture. We find that unlike co-utilizing communities, diauxic community assembly repeatably and spontaneously leads to communities with complementary resource preferences, namely communities where species prefer different resources as their top choice. Simulations and theory explain that the emergence of complementarity is driven by the disproportionate contribution of the top choice resource to the growth of a diauxic species. Additionally, we develop a geometric approach for analyzing serially diluted communities, with or without diauxie, which intuitively explains several additional emergent community properties, such as the apparent lack of species which grow fastest on a resource other than their most preferred resource. Overall, our work provides testable predictions for the assembly of natural as well as synthetic communities of diauxically shifting microbes.

Temporal assortment of cooperators in the spatial prisoner's dilemma

We study a spatial, one-shot prisoner's dilemma (PD) model in which selection operates on both an organism's behavioral strategy (cooperate or defect) and its decision of when to implement that strategy, which we depict as an organism's choice of one point in time, out of a set of discrete time slots, at which to carry out its PD strategy. Results indicate selection for cooperators across various time slots and parameter settings, including parameter settings in which cooperation would not evolve in an exclusively spatial model-as in work investigating exogenously imposed temporal networks. Moreover, in the presence of time slots, cooperators' portion of the population grows even under different combinations of spatial structure, transition rules, and update dynamics, though rates of cooperator fixation decline under pairwise comparison and synchronous updating. These findings indicate that, under certain evolutionary processes, merely existing in time and space promotes the evolution of cooperation.

Using predictive machine learning models for drug response simulation by calibrating patient-specific pathway signatures

The utility of pathway signatures lies in their capability to determine whether a specific pathway or biological process is dysregulated in a given patient. These signatures have been widely used in machine learning (ML) methods for a variety of applications including precision medicine, drug repurposing, and drug discovery. In this work, we leverage highly predictive ML models for drug response simulation in individual patients by calibrating the pathway activity scores of disease samples. Using these ML models and an intuitive scoring algorithm to modify the signatures of patients, we evaluate whether a given sample that was formerly classified as diseased, could be predicted as normal following drug treatment simulation. We then use this technique as a proxy for the identification of potential drug candidates. Furthermore, we demonstrate the ability of our methodology to successfully identify approved and clinically investigated drugs for four different cancers, outperforming six comparable state-of-the-art methods. We also show how this approach can deconvolute a drugs' mechanism of action and propose combination therapies. Taken together, our methodology could be promising to support clinical decision-making in personalized medicine by simulating a drugs' effect on a given patient.

Olive leaf-derived PPAR agonist complex induces collagen IV synthesis in human skin models

Introduction

Peroxisome proliferator‐activated receptor (PPAR) agonists are known to modulate the synthesis of dermal lipids and proteins including collagens. Olive (Olea europaea) leaves have been reported to contain PPAR‐binding ligands. Collagen IV, a major dermal‐epidermal junction (DEJ) protein, degrades with both age and disease. Here, we report the formulation of a novel multi‐ligand complex, Linefade, and its effects on collagen IV synthesis.

Methods

Linefade prepared from the leaves of Olea europaea contains 2% w/w plant extract solids dissolved in a mixture of glyceryl monoricinoleate and dimethyl isosorbide. In silico docking was performed with PPAR‐α (PDB ID: 2P54). Linefade was evaluated for PPAR‐α‐dependent transcription in a luciferase reporter assay system. Cell viability and collagen IV levels in human dermal fibroblast cultures were measured using the MTT method and ELISA assay, respectively. Transcriptome analysis was conducted on a full‐thickness reconstituted human skin (EpiDermFT) model. Ex vivo cell viability and collagen IV immunostaining were performed on human skin explants.

Results

In silico docking model of the major constituents (oleanolic acid and glyceryl monoricinoleate) produced a co‐binding affinity of −6.7 Kcal/mole. Linefade significantly increased PPAR‐α transcriptional activity in CHO cells and collagen IV synthesis in adult human dermal fibroblasts. Transcriptome analysis revealed that 1% Linefade modulated the expression of 280 genes with some related to epidermal differentiation, DEJ, PPAR, Nrf2 and retinoid pathways. An ex vivo human explant study showed that 1% Linefade, delivered via a triglycerides excipient, increased collagen IV levels along the dermal–epidermal junction by 52%.

Conclusion

In silico modelling and in vitro and ex vivo analyses confirmed Linefade‐mediated activation of PPAR‐α and stimulation of collagen IV synthesis.

Humans use minimum cost movements in a whole-body task

Humans have elegant bodies that allow gymnastics, piano playing, and tool use, but understanding how they do this in detail is difficult because their musculoskeletal systems are extraordinarily complicated. Nonetheless, common movements like walking and reaching can be stereotypical, and a very large number of studies have shown their energetic cost to be a major factor. In contrast, one might think that general movements are very individuated and intractable, but our previous study has shown that in an arbitrary set of whole-body movements used to trace large-scale closed curves, near-identical posture sequences were chosen across different subjects, both in the average trajectories of the body's limbs and in the variance within trajectories. The commonalities in that result motivate explanations for its generality. One explanation could be that humans also choose trajectories that are economical in cost. To test this hypothesis, we situate the tracing data within a forty eight degree of freedom human dynamic model that allows the computation of movement cost. Using the model to compare movement cost data from nominal tracings against various perturbed tracings shows that the latter are more energetically expensive, inferring that the original traces were chosen on the basis of minimum cost.

Commentary on Robinson et al. (2021): Evaluating theories of change for public health policies using computer model discovery methods

Recent developments in computer modelling—known as model discovery—could help to confirm the mechanisms underpinning Robinson and colleagues’ important early findings for the effectiveness of minimum unit pricing, and to test the complete theory of change underpinning this crucial evaluation.

Design principles of collateral sensitivity-based dosing strategies

Collateral sensitivity (CS)-based antibiotic treatments, where increased resistance to one antibiotic leads to increased sensitivity to a second antibiotic, may have the potential to limit the emergence of antimicrobial resistance. However, it remains unclear how to best design CS-based treatment schedules. To address this problem, we use mathematical modelling to study the effects of pathogen- and drug-specific characteristics for different treatment designs on bacterial population dynamics and resistance evolution. We confirm that simultaneous and one-day cycling treatments could supress resistance in the presence of CS. We show that the efficacy of CS-based cycling therapies depends critically on the order of drug administration. Finally, we find that reciprocal CS is not essential to suppress resistance, a result that significantly broadens treatment options given the ubiquity of one-way CS in pathogens. Overall, our analyses identify key design principles of CS-based treatment strategies and provide guidance to develop treatment schedules to suppress resistance.

Localization of Na+ channel clusters in narrowed perinexi of gap junctions enhances cardiac impulse transmission via ephaptic coupling: a model study

Abstract

It has been proposed that when gap junctional coupling is reduced in cardiac tissue, action potential propagation can be supported via ephaptic coupling, a mechanism mediated by negative electric potentials occurring in narrow intercellular clefts of intercalated discs (IDs). Recent studies showed that sodium (Na⁺) channels form clusters near gap junction plaques in nanodomains called perinexi, where the ID cleft is even narrower. To examine the electrophysiological relevance of Na⁺ channel clusters being located in perinexi, we developed a 3D finite element model of two longitudinally abutting cardiomyocytes, with a central Na⁺ channel cluster on the ID membranes. When this cluster was located in the perinexus of a closely positioned gap junction plaque, varying perinexal width greatly modulated impulse transmission from one cell to the other, with narrow perinexi potentiating ephaptic coupling. This modulation occurred via the interplay of Na⁺ currents, extracellular potentials in the cleft and patterns of current flow within the cleft. In contrast, when the Na⁺ channel cluster was located remotely from the gap junction plaque, this modulation by perinexus width largely disappeared. Interestingly, the Na⁺ current in the ID membrane of the pre‐junctional cell switched from inward to outward during excitation, thus contributing ions to the activating channels on the post‐junctional ID membrane. In conclusion, these results indicate that the localization of Na⁺ channel clusters in the perinexi of gap junction plaques is crucial for ephaptic coupling, which is furthermore greatly modulated by perinexal width. These findings are relevant for a comprehensive understanding of cardiac excitation.

Key points

Ephaptic coupling is a cardiac conduction mechanism involving nanoscale‐level interactions between the sodium (Na⁺) current and the extracellular potential in narrow intercalated disc clefts.

When gap junctional coupling is reduced, ephaptic coupling acts in conjunction with the classical cardiac conduction mechanism based on gap junctional current flow.

In intercalated discs, Na⁺ channels form clusters that are preferentially located in the periphery of gap junction plaques, in nanodomains known as perinexi, but the electrophysiological role of these perinexi has never been examined.

In our new 3D finite element model of two cardiac cells abutting each other with their intercalated discs, a Na⁺ channel cluster located inside a narrowed perinexus facilitated impulse transmission via ephaptic coupling.

Our simulations demonstrate the role of narrowed perinexi as privileged sites for ephaptic coupling in pathological situations when gap junctional coupling is decreased.

Development and qualification of a machine learning algorithm for automated hair counting

OBJECTIVE

Determining the amount of hair on the scalp has always been an important metric of patient satisfaction for hair growth and hair retention technologies. While simple in concept, this measurement is a difficult, resource intensive task for the dermatologist and the research scientist. Specifically, counting and measuring hair in phototrichogram images is very time consuming and labour intensive. Due to cost, often only a fraction of available images is manually analysed. There is a need for an automated method that can significantly increase speed and throughput while reducing the cost of counting and measuring hair in phototrichogram images.

METHODS

Recent advances in machine learning and deep convolutional neural networks (deep learning) have led to a revolution in the analysis of image, video, speech, text and other sensor data. Image diagnostics have seen remarkable improvements with completely automated methods outperforming both human experts and human-engineered analysis methods. Deep learning methods can also provide speed and cost benefits. To enable use of a deep learning, we created a data set of 288 manually annotated phototrichogram images with marked location and length of each hair (the training dataset). We designed a custom neural network architecture and custom image processing algorithms to best utilize the available training data and to maximize performance for hair counting and length measurement. The performance of the algorithm was qualified by comparing hair count and length measurements to an independent ground truth method, the semi-manual Canfield's Hair Metrix method.

RESULTS

Leveraging deep neural networks, we have developed capability to apply machine learning to reduce the time needed to acquire data from phototrichograms of patients' scalp from months to seconds. Our algorithm enables fast and fully automated hair counting and length measurement. The algorithm shows high agreement with human manually assisted analysis (ground truth).

CONCLUSIONS

We have trained and deployed an algorithm utilizing this technology and have demonstrated the reproducibility, accuracy and speed of this algorithm that, once deployed, requires little to no recurring cost or manual intervention for its operation. The method allows fast analysis of large number of images, reducing study cost and significantly reducing study analysis time.

Risk of proximal femoral nail antirotation (PFNA) implant failure upon different lateral femoral wall thickness in intertrochanteric fracture: a finite element analysis

Proximal Femoral Nail Antirotation (PFNA) has been commonly used to treat intertrochanteric fractures, despite the risk of implant failure. The integrity of the femur could influence the risk of implant failure. This study evaluated the influence of lateral femoral wall thickness on the potential of implant failure. A finite element model of the hip was reconstructed from the Computed Tomography of a female patient. Five intertrochanteric fracture models at different lateral femoral wall thickness (T1 = 27.6 mm, T2 = 25.4 mm, T3 = 23.4 mm, T4 = 21.4 mm, and T5 = 19.3 mm) were then created and fixed with PFNA. We simulated a critical loading condition by a high loading case during walking. Elastoplastic material models with yield stress and failure strain were applied to the bone and implant in which breakage can be simulated using the element deletion function. In addition, the stress and displacement of the implant and femur were analysed. Implant breakage occurred at the sides of the proximal nail canal in cases of T4 and T5 which was further supported by the higher maximum von Mises stress and nail displacement. The increased stress and displacement of the implant may implicate a reduction of stability and risk of implant failure. We suggested that precaution shall be taken when the wall thickness was less than 21.4 mm.

A genome-scale metabolic model of Saccharomyces cerevisiae that integrates expression constraints and reaction thermodynamics

Eukaryotic organisms play an important role in industrial biotechnology, from the production of fuels and commodity chemicals to therapeutic proteins. To optimize these industrial systems, a mathematical approach can be used to integrate the description of multiple biological networks into a single model for cell analysis and engineering. One of the most accurate models of biological systems include Expression and Thermodynamics FLux (ETFL), which efficiently integrates RNA and protein synthesis with traditional genome-scale metabolic models. However, ETFL is so far only applicable for E. coli. To adapt this model for Saccharomyces cerevisiae, we developed yETFL, in which we augmented the original formulation with additional considerations for biomass composition, the compartmentalized cellular expression system, and the energetic costs of biological processes. We demonstrated the ability of yETFL to predict maximum growth rate, essential genes, and the phenotype of overflow metabolism. We envision that the presented formulation can be extended to a wide range of eukaryotic organisms to the benefit of academic and industrial research.