A Biophysical Perspective on Enzyme Catalysis

Polyethylene Glycol Impacts Conformation and Dynamics of Escherichia coli Prolyl-tRNA Synthetase Via Crowding and Confinement Effects

Polyethylene glycol (PEG) is a flexible, nontoxic polymer commonly used in biological and medical research, and it is generally regarded as biologically inert. PEG molecules of variable sizes are also used as crowding agents to mimic intracellular environments. A recent study with PEG crowders revealed decreased catalytic activity of Escherichia coli prolyl-tRNA synthetase (Ec ProRS), where the smaller molecular weight PEGs had the maximum impact. The molecular mechanism of the crowding effects of PEGs is not clearly understood. PEG may impact protein conformation and dynamics, thus its function. In the present study, the effects of PEG molecules of various molecular weights and concentrations on the conformation and dynamics of Ec ProRS were investigated using a combined experimental and computational approach including intrinsic tryptophan fluorescence spectroscopy, atomic force microscopy, and atomistic molecular dynamic simulations. Results of the present study suggest that lower molecular weight PEGs in the dilute regime have modest effects on the conformational dynamics of Ec ProRS but impact the catalytic function primarily via the excluded volume effect; they form large clusters blocking the active site pocket. In contrast, the larger molecular weight PEGs in dilute to semidilute regimes have a significant impact on the protein's conformational dynamics; they wrap on the protein surface through noncovalent interactions. Thus, lower-molecular-weight PEG molecules impact protein dynamics and function via crowding effects, whereas larger PEGs induce confinement effects. These results have implications for the development of inhibitors for protein targets in a crowded cellular environment.

Four Parallel Pathways in T4 Ligase-Catalyzed Repair of Nicked DNA with Diverse Bending Angles

The structural diversity of biological macromolecules in different environments contributes complexity to enzymological processes vital for cellular functions. Fluorescence resonance energy transfer and electron microscopy are used to investigate the enzymatic reaction of T4 DNA ligase catalyzing the ligation of nicked DNA. The data show that both the ligase-AMP complex and the ligase-AMP-DNA complex can have four conformations. This finding suggests the parallel occurrence of four ligation reaction pathways, each characterized by specific conformations of the ligase-AMP complex that persist in the ligase-AMP-DNA complex. Notably, these complexes have DNA bending angles of ≈0°, 20°, 60°, or 100°. The mechanism of parallel reactions challenges the conventional notion of simple sequential reaction steps occurring among multiple conformations. The results provide insights into the dynamic conformational changes and the versatile attributes of T4 DNA ligase and suggest that the parallel multiple reaction pathways may correspond to diverse T4 DNA ligase functions. This mechanism may potentially have evolved as an adaptive strategy across evolutionary history to navigate complex environments.

Deeper Insights into Mixed Crowding through Enzyme Activity, Dynamics, and Crowder Diffusion

Given the fact that the cellular interior is crowded by many different kinds of macromolecules, it is important that in vitro studies be carried out in the presence of mixed crowder systems. In this regard, we have used binary crowders formed by the combination of some of the commonly used crowding agents, namely, Ficoll 70, Dextran 70, Dextran 40, and PEG 8000 (PEG 8), to study how these affect enzyme activity, dynamics, and crowder diffusion. The enzyme chosen is AK3L1, an isoform of adenylate kinase. To investigate its dynamics, we have carried out three single point mutations (A74C, A132C, and A209C) with the cysteine residues being labeled with a coumarin-based solvatochromic probe [CPM: (7-diethylamino-3-(4-maleimido-phenyl)-4-methylcoumarin)]. Both enzyme activity and dynamics decreased in the binary mixtures as compared with the sum of the individual crowders, suggesting a reduction in excluded volume (in the mixture). To gain deeper insights into the binary mixtures, fluorescence correlation spectroscopy studies were carried out using fluorescein isothiocyanate-labeled Dextran 70 and tetramethylrhodamine-labeled AK3L1 as the diffusion probes. Diffusion in binary mixtures was observed to be much more constrained (relative to the sum of the individual crowders) for the labeled enzyme as compared to the labeled crowder showing different environments being faced by the two species. This was further confirmed during imaging of the phase-separated droplets formed in the binary mixtures having PEG as one of the crowding agents. The interior of these droplets was found to be rich in crowders and densely packed, as shown by confocal and digital holographic microscopy images, with the enzymes predominantly residing outside these droplets, that is, in the relatively less crowded regions. Taken together, our data provide important insights into various aspects of the simplest form of mixed crowding, that is, composed of just two components, and also hint at the enhanced complexity that the cellular interior presents toward having a detailed and comprehensive understanding of the same.

An active site mutation induces oxygen reactivity in D-arginine dehydrogenase: A case of superoxide diverting protons

