Search

Interactions between biomolecules, electrons and protons are essential to many fundamental processes sustaining life. It is therefore of interest to build mathematical models of these bioelectrical processes not only to enhance understanding but also to enable computer models to complement in vitro and in vivo experiments.Such models can never be entirely accurate; it is nevertheless important that the models are compatible with physical principles. Network Thermodynamics, as implemented with bond graphs, provide one approach to creating physically compatible mathematical models of bioelectrical systems. This is illustrated using simple models of ion channels, redox reactions, proton pumps and electrogenic membrane transporters thus demonstrating that the approach can be used to build mathematical and computer models of a wide range of bioelectrical systems.

Whole-cell modelling is constrained by the laws of nature in general and the laws of thermodynamics in particular. This paper shows how one prolific source of information, stoichiometric models of biomolecular systems, can be integrated with thermodynamic principles using the bond graph approach to network thermodynamics.

Energy-based modelling brings engineering insight to the understanding of biomolecular systems. It is shown how well-established control engineering concepts, such as loop-gain, arise from energy feedback loops and are therefore amenable to control engineering insight. In particular, a novel method is introduced to allow the transfer function based approach of classical linear control to be utilised in the analysis of feedback systems modelled by network thermodynamics and thus amalgamate energy-based modelling with control systems analysis. The approach is illustrated using a class of metabolic cycles with activation and inhibition leading the concept of Cyclic Flow Modulation.

BondGraphTools is a Python library for scripted modelling of complex multi-physics systems. In contrast to existing modelling solutions, BondGraphTools is based upon the well established bond graph methodology, provides a programming interface for symbolic model composition, and is intended to be used in conjunction with the existing scientific Python toolchain. Here we discuss the design, implementation and use of BondGraphTools, demonstrate how it can be used to accelerate systems modelling with an example from optomechanics, and comment on current and future applications in cross-domain modelling, particularly in systems biology.

Systems biology and whole-cell modelling are demanding increasingly comprehensive mathematical models of cellular biochemistry. These models require the development of simplified models of specific processes which capture essential biophysical features but without unnecessarily complexity. Recently there has been renewed interest in thermodynamically-based modelling of cellular processes. Here we present an approach to developing of simplified yet thermodynamically consistent (hence physically plausible) models which can readily be incorporated into large scale biochemical descriptions but which do not require full mechanistic detail of the underlying processes. We illustrate the approach through development of a simplified, physically plausible model of the mitochondrial electron transport chain and show that the simplified model behaves like the full system.

The Bond Graph approach and the Chemical Reaction Network approach to modelling biomolecular systems developed independently. This paper brings together the two approaches by providing a bond graph interpretation of the chemical reaction network concept of complexes. Both closed and open systems are discussed. The method is illustrated using a simple enzyme-catalysed reaction and a trans-membrane transporter.

Membrane transporters contribute to the regulation of the internal environment of cells by translocating substrates across cell membranes. Like all physical systems, the behaviour of membrane transporters is constrained by the laws of thermodynamics. However, many mathematical models of transporters, especially those incorporated into whole-cell models, are not thermodynamically consistent, leading to unrealistic behaviour. In this paper we use a physics-based modelling framework, in which the transfer of energy is explicitly accounted for, to develop thermodynamically consistent models of transporters. We then apply this methodology to model two specific transporters: the cardiac sarcoplasmic/endoplasmic Ca$^{2+}$ ATPase (SERCA) and the cardiac Na$^+$/K$^+$ ATPase.

Efficient energy transduction is one driver of evolution; and thus understanding biomolecular energy transduction is crucial to understanding living organisms. As an energy-orientated modelling methodology, bond graphs provide a useful approach to describing and modelling the efficiency of living systems. This paper gives some new results on the efficiency of metabolism based on bond graph models of the key metabolic processes: glycolysis., Comment: Accepted for 13th International Conference on Bond Graph Modeling (ICBGM 18), July 9-12, 2018, Bordeaux, France

