The dynamical nature of enzymatic catalysis

CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

A Bioinspired Copper-Pair Catalyst in Metal-Organic Frameworks for Molecular Dioxygen Activation and Aerobic Oxidative C-N Coupling

Open Cu sites were loaded to the UiO-67 metal-organic framework (MOF) skeleton by introduction of flexible Cu-binding pyridylmethylamine (pyma) side chains to the biphenyldicarboxylate linkers. Distance between Cu centers in the MOF pores was tuned by controlling the density of metal-binding side chains. "Interacted" Cu-pair or "isolated" monomeric Cu sites were achieved with high and low (pyma)Cu side chain loading, respectively. Spectroscopic and theoretical studies indicate that "interacted" Cu pairs can effectively bind and activate molecular dioxygen to form Cu2O2 clusters, which showed high catalytic activity for aerobic oxidative C-N coupling. On the contrary, MOF catalyst bearing isolated monomeric Cu sites only showed modest catalytic activity. Enhancement in catalytic performance for the Cu-pair catalyst is attributed to the remote synergistic effect of the paired Cu site, which binds molecular dioxygen and cleaves the O═O bond in a collaborative manner. This work demonstrates that noncovalently interacted metal-pair sites can effectively activate inert small molecules and promote heterogeneous catalytic processes.

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Anisotropic dynamics of an interfacial enzyme active site observed using tethered substrate analogs and ultrafast 2D IR spectroscopy

Enzymes accelerate the rates of biomolecular reactions by many orders of magnitude compared to bulk solution, and it is widely understood that this catalytic effect arises from a combination of polar pre-organization and electrostatic transition state stabilization. A number of recent reports have also implicated ultrafast (femtosecond-picosecond) timescale motions in enzymatic activity. However, complications arising from spatially-distributed disorder, the occurrence of multiple substrate binding modes, and the influence of hydration dynamics on solvent-exposed active sites still confound many experimental studies. Here we use ultrafast two-dimensional infrared (2D IR) spectroscopy and covalently-tethered substrate analogs to examine dynamical properties of the promiscuous Pyrococcus horikoshii ene-reductase (PhENR) active site in two binding configurations mimicking proposed "inactive" and "reactive" Michaelis complexes. Spectral diffusion measurements of aryl-nitrile substrate analogs reveal an end-to-end tradeoff between fast (sub-ps) and slow (>5 ps) motions. Fermi resonant aryl-azide analogs that sense interactions of coupled oscillators are described. Lineshape and quantum beat analyses of these probes reveal characteristics that correlate with aryl-nitrile frequency fluctuation correlation functions parameters, demonstrating that this anisotropy is an intrinsic property of the water-exposed active site, where countervailing gradients of fast dynamics and disorder in the reactant ground state are maintained near the hydration interface. Our results suggest several plausible factors leading to state-selective rate enhancement and promiscuity in PhENR. This study also highlights a strategy to detect perturbations to vibrational modes outside the transparent window of the mid-IR spectrum, which may be extended to other macromolecular systems.

Deciphering Evolutionary Trajectories of Lactate Dehydrogenases Provides New Insights into Allostery

Lactate dehydrogenase (LDH, EC.1.1.127) is an important enzyme engaged in the anaerobic metabolism of cells, catalyzing the conversion of pyruvate to lactate and NADH to NAD+. LDH is a relevant enzyme to investigate structure-function relationships. The present work provides the missing link in our understanding of the evolution of LDHs. This allows to explain (i) the various evolutionary origins of LDHs in eukaryotic cells and their further diversification and (ii) subtle phenotypic modifications with respect to their regulation capacity. We identified a group of cyanobacterial LDHs displaying eukaryotic-like LDH sequence features. The biochemical and structural characterization of Cyanobacterium aponinum LDH, taken as representative, unexpectedly revealed that it displays homotropic and heterotropic activation, typical of an allosteric enzyme, whereas it harbors a long N-terminal extension, a structural feature considered responsible for the lack of allosteric capacity in eukaryotic LDHs. Its crystallographic structure was solved in 2 different configurations typical of the R-active and T-inactive states encountered in allosteric LDHs. Structural comparisons coupled with our evolutionary analyses helped to identify 2 amino acid positions that could have had a major role in the attenuation and extinction of the allosteric activation in eukaryotic LDHs rather than the presence of the N-terminal extension. We tested this hypothesis by site-directed mutagenesis. The resulting C. aponinum LDH mutants displayed reduced allosteric capacity mimicking those encountered in plants and human LDHs. This study provides a new evolutionary scenario of LDHs that unifies descriptions of regulatory properties with structural and mutational patterns of these important enzymes.

Characterization of the pH-dependent protein stability of 3α-hydroxysteroid dehydrogenase/carbonyl reductase by differential scanning fluorimetry

The characterization of protein stability is essential for understanding the functions of proteins. Hydroxysteroid dehydrogenase is involved in the biosynthesis of steroid hormones and the detoxification of xenobiotic carbonyl compounds. However, the stability of hydroxysteroid dehydrogenases has not yet been characterized in detail. Here, we determined the changes in Gibbs free energy, enthalpy, entropy, and heat capacity of unfolding for 3α-hydroxysteroid dehydrogenase/carbonyl reductase (3α-HSD/CR) by varying the pH and urea concentration through differential scanning fluorimetry and presented pH-dependent protein stability as a function of temperature. 3α-HSD/CR shows the maximum stability of 30.79 kJ mol-1 at 26.4°C, pH 7.6 and decreases to 7.74 kJ mol-1 at 25.7°C, pH 4.5. The change of heat capacity of 30.25 ± 1.38 kJ mol-1  K-1 is obtained from the enthalpy of denaturation as a function of melting temperature at varied pH. Two proton uptakes are linked to protein unfolding from residues with differential pKa of 4.0 and 6.5 in the native and denatured states, respectively. The large positive heat capacity change indicated that hydrophobic interactions played an important role in the folding of 3α-HSD/CR. These studies reveal the mechanism of protein unfolding in HSD and provide a convenient method to extract thermodynamic parameters for characterizing protein stability using differential scanning fluorimetry.

