Fluctuating enzymes: lessons from single-molecule studies

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

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.

Probing conformational kinetics of catalase with and without magnetic field by single-entity collision electrochemistry

The conformational motions of enzymes are crucial for their catalytic activities, but these fluctuations are usually spontaneous and unsynchronized and thus difficult to obtain from ensemble-averaged measurements. Here, we employ label-free single-entity electrochemical measurements to monitor in real time the fluctuating enzymatic behavior of single catalase molecules toward the degradation of hydrogen peroxide. By probing the electrochemical signals of single catalase molecules at a carbon nanoelectrode, we were able to observe three distinct current traces that could be attributed to conformational changes on the sub-millisecond timescale. Whereas, nearly uniform single long peaks were observed for single catalase molecules under a moderate magnetic field due to the restricted conformational changes of catalase. By combining high-resolution current signals with a multiphysics simulation model, we studied the catalytic kinetics of catalase with and without a magnetic field, and further estimated the maximum catalytic rate and conformational transition rate. This work introduces a new complementary approach to existing single-molecule enzymology, giving further insight into the enzymatic reaction mechanism.

Molecular noise-induced activator-inhibitor duality in enzyme inhibition kinetics

Classical theories of enzyme inhibition kinetics predict a monotonic decrease in the mean catalytic activity with the increase in inhibitor concentration. The steady-state result, derived from deterministic mass action kinetics, ignores molecular noise in enzyme-inhibition mechanisms. Here, we present a stochastic generalization of enzyme inhibition kinetics to mesoscopic enzyme concentrations by systematically accounting for molecular noise in competitive and uncompetitive mechanisms of enzyme inhibition. Our work reveals an activator-inhibitor duality as a non-classical effect in the transient regime in which inhibitors tend to enhance enzymatic activity. We introduce statistical measures that quantify this counterintuitive response through the stochastic analog of the Lineweaver-Burk plot that shows a merging of the inhibitor-dependent velocity with the Michaelis-Menten velocity. The statistical measures of mean and temporal fluctuations - fractional enzyme activity and waiting time correlations - show a non-monotonic rise with the increase in inhibitors before subsiding to their baseline value. The inhibitor and substrate dependence of the fractional enzyme activity yields kinetic phase diagrams for non-classical activator-inhibitor duality. Our work links this duality to a molecular memory effect in the transient regime, arising from positive correlations between consecutive product turnover times. The vanishing of memory in the steady state recovers all the classical results.

Accurate quantification of DNA using on-site PCR (osPCR) by characterizing DNA amplification at single-molecule resolution

Despite the need in various applications, accurate quantification of nucleic acids still remains a challenge. The widely-used qPCR has reduced accuracy at ultralow template concentration and is susceptible to nonspecific amplifications. The more recently developed dPCR is costly and cannot handle high-concentration samples. We combine the strengths of qPCR and dPCR by performing PCR in silicon-based microfluidic chips and demonstrate high quantification accuracy in a large concentration range. Importantly, at low template concentration, we observe on-site PCR (osPCR), where only certain sites of the channel show amplification. The sites have almost identical ct values, showing osPCR is a quasi-single molecule phenomenon. Using osPCR, we can measure both the ct values and the absolute concentration of templates in the same reaction. Additionally, osPCR enables identification of each template molecule, allowing removal of nonspecific amplification during quantification and greatly improving quantification accuracy. We develop sectioning algorithm that improves the signal amplitude and demonstrate improved detection of COVID in patient samples.

Estimating and Assessing Differential Equation Models with Time-Course Data

Ordinary differential equation (ODE) models are widely used to describe chemical or biological processes. This Article considers the estimation and assessment of such models on the basis of time-course data. Due to experimental limitations, time-course data are often noisy, and some components of the system may not be observed. Furthermore, the computational demands of numerical integration have hindered the widespread adoption of time-course analysis using ODEs. To address these challenges, we explore the efficacy of the recently developed MAGI (MAnifold-constrained Gaussian process Inference) method for ODE inference. First, via a range of examples we show that MAGI is capable of inferring the parameters and system trajectories, including unobserved components, with appropriate uncertainty quantification. Second, we illustrate how MAGI can be used to assess and select different ODE models with time-course data based on MAGI's efficient computation of model predictions. Overall, we believe MAGI is a useful method for the analysis of time-course data in the context of ODE models, which bypasses the need for any numerical integration.

Amalgamation of DNAzymes and Nanozymes in a Coronazyme

