Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray-derived conformational ensembles

A structural perspective on the temperature-dependent activity of enzymes

Enzymes are biomolecular catalysts whose activity varies with temperature. Unlike for small-molecule catalysts, the structural ensembles of enzymes can vary substantially with temperature, and it is in general unclear how this modulates the temperature dependence of activity. Here multi-temperature X-ray crystallography was used to record structural changes from -20°C to 40°C for a mesophilic enzyme in complex with inhibitors mimicking substrate-, intermediate-, and product-bound states, representative of major complexes underlying the kinetic constantk c a t . Both inhibitors, substrates and catalytically relevant loop motifs increasingly populate catalytically competent conformations as temperature increases. These changes occur even in temperature ranges where kinetic measurements show roughly linear Arrhenius/Eyring behavior where parameters characterizing the system are assumed to be temperature independent. Simple analysis shows that linear Arrhenius/Eyring behavior can still be observed when the underlying activation energy / enthalpy values vary with temperature, e.g., due to structural changes, and that the underlying thermodynamic parameters can be far from values derived from Arrhenius/Eyring model fits. Our results indicate a critical role for temperature-dependent atomic-resolution structural data in interpreting temperature-dependent kinetic data from enzymatic systems.

Structures of human PTP1B variants reveal allosteric sites to target for weight loss therapy

Protein Tyrosine Phosphatase 1B (PTP1B) is a negative regulator of leptin signaling whose disruption protects against diet-induced obesity in mice. We investigated whether structural characterization of human PTP1B variant proteins might reveal precise mechanisms to target for weight loss therapy. We selected 12 rare variants for functional characterization from exomes from 997 people with persistent thinness and 200,000 people from UK Biobank. Seven of 12 variants impaired PTP1B function by increasing leptin-stimulated STAT3 phosphorylation in cells. Using room-temperature X-ray crystallography, hydrogen-deuterium exchange mass spectrometry, and computational modeling, we determined that human variants modulate the 3-dimensional structure of PTP1B through distinct allosteric conduits that energetically link distal, highly ligandable structural regions to the active site. These studies inform the design of allosteric PTP1B inhibitors for the treatment of obesity.

Native dynamics and allosteric responses in PTP1B probed by high-resolution HDX-MS

Protein tyrosine phosphatase 1B (PTP1B) is a validated therapeutic target for obesity, diabetes, and certain types of cancer. In particular, allosteric inhibitors hold potential for therapeutic use, but an incomplete understanding of conformational dynamics and allostery in this protein has hindered their development. Here, we interrogate solution dynamics and allosteric responses in PTP1B using high-resolution hydrogen-deuterium exchange mass spectrometry (HDX-MS), an emerging and powerful biophysical technique. Using HDX-MS, we obtain a detailed map of backbone amide exchange that serves as a proxy for the solution dynamics of apo PTP1B, revealing several flexible loops interspersed among more constrained and rigid regions within the protein structure, as well as local regions that exchange faster than expected from their secondary structure and solvent accessibility. We demonstrate that our HDX rate data obtained in solution adds value to estimates of conformational heterogeneity derived from a pseudo-ensemble constructed from ~200 crystal structures of PTP1B. Furthermore, we report HDX-MS maps for PTP1B with active-site versus allosteric small-molecule inhibitors. These maps suggest distinct and widespread effects on protein dynamics relative to the apo form, including changes in locations distal (>35 Å) from the respective ligand binding sites. These results illuminate that allosteric inhibitors of PTP1B can induce unexpected changes in dynamics that extend beyond the previously understood allosteric network. Together, our data suggest a model of BB3 allostery in PTP1B that combines conformational restriction of active-site residues with compensatory liberation of distal residues that aid in entropic balancing. Overall, our work showcases the potential of HDX-MS for elucidating aspects of protein conformational dynamics and allosteric effects of small-molecule ligands and highlights the potential of integrating HDX-MS alongside other complementary methods, such as room-temperature X-ray crystallography, NMR spectroscopy, and molecular dynamics simulations, to guide the development of new therapeutics.

Opportunities and challenges for plastic depolymerization by biomimetic catalysis

