Relating Protein Motion to Catalysis

A Biodegradable Secondary Battery and its Biodegradation Mechanism for Eco-Friendly Energy-Storage Systems

The production of rechargeable batteries is rapidly expanding, and there are going to be new challenges in the near future about how the potential environmental impact caused by the disposal of the large volume of the used batteries can be minimized. Herein, a novel strategy is proposed to address these concerns by applying biodegradable device technology. An eco-friendly and biodegradable sodium-ion secondary battery (SIB) is developed through extensive material screening followed by the synthesis of biodegradable electrodes and their seamless assembly with an unconventional biodegradable separator, electrolyte, and package. Each battery component decomposes in nature into non-toxic compounds or elements via hydrolysis and/or fungal degradation, with all of the biodegradation products naturally abundant and eco-friendly. Detailed biodegradation mechanisms and toxicity influence of each component on living organisms are determined. In addition, this new SIB delivers performance comparable to that of conventional non-degradable SIBs. The strategy and findings suggest a novel eco-friendly biodegradable paradigm for large-scale rechargeable battery systems.

A Single-Point Mutation in d-Arginine Dehydrogenase Unlocks a Transient Conformational State Resulting in Altered Cofactor Reactivity

Proteins are inherently dynamic, and proper enzyme function relies on conformational flexibility. In this study, we demonstrated how an active site residue changes an enzyme's reactivity by modulating fluctuations between conformational states. Replacement of tyrosine 249 (Y249) with phenylalanine in the active site of the flavin-dependent d-arginine dehydrogenase yielded an enzyme with both an active yellow FAD (Y249F-y) and an inactive chemically modified green FAD, identified as 6-OH-FAD (Y249F-g) through various spectroscopic techniques. Structural investigation of Y249F-g and Y249F-y variants by comparison to the wild-type enzyme showed no differences in the overall protein structure and fold. A closer observation of the active site of the Y249F-y enzyme revealed an alternative conformation for some active site residues and the flavin cofactor. Molecular dynamics simulations probed the alternate conformations observed in the Y249F-y enzyme structure and showed that the enzyme variant with FAD samples a metastable conformational state, not available to the wild-type enzyme. Hybrid quantum/molecular mechanical calculations identified differences in flavin electronics between the wild type and the alternate conformation of the Y249F-y enzyme. The computational studies further indicated that the alternate conformation in the Y249F-y enzyme is responsible for the higher spin density at the C6 atom of flavin, which is consistent with the formation of 6-OH-FAD in the variant enzyme. The observations in this study are consistent with an alternate conformational space that results in fine-tuning the microenvironment around a versatile cofactor playing a critical role in enzyme function.

Flexibility of enzymatic transitions as a hallmark of optimized enzyme steady-state kinetics and thermodynamics

We investigate the relations between the enzyme kinetic flexibility, the rate of entropy production, and the Shannon information entropy in a steady-state enzyme reaction. All these quantities are maximized with respect to enzyme rate constants. We show that the steady-state, which is characterized by the most flexible enzymatic transitions between the enzyme conformational states, coincides with the global maxima of the Shannon information entropy and the rate of entropy production. This steady-state of an enzyme is referred to as globally optimal. This theoretical approach is then used for the analysis of the kinetic and the thermodynamic performance of the enzyme triose-phosphate isomerase. The analysis reveals that there exist well-defined maxima of the kinetic flexibility, the rate of entropy production, and the Shannon information entropy with respect to any arbitrarily chosen rate constant of the enzyme and that these maxima, calculated from the measured kinetic rate constants for the triose-phosphate isomerase are lower, however of the same order of magnitude, as the maxima of the globally optimal state of the enzyme. This suggests that the triose-phosphate isomerase could be a well, but not fully evolved enzyme, as it was previously claimed. Herein presented theoretical investigations also provide clear evidence that the flexibility of enzymatic transitions between the enzyme conformational states is a requirement for the maximal Shannon information entropy and the maximal rate of entropy production.

Towards resolving the complex paramagnetic nuclear magnetic resonance (NMR) spectrum of small laccase: assignments of resonances to residue-specific nuclei

Laccases efficiently reduce dioxygen to water in an active site containing a tri-nuclear copper centre (TNC). The dynamics of the protein matrix is a determining factor in the efficiency in catalysis. To probe mobility, nuclear magnetic resonance (NMR) spectroscopy is highly suitable. However, several factors complicate the assignment of resonances to active site nuclei in laccases. The paramagnetic nature causes large shifts and line broadening. Furthermore, the presence of slow chemical exchange processes of the imidazole rings of copper ligand results in peak doubling. A third complicating factor is that the enzyme occurs in two states, the native intermediate (NI) and resting oxidized (RO) states, with different paramagnetic properties. The present study aims at resolving the complex paramagnetic NMR spectra of the TNC of Streptomyces coelicolor small laccase (SLAC). With a combination of paramagnetically tailored NMR experiments, all eight His Nδ 1 and Hδ 1 resonances for the NI state are identified, as well as His Hβ protons for the RO state. With the help of second-shell mutagenesis, selective resonances are tentatively assigned to the histidine ligands of the copper in the type-2 site. This study demonstrates the utility of the approaches used for the sequence-specific assignment of the paramagnetic NMR spectra of ligands in the TNC that ultimately may lead to a description of the underlying motion.

Acceleration of catalysis in dihydrofolate reductase by transient, site-specific photothermal excitation

