Perspective: Defining and quantifying the role of dynamics in enzyme catalysis

Human indoleamine-2,3-dioxygenase 2 cofactor lability and low substrate affinity explained by homology modeling, molecular dynamics and molecular docking

The human indoleamine-2,3-dioxygenase 2 (hIDO2) protein is growing of interest as it is increasingly implicated in multiple diseases (cancer, autoimmune diseases, COVID-19). However, it is only poorly reported in the literature. Its mode of action remains unknown because it does not seem to catalyze the reaction for which it is attributed: the degradation of the L-Tryptophan into N-formyl-kynurenine. This contrasts with its paralog, the human indoleamine-2,3-dioxygenase 1 (hIDO1), which has been extensively studied in the literature and for which several inhibitors are already in clinical trials. Yet, the recent failure of one of the most advanced hIDO1 inhibitors, the Epacadostat, could be caused by a still unknown interaction between hIDO1 and hIDO2. In order to better understand the mechanism of hIDO2, and in the absence of experimental structural data, a computational study mixing homology modeling, Molecular Dynamics, and molecular docking was conducted. The present article highlights an exacerbated lability of the cofactor as well as an inadequate positioning of the substrate in the active site of hIDO2, which might bring part of an answer to its lack of activity.Communicated by Ramaswamy H. Sarma.

Cryo-EM and Molecular Dynamics Simulations Reveal Hidden Conformational Dynamics Controlling Ammonia Transport in Human Asparagine Synthetase

How motions in enzymes might be linked to catalytic function is of considerable general interest. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating functionally relevant motions in proteins that are not easily amenable to study by other biophysical methods. Here we use 3D variability analysis (3DVA) on cryo-EM maps for wild type (WT) human asparagine synthetase (ASNS) and the R142I ASNS variant to identify conformational changes in the Arg-142 side chain, which mediates the formation of a catalytically relevant intramolecular tunnel. Our 3DVA results for WT ASNS are consistent with independent molecular dynamics (MD) simulations on a model generated from the X-ray structure of human ASNS. Moreover, MD simulations of computational models for the ASNS/β-aspartyl-AMP/MgPPi and R142I/β-aspartyl-AMP/MgPPi ternary complexes, suggest that the structural integrity of the tunnel is impaired in the R142I variant when β-aspartyl-AMP is present in the synthetase active site. The kinetic properties of the R142I ASNS variant support the proposed function of Arg-142. These studies illustrate the power of cryo-EM to identify localized motions and dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics might regulate function in multi-domain enzymes possessing multiple active sites.

Allosteric activation unveils protein-mass modulation of ATP phosphoribosyltransferase product release

Heavy-isotope substitution into enzymes slows down bond vibrations and may alter transition-state barrier crossing probability if this is coupled to fast protein motions. ATP phosphoribosyltransferase from Acinetobacter baumannii is a multi-protein complex where the regulatory protein HisZ allosterically enhances catalysis by the catalytic protein HisGS. This is accompanied by a shift in rate-limiting step from chemistry to product release. Here we report that isotope-labelling of HisGS has no effect on the nonactivated reaction, which involves negative activation heat capacity, while HisZ-activated HisGS catalytic rate decreases in a strictly mass-dependent fashion across five different HisGS masses, at low temperatures. Surprisingly, the effect is not linked to the chemical step, but to fast motions governing product release in the activated enzyme. Disruption of a specific enzyme-product interaction abolishes the isotope effects. Results highlight how altered protein mass perturbs allosterically modulated thermal motions relevant to the catalytic cycle beyond the chemical step.

Structure and dynamics of liquid water from ab initio simulations: adding Minnesota density functionals to Jacob's ladder

The accurate representation of the structural and dynamical properties of water is essential for simulating the unique behavior of this ubiquitous solvent. Here we assess the current status of describing liquid water using ab initio molecular dynamics, with a special focus on the performance of all the later generation Minnesota functionals. Findings are contextualized within the current knowledge on DFT for describing bulk water under ambient conditions and compared to experimental data. We find that, contrary to the prevalent idea that local and semilocal functionals overstructure water and underestimate dynamical properties, M06-L, revM06-L, and M11-L understructure water, while MN12-L and MN15-L overdistance water molecules due to weak cohesive effects. This can be attributed to a weakening of the hydrogen bond network, which leads to dynamical fingerprints that are over fast. While most of the hybrid Minnesota functionals (M06, M08-HX, M08-SO, M11, MN12-SX, and MN15) also yield understructured water, their dynamical properties generally improve over their semilocal counterparts. It emerges that exact exchange is a crucial component for accurately describing hydrogen bonds, which ultimately leads to corrections in both the dynamical and structural properties. However, an excessive amount of exact exchange strengthens hydrogen bonds and causes overstructuring and slow dynamics (M06-HF). As a compromise, M06-2X is the best performing Minnesota functional for water, and its D3 corrected variant shows very good structural agreement. From previous studies considering nuclear quantum effects (NQEs), the hybrid revPBE0-D3, and the rung-5 RPA (RPA@PBE) have been identified as the only two approximations that closely agree with experiments. Our results suggest that the M06-2X(-D3) functionals have the potential to further improve the reproduction of experimental properties when incorporating NQEs through path integral approaches. This work provides further proof that accurate modeling of water interactions requires the inclusion of both exact exchange and balanced (non-local) correlation, highlighting the need for higher rungs on Jacob's ladder to achieve predictive simulations of complex biological systems in aqueous environments.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study, we employ this paradigm to answer a basic question: in enzyme superfamilies, where the catalytic mechanism, active sites, and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as representatives of the conserved protein tyrosine phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of the catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

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

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

Ratcheting synthesis

