Michaelis-Menten equation and detailed balance in enzymatic networks

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.

The linear framework II: using graph theory to analyse the transient regime of Markov processes

The linear framework uses finite, directed graphs with labelled edges to model biomolecular systems. Graph vertices represent chemical species or molecular states, edges represent reactions or transitions and edge labels represent rates that also describe how the system is interacting with its environment. The present paper is a sequel to a recent review of the framework that focussed on how graph-theoretic methods give insight into steady states as rational algebraic functions of the edge labels. Here, we focus on the transient regime for systems that correspond to continuous-time Markov processes. In this case, the graph specifies the infinitesimal generator of the process. We show how the moments of the first-passage time distribution, and related quantities, such as splitting probabilities and conditional first-passage times, can also be expressed as rational algebraic functions of the labels. This capability is timely, as new experimental methods are finally giving access to the transient dynamic regime and revealing the computations and information processing that occur before a steady state is reached. We illustrate the concepts, methods and formulas through examples and show how the results may be used to illuminate previous findings in the literature.

Identifying the active sites in unequal iron-nitrogen single-atom catalysts

Single-atom catalysts (SACs) have become one of the most attractive frontier research fields in catalysis and energy conversion. However, due to the atomic heterogeneity of SACs and limitations of ensemble-averaged measurements, the essential active sites responsible for governing specific catalytic properties and mechanisms remain largely concealed. In this study, we develop a quantitative method of single-atom catalysis-fluorescence correlation spectroscopy (SAC-FCS), leveraging the atomic structure-dependent catalysis kinetics and single-turnover resolution of single-molecule fluorescence microscopy. This method enables us to investigate the oxidase-like single-molecule catalysis on unidentical iron-nitrogen (Fe-N) coordinated SACs, quantifying the active sites and their kinetic parameters. The findings reveal the significant differences of single sites from the average behaviors and corroborate the oxidase-like catalytic mechanism of the Fe-N active sites. We anticipate that the method will give essential insights into the rational design and application of SACs.

Impact of Solvent on State-to-State Population Transport in Multistate Systems Using Coherences

Understanding the pathways taken by a quantum particle during a transport process is an enormous challenge. There are broadly two different aspects of the problem that affect the route taken. First is obviously the couplings between the various sites, which translates into the intrinsic "strength" of a state-to-state channel. Apart from these inter-state couplings, the relative coupling strengths and timescales of the solvent modes form the second factor. This impact of the dissipative environment is significantly more difficult to analyze. Building on the recently derived relations between coherences and population derivatives, we present an analysis of the transport that allows us to account for both the effects in a rigorous manner. We demonstrate the richness hidden behind the transport even for a relatively simple system, a 4-site coarse-grained model of the Fenna-Matthews-Olson complex. The effect of the local dissipative media is highly nontrivial. We show that while the impact on the total site population may be small, there are noticeable changes to the pathway taken by the transport process. We also demonstrate how an analysis in a similar spirit can be done using the Förster approximation. The ability to untangle the dynamics at a greater granularity opens up possibilities in terms of design of novel systems with an eye toward quantum control.

Analytical Decomposition of Transition Flux to Cycle Durations via Integration of Transition Times

Rigorous methods of decomposing kinetic networks to cycles are available, but the solutions usually contain entangled transition rates, which are difficult to analyze. This study proposes a new method of decomposing net transition flux to cycle durations, and the duration of each cycle is an integration of the transition times along the cycle. The method provides a series of neat dependences from the basic kinetic variables to the final flux, which support direct analysis based on the formulas. An assisting transformation diagram from symmetric conductivity to asymmetric conductivity is provided, which largely simplifies the application of the method. The method is likely a useful analytical tool for many studies relevant to kinetics and networks. Applications of the method shall provide new kinetic and thermodynamic information for the studied system.

Computational Aspects of Single-Molecule Kinetics for Coupled Catalytic Cycles: A Spectral Analysis