Artificial enzymes such as nanozymes and DNAzymes are economical and stable alternatives to natural enzymes. By coating Au nanoparticles (AuNPs) with a DNA corona (AuNP@DNA), we amalgamated nanozymes and DNAzymes into a new artificial enzyme with catalytic efficiency 5 times higher than AuNP nanozymes, 10 times higher than other nanozymes, and significantly greater than most of the DNAzymes on the same oxidation reaction. The AuNP@DNA demonstrates excellent specificity as its reactivity on a reduction reaction does not change with respect to pristine AuNP. Single-molecule fluorescence and force spectroscopies and density functional theory (DFT) simulations indicate a long-range oxidation reaction initiated by radical production on the AuNP surface, followed by radical transport to the DNA corona, where the binding and turnover of substrates take place. The AuNP@DNA is named coronazyme because of its natural enzyme mimicking capability through the well-orchestrated structures and synergetic functions. By incorporating different nanocores and corona materials beyond DNAs, we anticipate that the coronazymes represent generic enzyme mimics to carry out versatile reactions in harsh environments.

Theoretical Tools to Quantify Stochastic Fluctuations in Single-Molecule Catalysis by Enzymes and Nanoparticles

Single-molecule microscopic techniques allow the counting of successive turnover events and the study of the time-dependent fluctuations of the catalytic activities of individual enzymes and different sites on a single heterogeneous nanocatalyst. It is important to establish theoretical methods to obtain the statistical measurements of such stochastic fluctuations that provide insight into the catalytic mechanism. In this review, we discuss a few theoretical frameworks for evaluating the first passage time distribution functions using a self-consistent pathway approach and chemical master equations, to establish a connection with experimental observables. The measurable probability distribution functions and their moments depend on the molecular details of the reaction and provide a way to quantify the molecular mechanisms of the reaction process. The statistical measurements of these fluctuations should provide insight into the enzymatic mechanism.

Structural Insights into Binding of Remdesivir Triphosphate within the Replication-Transcription Complex of SARS-CoV-2

Remdesivir is an adenosine analogue that has a cyano substitution in the C1' position of the ribosyl moiety and a modified base structure to stabilize the linkage of the base to the C1' atom with its strong electron-withdrawing cyano group. Within the replication-transcription complex (RTC) of SARS-CoV-2, the RNA-dependent RNA polymerase nsp12 selects remdesivir monophosphate (RMP) over adenosine monophosphate (AMP) for nucleotide incorporation but noticeably slows primer extension after the added RMP of the RNA duplex product is translocated by three base pairs. Cryo-EM structures have been determined for the RTC with RMP at the nucleotide-insertion (i) site or at the i + 1, i + 2, or i + 3 sites after product translocation to provide a structural basis for a delayed-inhibition mechanism by remdesivir. In this study, we applied molecular dynamics (MD) simulations to extend the resolution of structures to the measurable maximum that is intrinsically limited by MD properties of these complexes. Our MD simulations provide (i) a structural basis for nucleotide selectivity of the incoming substrates of remdesivir triphosphate over adenosine triphosphate and of ribonucleotide over deoxyribonucleotide, (ii) new detailed information on hydrogen atoms involved in H-bonding interactions between the enzyme and remdesivir, and (iii) direct information on the catalytically active complex that is not easily captured by experimental methods. Our improved resolution of interatomic interactions at the nucleotide-binding pocket between remedesivir and the polymerase could help to design a new class of anti-SARS-CoV-2 inhibitors.

Effects of Hydrodynamic Backflow on the Transmission Coefficient of a Barrier-Crossing Brownian Particle

The slow power law decay of the velocity autocorrelation function of a particle moving stochastically in a condensed-phase fluid is widely attributed to the momentum that fluid molecules displaced by the particle transfer back to it during the course of its motion. The forces created by this backflow effect are known as Basset forces, and they have been found in recent analytical work and numerical simulations to be implicated in a number of interesting dynamical phenomena, including boosted particle mobility in tilted washboard potentials. Motivated by these findings, the present paper is an investigation of the role of backflow in thermally activated barrier crossing, the governing process in essentially all condensed-phase chemical reactions. More specifically, it is an exact analytical calculation, carried out within the framework of the reactive-flux formalism, of the transmission coefficient κ(t) of a Brownian particle that crosses an inverted parabola under the influence of a colored noise process originating in the Basset force and a Markovian time-local friction. The calculation establishes that κ(t) is significantly enhanced over its backflow-free limit.

Structural Basis for Reduced Dynamics of Three Engineered HNH Endonuclease Lys-to-Ala Mutants for the Clustered Regularly Interspaced Short Palindromic Repeat (CRISPR)-Associated 9 (CRISPR/Cas9) Enzyme