Plastic waste has imposed significant burdens on the environment. Chemical recycling allows for repeated regeneration of plastics without deterioration in quality, but often requires harsh reaction conditions, thus being environmentally unfriendly. Enzymatic catalysis offers a promising solution for recycling under mild conditions, but it faces inherent limitations such as poor stability, high cost, and narrow substrate applicability. Biomimetic catalysis may provide a new avenue by combining high enzyme-like activity with the stability of inorganic materials. Biomimetic catalysis has demonstrated great potential in biomass conversion and has recently shown promising progress in plastic degradation. This perspective discusses biomimetic catalysis for plastic degradation from two perspectives: the imitation of the active centers and the imitation of the substrate-binding clefts. Given the chemical similarity between biomass and plastics, relevant work is also included in the discussion to draw inspiration. We conclude this perspective by highlighting the challenges and opportunities in achieving sustainable plastic recycling via a biomimetic approach.

Nanohybrid nanozyme based colourimetric immunosensor for porcine gelatin

Enzyme mimicking nanomaterials, nanozymes, have gained considerable interest in the scientific community because of their superior properties compared to natural enzymes, including their high stability at extreme conditions, cheaper availability, and ease of synthesis. Herein, we report novel colloidal gold nanoparticles - graphene nanoplatelets - chitosan (CS) with peroxidase mimicking properties used to carry out highly sensitive and selective immunoassay for porcine gelatin detection. The interaction between anti-gelatin antibody conjugated nanozyme with porcine gelatin proteins produced an ultrasensitive immunoassay response in the form of a colourimetric signal directly proportional to the porcine gelatin protein concentration. The nanozyme produced a colourimetric response in the presence of its substrate, 3,3',5,5'-tetramethylbenzidine (TMB) and hydrogen peroxide (H2O2), demonstrating its peroxidase mimicking properties. The results revealed that the nanozyme exhibited remarkable selectivity and sensitivity in the assay, detecting proteins at concentrations as low as 86.42 pg/mL. Additionally, the immunosensor demonstrated a broad linear detection range spanning from 200 pg/mL to 2 ng/mL.

Changes in an enzyme ensemble during catalysis observed by high-resolution XFEL crystallography

Enzymes populate ensembles of structures necessary for catalysis that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography at an x-ray free electron laser to observe catalysis in a designed mutant isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations, and formation of the thioimidate intermediate selects for catalytically competent substates. The influence of cysteine ionization on the ICH ensemble is validated by determining structures of the enzyme at multiple pH values. Large molecular dynamics simulations in crystallo and time-resolved electron density maps show that Asp17 ionizes during catalysis and causes conformational changes that propagate across the dimer, permitting water to enter the active site for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

Screening, cloning, immobilization and application prospects of a novel β-glucosidase from the soil metagenome

The soil environment for straw return is a rich and valuable library containing many microorganisms and proteins. In this study, we aimed to screen a high-quality β-glucosidase (BGL) from the soil metagenomic library and to overcome the limitation of the low extraction rate of resveratrol in Polygonum cuspidatum. This includes the construction of a soil metagenomic library, screening of BGL, bioinformatics analysis, cloning, expression, immobilization, enzymatic property analysis, and application for the transformation of polydatin. The results showed that the soil metagenomic library of straw return was successfully constructed, and a novel BGL was screened. The identified 1356 bp long BGL belonged to the glycoside hydrolase 1 (GH1) family and was named Bgl1356. After successful cloning and expression of Bgl1356, it was immobilized using chitosan. The optimum temperature of immobilized Bgl1356 was 50 °C, and the pH was 5. It exhibited good tolerance for various metal ions (CO2+, Ni2+, Cu2+, Mn2+, Na2+, Ca2+, and Ag+) and organic solvents (DMSO, Triton-X-10, and ethanol). Enzymatic kinetics assays showed that Bgl1356 had good affinity for the substrate, and the specific enzyme activity was 234.03 U/mg. The conversion rate of polydatin by immobilized Bgl1356 was 95.70 ± 1.08%, facilitating the production of high amounts of resveratrol. Thus, this paper reports a novel temperature-, organic solvent-, and metal ion-tolerant BGL that has good application prospects in the pharmaceutical industry.

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

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

High-resolution double vision of the allosteric phosphatase PTP1B