Catalysis from single active sites is analyzed using methods developed from single-molecule kinetics. Using a stochastic Markov-state description, the observable properties of general catalytic networks of reactions are expressed using an eigenvalue decomposition of the transition matrix for the Markov process. By the use of a sensitivity analysis, the necessary eigenvalues and eigenvectors are related to the energies of controlling barriers and wells located along the reaction routes. A generalization of the energetic span theory allows the eigenvalues to be computed from several activation energies corresponding to distinct barrier-well pairings. The formalism is demonstrated for model problems and for a physically realistic mechanism for an alkene hydrogenation reaction on a single-atom catalyst. The spectral analysis permits a hierarchy of timescales to be identified from the single-molecule signal, which correspond to specific relaxation modes in the network.

Nonlinear Analytics for Electrochemical Biosensor Design Using Enzyme Aggregates and Delayed Mass Action

The paper is devoted to the extension of Brown’s model of enzyme kinetics to the case with distributed delays. Firstly, we construct a multi-substrate multi-inhibitor model using discrete and distributed delays. Furthermore, we consider simplified models including one substrate and one inhibitor, for which an experimental study has been performed. The algorithm of parameter identifications was developed which was tested on the experimental data of solution conductivity. Both the model and Kohlrausch’s law parameters are obtained as a result of the optimization procedure. Comparison of plots constructed with the help of the estimated parameters has shown that in such case the model with distributed delays is more chemically adequate in comparison with the discrete one. The methods of generalization of the results to the multi-substrate multi-inhibitor cases are discussed.

Understanding the Optimal Cooperativity of Human Glucokinase: Kinetic Resonance in Nonequilibrium Conformational Fluctuations

The cooperativity of a monomeric enzyme arises from dynamic correlation instead of spatial correlation and is a consequence of nonequilibrium conformation fluctuations. We investigate the conformation-modulated kinetics of human glucokinase, a monomeric enzyme with important physiological functions, using a five-state kinetic model. We derive the non-Michealis-Menten (MM) correction term of the activity (i.e., turnover rate), predict its relationship to cooperativity, and reveal the violation of conformational detailed balance. Most importantly, we reproduce and explain the observed resonance effect in human glucokinase (i.e., maximal cooperativity when the conformational fluctuation rate is comparable to the catalytic rate). With the realistic parameters, our theoretical results are in quantitative agreement with the reported measurement by Miller and co-workers. The analysis can be extended to a general chemical network beyond the five-state model, suggesting the generality of kinetic cooperativity and resonance.

Nonequilibrium Thermodynamics in Cell Biology: Extending Equilibrium Formalism to Cover Living Systems

We discuss new developments in the nonequilibrium dynamics and thermodynamics of living systems, giving a few examples to demonstrate the importance of nonequilibrium thermodynamics for understanding biological dynamics and functions. We study single-molecule enzyme dynamics, in which the nonequilibrium thermodynamic and dynamic driving forces of chemical potential and flux are crucial for the emergence of non-Michaelis-Menten kinetics. We explore single-gene expression dynamics, in which nonequilibrium dissipation can suppress fluctuations. We investigate the cell cycle and identify the nutrition supply as the energy input that sustains the stability, speed, and coherence of cell cycle oscillation, from which the different vital phases of the cell cycle emerge. We examine neural decision-making processes and find the trade-offs among speed, accuracy, and thermodynamic costs that are important for neural function. Lastly, we consider the thermodynamic cost for specificity in cellular signaling and adaptation.

Shear viscosity of two-state enzyme solutions

We discuss the shear viscosity of a Newtonian solution of catalytic enzymes and substrate molecules. The enzyme is modeled as a two-state dimer consisting of two spherical domains connected with an elastic spring. The enzymatic conformational dynamics is induced by the substrate binding and such a process is represented by an additional elastic spring. Employing the Boltzmann distribution weighted by the waiting times of enzymatic species in each catalytic cycle, we obtain the shear viscosity of dilute enzyme solutions as a function of substrate concentration and its physical properties. The substrate affinity distinguishes between fast and slow enzymes, and the corresponding viscosity expressions are obtained. Furthermore, we connect the obtained viscosity with the diffusion coefficient of a tracer particle in enzyme solutions.