Many bacteria possess type-II immunity against invading phages or plasmids known as the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated 9 (Cas9) system to detect and degrade the foreign DNA sequences. The Cas9 protein has two endonucleases responsible for double-strand breaks (the HNH domain for cleaving the target strand of DNA duplexes and RuvC domain for the nontarget strand, respectively) and a single-guide RNA-binding domain where the RNA and target DNA strands are base-paired. Three engineered single Lys-to-Ala HNH mutants (K810A, K848A, and K855A) exhibit an enhanced substrate specificity for cleavage of the target DNA strand. We report in this study that in the wild-type (wt) enzyme, D835, Y836, and D837 within the Y836-containing loop (comprising E827-D837) adjacent to the catalytic site have uncharacterizable broadened ¹H¹⁵N nuclear magnetic resonance (NMR) features, whereas remaining residues in the loop have different extents of broadened NMR spectra. We find that this loop in the wt enzyme exhibits three distinct conformations over the duration of the molecular dynamics simulations, whereas the three Lys-to-Ala mutants retain only one conformation. The versatility of multiple alternate conformations of this loop in the wt enzyme could help to recruit noncognate DNA substrates into the HNH active site for cleavage, thereby reducing its substrate specificity relative to the three mutants. Our study provides further experimental and computational evidence that Lys-to-Ala substitutions reduce dynamics of proteins and thus increase their stability.

Electrochemical Zero-Mode Waveguide Potential-Dependent Fluorescence of Glutathione Reductase at Single-Molecule Occupancy

Understanding functional states of individual redox enzymes is important because electron-transfer reactions are fundamental to life, and single-enzyme molecules exhibit molecule-to-molecule heterogeneity in their properties, such as catalytic activity. Zero-mode waveguides (ZMW) constitute a powerful tool for single-molecule studies, enabling investigations of binding reactions up to the micromolar range due to the ability to trap electromagnetic radiation in zeptoliter-scale observation volumes. Here, we report the potential-dependent fluorescence dynamics of single glutathione reductase (GR) molecules using a bimodal electrochemical ZMW (E-ZMW), where a single-ring electrode embedded in each of the nanopores of an E-ZMW array simultaneously serves to control electrochemical potential and to confine optical radiation within the nanopores. Here, the redox state of GR is manipulated using an external potential control of the Au electrode in the presence of a redox mediator, methyl viologen (MV). Redox-state transitions in GR are monitored by correlating electrochemical and spectroscopic signals from freely diffusing MV/GR in 60 zL effective observation volumes at single GR molecule average pore occupancy, ⟨n⟩ ∼ 0.8. Fluorescence intensities decrease (increase) at reducing (oxidizing) potentials for MV due to the MV-mediated control of the GR redox state. The spectroelectrochemical response of GR to the enzyme substrate, i.e., glutathione disulfide (GSSG), shows that GSSG promotes GR oxidation via enzymatic reduction. The capabilities of E-ZMWs to probe spectroelectrochemical phenomena in zL-scale-confined environments show great promise for the study of single-enzyme reactions and can be extended to important technological applications, such as those in molecular diagnostics.

Particle dynamics in viscoelastic media: Effects of non-thermal white noise on barrier crossing rates

The growing interest in the dynamics of self-driven particle motion has brought increased attention to the effects of non-thermal noise on condensed phase diffusion. Thanks to data recently collected by Ferrer et al. on activated dynamics in the presence of memory [Phys. Rev. Lett. 126, 108001 (2021)], some of these effects can now be characterized quantitatively. In the present paper, the data collected by Ferrer et al. are used to calculate the extent to which non-thermal white noise alters the time taken by single micron-sized silica particles in a viscoelastic medium to cross the barrier separating the two wells of an optically created bistable potential. The calculation-based on a generalized version of Kramers's flux-over-population approach-indicates that the added noise causes the barrier crossing rate (compared to the noise-free case) to first increase as a function of the noise strength and then to plateau to a constant value. The precise degree of rate enhancement may depend on how the data from the experiments conducted by Ferrer et al. are used in the flux-over-population approach. As claimed by Ferrer et al., this approach predicts barrier crossing times for the original silica-fluid system that agree almost perfectly with their experimental counterparts. However, this near-perfect agreement between theory and experiment is only achieved if the theoretical crossing times are obtained from the most probable values of a crossing time distribution constructed from the distributions of various parameters in Kramers's rate expression. If the mean values of these parameters are used in the expression instead, as would be commonly done, the theoretical crossing times are found to be as much as 1.5 times higher than the experimental values. However, these times turn out to be consistent with an alternative model of viscoelastic barrier crossing based on a mean first passage time formalism, which also uses mean parameter values in its rate expression. The rate enhancements predicted for barrier crossing under non-thermal noise are based on these mean parameter values and are open to experimental verification.