We have studied the role of protein dynamics in chemical catalysis in the enzyme dihydrofolate reductase (DHFR), using a pump-probe method that employs pulsed-laser photothermal heating of a gold nanoparticle (AuNP) to directly excite a local region of the protein structure and transient absorbance to probe the effect on enzyme activity. Enzyme activity is accelerated by pulsed-laser excitation when the AuNP is attached close to a network of coupled motions in DHFR (on the FG loop, containing residues 116-132, or on a nearby alpha helix). No rate acceleration is observed when the AuNP is attached away from the network (distal mutant and His-tagged mutant) with pulsed excitation, or for any attachment site with continuous wave excitation. We interpret these results within an energy landscape model in which transient, site-specific addition of energy to the enzyme speeds up the search for reactive conformations by activating motions that facilitate this search.

Protein Mass-Modulated Effects in Alkaline Phosphatase

Recent experimental studies engaging isotopically substituted protein (heavy protein) have revealed that many, but not all, enzymatic systems exhibit altered chemical steps in response to an altered mass. The results have been interpreted as femtosecond protein dynamics at the active site being linked (or not) to transition-state barrier crossing. An altered enzyme mass can influence several kinetic parameters (kcat, Km, and kchem) in amounts of ≤30% relative to light enzymes. An early report on deuterium-labeled Escherichia coli alkaline phosphatase (AP) showed an unusually large enzyme kinetic isotope effect on kcat. We examined steady-state and chemical step properties of native AP, [2H]AP, and [2H,13C,15N]AP to characterize the role of heavy enzyme protein dynamics in reactions catalyzed by AP. Both [2H]- and [2H,13C,15N]APs showed unaltered steady-state and single-turnover rate constants. These findings characterize AP as one of the enzymes in which mass-dependent catalytic site dynamics is dominated by reactant-linked atomic motions. Two catalytic site zinc ions activate the oxygen nucleophiles in the catalytic site of AP. The mass of the zinc ions is unchanged in light and heavy APs. They are essentially linked to catalysis and provide a possible explanation for the loss of linkage between catalysis and protein mass in these enzymes.

Low-Frequency Protein Motions Coupled to Catalytic Sites

This review examines low-frequency vibrational modes of proteins and their coupling to enzyme catalytic sites. That protein motions are critical to enzyme function is clear, but the kinds of motions present in proteins and how they are involved in function remain unclear. Several models of enzyme-catalyzed reaction suggest that protein dynamics may be involved in the chemical step of the catalyzed reaction, but the evidence in support of such models is indirect. Spectroscopic studies of low-frequency protein vibrations consistently show that there are underdamped modes of the protein with frequencies in the tens of wavenumbers where overdamped behavior would be expected. Recent studies even show that such underdamped vibrations modulate enzyme active sites. These observations suggest that increasingly sophisticated spectroscopic methods will be able to unravel the link between low-frequency protein vibrations and enzyme function.

Rhodium catalysts with cofactor mimics for the biomimetic reduction of CN bonds

Bio-inspired reduction of CN bonds was successfully performed using rhodium catalysts containing cofactor mimics. The intramolecular cooperation between rhodium and cofactor mimics enabled the transformation with good selectivity. A plausible mechanism was also proposed.

How directed evolution reshapes the energy landscape in an enzyme to boost catalysis

The advent of biocatalysts designed computationally and optimized by laboratory evolution provides an opportunity to explore molecular strategies for augmenting catalytic function. Applying a suite of nuclear magnetic resonance, crystallography, and stopped-flow techniques to an enzyme designed for an elementary proton transfer reaction, we show how directed evolution gradually altered the conformational ensemble of the protein scaffold to populate a narrow, highly active conformational ensemble and accelerate this transformation by nearly nine orders of magnitude. Mutations acquired during optimization enabled global conformational changes, including high-energy backbone rearrangements, that cooperatively organized the catalytic base and oxyanion stabilizer, thus perfecting transition-state stabilization. The development of protein catalysts for many chemical transformations could be facilitated by explicitly sampling conformational substates during design and specifically stabilizing productive substates over all unproductive conformations.

Characterization of the ExoU activation mechanism using EPR and integrative modeling

ExoU, a type III secreted phospholipase effector of Pseudomonas aeruginosa, serves as a prototype to model large, dynamic, membrane-associated proteins. ExoU is synergistically activated by interactions with membrane lipids and ubiquitin. To dissect the activation mechanism, structural homology was used to identify an unstructured loop of approximately 20 residues in the ExoU amino acid sequence. Mutational analyses indicate the importance of specific loop amino acid residues in mediating catalytic activity. Engineered disulfide cross-links show that loop movement is required for activation. Site directed spin labeling EPR and DEER (double electron-electron resonance) studies of apo and holo states demonstrate local conformational changes at specific sites within the loop and a conformational shift of the loop during activation. These data are consistent with the formation of a substrate-binding pocket providing access to the catalytic site. DEER distance distributions were used as constraints in RosettaDEER to construct ensemble models of the loop in both apo and holo states, significantly extending the range for modeling a conformationally dynamic loop.

Enzyme dynamics: Looking beyond a single structure

Conventional understanding of how enzymes function strongly emphasizes the role of structure. However, increasing evidence clearly indicates that enzymes do not remain fixed or operate exclusively in or close to their native structure. Different parts of the enzyme (from individual residues to full domains) undergo concerted motions on a wide range of time-scales, including that of the catalyzed reaction. Information obtained on these internal motions and conformational fluctuations has so far uncovered and explained many aspects of enzyme mechanisms, which could not have been understood from a single structure alone. Although there is wide interest in understanding enzyme dynamics and its role in catalysis, several challenges remain. In addition to technical difficulties, the vast majority of investigations are performed in dilute aqueous solutions, where conditions are significantly different than the cellular milieu where a large number of enzymes operate. In this review, we discuss recent developments, several challenges as well as opportunities related to this topic. The benefits of considering dynamics as an integral part of the enzyme function can also enable new means of biocatalysis, engineering enzymes for industrial and medicinal applications.