Quantifying the flux as the driving force for nonequilibrium dynamics and thermodynamics in non-Michaelis-Menten enzyme kinetics

The driving force for active physical and biological systems is determined by both the underlying landscape and nonequilibrium curl flux. While landscape can be experimentally quantified from the histograms of the collected real-time trajectories of the observables, quantifying the experimental flux remains challenging. In this work, we studied the single-molecule enzyme dynamics of horseradish peroxidase with dihydrorhodamine 123 and hydrogen peroxide (H2O2) as substrates. Surprisingly, significant deviations in the kinetics from the conventional Michaelis-Menten reaction rate were observed. Instead of a linear relationship between the inverse of the enzyme kinetic rate and the inverse of substrate concentration, a nonlinear relationship between the two emerged. We identified nonequilibrium flux as the origin of such non-Michaelis-Menten enzyme rate behavior. Furthermore, we quantified the nonequilibrium flux from experimentally obtained fluorescence correlation spectroscopy data and showed this flux to led to the deviations from the Michaelis-Menten kinetics. We also identified and quantified the nonequilibrium thermodynamic driving forces as the chemical potential and entropy production for such non-Michaelis-Menten kinetics. Moreover, through isothermal titration calorimetry measurements, we identified and quantified the origin of both nonequilibrium dynamic and thermodynamic driving forces as the heat absorbed (energy input) into the enzyme reaction system. Furthermore, we showed that the nonequilibrium driving forces led to time irreversibility through the difference between the forward and backward directions in time and high-order correlations were associated with the deviations from Michaelis-Menten kinetics. This study provided a general framework for experimentally quantifying the dynamic and thermodynamic driving forces for nonequilibrium systems.

Chemical Kinetics Roots and Methods to Obtain the Probability Distribution Function Evolution of Reactants and Products in Chemical Networks Governed by a Master Equation

In this paper first, we review the physical root bases of chemical reaction networks as a Markov process in multidimensional vector space. Then we study the chemical reactions from a microscopic point of view, to obtain the expression for the propensities for the different reactions that can happen in the network. These chemical propensities, at a given time, depend on the system state at that time, and do not depend on the state at an earlier time indicating that we are dealing with Markov processes. Then the Chemical Master Equation (CME) is deduced for an arbitrary chemical network from a probability balance and it is expressed in terms of the reaction propensities. This CME governs the dynamics of the chemical system. Due to the difficulty to solve this equation two methods are studied, the first one is the probability generating function method or z-transform, which permits to obtain the evolution of the factorial moment of the system with time in an easiest way or after some manipulation the evolution of the polynomial moments. The second method studied is the expansion of the CME in terms of an order parameter (system volume). In this case we study first the expansion of the CME using the propensities obtained previously and splitting the molecular concentration into a deterministic part and a random part. An expression in terms of multinomial coefficients is obtained for the evolution of the probability of the random part. Then we study how to reconstruct the probability distribution from the moments using the maximum entropy principle. Finally, the previous methods are applied to simple chemical networks and the consistency of these methods is studied.

Charge Transfer between [4Fe4S] Proteins and DNA Is Unidirectional: Implications for Biomolecular Signaling

Recent experiments suggest that DNA-mediated charge transport might enable signaling between the [4Fe4S] clusters in the C-terminal domains of human DNA primase and polymerase α, as well as the signaling between other replication and repair high-potential [4Fe4S] proteins. Our theoretical study demonstrates that the redox signaling cannot be accomplished exclusively by DNA-mediated charge transport because part of the charge transfer chain has an unfavorable free energy profile. We show that hole or excess electron transfer between a [4Fe4S] cluster and a nucleic acid duplex through a protein medium can occur within microseconds in one direction, while it is kinetically hindered in the opposite direction. We present a set of signaling mechanisms that may occur with the assistance of oxidants or reductants, using the allowed charge transfer processes. These mechanisms would enable the coordinated action of [4Fe4S] proteins on DNA, engaging the [4Fe4S] oxidation state dependence of the protein-DNA binding affinity.

