Accuracy of Substrate Selection by Enzymes Is Controlled by Kinetic Discrimination

Kinetic Mechanism for Fidelity of CRISPR-Cas9 Variants

CRISPR-Cas9 is a nuclease creating DNA breaks at sites with sufficient complementarity to the RNA guide. Notably, Cas9 does not require exact RNA-DNA complementarity and can cleave off-target sequences. Various high-accuracy Cas9 variants have been developed, but the precise mechanism of how these variants achieve higher accuracy remains unclear. Here, we develop a kinetic model of Cas9 substrate selection and cleavage parametrized by data from the literature, including single-molecule Förster resonance energy transfer (FRET) measurements. Based on observed FRET transition statistics, we predict that the Cas9 substrate recognition and cleavage mechanism must allow for HNH domain transitions independent of substrate binding. Additionally, we show that the enhancement in Cas9 substrate specificity must be due to changes in kinetics rather than changes in substrate binding affinities. Finally, we use our model to identify kinetic parameters for HNH domain transitions that can be perturbed to enable high-accuracy cleavage while maintaining cleavage speeds.

Bayesian Inference of Binding Kinetics from Fluorescence Time Series

The study of binding kinetics via the analysis of fluorescence time traces is often confounded by measurement noise and photophysics. Although photoblinking can be mitigated by using labels less likely to photoswitch, photobleaching generally cannot be eliminated. Current methods for measuring binding and unbinding rates are, therefore, limited by concurrent photobleaching events. Here, we propose a method to infer binding and unbinding rates alongside photobleaching rates using fluorescence intensity traces. Our approach is a two-stage process involving analyzing individual regions of interest (ROIs) with a hidden Markov model to infer the fluorescence intensity levels of each trace. We then use the inferred intensity level state trajectory from all of the ROIs to infer kinetic rates. Our method has several advantages, including the ability to analyze noisy traces, account for the presence of photobleaching events, and provide uncertainties associated with the inferred binding kinetics. We demonstrate the effectiveness and reliability of our method through simulations and data from DNA origami binding experiments.

Bayesian Inference of Binding Kinetics from Fluorescence Time Series

The study of binding kinetics via the analysis of fluorescence time traces is often confounded by measurement noise and photophysics. Although photoblinking can be mitigated by using labels less likely to photoswitch, photobleaching generally cannot be eliminated. Current methods for measuring binding and unbinding rates are therefore limited by concurrent photobleaching events. Here, we propose a method to infer binding and unbinding rates alongside photobleaching rates using fluorescence intensity traces. Our approach is a two-stage process involving analyzing individual regions of interest (ROIs) with a Hidden Markov Model to infer the fluorescence intensity levels of each trace. We then use the inferred intensity level state trajectory from all ROIs to infer kinetic rates. Our method has several advantages, including the ability to analyze noisy traces, account for the presence of photobleaching events, and provide uncertainties associated with the inferred binding kinetics. We demonstrate the effectiveness and reliability of our method through simulations and data from DNA origami binding experiments.

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.

On Kinetic Constraints That Catalysis Imposes on Elementary Processes

Catalysis, the acceleration of product formation by a substance that is left unchanged, typically results from multiple elementary processes, including diffusion of the reactants toward the catalyst, chemical steps, and release of the products. While efforts to design catalysts are often focused on accelerating the chemical reaction on the catalyst, catalysis is a global property of the catalytic cycle that involves all processes. These are controlled by both intrinsic parameters such as the composition and shape of the catalyst and extrinsic parameters such as the concentration of the chemical species at play. We examine here the conditions that catalysis imposes on the different steps of a reaction cycle and the respective role of intrinsic and extrinsic parameters of the system on the emergence of catalysis by using an approach based on first-passage times. We illustrate this approach for various decompositions of a catalytic cycle into elementary steps, including non-Markovian decompositions, which are useful when the presence and nature of intermediate states are a priori unknown. Our examples cover different types of reactions and clarify the constraints on elementary steps and the impact of species concentrations on catalysis.

Reaction-path statistical mechanics of enzymatic kinetics