Kramers’ Theory and the Dependence of Enzyme Dynamics on Trehalose-Mediated Viscosity

The disaccharide trehalose is accumulated in the cytoplasm of some organisms in response to harsh environmental conditions. Trehalose biosynthesis and accumulation are important for the survival of such organisms by protecting the structure and function of proteins and membranes. Trehalose affects the dynamics of proteins and water molecules in the bulk and the protein hydration shell. Enzyme catalysis and other processes dependent on protein dynamics are affected by the viscosity generated by trehalose, as described by the Kramers’ theory of rate reactions. Enzyme/protein stabilization by trehalose against thermal inactivation/unfolding is also explained by the viscosity mediated hindering of the thermally generated structural dynamics, as described by Kramers’ theory. The analysis of the relationship of viscosity–protein dynamics, and its effects on enzyme/protein function and other processes (thermal inactivation and unfolding/folding), is the focus of the present work regarding the disaccharide trehalose as the viscosity generating solute. Finally, trehalose is widely used (alone or in combination with other compounds) in the stabilization of enzymes in the laboratory and in biotechnological applications; hence, considering the effect of viscosity on catalysis and stability of enzymes may help to improve the results of trehalose in its diverse uses/applications.

Identification of Alternative Allosteric Sites in Glycolytic Enzymes for Potential Use as Species-Specific Drug Targets

Three allosteric glycolytic enzymes, phosphofructokinase, glyceraldehyde-3 phosphate dehydrogenase and pyruvate kinase, associated with bacterial, parasitic and human species, were explored to identify potential allosteric sites that would be used as prime targets for species-specific drug design purposes using a newly developed approach which incorporates solvent mapping, elastic network modeling, sequence and structural alignments. The majority of binding sites detected by solvent mapping overlapped with the interface regions connecting the subunits, thus appeared as promising target sites for allosteric regulation. Each binding site was then evaluated by its ability to alter the global dynamics of the receptor defined by the percentage change in the frequencies of the lowest-frequency modes most significantly and as anticipated, the most effective ones were detected in the vicinity of the well-reported catalytic and allosteric sites. Furthermore, some of our proposed regions intersected with experimentally resolved sites which are known to be critical for activity regulation, which further validated our approach. Despite the high degree of structural conservation encountered between bacterial/parasitic and human glycolytic enzymes, the majority of the newly presented allosteric sites exhibited a low degree of sequence conservation which further increased their likelihood to be used as species-specific target regions for drug design studies.

Dissecting the Dynamics during Enzyme Catalysis: A Case Study of Pin1 Peptidyl-Prolyl Isomerase

Free energy surfaces have played a central role in studying protein conformational changes and enzymatic reactions over decades. Yet, free energy barriers and kinetics are highly dependent on the coordinates chosen to define the surface, and furthermore, the dynamics during the reactions are often overlooked. Our recent study on the Pin1-catalyzed isomerization reaction has indicated that the isomerization transition events remarkably deviate from the free energy path, highlighting the need to understand the reaction dynamics in more detail. To this end, here we investigate the reaction coordinates that describe the transition states of the free energy and transition pathways by minimizing the cross-entropy function. We show that the isomerization transition events can be expressed by the concerted changes in the improper torsion angle ζ and nearby backbone torsional angles of the ligand, whereas the transition state of the free energy surface involves changes in a broad range of coordinates including multiple protein-ligand interactions. The current result supports the previous finding that the isomerization transitions occur quickly from the conformational excited states, which is in sharp contrast to the slow and collective changes suggested from the free energy path. Our results further indicate that the coordinates derived from the transition trajectories are not sufficient for finding the transition states on the free energy surfaces due to the lack of information from conformational excited states.

Proteins as diverse, efficient, and evolvable scaffolds for artificial metalloenzymes

We have extracted and categorized the desirable properties of proteins that are adapted as the scaffolds for artificial metalloenzymes.

Molecular Oxygen Binding in the Mitochondrial Electron Transfer Flavoprotein

Reactive oxygen species such as superoxide are potentially harmful byproducts of the aerobic metabolism in the inner mitochondrial membrane, and complexes I, II, III of the electron transport chain have been identified as primary sources. The mitochondrial fatty acid b-oxidation pathway may also play a yet uncharacterized role in reactive oxygen species generation, apparently at the level of the electron transfer flavoprotein:ubiquinone oxidoreductase (ETF:QO) and/or its redox partner electron-transfer flavoprotein (ETF). These enzymes comprise a key pathway through which electrons are sequentially shuttled from several dehydrogenases to the respiratory chain. The exact mechanisms of superoxide production have not been fully established, but a crucial starting point would be the binding of molecular oxygen within one of the protein complexes. The present investigation offers a comprehensive computational approach for the determination of binding modes and characteristic binding times of small molecules inside proteins, which is then used to reveal several O₂ binding sites near the flavin adenine dinucleotide cofactor of the ETF enzyme. The binding sites are further characterized to extract the necessary parameters for further studies of possible electron transfer between flavin and O₂ leading to radical pair formation and possible superoxide production.

Evolution Conserves the Network of Coupled Residues in Dihydrofolate Reductase