Coacervate-Based Plexcitonic Assembly toward Peroxidase-like Activity and Ultrasensitive Glucose Sensing

Inbuilt catalytic centers anchored inside the confined architecture of artificial nanoreactors have gained tremendous attention owing to their vast applicability in various catalytic transformations. However, designing homogeneously distributed catalytic units with exposed surfaces in confined environment is a challenging task. Here, we have utilized quantum dot (QD)-embedded coacervate droplets (QD-Ds) as a confined compartment for the in situ synthesis of gold nanoparticles (Au NPs) without any additional reducing agent. High-resolution transmission electron microscopy images reveal homogeneous distribution of 5.6 ± 0.2 nm-sized Au NPs inside the QD-Ds (Au@QD-Ds). The in situ synthesized Au NPs are found to be stable over a period of 28 days without any agglomeration. Control experiments reveal that the free surface carboxylic acid groups of embedded QDs simultaneously act as reducing and stabilizing agents for Au NPs. Notably, these Au@QD-Ds exhibit superior peroxidase-like activity compared to bulk aqueous Au NPs and Au@QDs under similar experimental conditions. The observed peroxidase-like activity follows the classical Michaelis-Menten model inside the Au@QD-Ds via the fast electron-transfer pathway. The enhanced peroxidase-like activity has been explained by considering confinement, mass action, and the ligand-free surface of embedded Au NPs. The present plexcitonic nanocomposites exhibit excellent recyclability over several consecutive cycles without any compromise in their catalytic activity. Finally, a cascade reaction with glucose oxidase (GOx)-loaded Au@QD-Ds have been utilized for colorimetric detection of glucose with a limit of detection of 272 nM in solution as well as on filter paper. The present work highlights a facile and robust methodology for the fabrication of optically active functional hybrid plexcitonic assemblies and may find importance in various fields including bioanalytical chemistry and optoelectronics.

Bioinspired Active Site with a Coordination-Adaptive Organosulfonate Ligand for Catalytic Water Oxidation at Neutral pH

Many enzymes use adaptive frameworks to preorganize substrates, accommodate various structural and electronic demands of intermediates, and accelerate related catalysis. Inspired by biological systems, a Ru-based molecular water oxidation catalyst containing a configurationally labile ligand [2,2':6',2″-terpyridine]-6,6″-disulfonate was designed to mimic enzymatic framework, in which the sulfonate coordination is highly flexible and functions as both an electron donor to stabilize high-valent Ru and a proton acceptor to accelerate water dissociation, thus boosting the catalytic water oxidation performance thermodynamically and kinetically. The combination of single-crystal X-ray analysis, various temperature NMR, electrochemical techniques, and DFT calculations was utilized to investigate the fundamental role of the self-adaptive ligand, demonstrating that the on-demand configurational changes give rise to fast catalytic kinetics with a turnover frequency (TOF) over 2000 s^-1, which is compared to oxygen-evolving complex (OEC) in natural photosynthesis.

Allosteric communication between ligand binding domains modulates substrate inhibition in adenylate kinase

Enzymes play a vital role in life processes; they control chemical reactions and allow functional cycles to be synchronized. Many enzymes harness large-scale motions of their domains to achieve tremendous catalytic prowess and high selectivity for specific substrates. One outstanding example is provided by the three-domain enzyme adenylate kinase (AK), which catalyzes phosphotransfer between ATP to AMP. Here we study the phenomenon of substrate inhibition by AMP and its correlation with domain motions. Using single-molecule FRET spectroscopy, we show that AMP does not block access to the ATP binding site, neither by competitive binding to the ATP cognate site nor by directly closing the LID domain. Instead, inhibitory concentrations of AMP lead to a faster and more cooperative domain closure by ATP, leading in turn to an increased population of the closed state. The effect of AMP binding can be modulated through mutations throughout the structure of the enzyme, as shown by the screening of an extensive AK mutant library. The mutation of multiple conserved residues reduces substrate inhibition, suggesting that substrate inhibition is an evolutionary well conserved feature in AK. Combining these insights, we developed a model that explains the complex activity of AK, particularly substrate inhibition, based on the experimentally observed opening and closing rates. Notably, the model indicates that the catalytic power is affected by the microsecond balance between the open and closed states of the enzyme. Our findings highlight the crucial role of protein motions in enzymatic activity.

Minimalistic Artificial Catalysts with Esterase-Like Activity from Multivalent Nanofibers Formed by the Self-Assembly of Dipeptides

Imitating and incorporating the multiple key structural features observed in natural enzymes into a minimalistic molecule to develop an artificial catalyst with outstanding catalytic efficiency is an attractive topic for chemists. Herein, we designed and synthesized one class of minimalistic dipeptide molecules containing a terminal -SH group and a terminal His-Phe dipeptide head linked by a hydrophobic alkyl chain with different lengths, marked as HS-C _n+1-His-Phe (n = 4, 7, 11, 15, and 17; n + 1 represents the carbon atom number of the alkyl chain). The His (-imidazole), Phe (-CO₂ ^-) moieties, the terminal -SH group, and a long hydrophobic alkyl chain were found to have important contributions to achieve high binding ability leading to outstanding absolute catalytic efficiency (k _cat/K _M) toward the hydrolysis reactions of carboxylic ester substrates.

DHFR Mutants Modulate Their Synchronized Dynamics with the Substrate by Shifting Hydrogen Bond Occupancies