Enzymes are potent catalysts that increase biochemical reaction rates by several orders of magnitude. Flavoproteins are a class of enzymes whose classification relies on their ability to react with molecular oxygen (O2) during catalysis using ionizable active site residues. Pseudomonas aeruginosa D-arginine dehydrogenase (PaDADH) is a flavoprotein that oxidizes D-arginine for P. aeruginosa survival and biofilm formation. The crystal structure of PaDADH reveals the interaction of the glutamate 246 (E246) side chain with the substrate and at least three other active site residues, establishing a hydrogen bond network in the active site. Additionally, E246 likely ionizes to facilitate substrate binding during PaDADH catalysis. This study aimed to investigate how replacing the E246 residue with leucine affects PaDADH catalysis and its ability to react with O2 using steady-state kinetics coupled with pH profile studies. The data reveal a gain of O2 reactivity in the E246L variant, resulting in a reduced flavin semiquinone species and superoxide (O2•-) during substrate oxidation. The O2•- reacts with active site protons, resulting in an observed nonstoichiometric slope of 1.5 in the enzyme's log (kcat/Km) pH profile with D-arginine. Adding superoxide dismutase results in an observed correction of the slope to 1.0. This study demonstrates how O2•- can alter the slopes of limbs in the pH profiles of flavin-dependent enzymes and serves as a model for correcting nonstoichiometric slopes in elucidating reaction mechanisms of flavoproteins.

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.

Kinetic View of Enzyme Catalysis from Enhanced Sampling QM/MM Simulations

The rate constants of enzyme-catalyzed reactions (kcat) are often approximated from the barrier height of the reactive step. We introduce an enhanced sampling QM/MM approach that directly calculates the kinetics of enzymatic reactions, without introducing the transition-state theory assumptions, and takes into account the dynamical equilibrium between the reactive and non-reactive conformations of the enzyme/substrate complex. Our computed kcat values are in order-of-magnitude agreement with the experimental data for two representative enzymatic reactions.

Structural enzymology studies with the substrate 3S-hydroxybutanoyl-CoA: bifunctional MFE1 is a less efficient dehydrogenase than monofunctional HAD

Multifunctional enzyme, type-1 (MFE1) catalyzes the second and third step of the β-oxidation cycle, being, respectively, the 2E-enoyl-CoA hydratase (ECH) reaction (N-terminal part, crotonase fold) and the NAD+-dependent, 3S-hydroxyacyl-CoA dehydrogenase (HAD) reaction (C-terminal part, HAD fold). Structural enzymological properties of rat MFE1 (RnMFE1) as well as of two of its variants, namely the E123A variant (a glutamate of the ECH active site is mutated into alanine) and the BCDE variant (without domain A of the ECH part), were studied, using as substrate 3S-hydroxybutanoyl-CoA. Protein crystallographic binding studies show the hydrogen bond interactions of 3S-hydroxybutanoyl-CoA as well as of its 3-keto, oxidized form, acetoacetyl-CoA, with the catalytic glutamates in the ECH active site. Pre-steady state binding experiments with NAD+ and NADH show that the kon and koff rate constants of the HAD active site of monomeric RnMFE1 and the homologous human, dimeric 3S-hydroxyacyl-CoA dehydrogenase (HsHAD) for NAD+ and NADH are very similar, being the same as those observed for the E123A and BCDE variants. However, steady state and pre-steady state kinetic data concerning the HAD-catalyzed dehydrogenation reaction of the substrate 3S-hydroxybutanoyl-CoA show that, respectively, the kcat and kchem rate constants for conversion into acetoacetyl-CoA by RnMFE1 (and its two variants) are about 10 fold lower as when catalyzed by HsHAD. The dynamical properties of dehydrogenases are known to be important for their catalytic efficiency, and it is discussed that the greater complexity of the RnMFE1 fold correlates with the observation that RnMFE1 is a slower dehydrogenase than HsHAD.

Exploration of biochemical reactivity with a QM/MM growing string method

In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.

Theoretical Improvements in Enzyme Efficiency Associated with Noisy Rate Constants and Increased Dissipation

Previous studies have revealed the extraordinarily large catalytic efficiency of some enzymes. High catalytic proficiency is an essential accomplishment of biological evolution. Natural selection led to the increased turnover number, kcat, and enzyme efficiency, kcat/KM, of uni-uni enzymes, which convert a single substrate into a single product. We added or multiplied random noise with chosen rate constants to explore the correlation between dissipation and catalytic efficiency for ten enzymes: beta-galactosidase, glucose isomerase, β-lactamases from three bacterial strains, ketosteroid isomerase, triosephosphate isomerase, and carbonic anhydrase I, II, and T200H. Our results highlight the role of biological evolution in accelerating thermodynamic evolution. The catalytic performance of these enzymes is proportional to overall entropy production-the main parameter from irreversible thermodynamics. That parameter is also proportional to the evolutionary distance of β-lactamases PC1, RTEM, and Lac-1 when natural or artificial evolution produces the optimal or maximal possible catalytic efficiency. De novo enzyme design and attempts to speed up the rate-limiting catalytic steps may profit from the described connection between kinetics and thermodynamics.

Leveraging liquid-liquid phase separation and volume modulation to regulate the enzymatic activity of formate dehydrogenase

Engineering of reaction media is an exciting alternative for modulating kinetic properties of biocatalytic reactions. We addressed the combined effect of an aqueous two-phase system (ATPS) and high hydrostatic pressure on the kinetics of the Candida boidinii formate dehydrogenase-catalyzed oxidation of formate to CO2. Pressurization was found to lead to an increase of the binding affinity (decrease of KM, respectively) and a decrease of the turnover number, kcat. The experimental approach was supported using thermodynamic modeling with the electrolyte Perturbed-Chain Statistical Associating Fluid Theory (ePC-SAFT) equation of state to predict the liquid-liquid phase separation and the molecular crowding effect of the ATPS on the kinetic properties. The ePC-SAFT was able to quantitatively predict the KM-values of the substrate in both phases at 1 bar as well as up to a pressure of 1000 bar. The framework presented enables significant advances in bioprocess engineering, including the design of processes with significantly fewer experiments and trial-and-error approaches.