Understanding protein motions and their role in enzymatic reactions is an important and timely topic in enzymology. Protein motions that are involved in the chemical step of catalysis are particularly intriguing but difficult to identify. A global network of coupled residues in Escherichia coli dihydrofolate reductase (E. coli DHFR), which assists in catalyzing the chemical step, has previously been demonstrated through quantum mechanical/molecular mechanical and molecular dynamics simulations as well as bioinformatic analyses. A few specific residues (M42, G121, F125, and I14) were shown to function synergistically with measurements of single-turnover rates and the temperature dependence of intrinsic kinetic isotope effects (KIEsint) of site-directed mutants. This study hypothesizes that the global network of residues involved in the chemical step is evolutionarily conserved and probes homologous residues of the potential global network in human DHFR through measurements of the temperature dependence of KIEsint and computer simulations based on the empirical valence bond method. We study mutants M53W and S145V. Both of these remote residues are homologous to network residues in E. coli DHFR. Non-additive isotope effects on activation energy are observed between M53 and S145, indicating their synergistic effect on the chemical step in human DHFR, which suggests that both of these residues are part of a network affecting the chemical step in enzyme catalysis. This finding supports the hypothesis that human and E. coli DHFR share similar networks, consistent with evolutionary preservation of such networks.

Structural basis of glycogen metabolism in bacteria

The evolution of metabolic pathways is a major force behind natural selection. In the spotlight of such process lies the structural evolution of the enzymatic machinery responsible for the central energy metabolism. Specifically, glycogen metabolism has emerged to allow organisms to save available environmental surplus of carbon and energy, using dedicated glucose polymers as a storage compartment that can be mobilized at future demand. The origins of such adaptive advantage rely on the acquisition of an enzymatic system for the biosynthesis and degradation of glycogen, along with mechanisms to balance the assembly and disassembly rate of this polysaccharide, in order to store and recover glucose according to cell energy needs. The first step in the classical bacterial glycogen biosynthetic pathway is carried out by the adenosine 5'-diphosphate (ADP)-glucose pyrophosphorylase. This allosteric enzyme synthesizes ADP-glucose and acts as a point of regulation. The second step is carried out by the glycogen synthase, an enzyme that generates linear α-(1→4)-linked glucose chains, whereas the third step catalyzed by the branching enzyme produces α-(1→6)-linked glucan branches in the polymer. Two enzymes facilitate glycogen degradation: glycogen phosphorylase, which functions as an α-(1→4)-depolymerizing enzyme, and the debranching enzyme that catalyzes the removal of α-(1→6)-linked ramifications. In this work, we rationalize the structural basis of glycogen metabolism in bacteria to the light of the current knowledge. We describe and discuss the remarkable progress made in the understanding of the molecular mechanisms of substrate recognition and product release, allosteric regulation and catalysis of all those enzymes.

Formation of Unstable and very Reactive Chemical Species Catalyzed by Metalloenzymes: A Mechanistic Overview

Nature has tailored a wide range of metalloenzymes that play a vast array of functions in all living organisms and from which their survival and evolution depends on. These enzymes catalyze some of the most important biological processes in nature, such as photosynthesis, respiration, water oxidation, molecular oxygen reduction, and nitrogen fixation. They are also among the most proficient catalysts in terms of their activity, selectivity, and ability to operate at mild conditions of temperature, pH, and pressure. In the absence of these enzymes, these reactions would proceed very slowly, if at all, suggesting that these enzymes made the way for the emergence of life as we know today. In this review, the structure and catalytic mechanism of a selection of diverse metalloenzymes that are involved in the production of highly reactive and unstable species, such as hydroxide anions, hydrides, radical species, and superoxide molecules are analyzed. The formation of such reaction intermediates is very difficult to occur under biological conditions and only a rationalized selection of a particular metal ion, coordinated to a very specific group of ligands, and immersed in specific proteins allows these reactions to proceed. Interestingly, different metal coordination spheres can be used to produce the same reactive and unstable species, although through a different chemistry. A selection of hand-picked examples of different metalloenzymes illustrating this diversity is provided and the participation of different metal ions in similar reactions (but involving different mechanism) is discussed.

An extreme halophilic xylanase from camel rumen metagenome with elevated catalytic activity in high salt concentrations

An extreme halophilic xylanase, designated as XylCMS, was characterized by cloning and expression of the encoding gene from a camel rumen metagenome. XylCMS proved to be a GH11 xylanase with high identity to a hypothetical glycosyl hydrolase from Ruminococcus flavefaciens. XylCMS with a molecular weight of about 47 kDa showed maximum activity at pH 6 and 55 °C. The enzyme activity was significantly stimulated by NaCl in 1–5 M concentrations. Interestingly, the optimum temperature was not influenced by NaCl but the K_cat of the enzyme was enhanced by 2.7-folds at 37 °C and 1.2-folds at 55 °C. The K_m value was decreased with NaCl by 4.3-folds at 37 °C and 3.7-folds at 55 °C resulting in a significant increase in catalytic efficiency (K_cat/K_m) by 11.5-folds at 37 °C and 4.4-folds at 55 °C. Thermodynamic analysis indicated that the activation energy (E_a) and enthalpy (∆H) of the reaction were decreased with NaCl by 2.4 and threefold, respectively. From the observations and the results of fluorescence spectroscopy, it was concluded that NaCl at high concentrations improves both the flexibility and substrate affinity of XylCMS that are crucial for catalytic activity by influencing substrate binding, product release and the energy barriers of the reaction. XylCMS as an extreme halophilic xylanase with stimulated activity in artificial seawater and low water activity conditions has potentials for application in industrial biotechnology.

Discovery of processive catalysis by an exo-hydrolase with a pocket-shaped active site

