Fibromyalgia pathogenesis explained by a neuroendocrine multistable model.

Fibromyalgia (FM) is a central disorder characterized by chronic pain, fatigue, insomnia, depression, and other minor symptoms. Knowledge about pathogenesis is lacking, diagnosis difficult, clinical approach puzzling, and patient management disappointing. We conducted a theoretical study based on literature data and computational analysis, aimed at developing a comprehensive model of FM pathogenesis and addressing suitable therapeutic targets. We started from the evidence that FM must involve a dysregulation of central pain processing, is female prevalent, suggesting a role for the hypothalamus-pituitary-gonadal (HPG) axis, and is stress-related, suggesting a role for the HP-adrenocortical (HPA) axis. Central pathogenesis was supposed to involve a pain processing loop system including the thalamic ventroposterolateral nucleus (VPL), the primary somatosensory cortex (SSC), and the thalamic reticular nucleus (TRN). For decreasing GABAergic and/or increasing glutamatergic transmission, the loop system crosses a bifurcation point, switching from monostable to bistable, and converging on a high-firing-rate steady state supposed to be the pathogenic condition. Thereafter, we showed that GABAergic transmission is positively correlated with gonadal-hormone-derived neurosteroids, notably allopregnanolone, whereas glutamatergic transmission is positively correlated with stress-induced glucocorticoids, notably cortisol. Finally, we built a dynamic model describing a multistable, double-inhibitory loop between HPG and HPA axes. This system has a high-HPA/low-HPG steady state, allegedly reached in females under combined premenstrual/postpartum brain allopregnanolone withdrawal and stress condition, driving the thalamocortical loop to the high-firing-rate steady state, and explaining the connection between endocrine and neural mechanisms in FM pathogenesis. Our model accounts for FM female prevalence and stress correlation, suggesting the use of neurosteroid drugs as a possible solution to currently unsolved problems in the clinical treatment of the disease.

A closed-loop negative feedback model for the pancreas: A new paradigm and pathway to a cure.

BACKGROUND AND AIMS: To develop a model that describes how the pancreas functions, how the rate of synthesis of digestive enzymes is regulated, and finally what puts the pancreas to rest between meals. METHODS: We applied the principals of control theory to previously published canine data to develop a model for how the canine pancreas functions. Using this model, we then describe the steps needed to apply this model to the human pancreas. RESULTS: This new closed-loop negative feedback model describes what regulates digestive enzyme synthesis. This model is based on basolateral exocytosis of butyrylcholinesterase (BCHE) into the interstitial space. It is this level of BCHE * BCHE activity that controls the rate of canine pancreas digestive enzyme synthesis, and in the absence of stimulation from the vagus nerve, puts the pancreas to rest between meals. CONCLUSIONS: Finding secretagogue-specific inhibitory enzymes in the human pancreas that are analogous to BCHE in the canine, and blocking its associated receptors, may lead to a cure for human pancreatitis.

Iron Homeostatic Regulation in Saccharomyces cerevisiae: Introduction to a Computational Modeling Method.

Over the past 30 years, much has been learned regarding iron homeostatic regulation in budding yeast, S. cerevisiae, including the identity of many of the proteins and molecular-level regulatory mechanisms involved. Most advances have involved inferring such mechanisms based on the analysis of iron-dysregulation phenotypes arising in various genetic mutant strains. Still lacking is a cellular- or system-level understanding of iron homeostasis. These experimental advances are summarized in this review, and a method for developing cellular-level regulatory mechanisms in yeast is presented. The method employs the results of Mössbauer spectroscopy of whole cells and organelles, iron quantification of the same, and ordinary differential equation-based mathematical models. Current models are simplistic when compared to the complexity of iron homeostasis in real cells, yet they hold promise as a useful, perhaps even required, complement to the popular genetics-based approach. The fundamental problem in comprehending cellular regulatory mechanisms is that, given the complexities involved, different molecular-level mechanisms can often give rise to virtually indistinguishable cellular phenotypes. Mathematical models cannot eliminate this problem, but they can minimize it.