Synthetic chemistry has traditionally relied on reactions between reactants of high chemical potential and transformations that proceed energetically downhill to either a global or local minimum (thermodynamic or kinetic control). Catalysts can be used to manipulate kinetic control, lowering activation energies to influence reaction outcomes. However, such chemistry is still constrained by the shape of one-dimensional reaction coordinates. Coupling synthesis to an orthogonal energy input can allow ratcheting of chemical reaction outcomes, reminiscent of the ways that molecular machines ratchet random thermal motion to bias conformational dynamics. This fundamentally distinct approach to synthesis allows multi-dimensional potential energy surfaces to be navigated, enabling reaction outcomes that cannot be achieved under conventional kinetic or thermodynamic control. In this Review, we discuss how ratcheted synthesis is ubiquitous throughout biology and consider how chemists might harness ratchet mechanisms to accelerate catalysis, drive chemical reactions uphill and programme complex reaction sequences.

Internal charge-transfer in a metal-catalyzed oxidative addition reaction turns an inhibitive electric field stimulus to catalytic

In a metal-catalyzed oxidative addition, an oriented external electric field (EEF) catalyzes the reaction along one direction and inhibits it when applied in the opposite direction. Beyond a threshold value, the inhibitory direction becomes catalyzing by swapping the metal-to-ligand charge transfer (MLCT) to ligand-to-metal charge-transfer (LMCT) or vice versa. The change in direction of the charge-transfer mechanism triggers the inversion of the dipole moment along the reaction axis, that results in the resurgence of catalysis. The charge-transfer mechanism in metal-catalyzed oxidative addition is tunable by EEF.

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.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

Allostery in Protein Tyrosine Phosphatases is Enabled by Divergent Dynamics

Dynamics-driven allostery provides important insights into the working mechanics of proteins, especially enzymes. In this study we employ this paradigm to answer a basic question: in enzyme superfamilies where the catalytic mechanism, active sites and protein fold are conserved, what accounts for the difference in the catalytic prowess of the individual members? We show that when subtle changes in sequence do not translate to changes in structure, they do translate to changes in dynamics. We use sequentially diverse PTP1B, TbPTP1, and YopH as the representatives of the conserved Protein Tyrosine Phosphatase (PTP) superfamily. Using amino acid network analysis of group behavior (community analysis) and influential node dominance on networks (eigenvector centrality), we explain the dynamic basis of catalytic variations seen between the three proteins. Importantly, we explain how a dynamics-based blueprint makes PTP1B amenable to allosteric control and how the same is abstracted in TbPTP1 and YopH.

The Conformational Transitions and Dynamics of Burkholderia cepacia Lipase Regulated by Water-Oil Interfaces

Structural dynamics and conformational transitions are crucial for the activities of enzymes. As one of the most widely used industrial biocatalysts, lipase could be activated by the water-oil interfaces. The interface activations were believed to be dominated by the close-to-open transitions of the lid subdomains. However, the detailed mechanism and the roles of structure transitions are still under debate. In this study, the dynamic structures and conformational transitions of Burkholderia cepacia lipase (LipA) were investigated by combining all-atom molecular dynamics simulations, enhanced sampling simulation, and spectrophotometric assay experiments. The conformational transitions between the lid-open and lid-closed states of LipA in aqueous solution are directly observed by the computational simulation methods. The interactions between the hydrophobic residues on the two lid-subdomains are the driven forces for the LipA closing. Meanwhile, the hydrophobic environment provided by the oil interfaces would separate the interactions between the lid-subdomains and promote the structure opening of LipA. Moreover, our studies demonstrate the opening of the lids structure is insufficient to initiate the interfacial activation, providing explanations for the inability of interfacial activation of many lipases with lid structures.

Enzyme dynamics-a brief review

Debate has been simmering for some years regarding the importance of internal thermal motions of enzymes to catalysis. Recent developments in protein design may bring resolution of the more contentious points a little closer.

Perspective: How Fast Dynamics Affect Slow Function in Protein Machines

Internal motions in proteins take place on a broad range of time- and space-scales. The potential roles of these dynamics in the biochemical functions of proteins have intrigued biophysicists for many years, and multiple mechanisms to couple motions to function have been proposed. Some of these mechanisms have relied on equilibrium concepts. For example, the modulation of dynamics was proposed to change the entropy of a protein, hence affecting processes such as binding. This so-called dynamic allostery scenario has been demonstrated in several recent experiments. Perhaps even more intriguing may be models that involve out-of-equilibrium operation, which by necessity require the input of energy. We discuss several recent experimental studies that expose such potential mechanisms for coupling dynamics and function. In Brownian ratchets, for example, directional motion is promoted by switching a protein between two free energy surfaces. An additional example involves the effect of microsecond domain-closure dynamics of an enzyme on its much slower chemical cycle. These observations lead us to propose a novel two-time-scale paradigm for the activity of protein machines: fast equilibrium fluctuations take place on the microsecond-millisecond time scale, while on a slower time scale, free energy is invested in order to push the system out of equilibrium and drive functional transitions. Motions on the two time scales affect each other and are essential for the overall function of these machines.

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.

Protein Dynamics and Enzymatic Catalysis

This Perspective presents a review of our work and that of others in the highly controversial topic of the coupling of protein dynamics to reaction in enzymes. We have been involved in studying this topic for many years. Thus, this perspective will naturally present our own views, but it also is designed to present an overview of the variety of viewpoints of this topic, both experimental and theoretical. This is obviously a large and contentious topic.

Consolidation of temperature-dependent toxicity and thermoregulatory behavior into risk assessments of insecticides under thermal scenarios: A prospective study on Eremias argus

