Sensitivity analysis of biochemical systems using bond graphs.

The sensitivity of systems biology models to parameter variation can give insights into which parameters are most important for physiological function, and also direct efforts to estimate parameters. However, in general, kinetic models of biochemical systems do not remain thermodynamically consistent after perturbing parameters. To address this issue, we analyse the sensitivity of biological reaction networks in the context of a bond graph representation. We find that the parameter sensitivities can themselves be represented as bond graph components, mirroring potential mechanisms for controlling biochemistry. In particular, a sensitivity system is derived which re-expresses parameter variation as additional system inputs. The sensitivity system is then linearized with respect to these new inputs to derive a linear system which can be used to give local sensitivity to parameters in terms of linear system properties such as gain and time constant. This linear system can also be used to find so-called sloppy parameters in biological models. We verify our approach using a model of the Pentose Phosphate Pathway, confirming the reactions and metabolites most essential to maintaining the function of the pathway.

Network thermodynamics of biological systems: A bond graph approach.

Edmund Crampin (1973-2021) was at the forefront of Systems Biology research and his work will influence the field for years to come. This paper brings together and summarises the seminal work of his group in applying energy-based bond graph methods to biological systems. In particular, this paper: (a) motivates the need to consider energy in modelling biology; (b) introduces bond graphs as a methodology for achieving this; (c) describes extensions to modelling electrochemical transduction; (d) outlines how bond graph models can be constructed in a modular manner and (e) describes stoichiometric approaches to deriving fundamental properties of reaction networks. These concepts are illustrated using a new bond graph model of photosynthesis in chloroplasts.

Is Intermittent Control the Source of the Non-Linear Oscillatory Component (0.2-2Hz) in Human Balance Control?

OBJECTIVE: To explain the 0.2-2Hz oscillation in human balance. MOTIVATION: Oscillation (0.2-2 Hz) in the control signal (ankle moment) is sustained independently of external disturbances and exaggerated in Parkinson's disease. Does resonance or limit cycles in the neurophysiological feedback loop cause this oscillation? We investigate two linear (non-predictive, predictive) and one non-linear (intermittent-predictive) control model (NPC, PC, IPC). METHODS: Fourteen healthy participants, strapped to an actuated single segment robot with dynamics of upright standing, used natural haptic-visual feedback and myoelectric control signals from lower leg muscles to maintain balance. An input disturbance applied stepwise changes in external force. A linear time invariant model (ARX) extracted the delayed component of the control signal related linearly to the disturbance, leaving the remaining, larger, oscillatory non-linear component. We optimized model parameters and noise (observation, motor) to replicate concurrently (i) estimated-delay, (ii) time-series of the linear component, and (iii) magnitude-frequency spectrum and transient magnitude response of the non-linear component. Results (mean±S.D., p<0.05): NPC produced estimated delays (0.116±0.03s) significantly lower than experiment (0.145±0.04s). Overall fit (i)-(iii) was (79±7%, 83±7%, 84±6% for NPC, PC, IPC). IPC required little or no noise. Mean frequency of experimental oscillation (0.99±0.16 Hz) correlated trial by trial with closed loop resonant frequency (fres), not limit cycles, nor sampling rate. NPC (0.36±0.08Hz) and PC (0.86±0.4Hz) showed fres significantly lower than IPC (0.98±0.2Hz). CONCLUSION: Human balance control requires short-term prediction. SIGNIFICANCE: IPC mechanisms (prediction error, threshold related sampling, sequential re-initialization of open-loop predictive control) explain resonant gain without uncontrolled oscillation for healthy balance.

Analysing and simulating energy-based models in biology using BondGraphTools.

Like all physical systems, biological systems are constrained by the laws of physics. However, mathematical models of biochemistry frequently neglect the conservation of energy, leading to unrealistic behaviour. Energy-based models that are consistent with conservation of mass, charge and energy have the potential to aid the understanding of complex interactions between biological components, and are becoming easier to develop with recent advances in experimental measurements and databases. In this paper, we motivate the use of bond graphs (a modelling tool from engineering) for energy-based modelling and introduce, BondGraphTools, a Python library for constructing and analysing bond graph models. We use examples from biochemistry to illustrate how BondGraphTools can be used to automate model construction in systems biology while maintaining consistency with the laws of physics.

Modular assembly of dynamic models in systems biology.