No dose adjustment of metformin or substrates of organic cation transporters (OCT)1 and OCT2 and multidrug and toxin extrusion protein (MATE)1/2K with fostemsavir coadministration based on modeling approaches.

Fostemsavir is an approved gp120-directed attachment inhibitor and prodrug for the treatment of human immunodeficiency virus type 1 infection in combination with other antiretrovirals (ARVs) in heavily treatment-experienced adults with multi-drug resistance, intolerance, or safety concerns with their current ARV regimen. Initial in vitro studies indicated that temsavir, the active moiety of fostemsavir, and its metabolites, inhibited organic cation transporter (OCT)1, OCT2, and multidrug and toxin extrusion transporters (MATEs) at tested concentration of 100 uM, although risk assessment based on the current Food and Drug Administration in vitro drug-drug interaction (DDI) guidance using the mechanistic static model did not reveal any clinically relevant inhibition on OCTs and MATEs. However, a DDI risk was flagged with EMA static model predictions. Hence, a physiologically based pharmacokinetic (PBPK) model of fostemsavir/temsavir was developed to further assess the DDI risk potential of OCT and MATEs inhibition by temsavir and predict changes in metformin (a sensitive OCT and MATEs substrate) exposure. No clinically relevant impact on metformin concentrations across a wide range of temsavir concentrations was predicted; therefore, no dose adjustment is recommended for metformin when co-administered with fostemsavir.

Prey group defense and hunting cooperation among generalist-predators induce complex dynamics: a mathematical study.

Group defense in prey and hunting cooperation in predators are two important ecological phenomena and can occur concurrently. In this article, we consider cooperative hunting in generalist predators and group defense in prey under a mathematical framework to comprehend the enormous diversity the model could capture. To do so, we consider a modified Holling-Tanner model where we implement Holling type IV functional response to characterize grazing pattern of predators where prey species exhibit group defense. Additionally, we allow a modification in the attack rate of predators to quantify the hunting cooperation among them. The model admits three boundary equilibria and up to three coexistence equilibrium points. The geometry of the nontrivial prey and predator nullclines and thus the number of coexistence equilibria primarily depends on a specific threshold of the availability of alternative food for predators. We use linear stability analysis to determine the types of hyperbolic equilibrium points and characterize the non-hyperbolic equilibrium points through normal form and center manifold theory. Change in the model parameters leading to the occurrences of a series of local bifurcations from non-hyperbolic equilibrium points, namely, transcritical, saddle-node, Hopf, cusp and Bogdanov-Takens bifurcation; there are also occurrences of global bifurcations such as homoclinic bifurcation and saddle-node bifurcation of limit cycles. We observe two interesting closed 'bubble' form induced by global bifurcations due to change in the strength of hunting cooperation and the availability of alternative food for predators. A three dimensional bifurcation diagram, concerning the original system parameters, captures how the alternation in model formulation induces gradual changes in the bifurcation scenarios. Our model highlights the stabilizing effects of group or gregarious behaviour in both prey and predator, hence supporting the predator-herbivore regulation hypothesis. Additionally, our model highlights the occurrence of "saltatory equilibria" in ecological systems and capture the dynamics observed for lion-herbivore interactions.

Coordination, cooperation, competition, crowding and congestion of molecular motors: Theoretical models and computer simulations.

Cytoskeletal motor proteins are biological nanomachines that convert chemical energy into mechanical work to carry out various functions such as cell division, cell motility, cargo transport, muscle contraction, beating of cilia and flagella, and ciliogenesis. Most of these processes are driven by the collective operation of several motors in the crowded viscous intracellular environment. Imaging and manipulation of the motors with powerful experimental probes have been complemented by mathematical analysis and computer simulations of the corresponding theoretical models. In this article, we illustrate some of the key theoretical approaches used to understand how coordination, cooperation and competition of multiple motors in the crowded intra-cellular environment drive the processes that are essential for biological function of a cell. In spite of the focus on theory, experimentalists will also find this article as an useful summary of the progress made so far in understanding multiple motor systems.

Forecasting potential invaders to prevent future biological invasions worldwide.