Mathematical models of cardiac action potentials have become increasingly important in the study of heart disease and pharmacology, but concerns linger over their robustness during long periods of simulation, in particular due to issues such as model drift and non-unique steady states. Previous studies have linked these to violation of conservation laws, but only explored those issues with respect to charge conservation in specific models. Here, we propose a general and systematic method of identifying conservation laws hidden in models of cardiac electrophysiology by using bond graphs, and develop a bond graph model of the cardiac action potential to study long-term behaviour. Bond graphs provide an explicit energy-based framework for modelling physical systems, which makes them well-suited for examining conservation within electrophysiological models. We find that the charge conservation laws derived in previous studies are examples of the more general concept of a "conserved moiety". Conserved moieties explain model drift and non-unique steady states, generalising the results from previous studies. The bond graph approach provides a rigorous method to check for drift and non-unique steady states in a wide range of cardiac action potential models, and can be extended to examine behaviours of other excitable systems.

The Na$^+$/K$^+$ ATPase is an essential component of cardiac electrophysiology, maintaining physiological Na$^+$ and K$^+$ concentrations over successive heart beats. Terkildsen et al. (2007) developed a model of the ventricular myocyte Na$^+$/K$^+$ ATPase to study extracellular potassium accumulation during ischaemia, demonstrating the ability to recapitulate a wide range of experimental data, but unfortunately there was no archived code associated with the original manuscript. Here we detail an updated version of the model and provide CellML and MATLAB code to ensure reproducibility and reusability. We note some errors within the original formulation which have been corrected to ensure that the model is thermodynamically consistent, and although this required some reparameterisation, the resulting model still provides a good fit to experimental measurements that demonstrate the dependence of Na$^+$/K$^+$ ATPase pumping rate upon membrane voltage and metabolite concentrations. To demonstrate thermodynamic consistency we also developed a bond graph version of the model. We hope that these models will be useful for community efforts to assemble a whole-cell cardiomyocyte model which facilitates the investigation of cellular energetics.

A new approach to computing the equilibria and steady-states of biomolecular systems modelled by bond graphs is presented. The approach is illustrated using a model of a biomolecular cycle representing a membrane transporter and a model of the mitochondrial electron transport chain., Comment: Accepted by IEEE Transactions on NanoBioscience

The sensitivity properties of intermittent control are analysed and the conditions for a limit cycle derived theoretically and verified by simulation.

Engineering systems modelling and analysis based on the bond graph approach has been applied to biomolecular systems. In this context, the notion of a Faraday-equivalent chemical potential is introduced which allows chemical potential to be expressed in an analogous manner to electrical volts thus allowing engineering intuition to be applied to biomolecular systems. Redox reactions, and their representation by half-reactions, are key components of biological systems which involve both electrical and chemical domains. A bond graph interpretation of redox reactions is given which combines bond graphs with the Faraday- equivalent chemical potential. This approach is particularly relevant when the biomolecular system implements chemoelectrical transduction - for example chemiosmosis within the key metabolic pathway of mitochondria: oxidative phosphorylation. An alternative way of implementing computational modularity using bond graphs is introduced and used to give a physically based model of the mitochondrial electron transport chain. To illustrate the overall approach, this model is analysed using the Faraday-equivalent chemical potential approach and engineering intuition is used to guide affinity equalisation: a energy based analysis of the mitochondrial electron transport chain.

Decomposition of biomolecular reaction networks into pathways is a powerful approach to the analysis of metabolic and signalling networks. Current approaches based on analysis of the stoichiometric matrix reveal information about steady-state mass flows (reaction rates) through the network. In this work we show how pathway analysis of biomolecular networks can be extended using an energy-based approach to provide information about energy flows through the network. This energy-based approach is developed using the engineering-inspired bond graph methodology to represent biomolecular reaction networks. The approach is introduced using glycolysis as an exemplar; and is then applied to analyse the efficiency of free energy transduction in a biomolecular cycle model of a transporter protein (Sodium-Glucose Transport Protein 1, SGLT1). The overall aim of our work is to present a framework for modelling and analysis of biomolecular reactions and processes which considers energy flows and losses as well as mass transport.