Targeted Mutation of a Non-catalytic Gating Residue Increases the Rate of Pseudomonas aeruginosa d-Arginine Dehydrogenase Catalytic Turnover

Commercial food and l-amino acid industries rely on bioengineered d-amino acid oxidizing enzymes to detect and remove d-amino acid contaminants. However, the bioengineering of enzymes to generate faster biological catalysts has proven difficult as a result of the failure to target specific kinetic steps that limit enzyme turnover, kcat, and the poor understanding of loop dynamics critical for catalysis. Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) oxidizes most d-amino acids and is a good candidate for application in the l-amino acid and food industries. The side chain of the loop L2 E246 residue located at the entrance of the PaDADH active site pocket potentially favors the closed active site conformation and secures the substrate upon binding. This study used site-directed mutagenesis, steady-state, and rapid reaction kinetics to generate the glutamine, glycine, and leucine variants and investigate whether increasing the rate of product release could translate to an increased enzyme turnover rate. Upon E246 mutation to glycine, there was an increased rate of d-arginine turnover kcat from 122 to 500 s-1. Likewise, the kcat values increased 2-fold for the glutamine or leucine variants. Thus, we have engineered a faster biocatalyst for industrial applications by selectively increasing the rate of the PaDADH product release.

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.

How Flexibility Can Enhance Catalysis

Conformational changes are observed in many enzymes, but their role in catalysis is highly controversial. Here we present a theoretical model that illustrates how rigid catalysts can be fundamentally limited and how a conformational change induced by substrate binding can overcome this limitation, ultimately enabling barrier-free catalysis. The model is deliberately minimal, but the principle it illustrates is general and consistent with unique features of proteins as well as with previous informal proposals to explain the superiority of enzymes over other classes of catalysts. Implementing the discriminative switch suggested by the model could help overcome limitations currently encountered in the design of artificial catalysts.

New Insights into the Cooperativity and Dynamics of Dimeric Enzymes

A survey of protein databases indicates that the majority of enzymes exist in oligomeric forms, with about half of those found in the UniProt database being homodimeric. Understanding why many enzymes are in their dimeric form is imperative. Recent developments in experimental and computational techniques have allowed for a deeper comprehension of the cooperative interactions between the subunits of dimeric enzymes. This review aims to succinctly summarize these recent advancements by providing an overview of experimental and theoretical methods, as well as an understanding of cooperativity in substrate binding and the molecular mechanisms of cooperative catalysis within homodimeric enzymes. Focus is set upon the beneficial effects of dimerization and cooperative catalysis. These advancements not only provide essential case studies and theoretical support for comprehending dimeric enzyme catalysis but also serve as a foundation for designing highly efficient catalysts, such as dimeric organic catalysts. Moreover, these developments have significant implications for drug design, as exemplified by Paxlovid, which was designed for the homodimeric main protease of SARS-CoV-2.

Reactivity Tracking of an Enzyme Progress Coordinate

The reactivity of individual solvent-coupled protein configurations is used to track and resolve the progress coordinate for the core reaction sequence of substrate radical rearrangement and hydrogen atom transfer in the ethanolamine ammonia-lyase (EAL) enzyme from Salmonella enterica. The first-order decay of the substrate radical intermediate is the monitored reaction. Heterogeneous confinement from sucrose hydrates in the mesophase solvent surrounding the cryotrapped protein introduces distributed kinetics in the non-native decay of the substrate radical pair capture substate, which arise from an ensemble of configurational microstates. Reaction rates increase by >103-fold across the distribution to approach that for the native enabled substate for radical rearrangement, which reacts with monotonic kinetics. The native progress coordinate thus involves a collapse of the configuration space to generate optimized reactivity. Reactivity tracking reveals fundamental features of solvent-protein-reaction configurational coupling and leads to a model that refines the ensemble paradigm of enzyme catalysis for strongly adiabatic chemical steps.

Evidence for a short-lived resonance state in enzyme catalysis via rate-equation convolution

At the cellular level, all biological function relies on enzymes to provide catalytic acceleration of essential biochemical processes driving cellular metabolism. The enzyme is presumed to lower the activation energy barrier separating reactants from products, but the precise mechanism remains unresolved. Here we examine the temperature dependence of the enzyme-catalyzed dissociation of p-nitrophenyl-α-D-glucopyranoside (pNPG), a chromogenic analog for maltose, isomaltose, and sucrose disaccharide sugars, into p-nitrophenol (pNP) and glucose (monosaccharide). The enzymes of interest are the wild type and mutant forms of glucosidase MalL produced by the probiotic bacterium Bacillus subtilis. The per-enzyme production rates k(T) for the pNPG→ glucose reaction all show a characteristic temperature profile with an Arrhenius-like (approximately exponential) slow acceleration at low temperatures, rising through a point of inflexion to reach a maximum, then turning over to decline steeply towards zero production at high temperatures. This asymmetric profile is found to be well fitted by convolving an exponential growth function f(T) with a Gaussian temperature distribution g(T) to produce an exponentially modified Gaussian function h(T). To give a physical interpretation of the convolution components, we make the temperature mapping Θ≡T_{ref}-T where T_{ref} marks the temperature at which a given mutant becomes fully denatured (unfolded) and therefore inactive, then convert the convolution components to probability density functions which obey the convolution theorem of statistics. Working in Θ space, we identify f(Θ) as the density function for an Arrhenius-like transition from ground-state A to metastable-state B, and g(Θ) as the Gaussian distribution of offset-temperature fluctuations for the metastable state. By mapping the standard thermodynamic relations for temperature and energy fluctuations to the enzyme frame of reference, we are able to derive an expression for the lifetime for the metastable B state. For the 15 enzyme experiments, we obtain a mean value 〈Δt〉≳(29.0±1.3)×10^{-15}s, in remarkably good agreement with the ∼30-fs estimate for the period of glycosidic bond oscillations extracted from published infrared spectroscopy. We suggest that the metastable B state provides a low-energy target that has the effect of lowering the activation energy barrier by presenting an alternative axis for the reaction coordinate.

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.