Antibiotic resistance is a global health problem in which mutations occurring in functional proteins render drugs ineffective. The working mechanisms of the arising mutants are seldom apparent; a methodology to decipher these mechanisms systematically would render devising therapies to control the arising mutational pathways possible. Here we utilize C_α-C_β bond vector relaxations obtained from moderate length MD trajectories to determine conduits for functionality of the resistance conferring mutants of Escherichia coli dihydrofolate reductase. We find that the whole enzyme is synchronized to the motions of the substrate, irrespective of the mutation introducing gain-of-function or loss-of function. The total coordination of the motions suggests changes in the hydrogen bond dynamics with respect to the wild type as a possible route to determine and classify the mode-of-action of individual mutants. As a result, nine trimethoprim-resistant point mutations arising frequently in evolution experiments are categorized. One group of mutants that display the largest occurrence (L28R, W30G) work directly by modifying the dihydrofolate binding region. Conversely, W30R works indirectly by the formation of the E139-R30 salt bridge which releases energy resulting from tight binding by distorting the binding cavity. A third group (D27E, F153S, I94L) arising as single, resistance invoking mutants in evolution experiment trajectories allosterically and dynamically affects a hydrogen bonding motif formed at residues 59-69-71 which in turn modifies the binding site dynamics. The final group (I5F, A26T, R98P) consists of those mutants that have properties most similar to the wild type; these only appear after one of the other mutants is fixed on the protein structure and therefore display clear epistasis. Thus, we show that the binding event is governed by the entire enzyme dynamics while the binding site residues play gating roles. The adjustments made in the total enzyme in response to point mutations are what make quantifying and pinpointing their effect a hard problem. Here, we show that hydrogen bond dynamics recorded on sub-μs time scales provide the necessary fingerprints to decipher the various mechanisms at play.

Temperature-sensitive contacts in disordered loops tune enzyme I activity

Homologous enzymes with identical folds often exhibit different thermal and kinetic behaviors. Understanding how an enzyme sequence encodes catalytic activity at functionally optimal temperatures is a fundamental problem in biophysics. Recently it was shown that the residues that tune catalytic activities of thermophilic/mesophilic variants of the C-terminal domain of bacterial enzyme I (EIC) are largely localized within disordered loops, offering a model system with which to investigate this phenomenon. In this work, we use molecular dynamics simulations and mutagenesis experiments to reveal a mechanism of sequence-dependent activity tuning of EIC homologs. We find that a network of contacts in the catalytic loops is particularly sensitive to changes in temperature, with some contacts exhibiting distinct linear or nonlinear temperature-dependent trends. Moreover, these trends define structurally clustered dynamical modes and can distinguish regions that tend toward order or disorder at higher temperatures. Assaying several thermophilic EIC mutants, we show that complementary mesophilic mutations to the most temperature-sensitive positions exhibit the most enhanced activity, while mutations to relatively temperature insensitive positions exhibit the least enhanced activities. These results provide a mechanistic explanation of sequence-dependent temperature tuning and offer a computational method for rational enzyme modification.

Role of Enzyme and Active Site Conformational Dynamics in the Catalysis by α-Amylase Explored with QM/MM Molecular Dynamics

We assessed enzyme:substrate conformational dynamics and the rate-limiting glycosylation step of a human pancreatic α-amylase:maltopentose complex. Microsecond molecular dynamics simulations suggested that the distance of the catalytic Asp197 nucleophile to the anomeric carbon of the buried glucoside is responsible for most of the enzyme active site fluctuations and that both Asp197 and Asp300 interact the most with the buried glucoside unit. The buried glucoside binds either in a ⁴C₁ chair or ²S_O skew conformations, both of which can change to TS-like conformations characteristic of retaining glucosidases. Starting from four distinct enzyme:substrate complexes, umbrella sampling quantum mechanics/molecular mechanics simulations (converged within less than 1 kcal·mol^-1 within a total simulation time of 1.6 ns) indicated that the reaction occurrs with a Gibbs barrier of 13.9 kcal·mol ^-1, in one asynchronous concerted step encompassing an acid-base reaction with Glu233 followed by a loose S_N2-like nucleophilic substitution by the Asp197. The transition state is characterized by a ²H₃ half-chair conformation of the buried glucoside that quickly changes to the E₃ envelope conformation preceding the attack of the anomeric carbon by the Asp197 nucleophile. Thermodynamic analysis of the reaction supported that a water molecule tightly hydrogen bonded to the glycosidic oxygen of the substrate at the reactant state (∼1.6 Å) forms a short hydrogen bond with Glu233 at the transition state (∼1.7 Å) and lowers the Gibbs barrier in over 5 kcal·mol^-1. The resulting Asp197-glycosyl was mostly found in the ⁴C₁ conformation, although the more endergonic B_3,O conformation was also observed. Altogether, the combination of short distances for the acid-base reaction with the Glu233 and for the nucleophilic attack by the Asp197 nucleophile and the availability of water within hydrogen bonding distance of the glycosidic oxygen provides a reliable criteria to identify reactive conformations of α-amylase complexes.

Some Clues about Enzymes from Psychrophilic Microorganisms

Enzymes purified from psychrophilic microorganisms prove to be efficient catalysts at low temperatures and possess a great potential for biotechnological applications. The low-temperature catalytic activity has to come from specific structural fluctuations involving the active site region, however, the relationship between protein conformational stability and enzymatic activity is subtle. We provide a survey of the thermodynamic stability of globular proteins and their rationalization grounded in a theoretical approach devised by one of us. Furthermore, we provide a link between marginal conformational stability and protein flexibility grounded in the harmonic approximation of the vibrational degrees of freedom, emphasizing the occurrence of long-wavelength and excited vibrations in all globular proteins. Finally, we offer a close view of three enzymes: chloride-dependent α-amylase, citrate synthase, and β-galactosidase.

Impact of protein conformational diversity on AlphaFold predictions

MOTIVATION

After the outstanding breakthrough of AlphaFold in predicting protein 3D models, new questions appeared and remain unanswered. The ensemble nature of proteins, for example, challenges the structural prediction methods because the models should represent a set of conformers instead of single structures. The evolutionary and structural features captured by effective deep learning techniques may unveil the information to generate several diverse conformations from a single sequence. Here, we address the performance of AlphaFold2 predictions obtained through ColabFold under this ensemble paradigm.

RESULTS