Substrates associate and products dissociate from enzyme catalytic sites rapidly, which hampers investigations of their trajectories. The high-resolution structure of the native Hordeum exo-hydrolase HvExoI isolated from seedlings reveals that non-covalently trapped glucose forms a stable enzyme-product complex. Here, we report that the alkyl β-d-glucoside and methyl 6-thio-β-gentiobioside substrate analogues perfused in crystalline HvExoI bind across the catalytic site after they displace glucose, while methyl 2-thio-β-sophoroside attaches nearby. Structural analyses and multi-scale molecular modelling of nanoscale reactant movements in HvExoI reveal that upon productive binding of incoming substrates, the glucose product modifies its binding patterns and evokes the formation of a transient lateral cavity, which serves as a conduit for glucose departure to allow for the next catalytic round. This path enables substrate-product assisted processive catalysis through multiple hydrolytic events without HvExoI losing contact with oligo- or polymeric substrates. We anticipate that such enzyme plasticity could be prevalent among exo-hydrolases.

Enzyme substrates and products often diffuse too rapidly to assess the catalytic implications of these movements. Here, the authors characterise the structural basis of product and substrate diffusion for an exo-hydrolase and discover a substrate-product assisted processive catalytic mechanism.

The Effect of Cofactor Binding on the Conformational Plasticity of the Biological Receptors in Artificial Metalloenzymes: The Case Study of LmrR

The design of Artificial Metalloenzymes (ArMs), which result from the incorporation of organometallic cofactors into biological structures, has grown steadily in the last two decades and important new-to-Nature reactions have been reached. These type of exercises could greatly benefit from an understanding of the structural impact that the inclusion of organometallic moieties may have on the biological host. To date though, our understanding of this phenomenon is highly partial. This lack of knowledge is one of the elements that condition that first-generation ArMs generally display relatively poor catalytic profiles. In this work, we approach this matter by assessing the dynamics and stability of a series of ArMs resulting from the inclusion, via different anchoring strategies, of a variety of organometallic cofactors into the Lactococcal multidrug resistance regulator (LmrR) protein. To this aim, we coupled standard force field-based techniques such as Protein-Ligand Docking and Molecular Dynamics simulations with a variety of trajectory convergence analyses, capable of assessing both the stability and flexibility of the different systems under study upon the binding of cofactors. Together with the experimental evidence obtained in other studies, we provide an overview on how these changes can affect the catalytic outcomes obtained from the different ArMs. Fundamentally, our results show that the convergence analysis used in this work can assess how the inclusion of synthetic metallic cofactors in proteins can condition different structural modulations of their host. Those conformational modifications are key to the success of the desired catalytic activity and their proper identification can be wisely used to improve the quality and the rate of success of the ArMs.

The structure of an iron-containing alcohol dehydrogenase from a hyperthermophilic archaeon in two chemical states

An iron-containing alcohol dehydrogenase (FeADH) from the hyperthermophilic archaeon Thermococcus thioreducens was crystallized in unit cells belonging to space groups P21, P212121 and P43212, and the crystal structures were solved at 2.4, 2.1 and 1.9 Å resolution, respectively, by molecular replacement using the FeADH from Thermotoga maritima (Schwarzenbacher et al., 2004) as a model. In the monoclinic and orthorhombic crystals the dehydrogenase (molecular mass 41.5 kDa) existed as a dimer containing a twofold noncrystallographic symmetry axis, which was crystallographic in the tetragonal crystals. In the monoclinic and orthorhombic asymmetric units one molecule contained iron and an NADP molecule, while the other did not. The tetragonal crystals lacked both iron and NADP. The structure is very similar to that of the FeADH from T. maritima (average r.m.s. difference for Cα atoms of 1.8 Å for 341 aligned atoms). The iron, which is internally sequestered, is bound entirely by amino acids from one domain: three histidines and one aspartic acid. The coenzyme is in an extended conformation, a feature that is common to the large superfamily of NADH-dependent dehydrogenases that share a classical nucleotide-binding domain. A long broad tunnel passes entirely through the enzyme between the two domains, completely encapsulating the coenzyme.

Structure of glycerol dehydrogenase (GldA) from Escherichia coli

Escherichia coli (strain K-12, substrain MG1655) glycerol dehydrogenase (GldA) is required to catalyze the first step in fermentative glycerol metabolism. The protein was expressed and purified to homogeneity using a simple combination of heat-shock and chromatographic methods. The high yield of the protein (∼250 mg per litre of culture) allows large-scale production for potential industrial applications. Purified GldA exhibited a homogeneous tetrameric state (∼161 kDa) in solution and relatively high thermostability (T_m = 65.6°C). Sitting-drop sparse-matrix screens were used for protein crystallization. An optimized condition with ammonium sulfate (2 M) provided crystals suitable for diffraction, and a binary structure containing glycerol in the active site was solved at 2.8 Å resolution. Each GldA monomer consists of nine β-strands, thirteen α-helices, two 3₁₀-helices and several loops organized into two domains, the N- and C-terminal domains; the active site is located in a deep cleft between the two domains. The N-terminal domain contains a classic Rossmann fold for NAD⁺ binding. The O₁ and O₂ atoms of glycerol serve as ligands for the tetrahedrally coordinated Zn²⁺ ion. The orientation of the glycerol within the active site is mainly stabilized by van der Waals and electrostatic interactions with the benzyl ring of Phe245. Computer modeling suggests that the glycerol molecule is sandwiched by the Zn²⁺ and NAD⁺ ions. Based on this, the mechanism for the relaxed substrate specificity of this enzyme is also discussed.

Conformational Excitation and Nonequilibrium Transition Facilitate Enzymatic Reactions: Application to Pin1 Peptidyl-Prolyl Isomerase