Structural conditions on complex networks for the Michaelis-Menten input-output response

The Michaelis-Menten (MM) fundamental formula describes how the rate of enzyme catalysis depends on substrate concentration. The familiar hyperbolic relationship was derived by timescale separation for a network of three reactions. The same formula has subsequently been found to describe steady-state input-output responses in many biological contexts, including single-molecule enzyme kinetics, gene regulation, transcription, translation, and force generation. Previous attempts to explain its ubiquity have been limited to networks with regular structure or simplifying parametric assumptions. Here, we exploit the graph-based linear framework for timescale separation to derive general structural conditions under which the MM formula arises. The conditions require a partition of the graph into two parts, akin to a "coarse graining" into the original MM graph, and constraints on where and how the input variable occurs. Other features of the graph, including the numerical values of parameters, can remain arbitrary, thereby explaining the formula's ubiquity. For systems at thermodynamic equilibrium, we derive a necessary and sufficient condition. For systems away from thermodynamic equilibrium, especially those with irreversible reactions, distinct structural conditions arise and a general characterization remains open. Nevertheless, our results accommodate, in much greater generality, all examples known to us in the literature.

Effect of Substrate Number Fluctuations in Stochastic Enzyme Kinetics

Conventional studies on enzyme kinetics assume that the substrate concentration remains constant. However, for catalytic reactions taking place in intracellular compartments, the substrate molecules are fed in and out of the compartment and are catalyzed into product molecules within the compartment. In this work, we use a chemical master equation approach to study the stochastic kinetics of a single enzyme for different reaction pathways with one or more intermediate states. We obtain velocity expressions that deviate from the Michaelis-Menten expression. We study the coefficient of variation, which is a measure of the noise strength as a function of the mean substrate concentration for systems where there is influx or/and outflux of substrate molecules. Our results show that the noise strength decreases with the increase in the substrate concentration and finally remains the same when the substrate is present in abundance.

Single-Molecule Kinetics of an Enzyme in the Presence of Multiple Substrates

Herein, the catalytic activity of a single enzyme in the presence of multiple substrates is studied. Three different mechanisms of bisubstrate binding, namely, ordered sequential, random sequential and ping-pong nonsequential pathway, are broadly discussed. By means of the chemical master equation approach, exact expressions for the waiting-time distributions, the mean turnover time and the randomness parameter as a function of the substrate concentration, such that both concentrations are fixed, but one of them is changed quasi-statically are obtained. The randomness parameter is not equal to unity at intermediate to high substrate concentrations, which indicates the presence of multiple rate-limiting steps in the reaction pathway in all three modes of bisubstrate binding. This arises due to transitions between the free enzyme and the enzyme-substrate complexes that occur on comparable timescales. Such turnover statistics of the single enzyme can also distinguish between the different types of bisubstrate binding mechanisms.

First-passage times in renewal and nonrenewal systems

Fluctuations in stochastic systems are usually characterized by full counting statistics, which analyzes the distribution of the number of events taking place in the fixed time interval. In an alternative approach, the distribution of the first-passage times, i.e., the time delays after which the counting variable reaches a certain threshold value, is studied. This paper presents the approach to calculate the first-passage time distribution in systems in which the analyzed current is associated with an arbitrary set of transitions within the Markovian network. Using this approach, it is shown that when the subsequent first-passage times are uncorrelated, there exist strict relations between the cumulants of the full counting statistics and the first-passage time distribution. On the other hand, when the correlations of the first-passage times are present, their distribution may provide additional information about the internal dynamics of the system in comparison to the full counting statistics; for example, it may reveal the switching between different dynamical states of the system. Additionally, I show that breaking of the fluctuation theorem for first-passage times may reveal the multicyclic nature of the Markovian network.