Protein tyrosine phosphatase 1B (PTP1B) plays important roles in cellular homeostasis and is a highly validated therapeutic target for multiple human ailments, including diabetes, obesity and breast cancer. However, much remains to be learned about how conformational changes may convey information through the structure of PTP1B to enable allosteric regulation by ligands or functional responses to mutations. High-resolution X-ray crystallography can offer unique windows into protein conformational ensembles, but comparison of even high-resolution structures is often complicated by differences between data sets, including non-isomorphism. Here, the highest resolution crystal structure of apo wild-type (WT) PTP1B to date is presented out of a total of ∼350 PTP1B structures in the PDB. This structure is in a crystal form that is rare for PTP1B, with two unique copies of the protein that exhibit distinct patterns of conformational heterogeneity, allowing a controlled comparison of local disorder across the two chains within the same asymmetric unit. The conformational differences between these chains are interrogated in the apo structure and between several recently reported high-resolution ligand-bound structures. Electron-density maps in a high-resolution structure of a recently reported activating double mutant are also examined, and unmodeled alternate conformations in the mutant structure are discovered that coincide with regions of enhanced conformational heterogeneity in the new WT structure. These results validate the notion that these mutations operate by enhancing local dynamics, and suggest a latent susceptibility to such changes in the WT enzyme. Together, these new data and analysis provide a detailed view of the conformational ensemble of PTP1B and highlight the utility of high-resolution crystallography for elucidating conformational heterogeneity with potential relevance for function.

Substrate Positioning Dynamics Involves a Non-Electrostatic Component to Mediate Catalysis

Substrate positioning dynamics (SPD) orients the substrate in the active site, thereby influencing catalytic efficiency. However, it remains unknown whether SPD effects originate primarily from electrostatic perturbation inside the enzyme or can independently mediate catalysis with a significant non-electrostatic component. In this work, we investigated how the non-electrostatic component of SPD affects transition state (TS) stabilization. Using high-throughput enzyme modeling, we selected Kemp eliminase variants with similar electrostatics inside the enzyme but significantly different SPD. The kinetic parameters of these mutants were experimentally characterized. We observed a valley-shaped, two-segment linear correlation between the TS stabilization free energy (converted from kinetic parameters) and substrate positioning index (a metric to quantify SPD). The energy varies by approximately 2 kcal/mol. Favorable SPD was observed for the distal mutant R154W, increasing the proportion of reactive conformations and leading to the lowest activation free energy. These results indicate the substantial contribution of the non-electrostatic component of SPD to enzyme catalytic efficiency.

Mutexa: A Computational Ecosystem for Intelligent Protein Engineering

Protein engineering holds immense promise in shaping the future of biomedicine and biotechnology. This Review focuses on our ongoing development of Mutexa, a computational ecosystem designed to enable "intelligent protein engineering". In this vision, researchers will seamlessly acquire sequences of protein variants with desired functions as biocatalysts, therapeutic peptides, and diagnostic proteins through a finely-tuned computational machine, akin to Amazon Alexa's role as a versatile virtual assistant. The technical foundation of Mutexa has been established through the development of a database that combines and relates enzyme structures and their respective functions (e.g., IntEnzyDB), workflow software packages that enable high-throughput protein modeling (e.g., EnzyHTP and LassoHTP), and scoring functions that map the sequence-structure-function relationship of proteins (e.g., EnzyKR and DeepLasso). We will showcase the applications of these tools in benchmarking the convergence conditions of enzyme functional descriptors across mutants, investigating protein electrostatics and cavity distributions in SAM-dependent methyltransferases, and understanding the role of nonelectrostatic dynamic effects in enzyme catalysis. Finally, we will conclude by addressing the future steps and fundamental challenges in our endeavor to develop new Mutexa applications that assist the identification of beneficial mutants in protein engineering.

Native dynamics and allosteric responses in PTP1B probed by high-resolution HDX-MS

Protein tyrosine phosphatase 1B (PTP1B) is a validated therapeutic target for obesity, diabetes, and certain types of cancer. In particular, allosteric inhibitors hold potential for therapeutic use, but an incomplete understanding of conformational dynamics and allostery in this protein has hindered their development. Here, we interrogate solution dynamics and allosteric responses in PTP1B using high-resolution hydrogen-deuterium exchange mass spectrometry (HDX-MS), an emerging and powerful biophysical technique. Using HDX-MS, we obtain a detailed map of the solution dynamics of apo PTP1B, revealing several flexible loops interspersed among more constrained and rigid regions within the protein structure, as well as local regions that exchange faster than expected from their secondary structure and buriedness. We demonstrate that our HDX rate data obtained in solution adds value to predictions of dynamics derived from a pseudo-ensemble constructed from ~200 crystal structures of PTP1B. Furthermore, we report HDX-MS maps for PTP1B with active-site vs. allosteric small-molecule inhibitors. These maps reveal distinct, dramatic, and widespread effects on protein dynamics relative to the apo form, including changes to dynamics in locations distal (>35 Å) from the respective ligand binding sites. These results help shed light on the allosteric nature of PTP1B and the surprisingly far-reaching consequences of inhibitor binding in this important protein. Overall, our work showcases the potential of HDX-MS for elucidating protein conformational dynamics and allosteric effects of small-molecule ligands, and highlights the potential of integrating HDX-MS alongside other complementary methods to guide the development of new therapeutics.