Energy-based bond graph modelling of biomolecular systems is extended to include chemoelectrical trans- duction thus enabling integrated thermodynamically-compliant modelling of chemoelectrical systems in general and excitable membranes in particular. Our general approach is illustrated by recreating a well-known model of an excitable membrane. This model is used to investigate the energy consumed during a membrane action potential thus contributing to the current debate on the trade-off between the speed of an action potential event and energy consumption. The influx of Na+ is often taken as a proxy for energy consumption; in contrast, this paper presents an energy based model of action potentials. As the energy based approach avoids the assumptions underlying the proxy approach it can be directly used to compute energy consumption in both healthy and diseased neurons. These results are illustrated by comparing the energy consumption of healthy and degenerative retinal ganglion cells using both simulated and in vitro data.

Bond graphs can be used to build thermodynamically-compliant hierarchical models of biomolecular systems. As bond graphs have been widely used to model, analyse and synthesise engineering systems, this paper suggests that they can play the same role in the modelling, analysis and synthesis of biomolecular systems. The particular structure of bond graphs arising from biomolecular systems is established and used to elucidate the relation between thermodynamically closed and open systems. Block diagram representations of the dynamics implied by these bond graphs are used to reveal implicit feedback structures and are linearised to allow the application of control-theoretical methods. Two concepts of modularity are examined: computational modularity where physical correctness is retained and behavioural modularity where module behaviour (such as ultrasensitivity) is retained. As well as providing computational modularity, bond graphs provide a natural formulation of behavioural modularity and reveal the sources of retroactivity. A bond graph approach to reducing retroactivity, and thus inter-module interaction, is shown to require a power supply such as that provided by the ATP = ADP + Pi reaction. The MAPK cascade (Raf-MEK-ERK pathway) is used as an illustrative example.

The bond graph approach to modelling biochemical networks is extended to allow hierarchical construction of complex models from simpler components. This is made possible by representing the simpler components as thermodynamically open systems exchanging mass and energy via ports. A key feature of this approach is that the resultant models are robustly thermodynamically compliant: the thermodynamic compliance is not dependent on precise numerical values of parameters. Moreover, the models are reusable due to the well-defined interface provided by the energy ports. To extract bond graph model parameters from parameters found in the literature, general and compact formulae are developed to relate free-energy constants and equilibrium constants. The existence and uniqueness of solutions is considered in terms of fundamental properties of stoichiometric matrices. The approach is illustrated by building a hierarchical bond graph model of glycogenolysis in skeletal muscle., Comment: 35 pages, 12 figures

Intermittent control has a long history in the physiological literature and there is strong experimental evidence that some human control systems are intermittent. Intermittent control has also appeared in various forms in the engineering literature. This article discusses a particular mathematical model of Event-driven Intermittent Control which brings together engineering and physiological insights and builds on and extends previous work in this area. Illustrative examples of the properties of Intermittent Control in a physiological context are given together with suggestions for future research directions in both physiology and engineering.

Thermodynamic aspects of chemical reactions have a long history in the Physical Chemistry literature. In particular, biochemical cycles - the building-blocks of biochemical systems - require a source of energy to function. However, although fundamental, the role of chemical potential and Gibb's free energy in the analysis of biochemical systems is often overlooked leading to models which are physically impossible. The bond graph approach was developed for modelling engineering systems where energy generation, storage and transmission are fundamental. The method focuses on how power flows between components and how energy is stored, transmitted or dissipated within components. Based on early ideas of network thermodynamics, we have applied this approach to biochemical systems to generate models which automatically obey the laws of thermodynamics. We illustrate the method with examples of biochemical cycles. We have found that thermodynamically compliant models of simple biochemical cycles can easily be developed using this approach. In particular, both stoichiometric information and simulation models can be developed directly from the bond graph. Furthermore, model reduction and approximation while retaining structural and thermodynamic properties is facilitated. Because the bond graph approach is also modular and scaleable, we believe that it provides a secure foundation for building thermodynamically compliant models of large biochemical networks.