Synchronization and Enhanced Catalysis of Mechanically Coupled Enzymes

We examine the stochastic dynamics of two enzymes that are mechanically coupled to each other, e.g., through an elastic substrate or a fluid medium. The enzymes undergo conformational changes during their catalytic cycle, which itself is driven by stochastic steps along a biased chemical free energy landscape. We find conditions under which the enzymes can synchronize their catalytic steps, and discover that the coupling can lead to a significant enhancement in their overall catalytic rate. Both effects can be understood as arising from a global bifurcation in the underlying dynamical system at sufficiently strong coupling. Our findings suggest that, despite their molecular scale, enzymes can be cooperative and improve their performance in metabolic clusters.

Mechanical Bond Approach to Introducing Self-Adaptive Active Sites in Covalent Organic Frameworks for Zinc-Catalyzed Organophosphorus Degradation

Mechanically interlocked molecules (MIMs) with discrete molecular components linked through a mechanical bond in space can be harnessed for the operation of molecular switches and machines, which shows huge potential to imitate the dynamic response of natural enzymes. In this work, rotaxane compounds were adopted as building monomers for the synthesis of a crown-ether ring mechanically intercalated covalence organic framework (COF). This incorporation of MIMs into open architecture implemented large amplitude motions, whose wheel slid along the axle in response to external stimulation. After impregnation with Zn²⁺ ions, the relative locations of two zinc active sites (crown-ether coordinated Zn(II) and bipyridine coordinated Zn(II)) are endowed with great flexibility to fit the conformational transformation of an organophosphorus agent during the hydrolytic process. Notably, the resulting self-adaptive binuclear zinc center in a crown-ether-threaded COF network is endowed with a record catalytic ability, with a rate over 85.5 μM min^-1 for organophosphorus degradation. The strategy of synthesis for porous artificial enzymes through the introduction of mechanically bound crown ether will enable significant breakthroughs and new synthetic concepts for the development of advanced biomimetic catalysts.

Accounting for the instantaneous disorder in the enzyme-substrate Michaelis complex to calculate the Gibbs free energy barrier of an enzyme reaction

Many enzyme reactions present instantaneous disorder. These dynamic fluctuations in the enzyme-substrate Michaelis complexes generate a wide range of energy barriers that cannot be experimentally observed, but that determine the measured kinetics of the reaction. These individual energy barriers can be calculated using QM/MM methods, but then the problem is how to deal with this dispersion of energy barriers to provide kinetic information. So far, the most usual procedure has implied the so-called exponential average of the energy barriers. In this paper, we discuss the foundations of this method, and we use the free energy perturbation theory to derive an alternative equation to get the Gibbs free energy barrier of the enzyme reaction. In addition, we propose a practical way to implement it. We have chosen four enzyme reactions as examples. In particular, we have studied the hydrolysis of a glycosidic bond catalyzed by the enzyme Thermus thermophilus β-glycosidase, and the mutant Y284P Ttb-gly, and the hydrogen abstraction reactions from C13 and C7 of arachidonic acid catalyzed by the enzyme rabbit 15-lipoxygenase-1.

Sensing Enzyme Activation Heat Capacity at the Single-Molecule Level Using Gold-Nanorod-Based Optical Whispering Gallery Modes

Here, we report a label-free gold nanoparticle-based single-molecule optical platform to study the immobilization, activity, and thermodynamics of single enzymes. The sensor uses plasmonic gold nanoparticles coupled to optical whispering gallery modes (WGMs) to probe enzyme conformational dynamics during turnover at a microsecond time resolution. Using a glucosidase enzyme as the model system, we explore the temperature dependence of the enzyme turnover at the single-molecule (SM) level. A recent physical model for understanding enzyme temperature dependencies (macromolecular rate theory; MMRT) has emerged as a powerful tool to study the relationship between enzyme turnover and thermodynamics. Using WGMs, SM enzyme measurements enable us to accurately track turnover as a function of conformational changes and therefore to quantitatively probe the key feature of the MMRT model, the activation heat capacity, at the ultimate level of SM. Our data shows that WGMs are extraordinarily sensitive to protein conformational change and can discern both multiple steps with turnover as well as microscopic conformational substates within those steps. The temperature dependence studies show that the MMRT model can be applied to a range of steps within turnover at the SM scale that is associated with conformational change. Our study validates the notion that MMRT captures differences in dynamics between states. The WGM sensors provide a platform for the quantitative analysis of SM activation heat capacity, applying MMRT to the label-free sensing of microsecond substates of active enzymes.