"Stereoelectronic Deprotection of Nitrogen": Recovering Nucleophilicity with a Conformational Change

Ureas are often thought of as "double amides" due to the obvious structural similarity of these functional groups. The main structural feature of an amide is its planarity, which is responsible for the conjugation between the nitrogen atom and carbonyl moiety and the decrease of amide nucleophilicity. Consequently, since amides are poor nucleophiles, ureas are often thought of as poor nucleophiles as well. Herein, we demonstrate that ureas can be distinctly different from amides. These differences can be amplified by rotation around one of the ureas' C-N bonds, which switches off the amide resonance and recovers the nucleophilicity of one of the nitrogen atoms. This conformational change can be further facilitated by the judicious introduction of steric bulk to disfavor the planar conformation. This change in reactivity is an example of "stereoelectronic deprotection," a concept when the desired reactivity of a functional group is produced by a conformational change rather than a chemical modification. This concept may be used complementarily to the traditional protecting groups. We also demonstrate both the viability and the utility of this concept by the synthesis of unusual 2-oxoimidazolium salts possessing quaternary nitrogen atoms at the urea moiety.

MAHOMES II: A webserver for predicting if a metal binding site is enzymatic

Recent advances have enabled high-quality computationally generated structures for proteins with no solved crystal structures. However, protein function data remains largely limited to experimental methods and homology mapping. Since structure determines function, it is natural that methods capable of using computationally generated structures for functional annotations need to be advanced. Our laboratory recently developed a method to distinguish between metalloenzyme and nonenzyme sites. Here we report improvements to this method by upgrading our physicochemical features to alleviate the need for structures with sub-angstrom precision and using machine learning to reduce training data labeling error. Our improved classifier identifies protein bound metal sites as enzymatic or nonenzymatic with 94% precision and 92% recall. We demonstrate that both adjustments increased predictive performance and reliability on sites with sub-angstrom variations. We constructed a set of predicted metalloprotein structures with no solved crystal structures and no detectable homology to our training data. Our model had an accuracy of 90%-97.5% depending on the quality of the predicted structures included in our test. Finally, we found the physicochemical trends that drove this model's successful performance were local protein density, second shell ionizable residue burial, and the pocket's accessibility to the site. We anticipate that our model's ability to correctly identify catalytic metal sites could enable identification of new enzymatic mechanisms and improve de novo metalloenzyme design success rates.

Engineering the Active Site Lid Dynamics to Improve the Catalytic Efficiency of Yeast Cytosine Deaminase

Conformational dynamics is important for enzyme catalysis. However, engineering dynamics to achieve a higher catalytic efficiency is still challenging. In this work, we develop a new strategy to improve the activity of yeast cytosine deaminase (yCD) by engineering its conformational dynamics. Specifically, we increase the dynamics of the yCD C-terminal helix, an active site lid that controls the product release. The C-terminal is extended by a dynamical single α-helix (SAH), which improves the product release rate by up to ~8-fold, and the overall catalytic rate kcat by up to ~2-fold. It is also shown that the kcat increase is due to the favorable activation entropy change. The NMR H/D exchange data indicate that the conformational dynamics of the transition state analog complex increases as the helix is extended, elucidating the origin of the enhanced catalytic entropy. This study highlights a novel dynamics engineering strategy that can accelerate the overall catalysis through the entropy-driven mechanism.

MAHOMES II: A webserver for predicting if a metal binding site is enzymatic

Recent advances have enabled high-quality computationally generated structures for proteins with no solved crystal structures. However, protein function data remains largely limited to experimental methods and homology mapping. Since structure determines function, it is natural that methods capable of using computationally generated structures for functional annotations need to be advanced. Our laboratory recently developed a method to distinguish between metalloenzyme and non-enzyme sites. Here we report improvements to this method by upgrading our physicochemical features to alleviate the need for structures with sub-angstrom precision and using machine learning to reduce training data labeling error. Our improved classifier identifies protein bound metal sites as enzymatic or non-enzymatic with 94% precision and 92% recall. We demonstrate that both adjustments increased predictive performance and reliability on sites with sub-angstrom variations. We constructed a set of predicted metalloprotein structures with no solved crystal structures and no detectable homology to our training data. Our model had an accuracy of 90 - 97.5% depending on the quality of the predicted structures included in our test. Finally, we found the physicochemical trends that drove this model's successful performance were local protein density, second shell ionizable residue burial, and the pocket's accessibility to the site. We anticipate that our model's ability to correctly identify catalytic metal sites could enable identification of new enzymatic mechanisms and improve de novo metalloenzyme design success rates.

Significance statement

Identification of enzyme active sites on proteins with unsolved crystallographic structures can accelerate discovery of novel biochemical reactions, which can impact healthcare, industrial processes, and environmental remediation. Our lab has developed an ML tool for predicting sites on computationally generated protein structures as enzymatic and non-enzymatic. We have made our tool available on a webserver, allowing the scientific community to rapidly search previously unknown protein function space.