Single-Molecule Plasmon Sensing: Current Status and Future Prospects

Single-molecule detection has long relied on fluorescent labeling with high quantum-yield fluorophores. Plasmon-enhanced detection circumvents the need for labeling by allowing direct optical detection of weakly emitting and completely nonfluorescent species. This review focuses on recent advances in single molecule detection using plasmonic metal nanostructures as a sensing platform, particularly using a single particle-single molecule approach. In the past decade two mechanisms for plasmon-enhanced single-molecule detection have been demonstrated: (1) by plasmonically enhancing the emission of weakly fluorescent biomolecules, or (2) by monitoring shifts of the plasmon resonance induced by single-molecule interactions. We begin with a motivation regarding the importance of single molecule detection, and advantages plasmonic detection offers. We describe both detection mechanisms and discuss challenges and potential solutions. We finalize by highlighting the exciting possibilities in analytical chemistry and medical diagnostics.

Conformational Nonequilibrium Enzyme Kinetics: Generalized Michaelis-Menten Equation

In a conformational nonequilibrium steady state (cNESS), enzyme turnover is modulated by the underlying conformational dynamics. On the basis of a discrete kinetic network model, we use an integrated probability flux balance method to derive the cNESS turnover rate for a conformation-modulated enzymatic reaction. The traditional Michaelis-Menten (MM) rate equation is extended to a generalized form, which includes non-MM corrections induced by conformational population currents within combined cyclic kinetic loops. When conformational detailed balance is satisfied, the turnover rate reduces to the MM functional form, explaining its general validity. For the first time, a one-to-one correspondence is established between non-MM terms and combined cyclic loops with unbalanced conformational currents. Cooperativity resulting from nonequilibrium conformational dynamics can be achieved in enzymatic reactions, and we provide a novel, rigorous means of predicting and characterizing such behavior. Our generalized MM equation affords a systematic approach for exploring cNESS enzyme kinetics.

Nonclassical Kinetics of Clonal yet Heterogeneous Enzymes

Enzyme-to-enzyme variation in the catalytic rate is ubiquitous among single enzymes created from the same genetic information, which persists over the lifetimes of living cells. Despite advances in single-enzyme technologies, the lack of an enzyme reaction model accounting for the heterogeneous activity of single enzymes has hindered a quantitative understanding of the nonclassical stochastic outcome of single enzyme systems. Here we present a new statistical kinetics and exactly solvable models for clonal yet heterogeneous enzymes with possibly nonergodic state dynamics and state-dependent reactivity, which enable a quantitative understanding of modern single-enzyme experimental results for the mean and fluctuation in the number of product molecules created by single enzymes. We also propose a new experimental measure of the heterogeneity and nonergodicity for a system of enzymes.

On the generality of Michaelian kinetics

The reversible Michaelis-Menten equation is shown to follow from a very broad class of steady-state kinetic models involving enzymes that adopt a unique free (i.e., not complexed to substrate/product) state in solution. In the case of enzymes with multiple free states/conformations (e.g., fluctuating, hysteretic, or co-operative monomeric enzymes), Michaelian behavior is still assured if the relative steady-state populations of free enzyme states are independent of substrate and product concentration. Prior models for Michaelian behavior in multiple conformer enzymes are shown to be special cases of this single condition.

Probing conformational dynamics of an enzymatic active site by an in situ single fluorogenic probe under piconewton force manipulation

Unraveling the conformational details of an enzyme during the essential steps of a catalytic reaction (i.e., enzyme-substrate interaction, enzyme-substrate active complex formation, nascent product formation, and product release) is challenging due to the transient nature of intermediate conformational states, conformational fluctuations, and the associated complex dynamics. Here we report our study on the conformational dynamics of horseradish peroxidase using single-molecule multiparameter photon time-stamping spectroscopy with mechanical force manipulation, a newly developed single-molecule fluorescence imaging magnetic tweezers nanoscopic approach. A nascent-formed fluorogenic product molecule serves as a probe, perfectly fitting in the enzymatic reaction active site for probing the enzymatic conformational dynamics. Interestingly, the product releasing dynamics shows the complex conformational behavior with multiple product releasing pathways. However, under magnetic force manipulation, the complex nature of the multiple product releasing pathways disappears and more simplistic conformations of the active site are populated.