Single enzyme electroanalysis

Traditional studies of enzymatic activity rely on the combined kinetics of millions of enzyme molecules to produce a product, an experimental approach that may wash out heterogeneities that exist between individual enzymes. Evaluating these properties on an enzyme-by-enzyme basis represents an unambiguous means of elucidating heterogeneities; however, the quantification of enzymatic activity at the single-enzyme level is fundamentally limited by the maximum catalytic rate, k_cat, inherent to a given enzyme. For electrochemical methods measuring current, single enzymes must turn over greater than 10⁷ molecules per second to produce a measurable signal on the order of 10^-12 A. Enzymes with this capability are extremely rare in nature, with typical k_cat values for biologically relevant enzymes falling between 1 and 10 000 s^-1. Thus, clever amplification strategies are necessary to electrochemically detect the vast majority of enzymes. This review details the progress toward the electroanalytical detection and evaluation of single enzyme kinetics largely focused on the nanoimpact method, a chronoamperometric detection strategy that monitors the change in the current-time profile associated with stochastic collisions of freely diffusing entities (e.g., enzymes) onto a microelectrode or nanoelectrode surface. We discuss the experimental setups and methods developed in the last decade toward the quantification of single molecule enzymatic rates. Special emphasis is given to the limitations of measurement science in the observation of single enzyme activity and feasible methods of signal amplification with reasonable bandwidth.

Multiplex Single-Molecule Kinetics of Nanopore-Coupled Polymerases

DNA polymerases have revolutionized the biotechnology field due to their ability to precisely replicate stored genetic information. Screening variants of these enzymes for specific properties gives the opportunity to identify polymerases with different features. We have previously developed a single-molecule DNA sequencing platform by coupling a DNA polymerase to an α-hemolysin pore on a nanopore array. Here, we use this approach to demonstrate a single-molecule method that enables rapid screening of polymerase variants in a multiplex manner. In this approach, barcoded DNA strands are complexed with polymerase variants and serve as templates for nanopore sequencing. Nanopore sequencing of the barcoded DNA reveals both the barcode identity and kinetic properties of the polymerase variant associated with the cognate barcode, allowing for multiplexed investigation of many polymerase variants in parallel on a single nanopore array. Further, we develop a robust classification algorithm that discriminates kinetic characteristics of the different polymerase mutants. As a proof of concept, we demonstrate the utility of our approach by screening a library of ∼100 polymerases to identify variants for potential applications of biotechnological interest. We anticipate our screening method to be broadly useful for applications that require polymerases with altered physical properties.

Drug discovery for enzymes

Enzymes are essential, physiological catalysts involved in all processes of life, including metabolism, cellular signaling and motility, as well as cell growth and division. They are attractive drug targets because of the presence of defined substrate-binding pockets, which can be exploited as binding sites for pharmaceutical enzyme inhibitors. Understanding the reaction mechanisms of enzymes and the molecular mode of action of enzyme inhibitors is indispensable for the discovery and development of potent, efficacious, and safe novel drugs. The combination of classical concepts of enzymology with new experimental and data analysis methods opens new routes for drug discovery.

Importance learning estimator for the site-averaged turnover frequency of a disordered solid catalyst

For disordered catalysts such as atomically dispersed "single-atom" metals on amorphous silica, the active sites inherit different properties from their quenched-disordered local environments. The observed kinetics are site-averages, typically dominated by a small fraction of highly active sites. Standard sampling methods require expensive ab initio calculations at an intractable number of sites to converge on the site-averaged kinetics. We present a new method that efficiently estimates the site-averaged turnover frequency (TOF). The new estimator uses the same importance learning algorithm [Vandervelden et al., React. Chem. Eng. 5, 77 (2020)] that we previously used to compute the site-averaged activation energy. We demonstrate the method by computing the site-averaged TOF for a simple disordered lattice model of an amorphous catalyst. The results show that with the importance learning algorithm, the site-averaged TOF and activation energy can now be obtained concurrently with orders of magnitude reduction in required ab initio calculations.

Fluoroacetate dehalogenase catalyzed dehalogenation of halogenated carboxylic acids: A QM/MM approach