Non-active Site Residue in Loop L4 Alters Substrate Capture and Product Release in d-Arginine Dehydrogenase

Numerous studies demonstrate that enzymes undergo multiple conformational changes during catalysis. The malleability of enzymes forms the basis for allosteric regulation: residues located far from the active site can exert long-range dynamical effects on the active site residues to modulate catalysis. The structure of Pseudomonas aeruginosa d-arginine dehydrogenase (PaDADH) shows four loops (L1, L2, L3, and L4) that span the substrate and the FAD-binding domains. Loop L4 comprises residues 329-336, spanning over the flavin cofactor. The I335 residue on loop L4 is ∼10 Å away from the active site and ∼3.8 Å from N(1)-C(2)═O atoms of the flavin. In this study, we used molecular dynamics and biochemical techniques to investigate the effect of the mutation of I335 to histidine on the catalytic function of PaDADH. Molecular dynamics showed that the conformational dynamics of PaDADH are shifted to a more closed conformation in the I335H variant. In agreement with an enzyme that samples more in a closed conformation, the kinetic data of the I335H variant showed a 40-fold decrease in the rate constant of substrate association (k₁), a 340-fold reduction in the rate constant of substrate dissociation from the enzyme-substrate complex (k₂), and a 24-fold decrease in the rate constant of product release (k₅), compared to that of the wild-type. Surprisingly, the kinetic data are consistent with the mutation having a negligible effect on the reactivity of the flavin. Altogether, the data indicate that the residue at position 335 has a long-range dynamical effect on the catalytic function in PaDADH.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

The origin of genetic and metabolic systems: Evolutionary structuralinsights

DNA is derived from reverse transcription and its origin is related to reverse transcriptase, DNA polymerase and integrase. The gene structure originated from the evolution of the first RNA polymerase. Thus, an explanation of the origin of the genetic system must also explain the evolution of these enzymes. This paper proposes a polymer structure model, termed the stable complex evolution model, which explains the evolution of enzymes and functional molecules. Enzymes evolved their functions by forming locally tightly packed complexes with specific substrates. A metabolic reaction can therefore be considered to be the result of adaptive evolution in this way when a certain essential molecule is lacking in a cell. The evolution of the primitive genetic and metabolic systems was thus coordinated and synchronized. According to the stable complex model, almost all functional molecules establish binding affinity and specific recognition through complementary interactions, and functional molecules therefore have the nature of being auto-reactive. This is thermodynamically favorable and leads to functional duplication and self-organization. Therefore, it can be speculated that biological systems have a certain tendency to maintain functional stability or are influenced by an inherent selective power. The evolution of dormant bacteria may support this hypothesis, and inherent selectivity can be unified with natural selection at the molecular level.

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.

Endohedral Functionalization for Structural Transformation of Polyoxovanadate-Based Metal-Organic Cube

Functionalized internal modifications of metal-organic polyhedra (MOPs) can endow properties and functions different from the original ones. Until now, there have been only a few examples of endohedral modifications of polyoxovanadate-based MOPs. Herein, an efficient coordination-driven strategy was chosen for the inner modification of two metal-organic cubes (MOCs) with different sizes, VMOC-1 and VMOC-4, constructed from polyoxovanadate clusters [V₆O₆(OCH₃)₉(SO₄)(CO₂)₃]^2- SBU and tetradentate ligands. Pyridinophosphonic acid with potential coordination capability was introduced to replace the sulfate of the hexavanadate cluster and graft the pyridine functional group inside the cage. The introduction of pyridylphosphate in the VMOC-4 system gave a cubic cage with a pyridyl endo-modified isomer. Interestingly, the smaller cubic cage VMOC-1 was induced to undergo structural transformation to obtain VMOC-py-1. The organic dyes adsorption of VMOC-py-1 and VMOC-1 showed that the endomodified structure could adsorb larger and more dyes, compared to the original cube.

Modeling Ensembles of Enzyme Reaction Pathways with Hi-MSM Reveals the Importance of Accounting for Pathway Diversity

Conventional quantum mechanical-molecular mechanics (QM/MM) simulation approaches for modeling enzyme reactions often assume that there is one dominant reaction pathway and that this pathway can be sampled starting from an X-ray structure of the enzyme. These assumptions reduce computational cost; however, their validity has not been extensively tested. This is due in part to the lack of a rigorous formalism for integrating disparate pathway information from dynamical QM/MM calculations. Here, we present a way to model ensembles of reaction pathways efficiently using a divide-and-conquer strategy through Hierarchical Markov State Modeling (Hi-MSM). This approach allows information on multiple, distinct pathways to be incorporated into a chemical kinetic model, and it allows us to test these two assumptions. Applying Hi-MSM to the reaction carried out by dihydrofolate reductase (DHFR) we find (i) there are multiple, distinct pathways significantly contributing to the overall flux of the reaction that the conventional approach does not identify and (ii) that the conventional approach does not identify the dominant reaction pathway. Thus, both assumptions underpinning the conventional approach are violated. Since DHFR is a relatively small enzyme, and configuration space scales exponentially with protein size, accounting for multiple reaction pathways is likely to be necessary for most enzymes.

Synchronizing Two Distinct Nano-Circular Sliding Motions in Six-Component Machinery for Double Catalysis

