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 featur
es 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 conc
ept 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 de
scribing 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.
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.
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 fl
ows (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.