Using a curated collection of apo-holo pairs of conformers, we found that AlphaFold2 predicts the holo form of a protein in ∼70% of the cases, being unable to reproduce the observed conformational diversity with the same error for both conformers. More importantly, we found that AlphaFold2's performance worsens with the increasing conformational diversity of the studied protein. This impairment is related to the heterogeneity in the degree of conformational diversity found between different members of the homologous family of the protein under study. Finally, we found that main-chain flexibility associated with apo-holo pairs of conformers negatively correlates with the predicted local model quality score plDDT, indicating that plDDT values in a single 3D model could be used to infer local conformational changes linked to ligand binding transitions.

AVAILABILITY AND IMPLEMENTATION

Data and code used in this manuscript are publicly available at https://gitlab.com/sbgunq/publications/af2confdiv-oct2021.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Towards the Accurate Thermodynamic Characterization of Enzyme Reaction Mechanisms

We employed QM/MM molecular dynamics (MD) simulations to characterize the rate-limiting step of the glycosylation reaction of pancreatic α-amylase with combined DFT/molecular dynamics methods (PBE/def2-SVP:AMBER). Upon careful choice of four starting active site conformations based on thorough reactivity criteria, Gibbs energy profiles were calculated with umbrella sampling simulations within a statistical convergence of 1-2 kcal⋅mol -1 . Nevertheless, Gibbs activation barriers and reaction energies still varied from 11.0 to 16.8 kcal⋅mol -1 and -6.3 to +3.8 kcal⋅mol -1 depending on the starting conformations, showing that despite significant state-of-the-art QM/MM MD sampling (0.5 ns/profile) the result still depends on the starting structure. The results supported the one step dissociative mechanism of Asp197 glycosylation preceded by an acid-base reaction by the Glu233, which are qualitatively similar to those from multi-PES QM/MM studies, and thus support the use of the latter to determine enzyme reaction mechanisms.

Application of Coarse-Grained (CG) Models to Explore Conformational Pathway of Large-Scale Protein Machines

Protein machines are clusters of protein assemblies that function in order to control the transfer of matter and energy in cells. For a specific protein machine, its working mechanisms are not only determined by the static crystal structures, but also related to the conformational transition dynamics and the corresponding energy profiles. With the rapid development of crystallographic techniques, the spatial scale of resolved structures is reaching up to thousands of residues, and the concomitant conformational changes become more and more complicated, posing a great challenge for computational biology research. Previously, a coarse-grained (CG) model aiming at conformational free energy evaluation was developed and showed excellent ability to reproduce the energy profiles by accurate electrostatic interaction calculations. In this study, we extended the application of the CG model to a series of large-scale protein machine systems. The spike protein trimer of SARS-CoV-2, ATP citrate lyase (ACLY) tetramer, and P4-ATPases systems were carefully studied and discussed as examples. It is indicated that the CG model is effective to depict the energy profiles of the conformational pathway between two endpoint structures, especially for large-scale systems. Both the energy change and energy barrier between endpoint structures provide reasonable mechanism explanations for the associated biological processes, including the opening of receptor binding domain (RBD) of spike protein, the phospholipid transportation of P4-ATPase, and the loop translocation of ACLY. Taken together, the CG model provides a suitable alternative in mechanistic studies related to conformational change in large-scale protein machines.

Correlating the Local and Global Dynamics of an Enzyme in the Crowded Milieu

Enzymes are dynamic biological macromolecules, with their catalytic function(s) being largely influenced by the changes in local fluctuations of amino acid side chains as well as global structural modulations that the enzyme undergoes. Such local and global motions can be highly affected inside the crowded physiological interior of the cell. Here, we have addressed the role of dynamic structural flexibility in affecting the activation energy barrier of a flexible multidomain enzyme adenylate kinase (AK3L1, UniProtKB: Q9UIJ7). Activation energy profiles of both local (at three different sites along the polypeptide backbone) and global dynamics of the enzyme have been monitored using solvation studies on the subnanosecond time scale and tryptophan quenching studies over the temperature range of 278-323 K, respectively, under crowded conditions (Ficoll 70, Dextran 40, Dextran 70, and PEG 8). This study not only provides the site-specific mapping of dynamics but reveals that the activation energies associated with these local motions undergo a significant decrease in the presence of macromolecular crowders, providing new insights into how crowding affects internal protein dynamics. The crowded scenario also aids in enhancing the coupling between the local and global motions of the enzyme. Moreover, select portions/regions of the enzyme when taken together can well mirror the overall dynamics of the biomolecule, showing possible energy hotspots along the polypeptide backbone.

Dynamic Tuning in Synthetic Glycosidase for Selective Hydrolysis of Alkyl and Aryl Glycosides

Enzymes use sophisticated conformational control to optimize the dynamics of their protein framework for efficient catalysis. Although it is difficult to employ a similar strategy to improve catalysis in a synthetic enzyme, we here report that modulation of the dynamics of the substrate in the active site is readily achievable in a complex between a molecularly imprinted nanoparticle and its acid cofactor, through tuning of the size and shape of the imprinted site. As the alkyl glucoside substrate is bound with increasing strength and held in a more tightly fitted pocket, the acid-catalyzed glycan hydrolysis becomes more difficult. A larger, wider active site, although less able to bind the substrate, affords a higher catalytic activity, likely due to easier alignment of the substrate and the acid cofactor for a general acid catalysis. The substrate selectivity is controlled by both the tightness of the aglycon-binding site and the orientation of the glycan-binding boroxole group.

Designing Sites in Heterogeneous Catalysis: Are We Reaching Selectivities Competitive With Those of Homogeneous Catalysts?

A critical review of different prominent nanotechnologies adapted to catalysis is provided, with focus on how they contribute to the improvement of selectivity in heterogeneous catalysis. Ways to modify catalytic sites range from the use of the reversible or irreversible adsorption of molecular modifiers to the immobilization or tethering of homogeneous catalysts and the development of well-defined catalytic sites on solid surfaces. The latter covers methods for the dispersion of single-atom sites within solid supports as well as the use of complex nanostructures, and it includes the post-modification of materials via processes such as silylation and atomic layer deposition. All these methodologies exhibit both advantages and limitations, but all offer new avenues for the design of catalysts for specific applications. Because of the high cost of most nanotechnologies and the fact that the resulting materials may exhibit limited thermal or chemical stability, they may be best aimed at improving the selective synthesis of high value-added chemicals, to be incorporated in organic synthesis schemes, but other applications are being explored as well to address problems in energy production, for instance, and to design greener chemical processes. The details of each of these approaches are discussed, and representative examples are provided. We conclude with some general remarks on the future of this field.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Energy Bilocalization Effect and the Emergence of Molecular Functions in Proteins