Conformational flexibility of protein is essential for enzyme catalysis. Yet, how protein's conformational rearrangements and dynamics contribute to catalysis remains highly controversial. To unravel protein's role in catalysis, it is inevitable to understand the static and dynamic mechanisms simultaneously. To this end, here the Pin1-catalyzed isomerization reaction is studied from the two perspectives. The static view indicates that the hydrogen bonds involving Pin1 rearrange in a tightly coupled manner with isomerization. In sharp contrast, the isomerization dynamics are found to be very rapid; protein's slow conformational rearrangements thus cannot occur simultaneously with isomerization, and the reaction proceeds in a nonequilibrium manner. The distinctive protein conformations necessary to stabilize the transition state are prepared a priori, i.e., as conformational excited states. The present result suggests that enzymatic reaction is not a simple thermal activation from equilibrium directly to the transition state, thus adding a novel perspective to Pauling's view.

Targeting Metalloenzymes for Therapeutic Intervention

Metalloenzymes are central to a wide range of essential biological activities, including nucleic acid modification, protein degradation, and many others. The role of metalloenzymes in these processes also makes them central for the progression of many diseases and, as such, makes metalloenzymes attractive targets for therapeutic intervention. Increasing awareness of the role metalloenzymes play in disease and their importance as a class of targets has amplified interest in the development of new strategies to develop inhibitors and ultimately useful drugs. In this Review, we provide a broad overview of several drug discovery efforts focused on metalloenzymes and attempt to map out the current landscape of high-value metalloenzyme targets.

Bulk gold catalyzes hydride transfer in the Cannizzaro and related reactions

Disproportionation reactions of benzaldehyde and benzyl alcohol are catalyzed by bulk gold with hot water as the only other reagent.

Why the Energy Landscape of Barnase Is Hierarchical

We have used NMR and computational methods to characterize the dynamics of the ribonuclease barnase over a wide range of timescales in free and inhibitor-bound states. Using temperature- and denaturant-dependent measurements of chemical shift, we show that barnase undergoes frequent and highly populated hinge bending. Using relaxation dispersion, we characterize a slower and less populated motion with a rate of 750 ± 200 s⁻¹, involving residues around the lip of the active site, which occurs in both free and bound states and therefore suggests conformational selection. Normal mode calculations characterize correlated hinge bending motions on a very rapid timescale. These three measurements are combined with previous measurements and molecular dynamics calculations on barnase to characterize its dynamic landscape on timescales from picoseconds to milliseconds and length scales from 0.1 to 2.5 nm. We show that barnase has two different large-scale fluctuations: one on a timescale of 10⁻⁹−10⁻⁶ s that has no free energy barrier and is a hinge bending that is determined by the architecture of the protein; and one on a timescale of milliseconds (i.e., 750 s⁻¹) that has a significant free energy barrier and starts from a partially hinge-bent conformation. These two motions can be described as hierarchical, in that the more highly populated faster motion provides a platform for the slower (less probable) motion. The implications are discussed. The use of temperature and denaturant is suggested as a simple and general way to characterize motions on the intermediate ns-μs timescale.

Nature of the Pre-Chemistry Ensemble in Mitogen-Activated Protein Kinases

In spite of the availability of a significant amount of structural detail on docking interactions involving mitogen-activated protein kinases (MAPKs) and their substrates, the mechanism by which the disordered phospho-acceptor on the substrate transiently interacts with the kinase catalytic elements and is phosphorylated, often with high efficiency, remains poorly understood. Here, this dynamic interaction is analyzed in the context of available biophysical and biochemical data for ERK2, an archetypal MAPK. A hypothesis about the nature of the ternary complex involving a MAPK, its substrate, and ATP immediately prior to the chemical step (the pre-chemistry complex) is proposed. It is postulated that the solution ensemble (the pre-chemistry ensemble) representing the pre-chemistry complex comprises several conformations that are linked by dynamics on multiple timescales. These individual conformations possess different intrinsic abilities to proceed through the chemical step. The overall rate of chemistry is therefore related to the microscopic nature of the pre-chemistry ensemble, its constituent conformational microstates, and their intrinsic abilities to yield a phosphorylated product. While characterizing these microstates within the pre-chemistry ensemble in atomic or near-atomic detail is an extremely challenging proposition, recent developments in hybrid methodologies that employ computational approaches driven by experimental data appear to provide the most promising path forward toward achieving this goal.

Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis

The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP⁺ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP⁺ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.

Catalytic enzymes are active matter

Catalysis and mobility of reactants in fluid are normally thought to be decoupled. Violating this classical paradigm, this paper presents the catalyst laws of motion. Comparing experimental data to the theory presented here, we conclude that part of the free energy released by chemical reaction is channeled into driving catalysts to execute wormlike trajectories by piconewton forces performing work of a few k_BT against fluid viscosity, where the rotational diffusion rate dictates the trajectory persistence length. This active motion agitates the fluid medium and produces antichemotaxis, the migration of catalyst down the gradient of the reactant concentration. Alternative explanations of enhanced catalyst mobility are examined critically.

Using a microscopic theory to analyze experiments, we demonstrate that enzymes are active matter. Superresolution fluorescence measurements—performed across four orders of magnitude of substrate concentration, with emphasis on the biologically relevant regime around or below the Michaelis–Menten constant—show that catalysis boosts the motion of enzymes to be superdiffusive for a few microseconds, enhancing their effective diffusivity over longer timescales. Occurring at the catalytic turnover rate, these fast ballistic leaps maintain direction over a duration limited by rotational diffusion, driving enzymes to execute wormlike trajectories by piconewton forces performing work of a few k_BT against viscosity. The boosts are more frequent at high substrate concentrations, biasing the trajectories toward substrate-poor regions, thus exhibiting antichemotaxis, demonstrated here experimentally over a wide range of aqueous concentrations. Alternative noncatalytic, passive mechanisms that predict chemotaxis, cross-diffusion, and phoresis, are critically analyzed. We examine the physical interpretation of our findings, speculate on the underlying mechanism, and discuss the avenues they open with biological and technological implications. These findings violate the classical paradigm that chemical reaction and motility are distinct processes, and suggest reaction–motion coupling as a general principle of catalysis.