It is widely acknowledged that the construction of large-scale dynamic models in systems biology requires complex modelling problems to be broken up into more manageable pieces. To this end, both modelling and software frameworks are required to enable modular modelling. While there has been consistent progress in the development of software tools to enhance model reusability, there has been a relative lack of consideration for how underlying biophysical principles can be applied to this space. Bond graphs combine the aspects of both modularity and physics-based modelling. In this paper, we argue that bond graphs are compatible with recent developments in modularity and abstraction in systems biology, and are thus a desirable framework for constructing large-scale models. We use two examples to illustrate the utility of bond graphs in this context: a model of a mitogen-activated protein kinase (MAPK) cascade to illustrate the reusability of modules and a model of glycolysis to illustrate the ability to modify the model granularity.

Network Thermodynamical Modeling of Bioelectrical Systems: A Bond Graph Approach.

Interactions among 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.

Modular dynamic biomolecular modelling with bond graphs: the unification of stoichiometry, thermodynamics, kinetics and data.

Renewed interest in dynamic simulation models of biomolecular systems has arisen from advances in genome-wide measurement and applications of such models in biotechnology and synthetic biology. In particular, genome-scale models of cellular metabolism beyond the steady state are required in order to represent transient and dynamic regulatory properties of the system. Development of such whole-cell models requires new modelling approaches. Here, we propose the energy-based bond graph methodology, which integrates stoichiometric models with thermodynamic principles and kinetic modelling. We demonstrate how the bond graph approach intrinsically enforces thermodynamic constraints, provides a modular approach to modelling, and gives a basis for estimation of model parameters leading to dynamic models of biomolecular systems. The approach is illustrated using a well-established stoichiometric model of Escherichia coli and published experimental data.

Energy-Based Modeling of the Feedback Control of Biomolecular Systems With Cyclic Flow Modulation.

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 to the concept of Cyclic Flow Modulation.

Physically-plausible modelling of biomolecular systems: A simplified, energy-based model of the mitochondrial electron transport chain.

Advances in systems biology and whole-cell modelling demand increasingly comprehensive mathematical models of cellular biochemistry. Such models require the development of simplified representations 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.

A thermodynamic framework for modelling membrane transporters.

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 Ca2+ ATPase (SERCA) and the cardiac Na+/K+ ATPase.

Bond graph modelling of the cardiac action potential: implications for drift and non-unique steady states.

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, generalizing 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.

Visuo-manual tracking: does intermittent control with aperiodic sampling explain linear power and non-linear remnant without sensorimotor noise?

KEY POINTS: A human controlling an external system is described most easily and conventionally as linearly and continuously translating sensory input to motor output, with the inevitable output remnant, non-linearly related to the input, attributed to sensorimotor noise. Recent experiments show sustained manual tracking involves repeated refractoriness (insensitivity to sensory information for a certain duration), with the temporary 200-500 ms periods of irresponsiveness to sensory input making the control process intrinsically non-linear. This evidence calls for re-examination of the extent to which random sensorimotor noise is required to explain the non-linear remnant. This investigation of manual tracking shows how the full motor output (linear component and remnant) can be explained mechanistically by aperiodic sampling triggered by prediction error thresholds. Whereas broadband physiological noise is general to all processes, aperiodic sampling is associated with sensorimotor decision making within specific frontal, striatal and parietal networks; we conclude that manual tracking utilises such slow serial decision making pathways up to several times per second. ABSTRACT: The human operator is described adequately by linear translation of sensory input to motor output. Motor output also always includes a non-linear remnant resulting from random sensorimotor noise from multiple sources, and non-linear input transformations, for example thresholds or refractory periods. Recent evidence showed that manual tracking incurs substantial, serial, refractoriness (insensitivity to sensory information of 350 and 550 ms for 1st and 2nd order systems respectively). Our two questions are: (i) What are the comparative merits of explaining the non-linear remnant using noise or non-linear transformations? (ii) Can non-linear transformations represent serial motor decision making within the sensorimotor feedback loop intrinsic to tracking? Twelve participants (instructed to act in three prescribed ways) manually controlled two systems (1st and 2nd order) subject to a periodic multi-sine disturbance. Joystick power was analysed using three models, continuous-linear-control (CC), continuous-linear-control with calculated noise spectrum (CCN), and intermittent control with aperiodic sampling triggered by prediction error thresholds (IC). Unlike the linear mechanism, the intermittent control mechanism explained the majority of total power (linear and remnant) (77-87% vs. 8-48%, IC vs. CC). Between conditions, IC used thresholds and distributions of open loop intervals consistent with, respectively, instructions and previous measured, model independent values; whereas CCN required changes in noise spectrum deviating from broadband, signal dependent noise. We conclude that manual tracking uses open loop predictive control with aperiodic sampling. Because aperiodic sampling is inherent to serial decision making within previously identified, specific frontal, striatal and parietal networks we suggest that these structures are intimately involved in visuo-manual tracking.