The heteroleptic multi-component double slider-on-deck system DS3 exhibits tight coupling of motional speed of two distinct nano-circular sliders (k₂₉₈ =77 and 41 kHz) despite a 2.2 nm separation. In comparison, the single sliders in DS1 and DS2 move at vastly different speed (k₂₉₈ =1.1 vs. 350 kHz). Synchronization of the motions in DS3 remains even when one slows the movement of the faster slider using small molecular brake pads. In contrast to the individual DS1 and DS2 systems, DS3 is a powerful catalyst for a two-step reaction by using the motion of both sliders to drive two catalytic processes.

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Insight into Substrate Recognition by Urea-Based Helical Foldamer Catalysts Using a DFT Global Optimization Approach

Peptides and foldamers have recently gained increasing attention as chiral catalysts to achieve challenging (asymmetric) transformations. We previously reported that short helically folded aliphatic oligoureas in combination with achiral Brønsted bases are effective H-bonding catalysts for C-C bond-forming reactions─i.e., the conjugate addition of 1,3-dicarbonyl pronucleophiles to nitroalkenes─with high reactivity and selectivity and at remarkably low chiral catalyst/substrate molar ratios. This theoretical investigation at the density functional theory level of theory, aims to both analyze how the substrates of the reaction interact with the foldamer catalyst and rationalize a chain-length dependence effect on the catalytic properties. We confirm that the first two ureas are the only H-bond donors available to interact with external molecules. Moreover, each urea site interacts with one of the two reactants allowing a short distance between the two reacting carbons, thus facilitating the conjugated addition. Additionally, it was observed that the molecular recognition and catalyst-substrate interactions are mainly governed by electrostatic interactions but not orbital interactions (see from NBO if this is finally true). On these grounds, an electrostatic potential (ESP) analysis showed an important internal charge separation in the catalyst, the positive ESP region being concentrated around the first two ureas, with its area extending as the number of residues increases.

Functional roles of enzyme dynamics in accelerating active site chemistry: Emerging techniques and changing concepts

With the growing acceptance of the contribution of protein conformational ensembles to enzyme catalysis and regulation, research in the field of protein dynamics has shifted toward an understanding of the atomistic properties of protein dynamical networks and the mechanisms and time scales that control such behavior. A full description of an enzymatic reaction coordinate is expected to extend beyond the active site and include site-specific networks that communicate with the protein/water interface. Advances in experimental tools for the spatial resolution of thermal activation pathways are being complemented by biophysical methods for visualizing dynamics in real time. An emerging multidimensional model integrates the impacts of bound substrate/effector on the distribution of protein substates that are in rapid equilibration near room temperature with reaction-specific protein embedded heat transfer conduits.

Kinetic solvent viscosity effects uncover an internal isomerization of the enzyme-substrate complex in Pseudomonas aeruginosa PAO1 NADH:Quinone oxidoreductase

NAD(P)H:quinone oxidoreductases (NQOs) play an essential protective role as antioxidants in the detoxification of quinones in both Prokaryotes and Eukaryotes. NQO from Pseudomonas aeruginosa PAO1 uses FMN to catalyze the two-electron reduction of various quinones with NADH. In this study, steady-state kinetics, kinetic solvent viscosity effects, and rapid reaction kinetics were used to determine which kinetic steps control the overall turnover of the enzyme with benzoquinone or juglone. The rate constant for flavin reduction (k_red) at pH 6.0 was 12.9 ± 0.3 s^-1, and the K_d for NADH was at least an order of magnitude lower than 90 μM. With benzoquinone, the k_cat value was 11.7 ± 0.3 s^-1, consistent with flavin reduction being almost entirely rate-limiting for overall turnover. With juglone, a k_cat value of 10.0 ± 0.5 s^-1 was recorded. The normalized plot of the relative solvent viscosity effects on the k_cat values established that hydride transfer from NADH to the FMN and quinol product release, with a calculated rate constant (k_P-rel) of 52 s^-1, are partially rate-limiting for the overall turnover of NQO. Kinetic solvent viscosity effects with glucose or sucrose revealed a hyperbolic dependence on the k_cat and k_cat/K_m values with benzoquinone or juglone, respectively, consistent with the presence of a solvent-sensitive internal isomerization of the enzyme-substrate complex (ES). The data demonstrate opposing effects of benzoquinone and juglone on the equilibrium of the NQO ES isomerization with glucose or sucrose. Thus, our study demonstrates how quinol substrate properties alter the equilibrium of NQO ES isomerization.

Resolution and characterization of contributions of select protein and coupled solvent configurational fluctuations to radical rearrangement catalysis in coenzyme B12-dependent ethanolamine ammonia-lyase

Coenzyme B₁₂ (adenosylcobalamin) -dependent ethanolamine ammonia-lyase (EAL) is the signature enzyme in ethanolamine utilization metabolism associated with microbiome homeostasis and disease conditions in the human gut. The enzyme conducts a complex choreography of bond-making/bond-breaking steps that rearrange substrate to products through a radical mechanism, with themes common to other coenzyme B₁₂-dependent and radical enzymes. The methods presented are targeted to test the hypothesis that particular, select protein and coupled solvent configurational fluctuations contribute to enzyme function. The general approach is to correlate enzyme function with an introduced perturbation that alters the properties (for example, degree of concertedness, or collectiveness) of protein and coupled solvent dynamics. Methods for sample preparation and low-temperature kinetic measurements by using temperature-step reaction initiation and time-resolved, full-spectrum electron paramagnetic resonance spectroscopy are detailed. A framework for interpretation of results obtained in ensemble systems under conditions of statistical equilibrium within the reacting, globally unstable state is presented. The temperature-dependence of the first-order rate constants for decay of the cryotrapped paramagnetic substrate radical state in EAL, through the chemical step of radical rearrangement, displays a piecewise-continuous Arrhenius dependence from 203 to 295K, punctuated by a kinetic bifurcation over 219-220K. The results reveal the obligatory contribution of a class of select collective protein and coupled solvent fluctuations to the interconversion of two resolved, sequential configurational substates, on the decay time scale. The select class of collective fluctuations also contributes to the chemical step. The methods and analysis are generally applicable to other coenzyme B₁₂-dependent and related radical enzymes.