Conformational entropy of a single peptide controlled under force governs protease recognition and catalysis

An immense repertoire of protein chemical modifications catalyzed by enzymes is available as proteomics data. Quantifying the impact of the conformational dynamics of the modified peptide remains challenging to understand the decisive kinetics and amino acid sequence specificity of these enzymatic reactions in vivo, because the target peptide must be disordered to accommodate the specific enzyme-binding site. Here, we were able to control the conformation of a single-molecule peptide chain by applying mechanical force to activate and monitor its specific cleavage by a model protease. We found that the conformational entropy impacts the reaction in two distinct ways. First, the flexibility and accessibility of the substrate peptide greatly increase upon mechanical unfolding. Second, the conformational sampling of the disordered peptide drives the specific recognition, revealing force-dependent reaction kinetics. These results support a mechanism of peptide recognition based on conformational selection from an ensemble that we were able to quantify with a torsional free-energy model. Our approach can be used to predict how entropy affects site-specific modifications of proteins and prompts conformational and mechanical selectivity.

Asp271 is critical for substrate interaction with the surface binding site in β-agarase a from Zobellia galactanivorans

In the marine environment agar degradation is assured by bacteria that contain large agarolytic systems with enzymes acting in various endo- and exo-modes. Agarase A (AgaA) is an endo-glycoside hydrolase of family 16 considered to initiate degradation of agarose. Agaro-oligosaccharide binding at a unique surface binding site (SBS) in AgaA from Zobellia galactanivorans was investigated by computational methods in conjunction with a structure/sequence guided approach of site-directed mutagenesis probed by surface plasmon resonance binding analysis of agaro-oligosaccharides of DP 4-10. The crystal structure has shown that agaro-octaose interacts via H-bonds and aromatic stacking along 7 subsites (L through R) of the SBS in the inactive catalytic nucleophile mutant AgaA-E147S. D271 is centrally located in the extended SBS where it forms H-bonds to galactose and 3,6-anhydrogalactose residues of agaro-octaose at subsites O and P. We propose D271 is a key residue in ligand binding to the SBS. Thus AgaA-E147S/D271A gave slightly decreasing K_D values from 625 ± 118 to 468 ± 13 μM for agaro-hexaose, -octaose, and -decaose, which represent 3- to 4-fold reduced affinity compared with AgaA-E147S. Molecular dynamics simulations and interaction analyses of AgaA-E147S/D271A indicated disruption of an extended H-bond network supporting that D271 is critical for the functional SBS. Notably, neither AgaA-E147S/W87A nor AgaA-E147S/W277A, designed to eliminate stacking with galactose residues at subsites O and Q, respectively, were produced in soluble form. W87 and W277 may thus control correct folding and structural integrity of AgaA.

Structurally- and dynamically-driven allostery of the chymotrypsin-like proteases of SARS, Dengue and Zika viruses

Coronavirus 3C-like and Flavivirus NS2B-NS3 proteases utilize the chymotrypsin fold to harbor their catalytic machineries but also contain additional domains/co-factors. Over the past decade, we aimed to decipher how the extra domains/co-factors mediate the catalytic machineries of SARS 3C-like, Dengue and Zika NS2B-NS3 proteases by characterizing their folding, structures, dynamics and inhibition with NMR, X-ray crystallography and MD simulations, and the results revealed: 1) the chymotrypsin fold of the SARS 3C-like protease can independently fold, while, by contrast, those of Dengue and Zika proteases lack the intrinsic capacity to fold without co-factors. 2) Mutations on the extra domain of SARS 3C-like protease can transform the active catalytic machinery into the inactive collapsed state by structurally-driven allostery. 3) Amazingly, even without detectable structural changes, mutations on the extra domain are sufficient to either inactivate or enhance the catalytic machinery of SARS 3C-like protease by dynamically-driven allostery. 4) Global networks of correlated motions have been identified: for SARS 3C-like protease, N214A inactivates the catalytic machinery by decoupling the network, while STI/A and STIF/A enhance by altering the patterns of the network. The global networks of Dengue and Zika proteases are coordinated by their NS2B-cofactors. 5) Natural products were identified to allosterically inhibit Zika and Dengue proteases through binding a pocket on the back of the active site. Therefore, by introducing extra domains/cofactors, nature develops diverse strategies to regulate the catalytic machinery embedded on the chymotrypsin fold through folding, structurally- and dynamically-driven allostery, all of which might be exploited to develop antiviral drugs.

High-pressure protein crystal structure analysis of Escherichia coli dihydrofolate reductase complexed with folate and NADP. Acta Crystallogr D Struct Biol 2018;74:895-905. [PMID: 30198899 PMCID: PMC6130465 DOI: 10.1107/s2059798318009397]

A high-pressure crystallographic study was conducted on Escherichia coli dihydrofolate reductase (ecDHFR) complexed with folate and NADP+ in crystal forms containing both the open and closed conformations of the M20 loop under high-pressure conditions of up to 800 MPa. At pressures between 270 and 500 MPa the crystal form containing the open conformation exhibited a phase transition from P21 to C2. Several structural changes in ecDHFR were observed at high pressure that were also accompanied by structural changes in the NADP+ cofactor and the hydration structure. In the crystal form with the closed conformation the M20 loop moved as the pressure changed, with accompanying conformational changes around the active site, including NADP+ and folate. These movements were consistent with the suggested hypothesis that movement of the M20 loop was necessary for ecDHFR to catalyze the reaction. In the crystal form with the open conformation the nicotinamide ring of the NADP+ cofactor undergoes a large flip as an intermediate step in the reaction, despite being in a crystalline state. Furthermore, observation of the water molecules between Arg57 and folate elucidated an early step in the substrate-binding pathway. These results demonstrate the possibility of using high-pressure protein crystallography as a method to capture high-energy substates or transient structures related to the protein reaction cycle.

