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.

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.

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.

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.

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.

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 and experiments. Such models can never be entirely accurate; it is nevertheless important that the models are compatible with physical principles.

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.

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.

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.

Bond Graph Representation of Chemical Reaction Networks.

The Bond Graph approach and the Chemical Reaction Network approach for modeling the biomolecular systems were 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.

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.

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.

Computing Biomolecular System Steady-States.

A new approach to compute the equilibria and the steady-states of biomolecular systems modeled 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.

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.

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.

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.

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.

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.

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.

Intermittent control models of human standing: similarities and differences.

Two architectures of intermittent control are compared and contrasted in the context of the single inverted pendulum model often used for describing standing in humans. The architectures are similar insofar as they use periods of open-loop control punctuated by switching events when crossing a switching surface to keep the system state trajectories close to trajectories leading to equilibrium. The architectures differ in two significant ways.

Human stick balancing: an intermittent control explanation.

There are two issues in balancing a stick pivoting on a finger tip (or mechanically on a moving cart): maintaining the stick angle near to vertical and maintaining the horizontal position within the bounds of reach or cart track. The (linearised) dynamics of the angle are second order (although driven by pivot acceleration), and so, as in human standing, control of the angle is not, by itself very difficult. However, once the angle is under control, the position dynamics are, in general, fourth order.

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.

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.

Frequency-domain identification of the human controller.

System identification techniques applied to experimental human-in-the-loop data provide an objective test of three alternative control-theoretical models of the human control system: non-predictive control, predictive control, and intermittent predictive control. A two-stage approach to the identification of a single-input single-output control system is used: first, the closed-loop frequency response is derived using the periodic property of the experimental data, followed by the fitting of a parametric model. While this approach is well-established for non-predictive and predictive control, it is here used for the first time with intermittent predictive control.