In this study, the temperature-dependent chemical toxicity of three insecticides and the resulting thermoregulatory (TR) behavior of the lizard Eremias argus have been consolidated into the current risk assessment framework. According to acute dermal toxicity assays, an increase of ambient temperature from 15 °C to 35 °C decreased the acute dermal toxicity of beta-cyfluthrin (BC) but increased the toxicity of chlorpyrifos (CPF). The toxicity of avermectin (AVM) did not show significant temperature-dependent responses. Based on thermal preference trials, lizards changed their body temperature via TR behavior to adaptively reduce toxicity under sub-lethal doses, which can be understood as a "self-rescue" behavior attenuating lethal effects. However, the risk quotient indicated that the effectiveness of this "self-rescue" behavior is limited. Metabolomics analysis showed that six different metabolites (i.e., creatine, glutamate, succinate, N-acetylaspartate, acetylcholine, and lactate) contributed to TR behavior changes. Biochemical assays and insecticide residue results demonstrated that the temperature-dependent toxicity of BC, CPF, and AVM affected lizards in the three aspects of biotransformation, oxidative stress, and neurometabolic interference. This work clarifies the ecotoxicological impacts of representative insecticides on reptiles from toxicological understanding to risk relevance. This knowledge may improve ecological predictions of agrochemical applications in the context of global climate change.

Enzymatische C4-Epimerisierung von UDP-Glucuronsäure: präzise gesteuerte Rotation eines transienten 4-Ketointermediats für eine invertierende Reaktion ohne Decarboxylierung

UDP‐Glucuronsäure(UDP‐GlcA)‐4‐Epimerase repräsentiert eine wichtige Fragestellung in der Enzymkatalyse: die Balance zwischen konformativer Flexibilität und genauer Positionierung. Das Enzym koordiniert die C4‐Oxidation des Substrats durch NAD+ mit der Rotation eines leicht decarboxylierbaren β‐Ketosäure‐Intermediats im aktiven Zentrum zur Ermöglichung der stereoinvertierenden Reduktion der Ketogruppe durch NADH. Wir zeigen hier die nur schwer erfassbare Rotationskoordinate des 4‐Ketointermediats. Distorsion des Zuckerrings in eine Boot‐Konformation erzeugt torsionale Mobilität in der Bindungstasche des Enzyms. Die Endpunkte der Rotation zeigen den 4‐Ketozucker in einer unverformten 4C1‐Sesselkonformation. Die äquatorial positionierte Carboxylatgruppe ist ungünstig für die 4‐Ketozucker‐Decarboxylierung. Varianten der Epimerase zeigen Decarboxylierung, wenn sie die Bindung mit der Carboxylatgruppe im entgegengesetzten Rotationsisomer des Substrats entfernen. R185A/D‐Substitutionen wandeln die Epimerase in UDP‐Xylose‐Synthasen um, welche UDP‐GlcA in stereospezifischen, konfigurationserhaltenden Reaktionen decarboxylieren.

Enzymatic C4-Epimerization of UDP-Glucuronic Acid: Precisely Steered Rotation of a Transient 4-Keto Intermediate for an Inverted Reaction without Decarboxylation

UDP-glucuronic acid (UDP-GlcA) 4-epimerase illustrates an important problem regarding enzyme catalysis: balancing conformational flexibility with precise positioning. The enzyme coordinates the C4-oxidation of the substrate by NAD⁺ and rotation of a decarboxylation-prone β-keto acid intermediate in the active site, enabling stereoinverting reduction of the keto group by NADH. We reveal the elusive rotational landscape of the 4-keto intermediate. Distortion of the sugar ring into boat conformations induces torsional mobility in the enzyme's binding pocket. The rotational endpoints show that the 4-keto sugar has an undistorted ⁴ C₁ chair conformation. The equatorially placed carboxylate group disfavors decarboxylation of the 4-keto sugar. Epimerase variants lead to decarboxylation upon removal of the binding interactions with the carboxylate group in the opposite rotational isomer of the substrate. Substitutions R185A/D convert the epimerase into UDP-xylose synthases that decarboxylate UDP-GlcA in stereospecific, configuration-retaining reactions.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Allosteric rescue of catalytically impaired ATP phosphoribosyltransferase variants links protein dynamics to active-site electrostatic preorganisation

ATP phosphoribosyltransferase catalyses the first step of histidine biosynthesis and is controlled via a complex allosteric mechanism where the regulatory protein HisZ enhances catalysis by the catalytic protein HisG_S while mediating allosteric inhibition by histidine. Activation by HisZ was proposed to position HisG_S Arg56 to stabilise departure of the pyrophosphate leaving group. Here we report active-site mutants of HisG_S with impaired reaction chemistry which can be allosterically restored by HisZ despite the HisZ:HisG_S interface lying ~20 Å away from the active site. MD simulations indicate HisZ binding constrains the dynamics of HisG_S to favour a preorganised active site where both Arg56 and Arg32 are poised to stabilise leaving-group departure in WT-HisG_S. In the Arg56Ala-HisG_S mutant, HisZ modulates Arg32 dynamics so that it can partially compensate for the absence of Arg56. These results illustrate how remote protein-protein interactions translate into catalytic resilience by restoring damaged electrostatic preorganisation at the active site.

Diversity oriented clicking delivers β-substituted alkenyl sulfonyl fluorides as covalent human neutrophil elastase inhibitors