Dehalogenation is one of the most important reactions in environmental pollution control, for instance, the degradation of persistent organic pollutants (POPs). Recently, fluoroacetate dehalogenase (FAcD) has been reported to catalyze the dehalogenation reactions, which shows great potential in treating halogenated pollutants. Here the dehalogenation mechanism catalyzed by FAcD was fully deciphered with the aid of quantum mechanics/molecular mechanics method. The results show that FAcD catalyzed dehalogenation efficiency follows the order of defluorination > dechlorination > debromination. The corresponding Boltzmann-weighted average barriers are 10.1, 19.7, and 20.9 kcal mol^-1. Positive/negative correlations between activation barriers and structural parameters (e.g. distance and angle) for FAcD catalyzed dechlorination and debromination were established. Based on the structure-energy relationship, we propose that mutation of the binding pocket amino acids (e.g. His155, Trp156, Tyr219) to smaller proton donor amino acids (e.g. Serine, Threonine, Cysteine, Asparagine) may increase the efficiency for dechlorination and debromination. The results may of practical value for the efficient degradation of chlorined and bromined pollutants by harnessing FAcD.

Catalysis and Electron Transfer in De Novo Designed Helical Scaffolds

The relationship between protein structure and function is one of the greatest puzzles within biochemistry. De novo metalloprotein design is a way to wipe the board clean and determine what is required to build in function from the ground up in an unrelated structure. This Review focuses on protein design efforts to create de novo metalloproteins within alpha-helical scaffolds. Examples of successful designs include those with carbonic anhydrase or nitrite reductase activity by incorporating a ZnHis3 or CuHis3 site, or that recapitulate the spectroscopic properties of unique electron-transfer sites in cupredoxins (CuHis2 Cys) or rubredoxins (FeCys4 ). This work showcases the versatility of alpha helices as scaffolds for metalloprotein design and the progress that is possible through careful rational design. Our studies cover the invariance of carbonic anhydrase activity with different site positions and scaffolds, refinement of our cupredoxin models, and enhancement of nitrite reductase activity up to 1000-fold.

High-Fidelity Single Molecule Quantification in a Flow Cytometer Using Multiparametric Optical Analysis

Microfluidic techniques are widely used for high-throughput quantification and discrete analysis of micron-scale objects but are difficult to apply to molecular-scale targets. Instead, single-molecule methods primarily rely on low-throughput microscopic imaging of immobilized molecules. Here we report that commercial-grade flow cytometers can detect single nucleic acid targets following enzymatic extension and dense labeling with multiple distinct fluorophores. We focus on microRNAs, short nucleic acids that can be extended by rolling circle amplification (RCA). We labeled RCA-extended microRNAs with multicolor fluorophores to generate repetitive nucleic acid products with submicron sizes and tunable multispectral profiles. By cross-correlating the multiparametric optical features, signal-to-background ratios were amplified 1600-fold to allow single-molecule detection across 4 orders of magnitude of concentration. The limit of detection was measured to be 47 fM, which is 100-fold better than gold-standard methods based on polymerase chain reaction. Furthermore, multiparametric analysis allowed discrimination of different microRNA sequences in the same solution using distinguishable optical barcodes. Barcodes can apply both ratiometric and colorimetric signatures, which could facilitate high-dimensional multiplexing. Because of the wide availability of flow cytometers, we anticipate that this technology can provide immediate access to high-throughput multiparametric single-molecule measurements and can further be adapted to the diverse range of molecular amplification methods that are continually emerging.

Harnessing fluoroacetate dehalogenase for defluorination of fluorocarboxylic acids: in silico and in vitro approach

Widely distributed fluorocarboxylic acids have aroused worldwide environmental concerns due to its toxicity, persistence, and bioaccumulation. Enzyme-based eco-friendly biodegradation techniques have become increasingly important in treating fluorocarboxylic acids. Here we utilized in silico and in vitro approaches to investigate the defluorination mechanism of fluoroacetate dehalogenase (FAcD) toward monofluoropropionic acids at atomic-level. The experimentally determined k_cat and k_M for defluorination of 2-fluoropropionic acid are 330 ± 60 min^-1 and 6.12 ± 0.13 mM. The in silico results demonstrated positive/negative correlations between activation barriers and structural parameters (e.g. distance and angle) under different enzymatic conformations. We also screened computationally and tested in vitro (enzyme assay and kinetic study) the catalytic proficiency of FAcD toward polyfluoropropionic acids and perfluoropropionic acids which are known to be challenging for enzymatic degradation. The results revealed potential degradation activity of FAcD enzyme toward 2,3,3,3-tetrafluoropropionic acids. Our work will initiate the development of a new "integrated approach" for enzyme engineering to degrade environmentally persistent fluorocarboxylic acids.

Conformational diversity induces nanosecond-timescale chemical disorder in the HIV-1 protease reaction pathway