Protein conformational ensembles in function: roles and mechanisms

The sequence-structure-function paradigm has dominated twentieth century molecular biology. The paradigm tacitly stipulated that for each sequence there exists a single, well-organized protein structure. Yet, to sustain cell life, function requires (i) that there be more than a single structure, (ii) that there be switching between the structures, and (iii) that the structures be incompletely organized. These fundamental tenets called for an updated sequence-conformational ensemble-function paradigm. The powerful energy landscape idea, which is the foundation of modernized molecular biology, imported the conformational ensemble framework from physics and chemistry. This framework embraces the recognition that proteins are dynamic and are always interconverting between conformational states with varying energies. The more stable the conformation the more populated it is. The changes in the populations of the states are required for cell life. As an example, in vivo, under physiological conditions, wild type kinases commonly populate their more stable "closed", inactive, conformations. However, there are minor populations of the "open", ligand-free states. Upon their stabilization, e.g., by high affinity interactions or mutations, their ensembles shift to occupy the active states. Here we discuss the role of conformational propensities in function. We provide multiple examples of diverse systems, including protein kinases, lipid kinases, and Ras GTPases, discuss diverse conformational mechanisms, and provide a broad outlook on protein ensembles in the cell. We propose that the number of molecules in the active state (inactive for repressors), determine protein function, and that the dynamic, relative conformational propensities, rather than the rigid structures, are the hallmark of cell life.

Insight into the High-Efficiency Benzo(a)pyrene Degradation Ability of Pseudomonas benzopyrenica BaP3 and Its Application in the Complete Bioremediation of Benzo(a)pyrene

Polycyclic aromatic hydrocarbons (PAHs) are common carcinogens. Benzo(a)pyrene is one of the most difficult high-molecular-weight (HMW) PAHs to remove. Biodegradation has become an ideal method to eliminate PAH pollutants from the environment. The existing research is mostly limited to low-molecular-weight PAHs; there is little understanding of HMW PAHs, particularly benzo(a)pyrene. Research into the biodegradation of HMW PAHs contributes to the development of microbial metabolic mechanisms and also provides new systems for environmental treatments. Pseudomonas benzopyrenica BaP3 is a highly efficient benzo(a)pyrene-degrading strain that is isolated from soil samples, but its mechanism of degradation remains unknown. In this study, we aimed to clarify the high degradation efficiency mechanism of BaP3. The genes encoding Rhd1 and Rhd2 in strain BaP3 were characterized, and the results revealed that rhd1 was the critical factor for high degradation efficiency. Molecular docking and enzyme activity determinations confirmed this conclusion. A recombinant strain that could completely mineralize benzo(a)pyrene was also proposed for the first time. We explained the mechanism of the high-efficiency benzo(a)pyrene degradation ability of BaP3 to improve understanding of the degradation mechanism of highly toxic PAHs and to provide new solutions to practical applications via synthetic biology.

Native and nonnative reactions in ethanolamine ammonia-lyase are actuated by different dynamics

We address the contribution of select classes of solvent-coupled configurational fluctuations to the complex choreography involved in configurational and chemical steps in an enzyme by comparing native and nonnative reactions conducted at different protein internal sites. The low temperature, first-order kinetics of covalent bond rearrangement of the cryotrapped substrate radical in coenzyme B12-dependent ethanolamine ammonia-lyase (EAL) from Salmonella enterica display a kink, or increase in slope, of the Arrhenius plot with decreasing temperature. The event is associated with quenching of a select class of reaction-actuating collective fluctuations in the protein hydration layer. For comparison, a nonnative, radical reaction of the protein interior cysteine sulfhydryl group with hydrogen peroxide (H2O2) is introduced by cryotrapping EAL in an aqueous H2O2 eutectic system. The low-temperature aqueous H2O2 protein hydration and mesodomain solvent phases surrounding cryotrapped EAL are characterized by using TEMPOL spin probe electron paramagnetic resonance spectroscopy, including a freezing transition of the eutectic phase that orders the protein hydration layer. Kinetics of the cysteine-H2O2 reaction in the EAL protein interior are monitored by DEPMPO spin trapping of hydroxyl radical product. In contrast to the native reaction, the linear Arrhenius relation for the nonnative cysteine-H2O2 reaction is maintained through the solvent-protein ordering transition. The nonnative reaction is coupled to the generic local, incremental fluctuations that are intrinsic to globular proteins. The comparative approach supports the proposal that select coupled solvent-protein configurational fluctuations actuate the native reaction, and suggests that select dynamical coupling contributes to the degree of catalysis in enzymes.