Diversity Oriented Clicking (DOC) is a discovery method geared toward the rapid synthesis of functional libraries. It combines the best attributes of both classical and modern click chemistries. DOC strategies center upon the chemical diversification of core "SuFExable" hubs-exemplified by 2-Substituted-Alkynyl-1-Sulfonyl Fluorides (SASFs)-enabling the modular assembly of compounds through multiple reaction pathways. We report here a range of stereoselective Michael-type addition pathways from SASF hubs including reactions with secondary amines, carboxylates, 1H-1,2,3-triazole, and halides. These high yielding conjugate addition pathways deliver unprecedented β-substituted alkenyl sulfonyl fluorides as single isomers with minimal purification, greatly enriching the repertoire of DOC and holding true to the fundamentals of modular click chemistry. Further, we demonstrate the potential for biological function - a key objective of click chemistry - of this family of SASF-derived molecules as covalent inhibitors of human neutrophil elastase.

A method for quantifying how the activity of an enzyme is affected by the net charge of its nearest crowded neighbor

The electrostatic effects of protein crowding have not been systematically explored. Rather, protein crowding is generally studied with co‐solvents or crowders that are electrostatically neutral, with no methods to measure how the net charge (Z) of a crowder affects protein function. For example, can the activity of an enzyme be affected electrostatically by the net charge of its neighbor in crowded milieu? This paper reports a method for crowding proteins of different net charge to an enzyme via semi‐random chemical crosslinking. As a proof of concept, RNase A was crowded (at distances ≤ the Debye length) via crosslinking to different heme proteins with Z = +8.50 ± 0.04, Z = +6.39 ± 0.12, or Z = −10.30 ± 1.32. Crosslinking did not disrupt the structure of proteins, according to amide H/D exchange, and did not inhibit RNase A activity. For RNase A, we found that the electrostatic environment of each crowded neighbor had significant effects on rates of RNA hydrolysis. Crowding with cationic cytochrome c led to increases in activity, while crowding with anionic “supercharged” cytochrome c or myoglobin diminished activity. Surprisingly, electrostatic crowding effects were amplified at high ionic strength (I = 0.201 M) and attenuated at low ionic strength (I = 0.011 M). This salt dependence might be caused by a unique set of electric double layers at the dimer interspace (maximum distance of 8 Å, which cannot accommodate four layers). This new method of crowding via crosslinking can be used to search for electrostatic effects in protein crowding.

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.

Modulation of α-Chymotrypsin Conjugated to Magnetic Nanoparticles by the Non-Heating Low-Frequency Magnetic Field: Molecular Dynamics, Reaction Kinetics, and Spectroscopy Analysis

Enzymes conjugated to magnetic nanoparticles (MNPs) undergo changes in the catalytic activity of the non-heating low-frequency magnetic field (LFMF). We apply in silico simulations by molecular dynamics (MD) and in vitro spectroscopic analysis of the enzyme kinetics and secondary structure to study α-chymotrypsin (CT) conjugated to gold-coated iron oxide MNPs. The latter are functionalized by either carboxylic or amino group moieties to vary the points of enzyme attachment. The MD simulation suggests that application of the stretching force to the CT globule by its amino or carboxylic groups causes shrinkage of the substrate-binding site but little if any changes in the catalytic triad. Consistent with this, in CT conjugated to MNPs by either amino or carboxylic groups, LFMF alters the Michaelis-Menten constant but not the apparent catalytic constant k _cat (= V _max/[E]_o). Irrespective of the point of conjugation to MNPs, the CT secondary structure was affected with nearly complete loss of α-helices and increase in the random structures in LFMF, as shown by attenuated total reflection Fourier transformed infrared spectroscopy. Both the catalytic activity and the protein structure of MNP-CT conjugates restored 3 h after the field exposure. We believe that such remotely actuated systems can find applications in advanced manufacturing, nanomedicine, and other areas.

Pairing structural reconstruction with catalytic competence to evaluate the mechanisms of key enzymes in the folate-mediated one-carbon pathway

Mammalian metabolism comprises a series of interlinking pathways that include two major cycles: the folate and methionine cycles. The folate-mediated metabolic cycle uses several oxidation states of tetrahydrofolate to carry activated one-carbon units to be readily used and interconverted within the cell. They are required for nucleotide synthesis, methylation and metabolism, and particularly for proliferation of cancer cells. Based on the latest progress in genome-wide CRISPR loss-of-function viability screening of 789 cell lines, we focus on the most cancer-dependent enzymes in this pathway, especially those that are hyperactivated in cancer, to provide new insight into the chemical basis for cancer drug development. Since the complete 3D structure of several of these enzymes of the one-carbon pathway in their active form are not available, we used homology modelling integrated with the interpretation of the reaction mechanism. In addition, have reconstructed the most likely scenario for the reactions taking place paired with their catalytic competence that provides a testable framework for this pathway.

Why Monoamine Oxidase B Preferably Metabolizes N-Methylhistamine over Histamine: Evidence from the Multiscale Simulation of the Rate-Limiting Step

Histamine levels in the human brain are controlled by rather peculiar metabolic pathways. In the first step, histamine is enzymatically methylated at its imidazole N^τ atom, and the produced N-methylhistamine undergoes an oxidative deamination catalyzed by monoamine oxidase B (MAO-B), as is common with other monoaminergic neurotransmitters and neuromodulators of the central nervous system. The fact that histamine requires such a conversion prior to oxidative deamination is intriguing since MAO-B is known to be relatively promiscuous towards monoaminergic substrates; its in-vitro oxidation of N-methylhistamine is about 10 times faster than that for histamine, yet this rather subtle difference appears to be governing the decomposition pathway. This work clarifies the MAO-B selectivity toward histamine and N-methylhistamine by multiscale simulations of the rate-limiting hydride abstraction step for both compounds in the gas phase, in aqueous solution, and in the enzyme, using the established empirical valence bond methodology, assisted by gas-phase density functional theory (DFT) calculations. The computed barriers are in very good agreement with experimental kinetic data, especially for relative trends among systems, thereby reproducing the observed MAO-B selectivity. Simulations clearly demonstrate that solvation effects govern the reactivity, both in aqueous solution as well as in the enzyme although with an opposing effect on the free energy barrier. In the aqueous solution, the transition-state structure involving histamine is better solvated than its methylated analog, leading to a lower barrier for histamine oxidation. In the enzyme, the higher hydrophobicity of N-methylhistamine results in a decreased number of water molecules at the active side, leading to decreased dielectric shielding of the preorganized catalytic electrostatic environment provided by the enzyme. This renders the catalytic environment more efficient for N-methylhistamine, giving rise to a lower barrier relative to histamine. In addition, the transition state involving N-methylhistamine appears to be stabilized by the surrounding nonpolar residues to a larger extent than with unsubstituted histamine, contributing to a lower barrier with the former.