Proteins are among the most complex molecular structures, which have evolved to develop broad functions, such as energy conversion and transport, information storage and processing, communication, and regulation of chemical reactions. However, the mechanisms by which these dynamical entities coordinate themselves to perform biological tasks remain hotly debated. Here, a physical theory is presented to explain how functional dynamical behavior possibly emerge in complex/macro molecules, thanks to the effect that we term bilocalization of thermal vibrations. More specifically, our approach allows us to understand how structural irregularities lead to a partitioning of the energy of the vibrations into two distinct sets of molecular domains, corresponding to slow and fast motions. This shape-encoded spectral allocation, associated to the genetic sequence, provides a close access to a wide reservoir of dynamical patterns, and eventually allows the emergence of biological functions by natural selection. To illustrate our approach, the SPIKE protein structure of SARS-COV2 is considered.

Antibodies targeting enzyme inhibition as potential tools for research and drug development

Antibodies have transformed biomedical research and are now being used for different experimental applications. Generally, the interaction of enzymes with their specific antibodies can lead to a reduction in their enzymatic activity. The effect of the antibody is dependent on its narrow i.e. the regions of the enzyme to which it is directed. The mechanism of this inhibition is rarely a direct combination of the antibodies with the catalytic site, but is rather due to steric hindrance, barring the substrate access to the active site. In several systems, however, the interaction with the antibody induces conformational changes on the enzyme that can either inhibit or enhance its catalytic activity. The extent of enzyme inhibition or enhancement is, therefore, a reflection of the nature and distribution of the various antigenic determinants on the enzyme molecule. Currently, the mode of action of many enzymes has been elucidated at the molecular level. We here review the molecular mechanisms and recent trends by which antibodies inhibit the catalytic activity of enzymes and provide examples of how specific antibodies can be useful for the neutralization of biologically active molecules.

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.

Structure, Function, and Thermodynamics of Lactate Dehydrogenases from Humans and the Malaria Parasite P. falciparum

Temperature adaptation is ubiquitous among all living organisms, yet the molecular basis for this process remains poorly understood. It can be assumed that for parasite-host systems, the same enzymes found in both organisms respond to the same selection factor (human body temperature) with similar structural changes. Herein, we report the existence of a reversible temperature-dependent structural transition for the glycolytic enzyme lactate dehydrogenase (LDH) from the malaria parasite Plasmodium falciparum (pfLDH) and human heart (hhLDH) occurring in the temperature range of human fever. This transition is observed for LDHs from psychrophiles, mesophiles, and moderate thermophiles in their operating temperature range. Thermodynamic analysis reveals unique thermodynamic signatures of the LDH-substrate complexes defining a specific temperature range to which human LDH is adapted and parasite LDH is not, despite their common mesophilic nature. The results of spectroscopic analysis combined with the available crystallographic data reveal the existence of an active center within pfLDH that imparts psychrophilic structural properties to the enzyme. This center consists of two pockets, one formed by the five amino acids (5AA insert) within the substrate specificity loop and the other by the active site, that mutually regulate one another in response to temperature and induce structural and functional changes in the Michaelis complex. Our findings pave the way toward a new strategy for malaria treatments and drug design using therapeutic agents that inactivate malarial LDH selectively at a specific temperature range of the cyclic malaria paroxysm.

Understanding enzyme behavior in a crowded scenario through modulation in activity, conformation and dynamics

Macromolecular crowding, inside the physiological interior, modulates the energy landscape of biological macromolecules in multiple ways. Amongst these, enzymes occupy a special place and hence understanding the function of the same in the crowded interior is of utmost importance. In this study, we have investigated the manner in which the multidomain enzyme, AK3L1 (PDB ID: 1ZD8), an isoform of adenylate kinase, has its features affected in presence of commonly used crowders (PEG 8, Dextran 40, Dextran 70, and Ficoll 70). Michaelis Menten plots reveal that the crowders in general enhance the activity of the enzyme, with the K_m and V_max values showing significant variations. Ficoll 70, induced the maximum activity for AK3L1 at 100 g/L, beyond which the activity reduced. Ensemble FRET studies were performed to provide insights into the relative domain (LID and CORE) displacements in presence of the crowders. Solvation studies reveal that the protein matrix surrounding the probe CPM (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin) gets restricted in presence of the crowders, with Ficoll 70 providing the maximum rigidity, the same being linked to the decrease in the activity of the enzyme. Through our multipronged approach, we have observed a distinct correlation between domain displacement, enzyme activity and associated dynamics. Thus, keeping in mind the complex nature of enzyme activity and the surrounding bath of dense soup that the biological entity remains immersed in, indeed more such approaches need to be undertaken to have a better grasp of the "enzymes in the crowd".

Dynamic Connection between Enzymatic Catalysis and Collective Protein Motions

Enzymes employ a wide range of protein motions to achieve efficient catalysis of chemical reactions. While the role of collective protein motions in substrate binding, product release, and regulation of enzymatic activity is generally understood, their roles in catalytic steps per se remain uncertain. Here, molecular dynamics simulations, enzyme kinetics, X-ray crystallography, and nuclear magnetic resonance spectroscopy are combined to elucidate the catalytic mechanism of adenylate kinase and to delineate the roles of catalytic residues in catalysis and the conformational change in the enzyme. This study reveals that the motions in the active site, which occur on a time scale of picoseconds to nanoseconds, link the catalytic reaction to the slow conformational dynamics of the enzyme by modulating the free energy landscapes of subdomain motions. In particular, substantial conformational rearrangement occurs in the active site following the catalytic reaction. This rearrangement not only affects the reaction barrier but also promotes a more open conformation of the enzyme after the reaction, which then results in an accelerated opening of the enzyme compared to that of the reactant state. The results illustrate a linkage between enzymatic catalysis and collective protein motions, whereby the disparate time scales between the two processes are bridged by a cascade of intermediate-scale motion of catalytic residues modulating the free energy landscapes of the catalytic and conformational change processes.