The role of conformational diversity in enzyme catalysis has been a matter of analysis in recent studies. Pre-organization of the active site has been pointed out as the major source for enzymes' catalytic power. Following this line of thought, it is becoming clear that specific, instantaneous, non-rare enzyme conformations that make the active site perfectly pre-organized for the reaction lead to the lowest activation barriers that mostly contribute to the macroscopically observed reaction rate. The present work is focused on exploring the relationship between structure and catalysis in HIV-1 protease (PR) with an adiabatic mapping method, starting from different initial structures, collected from a classical MD simulation. The first, rate-limiting step of the HIV-1 PR catalytic mechanism was studied with the ONIOM QM/MM methodology (B3LYP/6-31G(d):ff99SB), with activation and reaction energies calculated at the M06-2X/6-311++G(2d,2p):ff99SB level of theory, in 19 different enzyme:substrate conformations. The results showed that the instantaneous enzyme conformations have two independent consequences on the enzyme's chemistry: they influence the barrier height, something also observed in the past in other enzymes, and they also influence the specific reaction pathway, which is something unusual and unexpected, challenging the "one enzyme-one substrate-one reaction mechanism" paradigm. Two different reaction mechanisms, with similar reactant probabilities and barrier heights, lead to the same gem-diol intermediate. Subtle nanosecond-timescale rearrangements in the active site hydrogen bonding network were shown to determine which reaction the enzyme follows. We named this phenomenon chemical disorder. The results make us realize the unexpected mechanistic consequences of conformational diversity in enzymatic reactivity.

Substrate binding versus escape dynamics in a pH-affected fungal beta-glucosidase revealed by molecular dynamics simulations

The cellulolytic ability of fungal species is important to both natural and engineered biocycling of plant matter. One essential step is the conversion of cellobiose into glucose catalyzed by beta-glucosidases. Mutagenesis studies have implicated altering the substrate binding pocket to influence the pH-activity profile of this enzyme. However, structural understanding of the pH-affected substrate binding environment is lacking. Here we conducted molecular dynamics simulations of fully hydrated TrBgl2, a beta-glucosidase of Trichoderma reesei, equilibrated at its optimal pH (pH 6) and two unfavorable pHs (pH 5 and pH 7.5). We identified structural arrangement of specific residues that facilitated substrate escape from the catalytic site at pH 5 but locked the bound substrate in an unfavorable orientation at pH 7.5. For comparative analysis, we also performed simulations of a mutated TrBgl2 with previously demonstrated improved catalysis as a function of pH. We captured the responsible conformational changes in the engineered substrate binding pocket.

Thermodynamics of Reactions Affected by Medium Reorganization

We present a thermodynamic analysis of the activation barrier for reactions which can be monitored through the difference in the energies of reactants and products defined as the reaction coordinate (electron and atom transfer, enzyme catalysis, etc.). The free-energy surfaces along the reaction coordinate are separated into the enthalpy and entropy surfaces. For the Gaussian statistics of the reaction coordinate, the free-energy surfaces are parabolas, and the entropy surface is an inverted parabola. Its maximum coincides with the transition state for reactions with zero value of the reaction free energy. Maximum entropic depression of the activation barrier, anticipated by the concept of transition-state ensembles, can be achieved for such reactions. From Onsager's reversibility, the entropy of equilibrium fluctuations encodes the entropic component of the activation barrier. The reorganization entropy thus becomes the critical parameter of the theory reducing the problem of activation entropy to the problem of reorganization entropy. Standard solvation theories do not allow reorganization entropy sufficient for the barrier depression. Complex media, characterized by many relaxation processes, need to be involved. Proteins provide several routes for achieving large entropic effects through incomplete (nonergodic) sampling of the complex energy landscape and by facilitating an active role of water in the reaction mechanism.

Insights into the catalytic mechanism of dehydrogenase BphB: A quantum mechanics/molecular mechanics study

The present study delineated the dehydrogenation mechanism of cis-2,3-dihydro-2,3-dihydroxybiphenyl (2,3-DDBPH) and cis-2,3-dihydro-2,3-dihydroxy-4,4'-dichlorobiphenyl (2,3-DD-4,4'-DBPH) by Pandoraea pnomenusa strain B-356 cis-2,3-dihydro-2,3-dihydroxybiphenyl dehydrogenase (BphB) in atomistic detail. The enzymatic process was investigated by a combined quantum mechanics/molecular mechanics (QM/MM) approach. Five different snapshots were extracted and calculated, which revealed that the Boltzmann-weighted average barriers of 2,3-DDBPH and 2,3-DD-4,4'-DBPH dehydrogenation processes are 10.7 and 11.5 kcal mol^-1, respectively. The established dehydrogenation mechanism provides new insight into the degradation processes of other chlorinated 2,3-DDBPH. In addition to Asn115, Ser142, and Lys149, the importance of Ile 89, Asn143, Pro184, Met 187, Thr189, and Lue 191 during the dehydrogenation process of 2,3-DDBPH and 2,3-DD-4,4'-DBPH were also highlighted to search for promising mutation targets for improving the catalytic efficiency of BphB.

