Abstract:Kinetic and thermodynamic models of biological systems are commonly used to connect microscopic features to system function in a bottom-up multiscale approach. The parameters of such models-free energy differences for equilibrium properties and in general rates for equilibrium and out-of-equilibrium observables-have to be measured by different experiments or calculated from multiple computer simulations. All such parameters necessarily come with uncertainties so that when they are naively combined in a full model of the process of interest, they will generally violate fundamental statistical mechanical equalities, namely detailed balance and an equality of forward/backward rate products in cycles due to Hill. If left uncorrected, such models can produce arbitrary outputs that are physically inconsistent. Here, we develop a maximum likelihood approach (named multibind) based on the so-called potential graph to combine kinetic or thermodynamic measurements to yield state-resolved models that are thermodynamically consistent while being most consistent with the provided data and their uncertainties. We demonstrate the approach with two theoretical models, a generic two-proton binding site and a simplified model of a sodium/proton antiporter. We also describe an algorithm to use the multibind approach to solve the inverse problem of determining microscopic quantities from macroscopic measurements and, as an example, we predict the microscopic pKa values and protonation states of a small organic molecule from 1D NMR data. The multibind approach is applicable to any thermodynamic or kinetic model that describes a system as transitions between well-defined states with associated free energy differences or rates between these states. A Python package multibind, which implements the approach described here, is made publicly available under the MIT Open Source license.

Elucidating Dynamical Behaviors in Kinetically Constrained Models via Energy-Activity Double-Biased Matrix Product State Analysis

Dynamical Phase Transition in Kinetically Constrained Models with Energy-Activity Double-Bias Trajectory Ensemble

Kinematic Analysis of Planar Deployable Structures with Angulated Beams Based on Equilibrium Matrix

Dynamic Modeling Method for Constrained System with Singular Mass Matrices

Efficient and systematic calculation of arbitrary observables for the matrix product state excitation ansatz

Data2Dynamics: a Modeling Environment Tailored to Parameter Estimation in Dynamical Systems.

Dynamic Event-Triggered Asynchronous MPC of Markovian Jump Systems with Disturbances

Accurate determination of low-energy eigenspectra with multi-target matrix product states

Dynamic Structure Factor from Real Time Evolution and Exact Correction Vectors with Matrix Product States

Modeling a Hidden Dynamical System Using Energy Minimization and Kernel Density Estimates

Inherent Structure Versus Geometric Metric For State Space Discretization

Modulated Energy Estimates for Singular Kernels and their Applications to Asymptotic Analyses for Kinetic Equations

Dynamic Canonical and Microcanonical Transition Matrix Analyses of Critical Behavior

Finite-size Scaling in Kinetics of Phase Separation in Certain Models of Aligning Active Particles

Free Energy Landscapes, Diffusion Coefficients, and Kinetic Rates from Transition Paths

A Computational Approach to Multistationarity of Power-Law Kinetic Systems

A Kinetically Constrained Model with a Thermodynamic Flavor

Energy-based dual-phase dynamics identification of clearance nonlinearities

Computing Long Timescale Biomolecular Dynamics using Quasi-Stationary Distribution Kinetic Monte Carlo (QSD-KMC)

An energy-momentum method for ordinary differential equations with an underlying $k$-polysymplectic manifold

Thermodynamically consistent determination of free energies and rates in kinetic cycle models