Biological innovation in the functional landscape of a model regulator, or the lactose operon repressor

The operon model was proposed six decades ago. And yet, despite all this time, the lactose operon repressor, LacI, remains a subject of major interest. While it is well established that LacI can exist in two functional forms, one that renders the operon inactive via binding of LacI to DNA and another, bound to an inducer that does not allow repression, how it switches from one to the other is still not well understood. The construction of a library of several tens of thousands of LacI mutants has revealed some unexpected features. In particular, the transition implemented in some of them reveals a new type of transcription regulation: band-pass (OFF/ON/OFF) and band-stop (ON/OFF/ON) filters. This makes it natural to think that it is the network of hydrogen bonds associated with the water bound to the molecule that allows the remote interconnection between the binding site to an inducer molecule and the one that binds it to the DNA.

Ile258Met mutation of Brucella melitensis 7α-hydroxysteroid dehydrogenase significantly enhances catalytic efficiency, cofactor affinity, and thermostability

NAD(H)-dependent 7α-hydroxysteroid dehydrogenase catalyzes the oxidation of chenodeoxycholic acid to 7-oxolithocholic acid. Here, we designed mutations of Ile258 adjacent to the catalytic pocket of Brucella melitensis 7α-hydroxysteroid dehydrogenase. The I258M variant gave a 4.7-fold higher k_cat, but 4.5-fold lower K_M, compared with the wild type, resulting in a 21.8-fold higher k_cat/K_M value for chenodeoxycholic acid oxidation. It presented a 2.0-fold lower K_M value with NAD⁺, suggesting stronger binding to the cofactor. I258M produced 7-oxolithocholic acid in the highest yield of 92.3% in 2 h, whereas the wild-type gave 88.4% in 12 h. The I258M mutation increased the half-life from 20.8 to 31.1 h at 30 °C. Molecular dynamics simulations indicated increased interactions and a modified tunnel improved the catalytic efficiency, and enhanced rigidity at three regions around the ligand-binding pocket increased the enzyme thermostability. This is the first report about significantly improved catalytic efficiency, cofactor affinity, and enzyme thermostability through single site-mutation of Brucella melitensis 7α-hydroxysteroid dehydrogenase. KEY POINTS: • Sequence and structure analysis guided the site mutation design. • Thermostability, catalytic efficiency and 7-oxo-LCA production were determined. • MD simulation was performed to indicate the improvement by I258M mutation.

Non-Coordinative Binding of O2 at the Active Center of a Copper-Dependent Enzyme

Molecular oxygen (O₂ ) is a sustainable oxidation reagent. O₂ is strongly oxidizing but kinetically stable and its final reaction product is water. For these reasons learning how to activate O₂ and how to steer its reactivity along desired reaction pathways is a longstanding challenge in chemical research.^[1] Activation of ground-state diradical O₂ can occur either via conversion to singlet oxygen or by one-electron reduction to superoxide. Many enzymes facilitate activation of O₂ by direct fomation of a metal-oxygen coordination complex concomitant with inner sphere electron transfer. The formylglycine generating enzyme (FGE) is an unusual mononuclear copper enzyme that appears to follow a different strategy. Atomic-resolution crystal structures of the precatalytic complex of FGE demonstrate that this enzyme binds O₂ juxtaposed, but not coordinated to the catalytic Cu^I . Isostructural complexes that contain Ag^I instead of Cu^I or nitric oxide instead of O₂ confirm that formation of the initial oxygenated complex of FGE does not depend on redox activity. A stepwise mechanism that decouples binding and activation of O₂ is unprecedented for metal-dependent oxidases, but is reminiscent of flavin-dependent enzymes.

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine-histidine-aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the "moving histidine" mechanism, thus calling into question the "ring flip" mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s^-1. Overall, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

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.

Recent advances in hydride transfer-involved C(sp3)–H activation reactions

This review summarizes the recent progresses (2016–2020) in the hydride transfer-enabled C(sp³)–H activation according to the reaction types, categorized into the intramolecular/intermolecular C(sp³)–H functionalization, and hydride reduction.

The physical origin of rate promoting vibrations in enzymes revealed by structural rigidity

Enzymes are the most efficient catalysts known to date. However, decades of research have failed to fully explain the catalytic power of enzymes, and most of the current attempts to uncloak the details of atomic motions at active sites remain incomplete. Here, a straightforward manner for understanding the interplay between the complex or irregular enzyme topology and dynamical effects at catalytic sites is introduced, by revealing how fast localized vibrations form spontaneously in the stiffest parts of the scaffold. While shedding light on a physical mechanism that allowed the selection of the picosecond (ps) timescale to increase the catalytic proficiency, this approach exposes the functional importance of localized motions as a by-product of the stability-function tradeoff in enzyme evolution. From this framework of analysis—directly accessible from available diffraction data—experimental strategies for engineering the catalytic rate in enzymatic proteins are proposed.

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.

Dissecting the Structural and Chemical Determinants of the "Open-to-Closed" Motion in the Mannosyltransferase PimA from Mycobacteria