Dependence of the Enzymatic Velocity on the Substrate Dissociation Rate

Enzymes are biological catalysts that play a fundamental role in all living systems by supporting the selectivity and the speed for almost all cellular processes. While the general principles of enzyme functioning are known, the specific details of how they work at the microscopic level are not always available. Simple Michaelis-Menten kinetics assumes that the enzyme-substrate complex has only one conformation that decays as a single exponential. As a consequence, the enzymatic velocity decreases as the dissociation (off) rate constant of the complex increases. Recently, Reuveni et al. [ Proc. Natl. Acad. Sci. USA 2014 , 111 , 4391 - 4396 ] showed that it is possible for the enzymatic velocity to increase when the off rate becomes higher, if the enzyme-substrate complex has many conformations which dissociate with the same off rate constant. This was done using formal mathematical arguments, without specifying the nature of the dynamics of the enzyme-substrate complex. In order to provide a physical basis for this unexpected result, we derive an analytical expression for the enzymatic velocity assuming that the enzyme-substrate complex has multiple states and its conformational dynamics is described by rate equations with arbitrary rate constants. By applying our formalism to a complex with two conformations, we show that the unexpected off rate dependence of the velocity can be readily understood: If one of the conformations is unproductive, the system can escape from this "trap" by dissociating, thereby giving the enzyme another chance to form the productive enzyme-substrate complex. We also demonstrate that the nonmonotonic off rate dependence of the enzymatic velocity is possible not only when all off rate constants are identical, but even when they are different. We show that for typical experimentally determined rate constants, the nonmonotonic off rate dependence can occur for micromolar substrate concentrations. Finally, we discuss the relation of this work to the problem of optimizing the flux through singly occupied membrane channels and transporters.

Michaelis-Menten reaction scheme as a unified approach towards the optimal restart problem

We study the effect of restart, and retry, on the mean completion time of a generic process. The need to do so arises in various branches of the sciences and we show that it can naturally be addressed by taking advantage of the classical reaction scheme of Michaelis and Menten. Stopping a process in its midst-only to start it all over again-may prolong, leave unchanged, or even shorten the time taken for its completion. Here we are interested in the optimal restart problem, i.e., in finding a restart rate which brings the mean completion time of a process to a minimum. We derive the governing equation for this problem and show that it is exactly solvable in cases of particular interest. We then continue to discover regimes at which solutions to the problem take on universal, details independent forms which further give rise to optimal scaling laws. The formalism we develop, and the results obtained, can be utilized when optimizing stochastic search processes and randomized computer algorithms. An immediate connection with kinetic proofreading is also noted and discussed.

Michaelis–Menten kinetics under non-isothermal conditions

We extend the celebrated Michaelis–Menten kinetics description of an enzymatic reaction taking into consideration the presence of a thermal driving force.

Single-molecule chemical reaction reveals molecular reaction kinetics and dynamics

Understanding the microscopic elementary process of chemical reactions, especially in condensed phase, is highly desirable for improvement of efficiencies in industrial chemical processes. Here we show an approach to gaining new insights into elementary reactions in condensed phase by combining quantum chemical calculations with a single-molecule analysis. Elementary chemical reactions in liquid-phase, revealed from quantum chemical calculations, are studied by tracking the fluorescence of single dye molecules undergoing a reversible redox process. Statistical analyses of single-molecule trajectories reveal molecular reaction kinetics and dynamics of elementary reactions. The reactivity dynamic fluctuations of single molecules are evidenced and probably arise from either or both of the low-frequency approach of the molecule to the internal surface of the SiO2 nanosphere or the molecule diffusion-induced memory effect. This new approach could be applied to other chemical reactions in liquid phase to gain more insight into their molecular reaction kinetics and the dynamics of elementary steps.

