Kinetic control of stationary flux ratios for a wide range of biochemical processes

Mutual Linearity of Nonequilibrium Network Currents

For continuous-time Markov chains and open unimolecular chemical reaction networks, we prove that any two stationary currents are linearly related upon perturbations of a single edge's transition rates, arbitrarily far from equilibrium. We extend the result to nonstationary currents in the frequency domain, provide and discuss an explicit expression for the current-current susceptibility in terms of the network topology, and discuss possible generalizations. In practical scenarios, the mutual linearity relation has predictive power and can be used as a tool for inference or model proof testing.

Steady-state probabilities for Markov jump processes in terms of powers of the transition rate matrix

Several types of dynamics at stationarity can be described in terms of a Markov jump process among a finite number N of representative sites. Before dealing with the dynamical aspects, one basic problem consists in expressing the a priori steady-state occupation probabilities of the sites. In particular, one wishes to go beyond the mere black-box computational tools and find expressions in which the jump rate constants appear explicitly, therefore allowing for a potential design/control of the network. For strongly connected networks admitting a unique stationary state with all sites populated, here we express the occupation probabilities in terms of a formula that involves powers of the transition rate matrix up to order N - 1. We also provide an expression of the derivatives with respect to the jump rate constants, possibly useful in sensitivity analysis frameworks. Although we refer to dynamics in (bio)chemical networks at thermal equilibrium or under nonequilibrium steady-state conditions, the results are valid for any Markov jump process under the same assumptions.

Insights into Error Control Mechanisms in Biological Processes: Copolymerization and Enzyme-Kinetics Revisited

The high fidelity observed in biological information processing ranging from replication to translation has stimulated significant research efforts to clarify the underlying microscopic picture. Theoretically, several approaches to analyze the error rates have been proposed. The copolymerization theory describes the addition and removal of monomers at the growing tip of a copolymer, leading to a closed set of nonlinear equations. On the other hand, enzyme-kinetics approaches formulate linear equations of biochemical networks, describing transitions between discrete chemical states. However, it is still unclear whether the error values computed by the two approaches agree. Moreover, there are conflicting interpretations on whether the error is under thermodynamic or kinetic discrimination control. In this work, we examine the error rate in persistent copying biochemical processes by specifically analyzing both theoretical approaches. The initial disagreement of the results between the two theories motivated us to rederive the formula for the error rate in the kinetic model. The error computed with the new method resulted in excellent agreement between both theoretical approaches and with Monte Carlo simulations. Furthermore, our theoretical analysis shows that the kinetic discrimination controls the error, even when the energy difference between adding the right and wrong products is very small. Our theoretical investigation gives important insights into the physical-chemical properties of complex biological processes by providing the quantitative framework to evaluate them.

Flexibility and sensitivity in gene regulation out of equilibrium

Cells adapt to environments and tune gene expression by controlling the concentrations of proteins and their kinetics in regulatory networks. In both eukaryotes and prokaryotes, experiments and theory increasingly attest that these networks can and do consume bio-chemical energy. How does this dissipation enable cellular behaviors unobtainable in equilibrium? This open question demands quantitative models that transcend thermodynamic equilibrium. Here we study the control of a simple, ubiquitous gene regulatory motif to explore the consequences of departing equilibrium in kinetic cycles. Employing graph theory, we find that dissipation unlocks nonmonotonicity and enhanced sensitivity of gene expression with respect to a transcription factor's concentration. These features allow a single transcription factor to act as both a repressor and activator at different levels or achieve outputs with multiple concentration regions of locally-enhanced sensitivity. We systematically dissect how energetically-driving individual transitions within regulatory networks, or pairs of transitions, generates more adjustable and sensitive phenotypic responses. Our findings quantify necessary conditions and detectable consequences of energy expenditure. These richer mathematical behaviors-feasibly accessed using biological energy budgets and rates-may empower cells to accomplish sophisticated regulation with simpler architectures than those required at equilibrium.

Significance Statement