Evolution of Optimized Hydride Transfer Reaction and Overall Enzyme Turnover in Human Dihydrofolate Reductase

Evolution of dihydrofolate reductase (DHFR) has been studied using the enzyme from Escherichia coli DHFR (ecDHFR) as a model, but less studies have used the enzyme from Homo sapiens DHFR (hsDHFR). Each enzyme maintains a short and narrow distribution of hydride donor-acceptor distances (DAD) at the tunneling ready state (TRS). Evolution of the enzyme was previously studied in ecDHFR where three key sites were identified as important to the catalyzed reaction. The corresponding sites in hsDHFR are F28, 62-PEKN, and 26-PPLR. Each of these sites was studied here through the creation of mutant variants of the enzyme and measurements of the temperature dependence of the intrinsic kinetic isotope effects (KIEs) on the reaction. F28 is mutated first to M (F28M) and then to the L of the bacterial enzyme (F28L). The KIEs of the F28M variant are larger and more temperature-dependent than wild-type (WT), suggesting a broader and longer average DAD at the TRS. To more fully mimic ecDHFR, we also study a triple mutant of the human enzyme (F32L-PP26N-PEKN62G). Remarkably, the intrinsic KIEs, while larger in magnitude, are temperature-independent like the WT enzymes. We also construct deletion mutations of hsDHFR removing both the 62-PEKN and 26-PPLR sequences. The results mirror those described previously for insertion mutants of ecDHFR. Taken together, these results suggest a balancing act during DHFR evolution between achieving an optimal TRS for hydride transfer and preventing product inhibition arising from the different intercellular pools of NADPH and NADP⁺ in prokaryotic and eukaryotic cells.

Mapping Enzyme Landscapes by Time-Resolved Crystallography with Synchrotron and X-Ray Free Electron Laser Light

Directly observing enzyme catalysis in real time at the molecular level has been a long-standing goal of structural enzymology. Time-resolved serial crystallography methods at synchrotron and X-ray free electron laser (XFEL) sources have enabled researchers to follow enzyme catalysis and other nonequilibrium events at ambient conditions with unprecedented time resolution. X-ray crystallography provides detailed information about conformational heterogeneity and protein dynamics, which is enhanced when time-resolved approaches are used. This review outlines the ways in which information about the underlying energy landscape of a protein can be extracted from X-ray crystallographic data, with an emphasis on new developments in XFEL and synchrotron time-resolved crystallography. The emerging view of enzyme catalysis afforded by these techniques can be interpreted as enzymes moving on a time-dependent energy landscape. Some consequences of this view are discussed, including the proposal that irreversible enzymes or enzymes that use covalent catalytic mechanisms may commonly exhibit catalysis-activated motions. Expected final online publication date for the Annual Review of Biophysics, Volume 51 is May 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Nanosecond structural dynamics of intrinsically disordered β-casein micelles by neutron spectroscopy

β-casein undergoes a reversible endothermic self-association, forming protein micelles of limited size. In its functional state, a single β-casein monomer is unfolded, which creates a high structural flexibility, which is supposed to play a major role in preventing the precipitation of calcium phosphate particles. We characterize the structural flexibility in terms of nanosecond molecular motions, depending on the temperature by quasielastic neutron scattering. Our major questions are: Does the self-association reduce the chain flexibility? How does the dynamic spectrum of disordered caseins differ from a compactly globular protein? How does the dynamic spectrum of β-casein in solution differ from that of a protein in hydrated powder states? We report on two relaxation processes on a nanosecond and a sub-nanosecond timescale for β-casein in solution. Both processes are analyzed by Brownian oscillator model, by which the spring constant can be defined in the isotropic parabolic potential. The slower process, which is analyzed by neutron spin echo, seems a characteristic feature of the unfolded structure. It requires bulk solvent and is not seen in hydrated protein powders. The faster process, which is analyzed by neutron backscattering, has a smaller amplitude and requires hydration water, which is also observed with folded proteins in the hydrated state. The self-association had no significant influence on internal relaxation, and thus, a β-casein protein monomer flexibility is preserved in the micelle. We derive spring constants of the faster and slower motions of β-caseins in solution and compared them with those of some proteins in various states (folded or hydrated powder).

An electrostatic duel: subtle differences in the catalytic performance of monoamine oxidase A and B isoenzymes elucidated at the residue level using quantum computations