Role of substrate unbinding in Michaelis-Menten enzymatic reactions

The Michaelis-Menten equation provides a hundred-year-old prediction by which any increase in the rate of substrate unbinding will decrease the rate of enzymatic turnover. Surprisingly, this prediction was never tested experimentally nor was it scrutinized using modern theoretical tools. Here we show that unbinding may also speed up enzymatic turnover--turning a spotlight to the fact that its actual role in enzymatic catalysis remains to be determined experimentally. Analytically constructing the unbinding phase space, we identify four distinct categories of unbinding: inhibitory, excitatory, superexcitatory, and restorative. A transition in which the effect of unbinding changes from inhibitory to excitatory as substrate concentrations increase, and an overlooked tradeoff between the speed and efficiency of enzymatic reactions, are naturally unveiled as a result. The theory presented herein motivates, and allows the interpretation of, groundbreaking experiments in which existing single-molecule manipulation techniques will be adapted for the purpose of measuring enzymatic turnover under a controlled variation of unbinding rates. As we hereby show, these experiments will not only shed first light on the role of unbinding but will also allow one to determine the time distribution required for the completion of the catalytic step in isolation from the rest of the enzymatic turnover cycle.

Automated maximum likelihood separation of signal from baseline in noisy quantal data

Data recordings often include high-frequency noise and baseline fluctuations that are not generated by the system under investigation, which need to be removed before analyzing the signal for the system's behavior. In the absence of an automated method, experimentalists fall back on manual procedures for removing these fluctuations, which can be laborious and prone to subjective bias. We introduce a maximum likelihood formalism for separating signal from a drifting baseline plus noise, when the signal takes on integer multiples of some value, as in ion channel patch-clamp current traces. Parameters such as the quantal step size (e.g., current passing through a single channel), noise amplitude, and baseline drift rate can all be optimized automatically using the expectation-maximization algorithm, taking the number of open channels (or molecules in the on-state) at each time point as a hidden variable. Our goal here is to reconstruct the signal, not model the (possibly highly complex) underlying system dynamics. Thus, our likelihood function is independent of those dynamics. This may be thought of as restricting to the simplest possible hidden Markov model for the underlying channel current, in which successive measurements of the state of the channel(s) are independent. The resulting method is comparable to an experienced human in terms of results, but much faster. FORTRAN 90, C, R, and JAVA codes that implement the algorithm are available for download from our website.

Exploration of the spontaneous fluctuating activity of single enzyme molecules

Single enzyme molecules display inevitable, stochastic fluctuations in their catalytic activity. In metabolism, for instance, the stochastic activity of individual enzymes is averaged out due to their high copy numbers per single cell. However, many processes inside cells rely on single enzyme activity, such as transcription, replication, translation, and histone modifications. Here we introduce the main theoretical concepts of stochastic single-enzyme activity starting from the Michaelis-Menten enzyme mechanism. Next, we discuss stochasticity of multi-substrate enzymes, of enzymes and receptors with multiple conformational states and finally, how fluctuations in receptor activity arise from fluctuations in signal concentration. This paper aims to introduce the exciting field of single-molecule enzyme kinetics and stochasticity to a wider audience of biochemists and systems biologists.

Aggregated markov model using time series of single molecule dwell times with minimum excessive information

Statistics of the dwell times, the stationary state distributions (SSDs), are often studied to infer the underlying kinetics from a single molecule finite-level time series. However, it is well known that the underlying kinetic scheme, a hidden Markov model (HMM), cannot be identified uniquely from the SSDs because some features of the underlying HMM are hidden by finite-level measurements. Here, we quantify the amount of excessive information in a given HMM that is not warranted by the measured SSDs and extract the HMM with minimum excessive information as the most objective representation of the data. The method is applied to a single molecule enzymatic turnover experiment, and the origin of dynamic disorder is discussed in terms of the network properties of the HMM.