Low Molecular Weight Inhibitors Targeting the RNA-Binding Protein HuR

The RNA-binding protein human antigen R (HuR) regulates stability, translation, and nucleus-to-cytoplasm shuttling of its target mRNAs. This protein has been progressively recognized as a relevant therapeutic target for several pathologies, like cancer, neurodegeneration, as well as inflammation. Inhibitors of mRNA binding to HuR might thus be beneficial against a variety of diseases. Here, we present the rational identification of structurally novel HuR inhibitors. In particular, by combining chemoinformatic approaches, high-throughput virtual screening, and RNA-protein pulldown assays, we demonstrate that the 4-(2-(2,4,6-trioxotetrahydropyrimidin-5(2H)-ylidene)hydrazineyl)benzoate ligand exhibits a dose-dependent HuR inhibition effect in binding experiments. Importantly, the chemical scaffold is new with respect to the currently known HuR inhibitors, opening up a new avenue for the design of pharmaceutical agents targeting this important protein.

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.

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and time-resolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Assessment of Enzyme Functionality at Metal-Organic Framework Interfaces Developed through Molecular Simulations

The catalytic efficiency and unrivaled selectivity with which enzymes convert substrates to products have been tapped for widespread chemical transformations within biomedical technology, biofuel production, gas sensing, and the upgrading of commodity chemicals, just to name a few. However, the feasibility of enzymes implementation is challenged by the lack of reusability and loss of native catalytic activity due to the irreversible biocatalyst denaturation at high temperatures and in the presence of industrial solvents. Enzyme immobilization, a prerequisite for enzyme reusability, offers controllable strategies for increased functional viability of the biocatalyst in a synthetic environment. Herein we used molecular dynamics (MD) simulations and probed the noncovalent interactions between model enzymes of technological interest, i.e., carbonic anhydrase (CA) and myeloperoxidase (MPO), with selected metal-organic frameworks (MOFs; MIL-160 and ZIF-8) of proven industrial implementation. We found that the CA and MPO can bind to MIL-160 at optimal binding energies of 201 and 501 kJ mol^-1, respectively, that are strongly influenced by the increased incidence of hydrogen bonding between enzymes and the frameworks. The free energy of binding of enzymes to ZIF-8, on the other hand, was found to be less strongly influenced by hydrogen bonding networks relative to the occurrence of hydrophobic-hydrophobic interactions that yielded 106 kJ mol^-1 for CA and 201 kJ mol^-1 for MPO.

Overcoming Cytosolic Delivery Barriers of Proteins Using Denatured Protein-Conjugated Mesoporous Silica Nanoparticles

Intracellular delivery of therapeutic proteins has increased advantages over current small-molecule drugs and gene therapies, especially in therapeutic efficacies for a broad spectrum of diseases. Hence, developing the protein therapeutics approach provides a needed alternative. Here, we designed a mesoporous silica nanoparticle (MSN)-mediated protein delivery approach and demonstrated effective intracellular delivery of the denatured superoxide dismutase (SOD) protein, overcoming the delivery challenges and achieving higher enzymatic activity than native SOD-conjugated MSNs. The denatured SOD-conjugated MSN delivery strategy provides benefits of reduced size and steric hindrance, increased protein flexibility without distorting its secondary structure, exposure of the cell-penetrating peptide transactivator of transcription for enhanced efficient delivery, and a change in the corona protein composition, enabling cytosolic delivery. After delivery, SOD displayed a specific activity around threefold higher than in our previous reports. Furthermore, the in vivo biosafety and therapeutic potential for neuron therapy were evaluated, demonstrating the biocompatibility and the effective antioxidant effect in Neuro-2a cells that protected neurite outgrowth from paraquat-induced reactive oxygen species attack. This study offers an opportunity to realize the druggable possibility of cytosolic proteins using MSNs.

Convergence in determining enzyme functional descriptors across Kemp eliminase variants