The origin of the immense catalytic power of enzymes remains one of the biggest unresolved questions in biochemistry, with electrostatics being one of the main contenders. Herein, we report results that not only confirm that electrostatics is the driving force behind enzyme catalysis, but also that it is capable of tuning subtle differences in the catalytic performance between structurally similar enzymes, as demonstrated using the example of isoenzymes, monoamine oxidases A and B. Using our own computationally efficient multiscale model [A. Prah, et al., ACS Catal., 2019, 9, 1231] we analyzed the rate-limiting step of the reaction between phenylethylamine and both isoenzymes and deduced that the electrostatic environment provided by isoenzyme B has a perceivably higher catalytic influence on all the considered parameters of the reaction (energy barrier, charge transfer, dipole moment, and HOMO-LUMO gap). This is in full agreement with the available experimental kinetic data and with our own simulations of the reaction in question. In-depth analysis of individual amino acid contributions of both isoenzymes to the barrier (based on the interaction between the electric field provided by the enzyme and the dipole moment of the reacting moiety) shows that the majority of the difference between the isoenzymes can be attributed to a small number of sizable differences between the aligned amino acid pairs, whereas in most of the pairs the difference in contribution to the barrier is vanishingly small. These results suggest that electrostatics largely controls the substrate selectivity of enzymes and validates our approach as being capable of discerning fine nuances in the selectivity of structurally related isoenzymes.

Unraveling the Nature of Extraframework Catalytic Ensembles in Zeolites: Flexibility and Dynamics of the Copper-Oxo Trimers in Mordenite

Extraframework cations define the chemical versatility of zeolite catalysts. Addressing their structural complexity and dynamic behavior represents one of the main fundamental challenges in the field. Herein, we present a computational approach for the identification and analysis of the accessible pool of intrazeolite extraframework complexes with a Cu/MOR catalyst as an industrially important model system. We employ ab initio molecular dynamics for capturing the ensemble of reactive isomers with the [Cu₃O₃]²⁺ stoichiometry confined in the mordenite channels. The high structural diversity of the generated isomers was ensured by concentrating the kinetic energy along the low-curvature directions of the potential energy surface (PES). Geometrically distinct [Cu₃O₃]²⁺ complexes were identified via a series of clustering procedures ensuring that one structure of each local minima is retained. The proposed procedure has resulted in a set of previously unknown peroxo-complexes, which are >50 kJ/mol more stable than the recently hypothesized chair-shaped structure. Our analysis demonstrates that the most stable peroxo-containing clusters can be formed under operando conditions from molecular oxygen and the Cu₃O unit, similar to that in methane monooxygenase (MMO) enzymes.

Temperature-Dependent Low-Frequency Modes in the Active Site of Bovine Carbonic Anhydrase II Probed by 2D-IR Spectroscopy

Enzyme catalysis achieves tremendous rate accelerations. Enzyme reaction centers provide a constraint geometry that preferentially binds an activated form of the substrate and thus lowers the energy barrier. However, this transition state picture neglects the flexibility of proteins and its role in enzymatic catalysis. Especially for proton transfer reactions, it has been suggested that motions of the protein modulate the donor-acceptor distance and prepare a tunneling-ready state. We report the detection of frequency fluctuations of an azide anion (N₃^-) bound in the active site of the protein carbonic anhydrase II, where a low-frequency mode of the protein has been proposed to facilitate proton transfer over two water molecules during the catalyzed reaction. 2D-IR spectroscopy resolves an underdamped low-frequency mode at about 1 THz (30 cm^-1). We find its frequency to be viscosity- and temperature-dependent and to decrease by 6 cm^-1 between 230 and 320 K, reporting the softening of the mode's potential.

Conformational dynamics of a multidomain protein by neutron scattering and computational analysis

The flexible conformations of a multidomain protein are responsible for its biological functions. Although MurD, a 47-kDa protein that consists of three domains, sequentially changes its domain conformation from an open form to a closed form through a semiclosed form in its enzymatic reaction, the domain dynamics in each conformation remains unclear. In this study, we verify the conformational dynamics of MurD in the corresponding three states (apo and ATP- and inhibitor-bound states) with a combination of small-angle x-ray and neutron scattering (SAXS and SANS), dynamic light scattering (DLS), neutron backscattering (NBS), neutron spin echo (NSE) spectroscopy, and molecular dynamics (MD) simulations. Applying principal component analysis of the MD trajectories, twisting and open-closed domain modes are identified as the major collective coordinates. The deviations of the experimental SAXS profiles from the theoretical calculations based on the known crystal structures become smaller in the ATP-bound state than in the apo state, and a further decrease is evident upon inhibitor binding. These results suggest that domain motions of the protein are suppressed step by step of each ligand binding. The DLS and NBS data yield collective and self-translational diffusion constants, respectively, and we used them to extract collective domain motions in nanometer and nanosecond scales from the NSE data. In the apo state, MurD shows both twisting and open-closed domain modes, whereas an ATP binding suppresses twisting domain motions, and a further reduction of open-closed mode is seen in the inhibitor-binding state. These observations are consistent with the structure modifications measured by the small-angle scattering as well as the MD simulations. Such changes in the domain dynamics associated with the sequential enzymatic reactions should be related to the affinity and reaction efficiency with a ligand that binds specifically to each reaction state.

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

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

Correlated Motions in Structural Biology

Correlated motions in proteins arising from the collective movements of residues have long been proposed to be fundamentally important to key properties of proteins, from allostery and catalysis to evolvability. Recent breakthroughs in structural biology have made it possible to capture proteins undergoing complex conformational changes, yet intrinsic correlated motions within a conformation remain one of the least understood facets of protein structure. For many decades, the analysis of total X-ray scattering held the promise of animating crystal structures with correlated motions. With recent advances in both X-ray detectors and data interpretation methods, this long-held promise can now be met. In this Perspective, we will introduce how correlated motions are captured in total scattering and provide guidelines for the collection, interpretation, and validation of data. As structural biology continues to push the boundaries, we see an opportunity to gain atomistic insight into correlated motions using total scattering as a bridge between theory and experiment.

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.