We introduce a reaction-path statistical mechanics formalism based on the principle of large deviations to quantify the kinetics of single-molecule enzymatic reaction processes under the Michaelis-Menten mechanism, which exemplifies an out-of-equilibrium process in the living system. Our theoretical approach begins with the principle of equal a priori probabilities and defines the reaction path entropy to construct a new nonequilibrium ensemble as a collection of possible chemical reaction paths. As a result, we evaluate a variety of path-based partition functions and free energies by using the formalism of statistical mechanics. They allow us to calculate the timescales of a given enzymatic reaction, even in the absence of an explicit boundary condition that is necessary for the equilibrium ensemble. We also consider the large deviation theory under a closed-boundary condition of the fixed observation time to quantify the enzyme-substrate unbinding rates. The result demonstrates the presence of a phase-separation-like, bimodal behavior in unbinding events at a finite timescale, and the behavior vanishes as its rate function converges to a single phase in the long-time limit.

Modularity in Biological Networks

Network modeling, from the ecological to the molecular scale has become an essential tool for studying the structure, dynamics and complex behavior of living systems. Graph representations of the relationships between biological components open up a wide variety of methods for discovering the mechanistic and functional properties of biological systems. Many biological networks are organized into a modular structure, so methods to discover such modules are essential if we are to understand the biological system as a whole. However, most of the methods used in biology to this end, have a limited applicability, as they are very specific to the system they were developed for. Conversely, from the statistical physics and network science perspective, graph modularity has been theoretically studied and several methods of a very general nature have been developed. It is our perspective that in particular for the modularity detection problem, biology and theoretical physics/network science are less connected than they should. The central goal of this review is to provide the necessary background and present the most applicable and pertinent methods for community detection in a way that motivates their further usage in biological research.

The role of energy cost on accuracy, sensitivity, specificity, speed and adaptation of T cell foreign and self recognition

The critical role of energy consumption in biological systems including T cell discrimination process has been investigated in various ways. The kinetic proofreading (KPR) in T cell recognition involving different levels of energy dissipation influences functional outcomes such as error rates and specificity. In this work, we study quantitatively how the energy cost influences error fractions, sensitivity, specificity, kinetic speed in terms of Mean First Passage Time (MFPT) and adaption errors. These provide the background to adequately understand T cell dynamics. It is found that energy plays a central role in the system that aims to achieve minimum error fractions and maximum sensitivity and specificity with the fastest speed under our kinetic scheme for which numerical values of kinetic parameters are specially chosen, but such a condition can be broken with varying data. Starting with the application of steady state approximation (SSA) to the evaluation of the concentration of each complex produced associated with KPR, which is used to quantify various observables, we present both analytical and numerical results in detail.

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.

Molecular mechanism of selective substrate engagement and inhibitor disengagement of cysteine synthase

O-acetyl serine sulfhydrylase (OASS), referred to as cysteine synthase (CS), synthesizes cysteine from O-acetyl serine (OAS) and sulfur in bacteria and plants. The inherent challenge for CS is to overcome 4 to 6 log-folds stronger affinity for its natural inhibitor, serine acetyltransferase (SAT), as compared with its affinity for substrate, OAS. Our recent study showed that CS employs a novel competitive-allosteric mechanism to selectively recruit its substrate in the presence of natural inhibitor. In this study, we trace the molecular features that control selective substrate recruitment. To generalize our findings, we used CS from three different bacteria (Haemophilus, Salmonella, and Mycobacterium) as our model systems and analyzed structural and substrate-binding features of wild-type CS and its ∼13 mutants. Results show that CS uses a noncatalytic residue, M120, located 20 Å away from the reaction center, to discriminate in favor of substrate. M120A and background mutants display significantly reduced substrate binding, catalytic efficiency, and inhibitor binding. Results shows that M120 favors the substrate binding by selectively enhancing the affinity for the substrate and disengaging the inhibitor by 20 to 286 and 5- to 3-folds, respectively. Together, M120 confers a net discriminative force in favor of substrate by 100- to 858-folds.

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.

Different time scales in dynamic systems with multiple outcomes