Molecular simulations have been extensively employed to accelerate biocatalytic discoveries. Enzyme functional descriptors derived from molecular simulations have been leveraged to guide the search for beneficial enzyme mutants. However, the ideal active-site region size for computing the descriptors over multiple enzyme variants remains untested. Here, we conducted convergence tests for dynamics-derived and electrostatic descriptors on 18 Kemp eliminase variants across six active-site regions with various boundary distances to the substrate. The tested descriptors include the root-mean-square deviation of the active-site region, the solvent accessible surface area ratio between the substrate and active site, and the projection of the electric field (EF) on the breaking C-H bond. All descriptors were evaluated using molecular mechanics methods. To understand the effects of electronic structure, the EF was also evaluated using quantum mechanics/molecular mechanics methods. The descriptor values were computed for 18 Kemp eliminase variants. Spearman correlation matrices were used to determine the region size condition under which further expansion of the region boundary does not substantially change the ranking of descriptor values. We observed that protein dynamics-derived descriptors, including RMSD_{active_site} and SASA_ratio, converge at a distance cutoff of 5 Å from the substrate. The electrostatic descriptor, EF_C-H, converges at 6 Å using molecular mechanics methods with truncated enzyme models and 4 Å using quantum mechanics/molecular mechanics methods with whole enzyme model. This study serves as a future reference to determine descriptors for predictive modeling of enzyme engineering.

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.

Large-Scale Ligand Perturbations of the Protein Conformational Landscape Reveal State-Specific Interaction Hotspots

Protein flexibility is important for ligand binding but often ignored in drug design. Considering proteins as ensembles rather than static snapshots creates opportunities to target dynamic proteins that lack FDA-approved drugs, such as the human chaperone, heat shock protein 90 (Hsp90). Hsp90α accommodates ligands with a dynamic lid domain, yet no comprehensive analysis relating lid conformations to ligand properties is available. To date, ∼300 ligand-bound Hsp90α crystal structures are deposited in the Protein Data Bank, which enables us to consider ligand binding as a perturbation of the protein conformational landscape. By estimating binding site volumes, we classified structures into distinct major and minor lid conformations. Supported by retrospective docking, each conformation creates unique hotspots that bind chemically distinguishable ligands. Clustering revealed insightful exceptions and the impact of crystal packing. Overall, Hsp90α’s plasticity provides a cautionary tale of overinterpreting individual crystal structures and motivates an ensemble-based view of drug design.

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.