Growing theoretical and experimental evidence demonstrates that cells can (and do) spend biochemical energy while regulating their genes. Here we explore the impact of departing from equilibrium in simple regulatory cycles, and learn that beyond increasing sensitivity, dissipation can unlock more flexible input-output behaviors that are otherwise forbidden without spending energy. These more complex behaviors could enable cells to perform more sophisticated functions using simpler systems than those needed at equilibrium.

The energy cost and optimal design of networks for biological discrimination

Many biological processes discriminate between correct and incorrect substrates through the kinetic proofreading mechanism that enables lower error at the cost of higher energy dissipation. Elucidating physico-chemical constraints for global minimization of dissipation and error is important for understanding enzyme evolution. Here, we identify theoretically a fundamental error-cost bound that tightly constrains the performance of proofreading networks under any parameter variations preserving the rate discrimination between substrates. The bound is kinetically controlled, i.e. completely determined by the difference between the transition state energies on the underlying free energy landscape. The importance of the bound is analysed for three biological processes. DNA replication by T7 DNA polymerase is shown to be nearly optimized, i.e. its kinetic parameters place it in the immediate proximity of the error-cost bound. The isoleucyl-tRNA synthetase (IleRS) of E. coli also operates close to the bound, but further optimization is prevented by the need for reaction speed. In contrast, E. coli ribosome operates in a high-dissipation regime, potentially in order to speed up protein production. Together, these findings establish a fundamental error-dissipation relation in biological proofreading networks and provide a theoretical framework for studying error-dissipation trade-off in other systems with biological discrimination.

Multiscale kinetic analysis of proteins

The stochasticity of molecular motion results in the existence of multiple kinetically relevant pathways in many biomolecular mechanisms. Because it is highly demanding to characterize them for complex systems, mechanisms are often described with a single-pathway perspective. However, kinetic network analysis and sub-ensemble experimental insight are increasingly demonstrating not only the existence of competing pathways but also the importance of kinetic selection in biology. This review focuses on advances in multiscale kinetic analysis of proteins, which connects molecular level information from simulations to macroscopic data to characterize mechanistic reaction networks and the reactive flux through them. We describe a range of methods used and highlight several examples where kinetic modeling has revealed functional importance of pathway heterogeneity.

Theoretical Analysis Reveals the Cost and Benefit of Proofreading in Coronavirus Genome Replication

Severe acute respiratory syndrome coronaviruses have unusually large RNA genomes replicated by a multiprotein complex containing an RNA-dependent RNA polymerase (RdRp). Exonuclease activity enables the RdRp complex to remove wrongly incorporated bases via proofreading, a process not utilized by other RNA viruses. However, it is unclear why the RdRp complex needs proofreading and what the associated trade-offs are. Here we investigate the interplay among the accuracy, speed, and energetic cost of proofreading in the RdRp complex using a kinetic model and bioinformatics analysis. We find that proofreading nearly optimizes the rate of functional virus production. However, we find that further optimization would lead to a significant increase in the proofreading cost. Unexpected importance of the cost minimization is further supported by other global analyses. We speculate that cost optimization could help avoid cell defense responses. Thus, proofreading is essential for the production of functional viruses, but its rate is limited by energy costs.

Do We Understand the Mechanisms Used by Biological Systems to Correct Their Errors?

Most cellular processes involved in biological information processing display a surprisingly low error rate despite the stochasticity of the underlying biochemical reactions and the presence of competing chemical species. Such high fidelity is the result of nonequilibrium kinetic proofreading mechanisms, i.e., the existence of dissipative pathways for correcting the reactions that went in the wrong direction. While proofreading was often studied from the perspective of error minimization, a number of recent studies have demonstrated that the underlying mechanisms need to consider the interplay of other characteristic properties such as speed, energy dissipation, and noise reduction. Here, we present current views and new insights on the mechanisms of error-correction phenomena and various trade-off scenarios in the optimization of the functionality of biological systems. Existing challenges and future directions are also discussed.