Optimization of enzyme parameters for fermentative production of biorenewable fuels and chemicals

Microbial biocatalysts such as Escherichia coli and Saccharomyces cerevisiae have been extensively subjected to Metabolic Engineering for the fermentative production of biorenewable fuels and chemicals. This often entails the introduction of new enzymes, deletion of unwanted enzymes and efforts to fine-tune enzyme abundance in order to attain the desired strain performance. Enzyme performance can be quantitatively described in terms of the Michaelis-Menten type parameters K_m, turnover number k_cat and K_i, which roughly describe the affinity of an enzyme for its substrate, the speed of a reaction and the enzyme sensitivity to inhibition by regulatory molecules. Here we describe examples of where knowledge of these parameters have been used to select, evolve or engineer enzymes for the desired performance and enabled increased production of biorenewable fuels and chemicals. Examples include production of ethanol, isobutanol, 1-butanol and tyrosine and furfural tolerance. The Michaelis-Menten parameters can also be used to judge the cofactor dependence of enzymes and quantify their preference for NADH or NADPH. Similarly, enzymes can be selected, evolved or engineered for the preferred cofactor preference. Examples of exporter engineering and selection are also discussed in the context of production of malate, valine and limonene.

Origins of concentration dependence of waiting times for single-molecule fluorescence binding

Binary fluorescence time series obtained from single-molecule imaging experiments can be used to infer protein binding kinetics, in particular, association and dissociation rate constants from waiting time statistics of fluorescence intensity changes. In many cases, rate constants inferred from fluorescence time series exhibit nonintuitive dependence on ligand concentration. Here, we examine several possible mechanistic and technical origins that may induce ligand dependence of rate constants. Using aggregated Markov models, we show under the condition of detailed balance that non-fluorescent bindings and missed events due to transient interactions, instead of conformation fluctuations, may underly the dependence of waiting times and thus apparent rate constants on ligand concentrations. In general, waiting times are rational functions of ligand concentration. The shape of concentration dependence is qualitatively affected by the number of binding sites in the single molecule and is quantitatively tuned by model parameters. We also show that ligand dependence can be caused by non-equilibrium conditions which result in violations of detailed balance and require an energy source. As to a different but significant mechanism, we examine the effect of ambient buffers that can substantially reduce the effective concentration of ligands that interact with the single molecules. To demonstrate the effects by these mechanisms, we applied our results to analyze the concentration dependence in a single-molecule experiment EGFR binding to fluorophore-labeled adaptor protein Grb2 by Morimatsu et al. [Proc. Natl. Acad. Sci. U.S.A. 104, 18013 (2007)].

Interpreting single turnover catalysis measurements with constrained mean dwell times

Observation of a chemical transformation at the single-molecule level yields a detailed view of kinetic pathways contributing to the averaged results obtained in a bulk measurement. Studies of a fluorogenic reaction catalyzed by gold nanoparticles have revealed heterogeneous reaction dynamics for these catalysts. Measurements on single nanoparticles yield binary trajectories with stochastic transitions between a dark state in which no product molecules are adsorbed and a fluorescent state in which one product molecule is present. The mean dwell time in either state gives information corresponding to a bulk measurement. Quantifying fluctuations from mean kinetics requires identifying properties of the fluorescence trajectory that are selective in emphasizing certain dynamic processes according to their time scales. We propose the use of constrained mean dwell times, defined as the mean dwell time in a state with the constraint that the immediately preceding dwell time in the other state is, for example, less than a variable time. Calculations of constrained mean dwell times for a kinetic model with dynamic disorder demonstrate that these quantities reveal correlations among dynamic fluctuations at different active sites on a multisite catalyst. Constrained mean dwell times are determined from measurements of single nanoparticle catalysis. The results indicate that dynamical fluctuations at different active sites are correlated, and that especially rapid reaction events produce particularly slowly desorbing product molecules.