Hydride Transfer Mechanism of Enzymatic Sugar Nucleotide C2 Epimerization Probed with a Loose-Fit CDP-Glucose Substrate

Transient oxidation–reduction through hydride transfer with tightly bound NAD coenzyme is used by a large class of sugar nucleotide epimerases to promote configurational inversion of carbon stereocenters in carbohydrate substrates. A requirement for the epimerases to coordinate hydride abstraction and re-addition with substrate rotation in the binding pocket poses a challenge for dynamical protein conformational selection linked to enzyme catalysis. Here, we studied the thermophilic C2 epimerase from Thermodesulfatator atlanticus (TaCPa2E) in combination with a slow CDP-glucose substrate (k_cat ≈ 1.0 min^–1; 60 °C) to explore the sensitivity of the enzymatic hydride transfer toward environmental fluctuations affected by temperature (20–80 °C). We determined noncompetitive primary kinetic isotope effects (KIE) due to ²H at the glucose C2 and showed that a normal KIE on the k_cat (^Dk_cat) reflects isotope sensitivity of the hydrogen abstraction to enzyme-NAD⁺ in a rate-limiting transient oxidation. The ^Dk_cat peaked at 40 °C was 6.1 and decreased to 2.1 at low (20 °C) and 3.3 at high temperature (80 °C). The temperature profiles for k_cat with the ¹H and ²H substrate showed a decrease in the rate below a dynamically important breakpoint (∼40 °C), suggesting an equilibrium shift to an impaired conformational landscape relevant for catalysis in the low-temperature region. Full Marcus-like model fits of the rate and KIE profiles provided evidence for a high-temperature reaction via low-frequency conformational sampling associated with a broad distribution of hydride donor–acceptor distances (long-distance population centered at 3.31 ± 0.02 Å), only poorly suitable for quantum mechanical tunneling. Collectively, dynamical characteristics of TaCPa2E-catalyzed hydride transfer during transient oxidation of CDP-glucose reveal important analogies to mechanistically simpler enzymes such as alcohol dehydrogenase and dihydrofolate reductase. A loose-fit substrate (in TaCPa2E) resembles structural variants of these enzymes by extensive dynamical sampling to balance conformational flexibility and catalytic efficiency.

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.

Electrostatic fingerprints of catalytically active amino acids in enzymes

The computed electrostatic and proton transfer properties are studied for 20 enzymes that represent all six major enzyme commission classes and a variety of different folds. The properties of aspartate, glutamate, and lysine residues that have been previously experimentally determined to be catalytically active are reported. The catalytic aspartate and glutamate residues studied here are strongly coupled to at least one other aspartate or glutamate residue and often to multiple other carboxylate residues with intrinsic pK_a differences less than 1 pH unit. Sometimes these catalytic acidic residues are also coupled to a histidine residue, such that the intrinsic pK_a of the acidic residue is higher than that of the histidine. All catalytic lysine residues studied here are strongly coupled to tyrosine or cysteine residues, wherein the intrinsic pK_a of the anion-forming residue is higher than that of the lysine. Some catalytic lysines are also coupled to other lysines with intrinsic pK_a differences within 1 pH unit. Some evidence of the possible types of interactions that facilitate nucleophilicity is discussed. The interactions reported here provide important clues about how side chain functional groups that are weak Brønsted acids or bases for the free amino acid in solution can achieve catalytic potency and become strong acids, bases or nucleophiles in the enzymatic environment.

Mechanisms of macromolecular reactions

During the past two decades, philosophers of biology have increasingly turned their attention to mechanisms of biological phenomena. Through analyses of mechanistic proposals advanced by biologists, the goal of these philosophers is to understand what a mechanism is and how mechanisms explain. These analyses have generally focused on mechanistic proposals for phenomenon that occur at the cellular or sub-cellular level, such as synapse firing, protein synthesis, or metabolic pathway operation. Little is said about the mechanisms of the macromolecular reactions that underpin these phenomena. These reactions comprise a diverse family of reaction types, and include protein folding, macromolecular complex formation, receptor-ligand interactions, and enzyme catalysis. In this paper, I develop an account of mechanism that focuses exclusively on macromolecular reactions. I begin by reviewing how mechanism is understood in enzymology, and how mechanistic concepts of enzymology apply to macromolecular reactions in general. We will see that the mechanism of a macromolecular reaction is most accurately described as a progression of reaction intermediates, where the evolution of intermediates, from one to the next, is characterized by an energetic coupling between chemistry and protein dynamics. I then make the case that this description necessitates a grounding in a process ontology. To describe the mechanism by which a macromolecular reaction occurs is to describe a process.