The ever-increasing and expanding globalisation of trade and transport underpins the escalating global problem of biological invasions. Developing biosecurity infrastructures is crucial to anticipate and prevent the transport and introduction of invasive alien species. Still, robust and defensible forecasts of potential invaders are rare, especially for species without known invasion history. Here, we aim to support decision-making by developing a quantitative invasion risk assessment tool based on invasion syndromes (i.e., generalising typical attributes of invasive alien species). We implemented a workflow based on 'Multiple Imputation with Chain Equation' to estimate invasion syndromes from imputed datasets of species' life-history and ecological traits and macroecological patterns. Importantly, our models disentangle the factors explaining (i) transport and introduction and (ii) establishment. We showcase our tool by modelling the invasion syndromes of 466 amphibians and reptile species with invasion history. Then, we project these models to amphibians and reptiles worldwide (16,236 species [c.76% global coverage]) to identify species with a risk of being unintentionally transported and introduced, and risk of establishing alien populations. Our invasion syndrome models showed high predictive accuracy with a good balance between specificity and generality. Unintentionally transported and introduced species tend to be common and thrive well in human-disturbed habitats. In contrast, those with established alien populations tend to be large-sized, are habitat generalists, thrive well in human-disturbed habitats, and have large native geographic ranges. We forecast that 160 amphibians and reptiles without known invasion history could be unintentionally transported and introduced in the future. Among them, 57 species have a high risk of establishing alien populations. Our reliable, reproducible, transferable, statistically robust and scientifically defensible quantitative invasion risk assessment tool is a significant new addition to the suite of decision-support tools needed for developing a future-proof preventative biosecurity globally.

A Mathematical Model for Fibrous Dysplasia: The Role of the Flow of Mutant Cells.

Fibrous dysplasia (FD) is a mosaic non-inheritable genetic disorder of the skeleton in which normal bone is replaced by structurally unsound fibro-osseous tissue. There is no curative treatment for FD, partly because its pathophysiology is not yet fully known. We present a simple mathematical model of the disease incorporating its basic known biology, to gain insight on the dynamics of the involved bone-cell populations, and shed light on its pathophysiology. We develop an analytical study of the model and study its basic properties. The existence and stability of steady states are studied, an analysis of sensitivity on the model parameters is done, and different numerical simulations provide findings in agreement with the analytical results. We discuss the model dynamics match with known facts on the disease, and how some open questions could be addressed using the model.

A strategy to detect metabolic changes induced by exposure to chemicals from large sets of condition-specific metabolic models computed with enumeration techniques.