Evolutionarily Related Dihydrofolate Reductases Perform Coequal Functions Yet Show Divergence in Their Trajectories

The enzyme dihydrofolate reductase (DHFR) catalyzes NADPH dependent reduction of dihydrofolate to tetrahydrofolate. It plays a crucial role in the DNA synthesis. The investigation of evolution of DHFR generates immense curiosity. It aids in predicting how the enzyme has adapted to the surroundings of various cell types. In spite of great similarity in the structure of E. coli DHFR and human DHFR, their primary sequences are divergent to a great extent, which is evident in variations in the kinetics mechanism of their catalysis. In presence of physiological levels of ligands, they possess distinct kinetics and different rate limiting steps. We have reviewed the process of their unfolding and refolding, their behaviour in denaturing conditions and in presence of various chaperones. Although there is structural similarity between these two homologous enzymes yet they have established distinct mechanisms to accomplish the coequal functions.

Role of conformational dynamics in the evolution of novel enzyme function

The free energy landscape concept that describes enzymes as an ensemble of differently populated conformational sub-states in dynamic equilibrium is key for evaluating enzyme activity, enantioselectivity, and specificity. Mutations introduced in the enzyme sequence can alter the populations of the pre-existing conformational states, thus strongly modifying the enzyme ability to accommodate alternative substrates, revert its enantiopreferences, and even increase the activity for some residual promiscuous reactions. In this feature article, we present an overview of the current experimental and computational strategies to explore the conformational free energy landscape of enzymes. We provide a series of recent publications that highlight the key role of conformational dynamics for the enzyme evolution towards new functions and substrates, and provide some perspectives on how conformational dynamism should be considered in future computational enzyme design protocols.

Green function of correlated genes in a minimal mechanical model of protein evolution

The function of proteins arises from cooperative interactions and rearrangements of their amino acids, which exhibit large-scale dynamical modes. Long-range correlations have also been revealed in protein sequences, and this has motivated the search for physical links between the observed genetic and dynamic cooperativity. We outline here a simplified theory of protein, which relates sequence correlations to physical interactions and to the emergence of mechanical function. Our protein is modeled as a strongly coupled amino acid network with interactions and motions that are captured by the mechanical propagator, the Green function. The propagator describes how the gene determines the connectivity of the amino acids and thereby, the transmission of forces. Mutations introduce localized perturbations to the propagator that scatter the force field. The emergence of function is manifested by a topological transition when a band of such perturbations divides the protein into subdomains. We find that epistasis-the interaction among mutations in the gene-is related to the nonlinearity of the Green function, which can be interpreted as a sum over multiple scattering paths. We apply this mechanical framework to simulations of protein evolution and observe long-range epistasis, which facilitates collective functional modes.

Cooperativity Principles in Self-Assembled Nanomedicine

Nanomedicine is a discipline that applies nanoscience and nanotechnology principles to the prevention, diagnosis, and treatment of human diseases. Self-assembly of molecular components is becoming a common strategy in the design and syntheses of nanomaterials for biomedical applications. In both natural and synthetic self-assembled nanostructures, molecular cooperativity is emerging as an important hallmark. In many cases, interplay of many types of noncovalent interactions leads to dynamic nanosystems with emergent properties where the whole is bigger than the sum of the parts. In this review, we provide a comprehensive analysis of the cooperativity principles in multiple self-assembled nanostructures. We discuss the molecular origin and quantitative modeling of cooperative behaviors. In selected systems, we describe the examples on how to leverage molecular cooperativity to design nanomedicine with improved diagnostic precision and therapeutic efficacy in medicine.

Direct observation of ultrafast large-scale dynamics of an enzyme under turnover conditions

The potential effect of conformational dynamics of enzymes on their chemical steps has been intensely debated recently. We use single-molecule FRET experiments on adenylate kinase (AK) to shed new light on this question. AK closes its domains to bring its two substrate close together for reaction. We show that domain closure takes only microseconds to complete, which is two orders of magnitude faster than the chemical reaction. Nevertheless, active-site mutants that reduce the rate of domain closure also reduce the reaction rate, suggesting a connection between the two phenomena. We propose that ultrafast domain closure is used by enzymes as a mechanism to optimize mutual orientation of substrates, a novel mode of coupling between conformational dynamics and catalysis.

The functional cycle of many proteins involves large-scale motions of domains and subunits. The relation between conformational dynamics and the chemical steps of enzymes remains under debate. Here we show that in the presence of substrates, domain motions of an enzyme can take place on the microsecond time scale, yet exert influence on the much-slower chemical step. We study the domain closure reaction of the enzyme adenylate kinase from Escherichia coli while in action (i.e., under turnover conditions), using single-molecule FRET spectroscopy. We find that substrate binding increases dramatically domain closing and opening times, making them as short as ∼15 and ∼45 µs, respectively. These large-scale conformational dynamics are likely the fastest measured to date, and are ∼100–200 times faster than the enzymatic turnover rate. Some active-site mutants are shown to fully or partially prevent the substrate-induced increase in domain closure times, while at the same time they also reduce enzymatic activity, establishing a clear connection between the two phenomena, despite their disparate time scales. Based on these surprising observations, we propose a paradigm for the mode of action of enzymes, in which numerous cycles of conformational rearrangement are required to find a mutual orientation of substrates that is optimal for the chemical reaction.