Reversible Regulating the Substrate Specificity of Enzymes in Microgels by a Phase Transition in Polymer Networks

Here, we report a distinct approach for regulating the substrate specificity of enzymes immobilized in microgels by a phase transition in polymer networks. The finding is demonstrated on glucose oxidase that is immobilized in thermoresponsive poly(N-isopropylacrylamide)-based microgels. Laser light scattering and enzymatic oxidation tests indicate that the broadened specificity appears at low temperatures, at which the gel matrix is in the relatively swollen state relative to its state at microgel synthesis temperature; upon heating to the relative higher temperatures, the gel matrix is not able to shrink further that offers a tight space in which the enzyme resides to retain high glucose specificity. It is proposed that polymer phase transition in the gel matrix mainly alter protein gates that control passage of substrates into active sites, making them open or close to a certain extent that enable reversible regulating the substrate specificity. The finding is also observed on bulk gels under a rational design, making it of potential interest in enzymatic biofuel cell applications.

Computational Enzyme Design at Zymvol

Directed evolution is the most recognized methodology for enzyme engineering. The main drawback resides in its random nature and in the limited sequence exploration; both require screening of thousands (if not millions) of variants to achieve a target function. Computer-driven approaches can limit laboratorial screening to a few hundred candidates, enabling and accelerating the development of industrial enzymes. In this book chapter, the technology adopted at Zymvol is described. An overview of the current development and future directions in the company is also provided.

A protein scaffold enables hydrogen evolution for a Ni-bisdiphosphine complex

An artificial metalloenzyme acting as a functional biomimic of hydrogenase enzymes was activated by assembly via covalent attachment of the molecular complex, [Ni(PN^glycineP)₂]^2-, within a structured protein scaffold. Electrocatalytic H₂ production was observed from pH 3.0 to 10.0 for the artificial enzyme, while no electrocatalytic activity was observed for similar [Ni(PNP)₂]²⁺ systems.

Correspondence insights into the role of genes in cell functionality. Comments on "The gene: An appraisal" by K. Baverstock

One of the most important goals of the post-genomic era is to understand the different sources of molecular information that regulate the functional and structural architecture of cells. In this regard, Prof. K. Baverstock underscores in his recent article "The gene: An appraisal" (Baverstock, 2021) that genes are not the leading elements in cellular functionality, inheritance and evolution. As a consequence, the theory of evolution based on the Neo-Darwinian synthesis, is inadequate for today's scientific evidence. Conversely, the author contends that life processes viewed on the basis of thermodynamics, complex system dynamics and self-organization provide a new framework for the foundations of Biology. I consider it necessary to comment on some essential aspects of this relevant work, and here I present a short overview of the main non-genetic sources of biomolecular order and complexity that underline the molecular dynamics and functionality of cells. These sources generate different processes of complexity, which encompasses from the most elementary levels of molecular activity to the emergence of systemic behaviors, and the information necessary to sustain them is not contained in the genome.

Quantification of Conformational Entropy Unravels Effect of Disordered Flanking Region in Coupled Folding and Binding

Intrinsic disorder (ID) constitutes a new dimension to the protein structure-function relationship. The ability to undergo conformational changes upon binding is a key property of intrinsically disordered proteins and remains challenging to study using conventional methods. A 1994 paper by R. S. Spolar and M. T. Record presented a thermodynamic approach for estimating changes in conformational entropy based on heat capacity changes, allowing quantification of residues folding upon binding. Here, we adapt the method for studies of intrinsically disordered proteins. We integrate additional data to provide a broader experimental foundation for the underlying relations and, based on >500 protein-protein complexes involving disordered proteins, reassess a key relation between polar and nonpolar surface area changes, previously determined using globular protein folding. We demonstrate the improved suitability of the adapted method to studies of the folded αα-hub domain RST from radical-induced cell death 1, whose interactome is characterized by ID. From extensive thermodynamic data, quantifying the conformational entropy changes upon binding, and comparison to the NMR structure, the adapted method improves accuracy for ID-based studies. Furthermore, we apply the method, in conjunction with NMR, to reveal hitherto undetected effects of interaction-motif context. Thus, inclusion of the disordered context of the DREB2A RST-binding motif induces structuring of the binding motif, resulting in major enthalpy-entropy compensation in the interaction interface. This study, also evaluating additional interactions, demonstrates the strength of the ID-adapted Spolar-Record thermodynamic approach for dissection of structural features of ID-based interactions, easily overlooked in traditional studies, and for translation of these into mechanistic knowledge.

Unifying Catalysis Framework to Dissect Proteasomal Degradation Paradigms

Diverging from traditional target inhibition, proteasomal protein degradation approaches have emerged as novel therapeutic modalities that embody distinct pharmacological profiles and can access previously undrugged targets. Small molecule degraders have the potential to catalytically destroy target proteins at substoichiometric concentrations, thus lowering administered doses and extending pharmacological effects. With this mechanistic premise, research efforts have advanced the development of small molecule degraders that benefit from stable and increased affinity ternary complexes. However, a holistic framework that evaluates different degradation modes from a catalytic perspective, including focusing on kinetically favored degradation mechanisms, is lacking. In this Outlook, we introduce the concept of an induced cooperativity spectrum as a unifying framework to mechanistically understand catalytic degradation profiles. This framework is bolstered by key examples of published molecular degraders extending from molecular glues to bivalent degraders. Critically, we discuss remaining challenges and future opportunities in drug discovery to rationally design and phenotypically screen for efficient degraders.

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".