The phosphatidyl-myo-inositol mannosyltransferase A (PimA) is an essential peripheral membrane glycosyltransferase that initiates the biosynthetic pathway of phosphatidyl-myo-inositol mannosides (PIMs), key structural elements and virulence factors of Mycobacterium tuberculosis. PimA undergoes functionally important conformational changes, including (i) α-helix-to-β-strand and β-strand-to-α-helix transitions and (ii) an "open-to-closed" motion between the two Rossmann-fold domains, a conformational change that is necessary to generate a catalytically competent active site. In previous work, we established that GDP-Man and GDP stabilize the enzyme and facilitate the switch to a more compact active state. To determine the structural contribution of the mannose ring in such an activation mechanism, we analyzed a series of chemical derivatives, including mannose phosphate (Man-P) and mannose pyrophosphate-ribose (Man-PP-RIB), and additional GDP derivatives, such as pyrophosphate ribose (PP-RIB) and GMP, by the combined use of X-ray crystallography, limited proteolysis, circular dichroism, isothermal titration calorimetry, and small angle X-ray scattering methods. Although the β-phosphate is present, we found that the mannose ring, covalently attached to neither phosphate (Man-P) nor PP-RIB (Man-PP-RIB), does promote the switch to the active compact form of the enzyme. Therefore, the nucleotide moiety of GDP-Man, and not the sugar ring, facilitates the "open-to-closed" motion, with the β-phosphate group providing the high-affinity binding to PimA. Altogether, the experimental data contribute to a better understanding of the structural determinants involved in the "open-to-closed" motion not only observed in PimA but also visualized and/or predicted in other glycosyltransfeases. In addition, the experimental data might prove to be useful for the discovery and/or development of PimA and/or glycosyltransferase inhibitors.

Driving Protein Conformational Cycles in Physiology and Disease: "Frustrated" Amino Acid Interaction Networks Define Dynamic Energy Landscapes: Amino Acid Interaction Networks Change Progressively Along Alpha Tryptophan Synthase's Catalytic Cycle

A general framework by which dynamic interactions within a protein will promote the necessary series of structural changes, or "conformational cycle," required for function is proposed. It is suggested that the free-energy landscape of a protein is biased toward this conformational cycle. Fluctuations into higher energy, although thermally accessible, conformations drive the conformational cycle forward. The amino acid interaction network is defined as those intraprotein interactions that contribute most to the free-energy landscape. Some network connections are consistent in every structural state, while others periodically change their interaction strength according to the conformational cycle. It is reviewed here that structural transitions change these periodic network connections, which then predisposes the protein toward the next set of network changes, and hence the next structural change. These concepts are illustrated by recent work on tryptophan synthase. Disruption of these dynamic connections may lead to aberrant protein function and disease states.

Study on effects of co-solvents on the structure of DhaA by molecular dynamics simulation

With the increasing application of enzymes in various research fields, the choices of co-solvents in enzymatic preparations which directly related to the catalytic activity have been attracted attention. Thus, researching on the stabilization or destabilization behaviors of enzymes in different solvents is extremely essential. In this study, the structural changes of DhaA in two typical aprotic co-solvents (acetonitrile and tetrahydrofuran) were firstly investigated by molecular dynamics (MD) simulation. The simulation results revealed the strong van der Waals force between co-solvents and DhaA which could induce the structural change of enzyme. Interestingly, the differences of molecular size and the electrostatic force with enzyme of two co-solvents led to quite different influences on DhaA. As for acetonitrile, solvent molecules could penetrate into the catalytic site of DhaA which promoted by the electrostatic interaction. On the contrary, tetrahydrofuran molecules were mainly distributed around the catalytic site due to the relative weak electrostatic interaction and steric resistance effect. It can be concluded that different co-solvent can affect the key domains, substrate pathway and catalytic pocket of DhaA.Communicated by Ramaswamy H. Sarma.

Enzyme-Substrate-Cofactor Dynamical Networks Revealed by High-Resolution Field Cycling Relaxometry

The remarkable power and specificity of enzyme catalysis rely on the dynamic alignment of the enzyme, substrates, and cofactors, yet the role of dynamics has usually been approached from the perspective of the protein. We have been using an underappreciated NMR technique, subtesla high-resolution field cycling 31P NMR relaxometry, to investigate the dynamics of the enzyme-bound substrates and cofactor on guanosine-5'-monophosphate reductase (GMPR). GMPR forms two dead end, yet catalytically competent, complexes that mimic distinct steps in the catalytic cycle: E·IMP·NADP+ undergoes a partial hydride transfer reaction, while E·GMP·NADP+ undergoes a partial deamination reaction. A different cofactor conformation is required for each partial reaction. Here we report the effects of mutations designed to perturb cofactor conformation and ammonia binding with the goal of identifying the structural features that contribute to the distinct dynamic signatures of the hydride transfer and deamination complexes. These experiments suggest that Asp129 is a central cog in a dynamic network required for both hydride transfer and deamination. In contrast, Lys77 modulates the conformation and mobility of substrates and cofactors in a reaction-specific manner. Thr105 and Tyr318 are part of a deamination-specific dynamic network that includes the 2'-OH of GMP. These residues have comparatively little effect on the dynamic properties of the hydride transfer complex. These results further illustrate the potential of high-resolution field cycling NMR relaxometry for the investigation of ligand dynamics. In addition, exchange experiments indicate that NH3/NH4+ has a high affinity for the deamination complex but a low affinity for the hydride transfer complex, suggesting that the movement of ammonia may gate the cofactor conformational change. Collectively, these experiments reinforce the view that the enzyme, substrates, and cofactor are linked in intricate, reaction-specific, dynamic networks and demonstrate that distal portions of the substrates and cofactors are critical features in these networks.

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.

Hydrogen deuterium exchange defines catalytically linked regions of protein flexibility in the catechol O-methyltransferase reaction

Human catechol O-methyltransferase (COMT) has emerged as a model for understanding enzyme-catalyzed methyl transfer from S-adenosylmethionine (AdoMet) to small-molecule catecholate acceptors. Mutation of a single residue (tyrosine 68) behind the methyl-bearing sulfonium of AdoMet was previously shown to impair COMT activity by interfering with methyl donor-acceptor compaction within the activated ground state of the wild type enzyme [J. Zhang, H. J. Kulik, T. J. Martinez, J. P. Klinman, Proc. Natl. Acad. Sci. U.S.A. 112, 7954-7959 (2015)]. This predicts the involvement of spatially defined protein dynamical effects that further tune the donor/acceptor distance and geometry as well as the electrostatics of the reactants. Here, we present a hydrogen/deuterium exchange (HDX)-mass spectrometric study of wild type and mutant COMT, comparing temperature dependences of HDX against corresponding kinetic and cofactor binding parameters. The data show that the impaired Tyr68Ala mutant displays similar breaks in Arrhenius plots of both kinetic and HDX properties that are absent in the wild type enzyme. The spatial resolution of HDX below a break point of 15-20 °C indicates changes in flexibility across ∼40% of the protein structure that is confined primarily to the periphery of the AdoMet binding site. Above 20 °C, Tyr68Ala behaves more like WT in HDX, but its rate and enthalpic barrier remain significantly altered. The impairment of catalysis by Tyr68Ala can be understood in the context of a mutationally induced alteration in protein motions that becomes manifest along and perpendicular to the primary group transfer coordinate.