Catalytic Descriptors to Investigate Catalytic Power in the Reaction of Haloalkane Dehalogenase Enzyme with 1,2-Dichloroethane

Enzymes play a fundamental role in many biological processes. We present a theoretical approach to investigate the catalytic power of the haloalkane dehalogenase reaction with 1,2-dichloroethane. By removing the three main active-site residues one by one from haloalkane dehalogenase, we found two reactive descriptors: one descriptor is the distance difference between the breaking bond and the forming bond, and the other is the charge difference between the transition state and the reactant complex. Both descriptors scale linearly with the reactive barriers, with the three-residue case having the smallest barrier and the zero-residue case having the largest. The results demonstrate that, as the number of residues increases, the catalytic power increases. The predicted free energy barriers using the two descriptors of this reaction in water are 23.1 and 24.2 kcal/mol, both larger than the ones with any residues, indicating that the water solvent hinders the reactivity. Both predicted barrier heights agree well with the calculated one at 25.2 kcal/mol using a quantum mechanics and molecular dynamics approach, and also agree well with the experimental result at 26.0 kcal/mol. This study shows that reactive descriptors can also be used to describe and predict the catalytic performance for enzyme catalysis.

Cryo-kinetics Reveal Dynamic Effects on the Chemistry of Human Dihydrofolate Reductase

Effects of isotopic substitution on the rate constants of human dihydrofolate reductase (HsDHFR), an important target for anti‐cancer drugs, have not previously been characterized due to its complex fast kinetics. Here, we report the results of cryo‐measurements of the kinetics of the HsDHFR catalyzed reaction and the effects of protein motion on catalysis. Isotopic enzyme labeling revealed an enzyme KIE (k_H^LE/k_H^HE) close to unity above 0 °C; however, the enzyme KIE was increased to 1.72±0.15 at −20 °C, indicating that the coupling of protein motions to the chemical step is minimized under optimal conditions but enhanced at non‐physiological temperatures. The presented cryogenic approach provides an opportunity to probe the kinetics of mammalian DHFRs, thereby laying the foundation for characterizing their transition state structure.

Reciprocally-Coupled Gating: Strange Loops in Bioenergetics, Genetics, and Catalysis

Bioenergetics, genetic coding, and catalysis are all difficult to imagine emerging without pre-existing historical context. That context is often posed as a "Chicken and Egg" problem; its resolution is concisely described by de Grasse Tyson: "The egg was laid by a bird that was not a chicken". The concision and generality of that answer furnish no details-only an appropriate framework from which to examine detailed paradigms that might illuminate paradoxes underlying these three life-defining biomolecular processes. We examine experimental aspects here of five examples that all conform to the same paradigm. In each example, a paradox is resolved by coupling "if, and only if" conditions for reciprocal transitions between levels, such that the consequent of the first test is the antecedent for the second. Each condition thus restricts fluxes through, or "gates" the other. Reciprocally-coupled gating, in which two gated processes constrain one another, is self-referential, hence maps onto the formal structure of "strange loops". That mapping uncovers two different kinds of forces that may help unite the axioms underlying three phenomena that distinguish biology from chemistry. As a physical analog for Gödel's logic, biomolecular strange-loops provide a natural metaphor around which to organize a large body of experimental data, linking biology to information, free energy, and the second law of thermodynamics.

Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts

As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.

Dynamic Structural Biology Experiments at XFEL or Synchrotron Sources

Macromolecular crystallography (MX) leverages the methods of physics and the language of chemistry to reveal fundamental insights into biology. Often beautifully artistic images present MX results to support profound functional hypotheses that are vital to entire life science research community. Over the past several decades, synchrotrons around the world have been the workhorses for X-ray diffraction data collection at many highly automated beamlines. The newest tools include X-ray-free electron lasers (XFELs) located at facilities in the USA, Japan, Korea, Switzerland, and Germany that deliver about nine orders of magnitude higher brightness in discrete femtosecond long pulses. At each of these facilities, new serial femtosecond crystallography (SFX) strategies exploit slurries of micron-size crystals by rapidly delivering individual crystals into the XFEL X-ray interaction region, from which one diffraction pattern is collected per crystal before it is destroyed by the intense X-ray pulse. Relatively simple adaptions to SFX methods produce time-resolved data collection strategies wherein reactions are triggered by visible light illumination or by chemical diffusion/mixing. Thus, XFELs provide new opportunities for high temporal and spatial resolution studies of systems engaged in function at physiological temperature. In this chapter, we summarize various issues related to microcrystal slurry preparation, sample delivery into the X-ray interaction region, and some emerging strategies for time-resolved SFX data collection.

GRaSP: a graph-based residue neighborhood strategy to predict binding sites

Motivation

The discovery of protein–ligand-binding sites is a major step for elucidating protein function and for investigating new functional roles. Detecting protein–ligand-binding sites experimentally is time-consuming and expensive. Thus, a variety of in silico methods to detect and predict binding sites was proposed as they can be scalable, fast and present low cost.

Results

We proposed Graph-based Residue neighborhood Strategy to Predict binding sites (GRaSP), a novel residue centric and scalable method to predict ligand-binding site residues. It is based on a supervised learning strategy that models the residue environment as a graph at the atomic level. Results show that GRaSP made compatible or superior predictions when compared with methods described in the literature. GRaSP outperformed six other residue-centric methods, including the one considered as state-of-the-art. Also, our method achieved better results than the method from CAMEO independent assessment. GRaSP ranked second when compared with five state-of-the-art pocket-centric methods, which we consider a significant result, as it was not devised to predict pockets. Finally, our method proved scalable as it took 10–20 s on average to predict the binding site for a protein complex whereas the state-of-the-art residue-centric method takes 2–5 h on average.