Probing ligand binding of endothiapepsin by `temperature-resolved' macromolecular crystallography

Continuous developments in cryogenic X-ray crystallography have provided most of our knowledge of 3D protein structures, which has recently been further augmented by revolutionary advances in cryoEM. However, a single structural conformation identified at cryogenic temperatures may introduce a fictitious structure as a result of cryogenic cooling artefacts, limiting the overview of inherent protein physiological dynamics, which play a critical role in the biological functions of proteins. Here, a room-temperature X-ray crystallographic method using temperature as a trigger to record movie-like structural snapshots has been developed. The method has been used to show how TL00150, a 175.15 Da fragment, undergoes binding-mode changes in endothiapepsin. A surprising fragment-binding discrepancy was observed between the cryo-cooled and physiological temperature structures, and multiple binding poses and their interplay with DMSO were captured. The observations here open up new promising prospects for structure determination and interpretation at physiological temperatures with implications for structure-based drug discovery.

Evaluating the impact of X-ray damage on conformational heterogeneity in room-temperature (277 K) and cryo-cooled protein crystals

Cryo-cooling has been nearly universally adopted to mitigate X-ray damage and facilitate crystal handling in protein X-ray crystallography. However, cryo X-ray crystallographic data provide an incomplete window into the ensemble of conformations that is at the heart of protein function and energetics. Room-temperature (RT) X-ray crystallography provides accurate ensemble information, and recent developments allow conformational heterogeneity (the experimental manifestation of ensembles) to be extracted from single-crystal data. Nevertheless, high sensitivity to X-ray damage at RT raises concerns about data reliability. To systematically address this critical issue, increasingly X-ray-damaged high-resolution data sets (1.02-1.52 Å resolution) were obtained from single proteinase K, thaumatin and lysozyme crystals at RT (277 K). In each case a modest increase in conformational heterogeneity with X-ray damage was observed. Merging data with different extents of damage (as is typically carried out) had negligible effects on conformational heterogeneity until the overall diffraction intensity decayed to ∼70% of its initial value. These effects were compared with X-ray damage effects in cryo-cooled crystals by carrying out an analogous analysis of increasingly damaged proteinase K cryo data sets (0.9-1.16 Å resolution). X-ray damage-associated heterogeneity changes were found that were not observed at RT. This property renders it difficult to distinguish real from artefactual conformations and to determine the conformational response to changes in temperature. The ability to acquire reliable heterogeneity information from single crystals at RT, together with recent advances in RT data collection at accessible synchrotron beamlines, provides a strong motivation for the widespread adoption of RT X-ray crystallography to obtain conformational ensemble information.

Microalgae as a potential sustainable solution to environment health

Cyanobacteria such as Spirulina platensis secretes numerous biomolecules while consuming CO₂ for photosynthesis which can reduce the environmental pollution as it can also be grown in wastewater. These biomolecules can be further processed in numerous pathways such as feed, fuel, pharmaceuticals, and nutraceuticals. This study aims to screen the potential molecular mechanisms of pigments from cyanobacteria as antidiabetic type-2 candidates through molecular docking. The activities of the test compounds were compared to commercial diabetic drugs, such as acarbose, linagliptin and polydatin. The results indicated that the binding affinity of pheophytin, β-carotene, and phycocyanobilin to α-amylase were 0.4, 2, and 2.6 kcal/mol higher than that of acarbose with α-amylase. Binding affinity between pheophytin, β-carotene, and phycocyanobilin with α-glucosidase were found to be comparable, which resulted 1.2, and 1.6 kcal/mol higher than that of acarbose with α-glucosidase. Meanwhile, binding activity of β-carotene and phycocyanobilin with DPP-IV were 0.5 and 0.3 kcal/mol higher than that of linagliptin with DPP-IV, whereas pheophytin, β-carotene, and phycocyanobilin with Glucose-6-phosphate dehydrogenase (G6PD) were 0.2, 1, and 1.4 kcal/mol higher from that of polydatin with G6PD. Moreover, pheophytin, β-carotene and phycocyanobilin were likely to inhibit α-amylase, α-glucosidase, and DPP-IV competitively, while uncompetitively for G6PD. Thus, the integration of molecular docking and experimental approach, such as in vitro and in vivo studies may greatly improve the discovery of true bioactive compounds in cyanobacteria for type 2 diabetes mellitus drugs and treatments.

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 B12 (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 B12-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 B12-dependent and related radical enzymes.

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.

High throughput and quantitative enzymology in the genomic era

Accurate predictions from models based on physical principles are the ultimate metric of our biophysical understanding. Although there has been stunning progress toward structure prediction, quantitative prediction of enzyme function has remained challenging. Realizing this goal will require large numbers of quantitative measurements of rate and binding constants and the use of these ground-truth data sets to guide the development and testing of these quantitative models. Ground truth data more closely linked to the underlying physical forces are also desired. Here, we describe technological advances that enable both types of ground truth measurements. These advances allow classic models to be tested, provide novel mechanistic insights, and place us on the path toward a predictive understanding of enzyme structure and function.

Glutathione Transferases as Efficient Ketosteroid Isomerases

In addition to their well-established role in detoxication, glutathione transferases (GSTs) have other biological functions. We are focusing on the ketosteroid isomerase activity, which appears to contribute to steroid hormone biosynthesis in mammalian tissues. A highly efficient GST A3-3 is present in some, but not all, mammals. The alpha class enzyme GST A3-3 in humans and the horse shows the highest catalytic efficiency with kcat/Km values of approximately 107 M-1s-1, ranking close to the most active enzymes known. The expression of GST A3-3 in steroidogenic tissues suggests that the enzyme has evolved to support the activity of 3β-hydroxysteroid dehydrogenase, which catalyzes the formation of 5-androsten-3,17-dione and 5-pregnen-3,20-dione that are substrates for the double-bond isomerization catalyzed by GST A3-3. The dehydrogenase also catalyzes the isomerization, but its kcat of approximately 1 s-1 is 200-fold lower than the kcat values of human and equine GST A3-3. Inhibition of GST A3-3 in progesterone-producing human cells suppress the formation of the hormone. Glutathione serves as a coenzyme contributing a thiolate as a base in the isomerase mechanism, which also involves the active-site Tyr9 and Arg15. These conserved residues are necessary but not sufficient for the ketosteroid isomerase activity. A proper assortment of H-site residues is crucial to efficient catalysis by forming the cavity binding the hydrophobic substrate. It remains to elucidate why some mammals, such as rats and mice, lack GSTs with the prominent ketosteroid isomerase activity found in certain other species. Remarkably, the fruit fly Drosophila melanogaster, expresses a GSTE14 with notable steroid isomerase activity, even though Ser14 has evolved as the active-site residue corresponding to Tyr9 in the mammalian alpha class.

Reverse Protonation of Buried Ion-Pairs in Staphylococcal Nuclease Mutants

Although buried titratable residues in protein cavities are often of major functional importance, it is generally challenging to understand their properties such as the ionization state and factors of stabilization based on experimental studies alone. A specific set of examples involve buried Glu-Lys pairs in a series of variants of Staphylococcal nuclease, for which recent structural and thermodynamic studies appeared to suggest that both the stability and the ionization state of the buried Glu-Lys pair are sensitive to its orientation (i.e., Glu23-Lys36 vs Lys23-Glu36). To further clarify the situation, especially ionization states of the buried Glu-Lys pairs, we have conducted extensive molecular dynamics simulations and free energy computations. Microsecond molecular dynamics simulations show that the hydration level of the cavity depends on the orientation of the buried ion-pair therein as well as its ionization state; free energy simulations recapitulate the relative stability of Glu23-Lys36 (EK) vs Lys23-Glu36 (KE) mutants measured experimentally, although the difference is similar in magnitude regardless of the ionization state of the Glu-Lys pair. A complementary set of free energy simulations strongly suggests that, in contrast to the original suggestion in the experimental analysis, the Glu and Lys residues prefer to adopt their charge-neutral rather than the ionized states. This result is consistent with the low dielectric constant computed for water in the cavity, which makes it difficult for the protein cavity to stabilize a pair of charged Glu-Lys residues, even with water penetration. The current study highlights the role of free energy simulations in understanding the ionization state of buried titratable residues and the relevant energetic contributions, forming the basis for the rational design of buried charge networks in proteins.

Fundamentals to function: Quantitative and scalable approaches for measuring protein stability

Folding a linear chain of amino acids into a three-dimensional protein is a complex physical process that ultimately confers an impressive range of diverse functions. Although recent advances have driven significant progress in predicting three-dimensional protein structures from sequence, proteins are not static molecules. Rather, they exist as complex conformational ensembles defined by energy landscapes spanning the space of sequence and conditions. Quantitatively mapping the physical parameters that dictate these landscapes and protein stability is therefore critical to develop models that are capable of predicting how mutations alter function of proteins in disease and informing the design of proteins with desired functions. Here, we review the approaches that are used to quantify protein stability at a variety of scales, from returning multiple thermodynamic and kinetic measurements for a single protein sequence to yielding indirect insights into folding across a vast sequence space. The physical parameters derived from these approaches will provide a foundation for models that extend beyond the structural prediction to capture the complexity of conformational ensembles and, ultimately, their function.

Discovery of a hidden transient state in all bromodomain families

Bromodomains (BDs) are small protein modules that interact with acetylated marks in histones. These posttranslational modifications are pivotal to regulate gene expression, making BDs promising targets to treat several diseases. While the general structure of BDs is well known, their dynamical features and their interplay with other macromolecules are poorly understood, hampering the rational design of potent and selective inhibitors. Here, we combine extensive molecular dynamics simulations, Markov state modeling, and available structural data to reveal a transiently formed state that is conserved across all BD families. It involves the breaking of two backbone hydrogen bonds that anchor the ZA-loop with the αA helix, opening a cryptic pocket that partially occludes the one associated to histone binding. By analyzing more than 1,900 experimental structures, we unveil just two adopting the hidden state, explaining why it has been previously unnoticed and providing direct structural evidence for its existence. Our results suggest that this state is an allosteric regulatory switch for BDs, potentially related to a recently unveiled BD-DNA-binding mode.

Modulating Enzyme Function via Dynamic Allostery within Biliverdin Reductase B

The biliverdin reductase B (BLVRB) class of enzymes catalyze the NADPH-dependent reduction of multiple flavin substrates and are emerging as critical players in cellular redox regulation. However, the role of dynamics and allostery have not been addressed, prompting studies here that have revealed a position 15 Å away from the active site within human BLVRB (T164) that is inherently dynamic and can be mutated to control global micro-millisecond motions and function. By comparing the inherent dynamics through nuclear magnetic resonance (NMR) relaxation approaches of evolutionarily distinct BLVRB homologues and by applying our previously developed Relaxation And Single Site Multiple Mutations (RASSMM) approach that monitors both the functional and dynamic effects of multiple mutations to the single T164 site, we have discovered that the most dramatic mutagenic effects coincide with evolutionary changes and these modulate coenzyme binding. Thus, evolutionarily changing sites distal to the active site serve as dynamic "dials" to globally modulate motions and function. Despite the distal dynamic and functional coupling modulated by this site, micro-millisecond motions span an order of magnitude in their apparent kinetic rates of motions. Thus, global dynamics within BLVRB are a collection of partially coupled motions tied to catalytic function.