Sensitivity of Functional Loop Conformations on Long-Range Electrostatics: Implications for M20 Loop Dynamics in E. coli Dihydrofolate Reductase

In E. coli dihydrofolate reductase, unusual conformational motions of a functional M20 loop that interacts with substrate and coenzyme have been construed as evidence for dynamical effects in enzyme catalysis. By computing this loop's conformational free energies in the apoenzyme, we show that it is sensitive to the treatment of long-range electrostatic interactions and the solvation box size in modeling/simulations. These results provide important guidelines for computing reaction/binding free energy profiles of proteins with functional loops.

Redox-Responsive Host-Guest Chemistry of a Flexible Cage with Naphthalene Walls

"Naphthocage", a naphthalene-based organic cage, reveals very strong binding (up to 10¹⁰ M^-1) to aromatic (di)cationic guests, i.e., the tetrathiafulvalene mono- and dication and methyl viologen. Intercalation of the guests between two naphthalene walls is mediated by C-H···O, C-H···π, and cation···π interactions. The guests can be switched into and out of the cage by redox processes with high binding selectivity. Oxidation of the flexible cage itself in the absence of a guest leads to a stable radical cation with the oxidized naphthalene intercalated between and stabilized by the other two. Encapsulated guest cations are released from the cavity upon cage oxidation, paving the way to future applications in redox-controlled guest release or novel stimuli-responsive materials.

Effect of Protein Isotope Labeling on the Catalytic Mechanism of Lactate Dehydrogenase

We investigate how isotopic labeling of the enzyme lactate dehydrogenase (LDH) affects its function. LDH is of special interest because there exists a line of residues spanning the protein that are involved in the transition state (TS) of the chemical reaction coordinate (so-called promoting vibration). Hence, studies have been carried out on this protein (as well as others) using labeled protein (so-called heavy protein) along with measurements of single turnover k_cat yielding a KIE (=k_cat^light/k_cat^heavy) aimed at understanding the effect of labeling generally and more specifically this line of residues. Here, it is shown that ¹³C, ¹⁵N, and ²H atom labeling of hhLDH (human heart) affects its internal structure which in turn affects its dynamics and catalytic mechanism. Spectral studies employing advanced FTIR difference spectroscopy show that the height of the electronic potential surface of the TS is lowered (probably by ground state destabilization) by labeling. Moreover, laser-induced T-jump relaxation kinetic spectroscopy shows that the microsecond to millisecond nuclear motions internal to the protein are affected by labeling. While the effects are small, they are sufficient to contribute to the observed KIE values as well or even more than promoting vibration effects.

Can Conformational Changes of Proteins Be Represented in Torsion Angle Space? A Study with Rescaled Ridge Regression

Torsion angles are the natural degrees of freedom of protein structures. The ability to determine torsional variations corresponding to observed changes in Cartesian coordinates is highly valuable, notably to investigate the mechanisms of functional conformational changes or to develop computational models of protein dynamics. This issue is far from trivial in practice since the impact of modifying one torsion angle strongly depends on all other angles, and the compounding effects of small variations in bond lengths and valence angles can completely disrupt a protein fold. We demonstrate that naive strategies, such as directly comparing torsion angles between structures without correcting for variations in bond lengths and valence angles or fitting torsional variations without a proper regularization scheme, fail at producing an adequate representation of conformational changes in internal coordinates. In contrast, rescaled ridge regression, a method recently introduced to regularize multidimensional regressions with correlated explanatory variables, is shown to consistently identify a minimal set of torsion angles variations that closely reproduce changes in Cartesian coordinates. This torsional representation of conformational changes is shown to be robust to the choice of experimental structures. It also provides a better agreement with theoretical models of protein dynamics than the Cartesian representation, regarding notably the predominance of low-frequency normal modes in functional motions and the presence, in predicted equilibrium dynamics, of hints of natural selection for specific functional motions. The software is available at https://github.com/ugobas/tnm .

Interaction vs Preorganization in Enzyme Catalysis. A Dispute That Calls for Resolution

This essay focuses on the debate between Warshel et al. (proponents of preorganization) and Menger and Nome (proponents of spatiotemporal effects) over the source of fast enzyme catalysis. The Warshel model proposes that the main function of enzymes is to push the solvent coordinate toward the transition state. Other physical-organic factors (e.g., desolvation, entropic effects, ground state destabilization, etc.) do not, ostensibly, contribute substantially to the rate. Indeed, physical organic chemistry in its entirety was claimed to be "irrelevant to an enzyme's active site". Preorganization had been applied by Warshel to his "flagship" enzyme, ketosteroid isomerase, but we discuss troubling issues with their ensuing analysis. For example, the concepts of "general acid" and "general base", known to play a role in this enzyme's mechanism, are ignored in the text. In contrast, the spatiotemporal theory postulates that enzyme-like rates (i.e., accelerations >10⁸) occur when two functionalities are held rigidly at contact distances less than ca. 3 Å. Numerous diverse organic systems are shown to bear this out experimentally. Many of these are intramolecular systems where distances between functionalities are known. Among them are fast intramolecular systems where strain is actually generated during the reaction, thereby excluding steric compression as a source of the observed enzyme-like rates. Finally, the account ends with structural data from four active sites of enzymes, obtained by others, all showing contact distances between substrate analogues and enzyme. To our knowledge, contact distances less than the diameter of water are found universally among enzymes, and it is to this fact that we attribute their extremely fast rates given the assumption that enzymes, whatever their particular mechanism, obey elementary chemical principles.

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.