Availability and implementation

The source code and datasets are available at https://github.com/charles-abreu/GRaSP.

Supplementary information

Supplementary data are available at Bioinformatics online.

Distinct conformational dynamics and allosteric networks in alpha tryptophan synthase during active catalysis

Experimental observations of enzymes under active turnover conditions have brought new insight into the role of protein motions and allosteric networks in catalysis. Many of these studies characterize enzymes under dynamic chemical equilibrium conditions, in which the enzyme is actively catalyzing both the forward and reverse reactions during data acquisition. We have previously analyzed conformational dynamics and allosteric networks of the alpha subunit of tryptophan synthase under such conditions using NMR. We have proposed that this working state represents a four to one ratio of the enzyme bound with the indole-3-glycerol phosphate substrate (E:IGP) to the enzyme bound with the products indole and glyceraldehyde-3-phosphate (E:indole:G3P). Here, we analyze the inactive D60N variant to deconvolute the contributions of the substrate- and products-bound states to the working state. While the D60N substitution itself induces small structural and dynamic changes, the D60N E:IGP and E:indole:G3P states cannot entirely account for the conformational dynamics and allosteric networks present in the working state. The act of chemical bond breakage and/or formation, or possibly the generation of an intermediate, may alter the structure and dynamics present in the working state. As the enzyme transitions from the substrate-bound to the products-bound state, millisecond conformational exchange processes are quenched and new allosteric connections are made between the alpha active site and the surface which interfaces with the beta subunit. The structural ordering of the enzyme and these new allosteric connections may be important in coordinating the channeling of the indole product into the beta subunit.

Recent results in time resolved serial femtosecond crystallography at XFELs

Time-resolved serial femtosecond crystallography (tr-SFX) methods exploit slurries of crystalline samples that range in size from hundreds of nanometers to a few tens of micrometers, at near-physiological temperature and pressure, to generate atomic resolution models and probe authentic function with the same experiment. 'Dynamic structural biology' is often used to encompass the research philosophy and techniques. Reaction cycles for tr-SFX studies are initiated by photons or ligand addition/mixing strategies, wherein the latter are potentially generalizable across enzymology. Thus, dynamic structural biology often creates stop-motion molecular movies of macromolecular function. In metal-dependent systems, complementary spectroscopic information can also be collected from the same samples and X-ray pulses, which provides even more detailed mechanistic insights. These types of experimental data also complement quantum mechanical and classical dynamics numerical calculations. Correlated structural-functional results will yield more detailed mechanistic insights and will likely translate into better drugs and treatments impacting human health, and better catalysis for clean energy and agriculture.

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.

How far the vibrational exciton interactions are responsible for the generation of the infrared spectra of oxindole crystals?

Oxindole (indolin-2-one, Ox) is a unique and a crucial molecular system in spectroscopic studies. Indole is the core structure of many substances found in the human body (tryptophan, serotonin) and the indole alkaloids have highly differentiated pharmacological properties such as analgesic, anti-fever and anti-inflammatory. The Ox's structural results given in the Cambridge Structural Database revealed the existence of only one crystalline form of Ox, referred to the α-form. However, we have experimentally noticed the existence of two polymorphic forms during the crystallization of Ox. Furthermore, the significant spectral differences that we have observed in the solid state infrared spectra of these two forms additionally confirm the existence of the polymorphism phenomenon. Of the four polymorphic forms of Ox, two of them - α - and β-forms - were of particular interest. In the crystalline lattices of both polymorphs, we observed a similar pattern of molecular arrangements giving rise to the supramolecular synthon according to the terminology of Etter. Moreover, hydrogen bonds in the dimer of the α-form are found to be non-equivalent (non-centrosymmetric dimers), having a length of 2797 Å and 2979 Å, respectively. Comparatively, in the most densely packed crystalline structure of Ox, the β-form, the dimer is formed by a pair of almost identical intermolecular hydrogen bonds and consequently the crystals of β-form exhibited spectral properties typical to centrosymmetric hydrogen bond dimers. In addition, the spectroscopic studies that we have conducted to polymorphic forms of Ox, isotopically diluted with deuterium, show the dramatic influence of isotopic substitution in the hydrogen bridge on the infrared spectra of hydrogen bonding. Thus, the main goal of this work is the proposition of a theoretical approach that can describe the main features of the crystalline infrared spectra of the Ox polymorphs. The proposed approach is based on the phenomenon of the exciton coupling results directly from intermolecular interactions in the vibrationally excited state which leads to the delocalization of the excitation over the molecules in the lattice and to the Davydov splitting effect in the crystalline spectra.

Substrate Binding Stiffens Aspartate Aminotransferase by Altering the Enzyme Picosecond Vibrational Dynamics

Protein dynamics on various time scales from femtoseconds to milliseconds impacts biological function by driving proteins to conformations conducive to ligand binding and creating functional states in enzyme catalysis. Neutron vibrational spectroscopy carried out by measuring inelastic neutron scattering from protein molecules in combination with molecular simulations has the unique ability of detecting and visualizing changes in the picosecond protein vibrational dynamics due to ligand binding. Here we present neutron vibrational spectra of a homodimeric pyridoxal 5'-phosphate-dependent enzyme, aspartate aminotransferase, obtained from the open internal aldimine and closed external aldimine conformational states. We observe that in the external aldimine state the protein structure stiffens relative to the internal aldimine state, indicating rigidified vibrational dynamics on the picosecond time scale in the low-frequency regime of 5-50 cm^-1. Our molecular dynamics simulations indicate substantial changes in the picosecond dynamics of the enzyme secondary structure elements upon substrate binding, with the largest contributions from just two helices and the β-sheet.