Energy-based analysis of biomolecular pathways.

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.

Bond Graph Modeling of Chemiosmotic Biomolecular Energy Transduction.

Engineering systems modeling 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 analyzed 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.

Modular bond-graph modelling and analysis of biomolecular systems.

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 study suggests that they can play the same rôle 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 mitogen-activated protein kinase cascade (Raf-MEK-ERK pathway) is used as an illustrative example.

Energy-based analysis of biochemical cycles using bond graphs.

Thermodynamic aspects of chemical reactions have a long history in the physical chemistry literature. In particular, biochemical cycles 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 the 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.

Does the motor system need intermittent control?

Explanation of motor control is dominated by continuous neurophysiological pathways (e.g., transcortical, spinal) and the continuous control paradigm. Using new theoretical development, methodology, and evidence, we propose intermittent control, which incorporates a serial ballistic process within the main feedback loop, provides a more general and more accurate paradigm necessary to explain attributes highly advantageous for competitive survival and performance.

Interfacing sensory input with motor output: does the control architecture converge to a serial process along a single channel?

Modular organization in control architecture may underlie the versatility of human motor control; but the nature of the interface relating sensory input through task-selection in the space of performance variables to control actions in the space of the elemental variables is currently unknown. Our central question is whether the control architecture converges to a serial process along a single channel? In discrete reaction time experiments, psychologists have firmly associated a serial single channel hypothesis with refractoriness and response selection [psychological refractory period (PRP)]. Recently, we developed a methodology and evidence identifying refractoriness in sustained control of an external single degree-of-freedom system. We hypothesize that multi-segmental whole-body control also shows refractoriness. Eight participants controlled their whole body to ensure a head marker tracked a target as fast and accurately as possible. Analysis showed enhanced delays in response to stimuli with close temporal proximity to the preceding stimulus. Consistent with our preceding work, this evidence is incompatible with control as a linear time invariant process. This evidence is consistent with a single-channel serial ballistic process within the intermittent control paradigm with an intermittent interval of around 0.5 s. A control architecture reproducing intentional human movement control must reproduce refractoriness. Intermittent control is designed to provide computational time for an online optimization process and is appropriate for flexible adaptive control. For human motor control we suggest that parallel sensory input converges to a serial, single channel process involving planning, selection, and temporal inhibition of alternative responses prior to low dimensional motor output. Such design could aid robots to reproduce the flexibility of human control.

Refractoriness in sustained visuo-manual control: is the refractory duration intrinsic or does it depend on external system properties?

Researchers have previously adopted the double stimulus paradigm to study refractoriness in human neuromotor control. Currently, refractoriness, such as the Psychological Refractory Period (PRP) has only been quantified in discrete movement conditions. Whether refractoriness and the associated serial ballistic hypothesis generalises to sustained control tasks has remained open for more than sixty years. Recently, a method of analysis has been presented that quantifies refractoriness in sustained control tasks and discriminates intermittent (serial ballistic) from continuous control. Following our recent demonstration that continuous control of an unstable second order system (i.e. balancing a 'virtual' inverted pendulum through a joystick interface) is unnecessary, we ask whether refractoriness of substantial duration (~200 ms) is evident in sustained visual-manual control of external systems. We ask whether the refractory duration (i) is physiologically intrinsic, (ii) depends upon system properties like the order (0, 1(st), and 2(nd)) or stability, (iii) depends upon target jump direction (reversal, same direction). Thirteen participants used discrete movements (zero order system) as well as more sustained control activity (1(st) and 2(nd) order systems) to track unpredictable step-sequence targets. Results show a substantial refractory duration that depends upon system order (250, 350 and 550 ms for 0, 1(st) and 2(nd) order respectively, n=13, p<0.05), but not stability. In sustained control refractoriness was only found when the target reverses direction. In the presence of time varying actuators, systems and constraints, we propose that central refractoriness is an appropriate control mechanism for accommodating online optimization delays within the neural circuitry including the more variable processing times of higher order (complex) input-output relations.