Real-time tracking of single-molecule collagenase on native collagen and partially structured collagen-mimic substrates

The dynamic interactions of an individual matrix metalloproteinase-1 were imaged and monitored in the presence of either triple-helical or non-triple-helical, partially structured collagen-mimic substrates. The enzyme exhibited ten-fold increased catalytic turnover rates with the structurally modified substrate by skipping the triple-helix unwinding step during the catalytic pathway.

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.

Correlating Structural and Functional Heterogeneity of Immobilized Enzymes

Many nanobiotechnology applications rely on stable and efficient integration of functional biomacromolecules with synthetic nanomaterials. Unfortunately, the reasons for the ubiquitous loss of activity of immobilized enzymes remain poorly understood due to the difficulty in distinguishing between distinct molecular-level mechanisms. Here, we employ complementary single-molecule fluorescence methods that independently measure the impact of immobilization on the structure and function ( i. e., substrate binding kinetics) of nitroreductase (NfsB). Stochastic statistical modeling methods were used to unambiguously quantify the effects of immobilized NfsB structural dynamics on function, allowing us to explicitly separate effects due to conformation and accessibility. Interestingly, we found that nonspecifically tethered NfsB exhibited enhanced stability compared to site-specifically tethered NfsB; however, the folded state of site-specifically tethered NfsB had faster substrate binding rates, suggesting improved active site accessibility. This demonstrated an unexpected intrinsic trade-off associated with competing bioconjugation methods, suggesting that it may be necessary to balance conformational stability versus active site accessibility. This nuanced view of the impact of immobilization will facilitate a rational approach to the integration of enzymes with synthetic nanomaterials.

Computational evidence for the degradation mechanism of haloalkane dehalogenase LinB and mutants of Leu248 to 1-chlorobutane

The catalytic degradation ability of the haloalkane dehalogenase LinB toward 1-chlorobutane (1-CB) was studied using a combined quantum mechanics/molecular mechanics (QM/MM) approach. Two major processes are involved in the LinB-catalyzed removal of halogens: dechlorination and hydrolyzation. The present study confirmed the experimentally proposed reaction path at the molecular level. Moreover, based on nucleophilic substitution mechanism (SN2 reaction), dechlorination was found to be the rate-determining step of the entire reaction process. In this study, the Boltzmann-weighted average barrier for dechlorination was determined to be 17.0 kcal mol-1, which is fairly close to the experimental value (17.4 kcal mol-1). The state of His107 and the influence of Leu248 on the dechlorination process were also explored. In addition, an intriguing phenomenon was discovered: the potential energy barrier decreased by 7.5 kcal mol-1 when the Leu248 residue was mutated into Phe248. This discovery might be of great help to design new mutant enzymes or novel biocatalysts.

Real-time monitoring of conformational transitions of single-molecule histone deacetylase 8 with nanocircuits

Using single-molecule approaches, we directly observed the dynamic interaction between HDAC8 and various ligands as well as conformational interconversions during the catalytic reaction. Statistical analysis identified key kinetic parameters, demonstrating that the enzymatic activity is highly sensitive to both minor variations in the ligand structures and small synthetic molecules.

Electrochemical Zero-Mode Waveguide Studies of Single Enzyme Reactions

Because electron transfer reactions are fundamental to life processes, such as respiration, vision, and energy catabolism, it is critically important to understand the relationship between functional states of individual redox enzymes and the macroscopically observed phenotype, which results from averaging over all copies of the same enzyme. To address this problem, we have developed a new technology, based on a bifunctional nanoelectrochemical-nanophotonic architecture - the electrochemical zero mode waveguide (E-ZMW) - that can couple biological electron transfer reactions to luminescence, making it possible to observe single electron transfer events in redox enzymes. Here we describe E-ZMW architectures capable of supporting potential-controlled redox reactions with single copies of the oxidoreductase enzyme, glutathione reductase, GR, and extend these capabilities to electron transfer events where reactive oxygen species are synthesized within the ~100 zL volume of the nanopore.

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.