BACKGROUND: The growing abundance of in vitro omics data, coupled with the necessity to reduce animal testing in the safety assessment of chemical compounds and even eliminate it in the evaluation of cosmetics, highlights the need for adequate computational methodologies. Data from omics technologies allow the exploration of a wide range of biological processes, therefore providing a better understanding of mechanisms of action (MoA) related to chemical exposure in biological systems. However, the analysis of these large datasets remains difficult due to the complexity of modulations spanning multiple biological processes. RESULTS: To address this, we propose a strategy to reduce information overload by computing, based on transcriptomics data, a comprehensive metabolic sub-network reflecting the metabolic impact of a chemical. The proposed strategy integrates transcriptomic data to a genome scale metabolic network through enumeration of condition-specific metabolic models hence translating transcriptomics data into reaction activity probabilities. Based on these results, a graph algorithm is applied to retrieve user readable sub-networks reflecting the possible metabolic MoA (mMoA) of chemicals. This strategy has been implemented as a three-step workflow. The first step consists in building cell condition-specific models reflecting the metabolic impact of each exposure condition while taking into account the diversity of possible optimal solutions with a partial enumeration algorithm. In a second step, we address the challenge of analyzing thousands of enumerated condition-specific networks by computing differentially activated reactions (DARs) between the two sets of enumerated possible condition-specific models. Finally, in the third step, DARs are grouped into clusters of functionally interconnected metabolic reactions, representing possible mMoA, using the distance-based clustering and subnetwork extraction method. The first part of the workflow was exemplified on eight molecules selected for their known human hepatotoxic outcomes associated with specific MoAs well described in the literature and for which we retrieved primary human hepatocytes transcriptomic data in Open TG-GATEs. Then, we further applied this strategy to more precisely model and visualize associated mMoA for two of these eight molecules (amiodarone and valproic acid). The approach proved to go beyond gene-based analysis by identifying mMoA when few genes are significantly differentially expressed (2 differentially expressed genes (DEGs) for amiodarone), bringing additional information from the network topology, or when very large number of genes were differentially expressed (5709 DEGs for valproic acid). In both cases, the results of our strategy well fitted evidence from the literature regarding known MoA. Beyond these confirmations, the workflow highlighted potential other unexplored mMoA. CONCLUSION: The proposed strategy allows toxicology experts to decipher which part of cellular metabolism is expected to be affected by the exposition to a given chemical. The approach originality resides in the combination of different metabolic modelling approaches (constraint based and graph modelling). The application to two model molecules shows the strong potential of the approach for interpretation and visual mining of complex omics in vitro data. The presented strategy is freely available as a python module ( https://pypi.org/project/manamodeller/ ) and jupyter notebooks ( https://github.com/LouisonF/MANA ).

Inter-Species Pharmacokinetic Modeling and Scaling for Drug Repurposing of Pyronaridine and Artesunate.

Even though several new targets (mostly viral infection) for drug repurposing of pyronaridine and artesunate have recently emerged in vitro and in vivo, inter-species pharmacokinetic (PK) data that can extend nonclinical efficacy to humans has not been reported over 30 years of usage. Since extrapolation of animal PK data to those of humans is essential to predict clinical outcomes for drug repurposing, this study aimed to investigate inter-species PK differences in three animal species (hamster, rat, and dog) and to support clinical translation of a fixed-dose combination of pyronaridine and artesunate. PK parameters (e.g., steady-state volume of distribution (Vss), clearance (CL), area under the concentration-time curve (AUC), mean residence time (MRT), etc.) of pyronaridine, artesunate, and dihydroartemisinin (an active metabolite of artesunate) were determined by non-compartmental analysis. In addition, one- or two-compartment PK modeling was performed to support inter-species scaling. The PK models appropriately described the blood concentrations of pyronaridine, artesunate, and dihydroartemisinin in all animal species, and the estimated PK parameters in three species were integrated for inter-species allometric scaling to predict human PKs. The simple allometric equation (Y = a × Wb) well explained the relationship between PK parameters and the actual body weight of animal species. The results from the study could be used as a basis for drug repurposing and support determining the effective dosage regimen for new indications based on in vitro/in vivo efficacy data and predicted human PKs in initial clinical trials.

Minimal Design of the Elephant Trunk as an Active Filament.

One of the key problems in active materials is the control of shape through actuation. A fascinating example of such control is the elephant trunk, a long, muscular, and extremely dexterous organ with multiple vital functions. The elephant trunk is an object of fascination for biologists, physicists, and children alike. Its versatility relies on the intricate interplay of multiple unique physical mechanisms and biological design principles. Here, we explore these principles using the theory of active filaments and build, theoretically, computationally, and experimentally, a minimal model that explains and accomplishes some of the spectacular features of the elephant trunk.

Poissonian Cellular Potts Models Reveal Nonequilibrium Kinetics of Cell Sorting.

Cellular Potts models are broadly applied across developmental biology and cancer research. We overcome limitations of the traditional approach, which reinterprets a modified Metropolis sampling as ad hoc dynamics, by introducing a physical timescale through Poissonian kinetics and by applying principles of stochastic thermodynamics to separate thermal and relaxation effects from athermal noise and nonconservative forces. Our method accurately describes cell-sorting dynamics in mouse-embryo development and identifies the distinct contributions of nonequilibrium processes, e.g., cell growth and active fluctuations.

Use of Individual-Based Mathematical Modelling to Understand More About Antibiotic Resistance Within-Host.

To model complex systems, individual-based models (IBMs), sometimes called "agent-based models" (ABMs), describe a simplification of the system through an adequate representation of the elements. IBMs simulate the actions and interaction of discrete individuals/agents within a system in order to discover the pattern of behavior that comes from these interactions. Examples of individuals/agents in biological systems are individual immune cells and bacteria that act independently with their own unique attributes defined by behavioral rules. In IBMs, each of these agents resides in a spatial environment and interactions are guided by predefined rules. These rules are often simple and can be easily implemented. It is expected that following the interaction guided by these rules we will have a better understanding of agent-agent interaction as well as agent-environment interaction. Stochasticity described by probability distributions must be accounted for. Events that seldom occur such as the accumulation of rare mutations can be easily modeled.Thus, IBMs are able to track the behavior of each individual/agent within the model while also obtaining information on the results of their collective behaviors. The influence of impact of one agent with another can be captured, thus allowing a full representation of both direct and indirect causation on the aggregate results. This means that important new insights can be gained and hypotheses tested.

Sample Size Determination and Study Design Impact on Dose-Scale Pharmacodynamic Bioequivalence: a Case Study Using Orlistat.

Dose-scale pharmacodynamic bioequivalence is recommended for evaluating the consistency of generic and innovator formulations of certain locally acting drugs, such as orlistat. This study aimed to investigate the standard methodology for sample size determination and the impact of study design on dose-scale pharmacodynamic bioequivalence using orlistat as the model drug. A population pharmacodynamic model of orlistat was developed using NONMEM 7.5.1 and utilized for subsequent simulations. Three different study designs were evaluated across various predefined relative bioavailability ratios of test/reference (T/R) formulations. These designs included Study Design 1 (2×1 crossover with T1 60 mg, R1 60 mg, and R2 120 mg), Study Design 2 (2×1 crossover with T2 120 mg, R1 60 mg, and R2 120 mg), and Study Design 3 (2×2 crossover with T1 60 mg, T2 120 mg, R1 60 mg, and R2 120 mg). Sample sizes were determined using a stochastic simulation and estimation approach. Under the same T/R ratio and power, Study Design 3 required the minimum sample size for bioequivalence, followed by Study Design 1, while Study Design 2 performed the worst. For Study Designs 1 and 3, a larger sample size was needed on the T/R ratio < 1.0 side for the same power compared to that on the T/R ratio > 1.0 side. The opposite asymmetry was observed for Study Design 2. We demonstrated that Study Design 3 is most effective for reducing the sample size for orlistat bioequivalence studies, and the impact of T/R ratio on sample size shows asymmetry.

Design Principles for Perfect Adaptation in Biological Networks with Nonlinear Dynamics.

Establishing a mapping between the emergent biological properties and the repository of network structures has been of great relevance in systems and synthetic biology. Adaptation is one such biological property of paramount importance that promotes regulation in the presence of environmental disturbances. This paper presents a nonlinear systems theory-driven framework to identify the design principles for perfect adaptation with respect to external disturbances of arbitrary magnitude. Based on the prior information about the network, we frame precise mathematical conditions for adaptation using nonlinear systems theory. We first deduce the mathematical conditions for perfect adaptation for constant input disturbances. Subsequently, we translate these conditions to specific necessary structural requirements for adaptation in networks of small size and then extend to argue that there exist only two classes of architectures for a network of any size that can provide local adaptation in the entire state space, namely, incoherent feed-forward (IFF) structure and negative feedback loop with buffer node (NFB). The additional positiveness constraints further narrow the admissible set of network structures. This also aids in establishing the global asymptotic stability for the steady state given a constant input disturbance. The proposed method does not assume any explicit knowledge of the underlying rate kinetics, barring some minimal assumptions. Finally, we also discuss the infeasibility of certain IFF networks in providing adaptation in the presence of downstream connections. Moreover, we propose a generic and novel algorithm based on non-linear systems theory to unravel the design principles for global adaptation. Detailed and extensive simulation studies corroborate the theoretical findings.

Identifying Bayesian optimal experiments for uncertain biochemical pathway models.

Pharmacodynamic (PD) models are mathematical models of cellular reaction networks that include drug mechanisms of action. These models are useful for studying predictive therapeutic outcomes of novel drug therapies in silico. However, PD models are known to possess significant uncertainty with respect to constituent parameter data, leading to uncertainty in the model predictions. Furthermore, experimental data to calibrate these models is often limited or unavailable for novel pathways. In this study, we present a Bayesian optimal experimental design approach for improving PD model prediction accuracy. We then apply our method using simulated experimental data to account for uncertainty in hypothetical laboratory measurements. This leads to a probabilistic prediction of drug performance and a quantitative measure of which prospective laboratory experiment will optimally reduce prediction uncertainty in the PD model. The methods proposed here provide a way forward for uncertainty quantification and guided experimental design for models of novel biological pathways.

Predicting emergent animal biodiversity patterns across multiple scales.

Restoring biodiversity-based resilience and ecosystem multi-functionality needs to be informed by more accurate predictions of animal biodiversity responses to environmental change. Ecological models make a substantial contribution to this understanding, especially when they encode the biological mechanisms and processes that give rise to emergent patterns (population, community, ecosystem properties and dynamics). Here, a distinction between 'mechanistic' and 'process-based' ecological models is established to review existing approaches. Mechanistic and process-based ecological models have made key advances to understanding the structure, function and dynamics of animal biodiversity, but are typically designed to account for specific levels of biological organisation and spatiotemporal scales. Cross-scale ecological models, which predict emergent co-occurring biodiversity patterns at interacting scales of space, time and biological organisation, is a critical next step in predictive ecology. A way forward is to first capitalise on existing models to systematically evaluate the ability of scale-explicit mechanisms and processes to predict emergent patterns at alternative scales. Such model intercomparisons will reveal mechanism to process transitions across fine to broad scales, overcome approach-specific barriers to model realism or tractability and identify gaps which necessitate the development of new fundamental principles. Key challenges surrounding model complexity and uncertainty would need to be addressed, and while opportunities from big data can streamline the integration of multiple scale-explicit biodiversity patterns, ambitious cross-scale field studies are also needed. Crucially, overcoming cross-scale ecological modelling challenges would unite disparate fields of ecology with the common goal of improving the evidence-base to safeguard biodiversity and ecosystems under novel environmental change.

Using Mathematical Optimization Models to Improve the Efficiency of Duckweeds (Wolffia arrhiza and Lemna minor) Micropropagation.

The method of plant micropropagation is widely used to obtain genetically homogeneous and infection-free plants for the needs of various industries and agriculture. Optimization of plant growth and development conditions plays a key role in economically successful micropropagation. Computer technologies have provided researchers with new approaches for modeling and a better understanding of the role of the factors involved in plant growth in vitro. To develop new models for optimizing growth conditions, we used plants with a high speed of vegetative in vitro reproduction, such as duckweed (Wolffia arrhiza and Lemna minor). Using the development of the optimal modeling of the biological processes, we have obtained the prescriptions for an individually balanced culture medium that enabled us to obtain 1.5-2.0 times more duckweed biomass with a 1.5 times higher protein concentration in the dry mass. Thus, we have demonstrated that the method of optimization modeling of the biological processes based on solving multinomial tasks from the series of quadratic equations can be used for the optimization of trophic needs of plants, specifically for micropropagation of duckweeds in vitro.

Robust trigger wave speed in Xenopus cytoplasmic extracts.

Self-regenerating trigger waves can spread rapidly through the crowded cytoplasm without diminishing in amplitude or speed, providing consistent, reliable, long-range communication. The macromolecular concentration of the cytoplasm varies in response to physiological and environmental fluctuations, raising the question of how or if trigger waves can robustly operate in the face of such fluctuations. Using Xenopus extracts, we find that mitotic and apoptotic trigger wave speeds are remarkably invariant. We derive a model that accounts for this robustness and for the eventual slowing at extremely high and low cytoplasmic concentrations. The model implies that the positive and negative effects of cytoplasmic concentration (increased reactant concentration vs. increased viscosity) are nearly precisely balanced. Accordingly, artificially maintaining a constant cytoplasmic viscosity during dilution abrogates this robustness. The robustness in trigger wave speeds may contribute to the reliability of the extremely rapid embryonic cell cycle.