Stochastic biochemical and transport processes have various final outcomes, and they can be viewed as dynamic systems with multiple exits. Many current theoretical studies, however, typically consider only a single time scale for each specific outcome, effectively corresponding to a single-exit process and assuming the independence of each exit process. However, the presence of other exits influences the statistical properties and dynamics measured at any specific exit. Here, we present theoretical arguments to explicitly show the existence of different time scales, such as mean exit times and inverse exit fluxes, for dynamic processes with multiple exits. This implies that the statistics of any specific exit dynamics cannot be considered without taking into account the presence of other exits. Several illustrative examples are described in detail using analytical calculations, mean-field estimates, and kinetic Monte Carlo computer simulations. The underlying microscopic mechanisms for the existence of different time scales are discussed. The results are relevant for understanding the mechanisms of various biological, chemical, and industrial processes, including transport through channels and pores.

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

One of the most intriguing features of biological systems is their ability to regulate the steady-state fluxes of the underlying biochemical reactions; however, the regulatory mechanisms and their physicochemical properties are not fully understood. Fundamentally, flux regulation can be explained with a chemical kinetic formalism describing the transitions between discrete states, with the reaction rates defined by an underlying free energy landscape. Which features of the energy landscape affect the flux distribution? Here we prove that the ratios of the steady-state fluxes of quasi-first-order biochemical processes are invariant to energy perturbations of the discrete states and are only affected by the energy barriers. In other words, the nonequilibrium flux distribution is under kinetic and not thermodynamic control. We illustrate the generality of this result for three biological processes. For the network describing protein folding along competing pathways, the probabilities of proceeding via these pathways are shown to be invariant to the stability of the intermediates or to the presence of additional misfolded states. For the network describing protein synthesis, the error rate and the energy expenditure per peptide bond is proven to be independent of the stability of the intermediate states. For molecular motors such as myosin-V, the ratio of forward to backward steps and the number of adenosine 5'-triphosphate (ATP) molecules hydrolyzed per step is demonstrated to be invariant to energy perturbations of the intermediate states. These findings place important constraints on the ability of mutations and drug perturbations to affect the steady-state flux distribution for a wide class of biological processes.

RIG-I Uses an ATPase-Powered Translocation-Throttling Mechanism for Kinetic Proofreading of RNAs and Oligomerization

RIG-I has a remarkable ability to specifically select viral 5'ppp dsRNAs for activation from a pool of cytosolic self-RNAs. The ATPase activity of RIG-I plays a role in RNA discrimination and activation, but the underlying mechanism was unclear. Using transient-state kinetics, we elucidated the ATPase-driven "kinetic proofreading" mechanism of RIG-I activation and RNA discrimination, akin to DNA polymerases, ribosomes, and T cell receptors. Even in the autoinhibited state of RIG-I, the C-terminal domain kinetically discriminates against self-RNAs by fast off rates. ATP binding facilitates dsRNA engagement but, interestingly, makes RIG-I promiscuous, explaining the constitutive signaling by Singleton-Merten syndrome-linked mutants that bind ATP without hydrolysis. ATP hydrolysis dissociates self-RNAs faster than 5'ppp dsRNA but, more importantly, drives RIG-I oligomerization through translocation, which we show to be regulated by helicase motif IVa. RIG-I translocates directionally from the dsRNA end into the stem region, and the 5'ppp end "throttles" translocation to provide a mechanism for threading and building a signaling-active oligomeric complex.

Energy-speed-accuracy relation in complex networks for biological discrimination

Discriminating between correct and incorrect substrates is a core process in biology, but how is energy apportioned between the conflicting demands of accuracy (μ), speed (σ), and total entropy production rate (P)? Previous studies have focused on biochemical networks with simple structure or relied on simplifying kinetic assumptions. Here, we use the linear framework for timescale separation to analytically examine steady-state probabilities away from thermodynamic equilibrium for networks of arbitrary complexity. We also introduce a method of scaling parameters that is inspired by Hopfield's treatment of kinetic proofreading. Scaling allows asymptotic exploration of high-dimensional parameter spaces. We identify in this way a broad class of complex networks and scalings for which the quantity σln(μ)/P remains asymptotically finite whenever accuracy improves from equilibrium, so that μ_{eq}/μ→0. Scalings exist, however, even for Hopfield's original network, for which σln(μ)/P is asymptotically infinite, illustrating the parametric complexity. Outside the asymptotic regime, numerical calculations suggest that, under more restrictive parametric assumptions, networks satisfy the bound, σln(μ/μ_{eq})/P<1, and we discuss the biological implications for discrimination by ribosomes and DNA polymerase. The methods introduced here may be more broadly useful for analyzing complex networks that implement other forms of cellular information processing.