Conformational Disorganization within the Active Site of a Recently Evolved Organophosphate Hydrolase Limits Its Catalytic Efficiency

Selection of a promiscuous minimalist cAMP phosphodiesterase from a library of de novo designed proteins

The ability of unevolved amino acid sequences to become biological catalysts was key to the emergence of life on Earth. However, billions of years of evolution separate complex modern enzymes from their simpler early ancestors. To probe how unevolved sequences can develop new functions, we use ultrahigh-throughput droplet microfluidics to screen for phosphoesterase activity amidst a library of more than one million sequences based on a de novo designed 4-helix bundle. Characterization of hits revealed that acquisition of function involved a large jump in sequence space enriching for truncations that removed >40% of the protein chain. Biophysical characterization of a catalytically active truncated protein revealed that it dimerizes into an α-helical structure, with the gain of function accompanied by increased structural dynamics. The identified phosphodiesterase is a manganese-dependent metalloenzyme that hydrolyses a range of phosphodiesters. It is most active towards cyclic AMP, with a rate acceleration of ~109 and a catalytic proficiency of >1014 M-1, comparable to larger enzymes shaped by billions of years of evolution.

Dynamic Roles of Insect Carboxyl/Cholinesterases in Chemical Adaptation

Insects have evolved several intricate defense mechanisms to adapt to their chemical environment. Due to their versatile capabilities in hydrolytic biotransformation, insect carboxyl/cholinesterases (CCEs) play vital roles in the development of pesticide resistance, facilitating the adaptation of insects to their host plants, and manipulating insect behaviors through the olfaction system. CCEs confer insecticide resistance through the mechanisms of qualitative or quantitative changes of CCE-mediated enhanced metabolism or target-site insensitivity, and may contribute to the host plant adaptation. CCEs represent the first odorant-degrading enzymes (ODEs) discovered to degrade insect pheromones and plant odors and remain the most promising ODE candidates. Here, we summarize insect CCE classification, currently characterized insect CCE protein structure characteristics, and the dynamic roles of insect CCEs in chemical adaptation.

Ultrahigh-Throughput Directed Evolution of a Metal-Free α/β-Hydrolase with a Cys-His-Asp Triad into an Efficient Phosphotriesterase

Finding new mechanistic solutions for biocatalytic challenges is key in the evolutionary adaptation of enzymes, as well as in devising new catalysts. The recent release of man-made substances into the environment provides a dynamic testing ground for observing biocatalytic innovation at play. Phosphate triesters, used as pesticides, have only recently been introduced into the environment, where they have no natural counterpart. Enzymes have rapidly evolved to hydrolyze phosphate triesters in response to this challenge, converging onto the same mechanistic solution, which requires bivalent cations as a cofactor for catalysis. In contrast, the previously identified metagenomic promiscuous hydrolase P91, a homologue of acetylcholinesterase, achieves slow phosphotriester hydrolysis mediated by a metal-independent Cys-His-Asp triad. Here, we probe the evolvability of this new catalytic motif by subjecting P91 to directed evolution. By combining a focused library approach with the ultrahigh throughput of droplet microfluidics, we increase P91's activity by a factor of ≈360 (to a k_cat/K_M of ≈7 × 10⁵ M^-1 s^-1) in only two rounds of evolution, rivaling the catalytic efficiencies of naturally evolved, metal-dependent phosphotriesterases. Unlike its homologue acetylcholinesterase, P91 does not suffer suicide inhibition; instead, fast dephosphorylation rates make the formation of the covalent adduct rather than its hydrolysis rate-limiting. This step is improved by directed evolution, with intermediate formation accelerated by 2 orders of magnitude. Combining focused, combinatorial libraries with the ultrahigh throughput of droplet microfluidics can be leveraged to identify and enhance mechanistic strategies that have not reached high efficiency in nature, resulting in alternative reagents with novel catalytic machineries.

Excavating phytochemicals from plants possessing antiviral activities for identifying SARS-CoV hemagglutinin-esterase inhibitors by diligent computational workflow

Coronaviruses (CoVs) belong to a group of RNA viruses that cause diseases in vertebrates including. Newer and deadlier than SARS CoV-2 are sought to appear in future for which the scientific community must be prepared with the strategies for their control. Spike protein (S-protein) of all the CoVs require angiotensin-converting enzyme2 (ACE2), while CoVs also require hemagglutinin-acetylesterase (HE) glycoprotein receptor to simultaneously interact with O-acetylated sialic acids on host cells, both these interactions enable viral particle to enter host cell leading to its infection. Target inhibition of viral S-protein and HE glycoprotein receptor can lead to a development of therapy against the SARS CoV-2. The proposition is to recognize molecules from the bundle of phytochemicals of medicinal plants known to possess antiviral potentials as a lead that could interact and mask the active site of, HE glycoprotein which would ideally bind to O-acetylated sialic acids on human host cells. Such molecules can be addressed as 'HE glycoprotein blockers'. A library of 110 phytochemicals from Withania somnifera, Asparagus racemosus, Zinziber officinalis, Allium sativum, Curcuma longa and Adhatoda vasica was constructed and was used under present study. In silico analysis was employed with plant-derived phytochemicals. The molecular docking, molecular dynamics simulations over the scale of 1000 ns (1 μs) and ADMET prediction revealed that the Withania somnifera (ashwagandha) and Asparagus racemosus (shatavari) plants possessed various steroidal saponins and alkaloids which could potentially inhibit the COVID-19 virus and even other CoVs targeted HE glycoprotein receptor.

Paraoxonase Mimic by a Nanoreactor Aggregate Containing Benzimidazolium Calix and l-Histidine: Demonstration of the Acetylcholine Esterase Activity

An anion-mediated preorganization approach was used to design and synthesize the benzimidazolium-based calix compound R1⋅2 ClO₄ ^- . X-ray crystallography analysis revealed that the hydrogen-bonding interactions between the benzimidazolium cations and N,N-dimethylformamide (DMF) helped R1⋅2 ClO₄ ^- encapsulate DMF molecule(s). A nanoreactor, with R1⋅2 ClO₄ ^- and l-histidine (l-His) as the components, was fabricated by using a neutralization method. The nanoreactor could detoxify paraoxon in 30 min. l-His played a vital role in this process. Paraoxonase is a well-known enzyme used for pesticide degradation. The Ellman's reagent was used to determine the percentage inhibition of the acetylcholinesterase (AChE) activity in the presence of the nanoreactor. The results indicated that the nanoreactor inhibited AChE inhibition.

The evolution and engineering of enzyme activity through tuning conformational landscapes

Proteins are dynamic molecules whose structures consist of an ensemble of conformational states. Dynamics contribute to protein function and a link to protein evolution has begun to emerge. This increased appreciation for the evolutionary impact of conformational sampling has grown from our developing structural biology capabilities and the exploration of directed evolution approaches, which have allowed evolutionary trajectories to be mapped. Recent studies have provided empirical examples of how proteins can evolve via conformational landscape alterations. Moreover, minor conformational substates have been shown to be involved in the emergence of new enzyme functions as they can become enriched through evolution. The role of remote mutations in stabilizing new active site geometries has also granted insight into the molecular basis underpinning poorly understood epistatic effects that guide protein evolution. Finally, we discuss how the growth of our understanding of remote mutations is beginning to refine our approach to engineering enzymes.

Overcoming insecticide resistance through computational inhibitor design

Insecticides allow control of agricultural pests and disease vectors and are vital for global food security and health. The evolution of resistance to insecticides, such as organophosphates (OPs), is a serious and growing concern. OP resistance often involves sequestration or hydrolysis of OPs by carboxylesterases. Inhibiting carboxylesterases could, therefore, restore the effectiveness of OPs for which resistance has evolved. Here, we use covalent virtual screening to produce nano-/picomolar boronic acid inhibitors of the carboxylesterase αE7 from the agricultural pest Lucilia cuprina as well as a common Gly137Asp αE7 mutant that confers OP resistance. These inhibitors, with high selectivity against human acetylcholinesterase and low to no toxicity in human cells and in mice, act synergistically with the OPs diazinon and malathion to reduce the amount of OP required to kill L. cuprina by up to 16-fold and abolish resistance. The compounds exhibit broad utility in significantly potentiating another OP, chlorpyrifos, against the common pest, the peach-potato aphid (Myzus persicae). These compounds represent a solution to OP resistance as well as to environmental concerns regarding overuse of OPs, allowing significant reduction of use without compromising efficacy.

Protein engineering: the potential of remote mutations

Engineered proteins, especially enzymes, are now commonly used in many industries owing to their catalytic power, specific binding of ligands, and properties as materials and food additives. As the number of potential uses for engineered proteins has increased, the interest in engineering or designing proteins to have greater stability, activity and specificity has increased in turn. With any rational engineering or design pursuit, the success of these endeavours relies on our fundamental understanding of the systems themselves; in the case of proteins, their structure-dynamics-function relationships. Proteins are most commonly rationally engineered by targeting the residues that we understand to be functionally important, such as enzyme active sites or ligand-binding sites. This means that the majority of the protein, i.e. regions remote from the active- or ligand-binding site, is often ignored. However, there is a growing body of literature that reports on, and rationalises, the successful engineering of proteins at remote sites. This minireview will discuss the current state of the art in protein engineering, with a particular focus on engineering regions that are remote from active- or ligand-binding sites. As the use of protein technologies expands, exploiting the potential improvements made possible through modifying remote regions will become vital if we are to realise the full potential of protein engineering and design.

The evolution of multiple active site configurations in a designed enzyme

Developments in computational chemistry, bioinformatics, and laboratory evolution have facilitated the de novo design and catalytic optimization of enzymes. Besides creating useful catalysts, the generation and iterative improvement of designed enzymes can provide valuable insight into the interplay between the many phenomena that have been suggested to contribute to catalysis. In this work, we follow changes in conformational sampling, electrostatic preorganization, and quantum tunneling along the evolutionary trajectory of a designed Kemp eliminase. We observe that in the Kemp Eliminase KE07, instability of the designed active site leads to the emergence of two additional active site configurations. Evolutionary conformational selection then gradually stabilizes the most efficient configuration, leading to an improved enzyme. This work exemplifies the link between conformational plasticity and evolvability and demonstrates that residues remote from the active sites of enzymes play crucial roles in controlling and shaping the active site for efficient catalysis.

Evolutionary repurposing of a sulfatase: A new Michaelis complex leads to efficient transition state charge offset

The recruitment and evolutionary optimization of promiscuous enzymes is key to the rapid adaptation of organisms to changing environments. Our understanding of the precise mechanisms underlying enzyme repurposing is, however, limited: What are the active-site features that enable the molecular recognition of multiple substrates with contrasting catalytic requirements? To gain insights into the molecular determinants of adaptation in promiscuous enzymes, we performed the laboratory evolution of an arylsulfatase to improve its initially weak phenylphosphonate hydrolase activity. The evolutionary trajectory led to a 100,000-fold enhancement of phenylphosphonate hydrolysis, while the native sulfate and promiscuous phosphate mono- and diester hydrolyses were only marginally affected (≤50-fold). Structural, kinetic, and in silico characterizations of the evolutionary intermediates revealed that two key mutations, T50A and M72V, locally reshaped the active site, improving access to the catalytic machinery for the phosphonate. Measured transition state (TS) charge changes along the trajectory suggest the creation of a new Michaelis complex (E•S, enzyme-substrate), with enhanced leaving group stabilization in the TS for the promiscuous phosphonate (β_leavinggroup from -1.08 to -0.42). Rather than altering the catalytic machinery, evolutionary repurposing was achieved by fine-tuning the molecular recognition of the phosphonate in the Michaelis complex, and by extension, also in the TS. This molecular scenario constitutes a mechanistic alternative to adaptation solely based on enzyme flexibility and conformational selection. Instead, rapid functional transitions between distinct chemical reactions rely on the high reactivity of permissive active-site architectures that allow multiple substrate binding modes.

Optimization of Cholinesterase-Based Catalytic Bioscavengers Against Organophosphorus Agents

Organophosphorus agents (OPs) are irreversible inhibitors of acetylcholinesterase (AChE). OP poisoning causes major cholinergic syndrome. Current medical counter-measures mitigate the acute effects but have limited action against OP-induced brain damage. Bioscavengers are appealing alternative therapeutic approach because they neutralize OPs in bloodstream before they reach physiological targets. First generation bioscavengers are stoichiometric bioscavengers. However, stoichiometric neutralization requires administration of huge doses of enzyme. Second generation bioscavengers are catalytic bioscavengers capable of detoxifying OPs with a turnover. High bimolecular rate constants (k_cat/K_m > 10⁶ M⁻¹min⁻¹) are required, so that low enzyme doses can be administered. Cholinesterases (ChE) are attractive candidates because OPs are hemi-substrates. Moderate OP hydrolase (OPase) activity has been observed for certain natural ChEs and for G117H-based human BChE mutants made by site-directed mutagenesis. However, before mutated ChEs can become operational catalytic bioscavengers their dephosphylation rate constant must be increased by several orders of magnitude. New strategies for converting ChEs into fast OPase are based either on combinational approaches or on computer redesign of enzyme. The keystone for rational conversion of ChEs into OPases is to understand the reaction mechanisms with OPs. In the present work we propose that efficient OP hydrolysis can be achieved by re-designing the configuration of enzyme active center residues and by creating specific routes for attack of water molecules and proton transfer. Four directions for nucleophilic attack of water on phosphorus atom were defined. Changes must lead to a novel enzyme, wherein OP hydrolysis wins over competing aging reactions. Kinetic, crystallographic, and computational data have been accumulated that describe mechanisms of reactions involving ChEs. From these studies, it appears that introducing new groups that create a stable H-bonded network susceptible to activate and orient water molecule, stabilize transition states (TS), and intermediates may determine whether dephosphylation is favored over aging. Mutations on key residues (L286, F329, F398) were considered. QM/MM calculations suggest that mutation L286H combined to other mutations favors water attack from apical position. However, the aging reaction is competing. Axial direction of water attack is not favorable to aging. QM/MM calculation shows that F329H+F398H-based multiple mutants display favorable energy barrier for fast reactivation without aging.

Pharmacophore-based virtual screening of catechol-o-methyltransferase (COMT) inhibitors to combat Alzheimer's disease

Alzheimer's disease (AD) is one of the most significant neurodegenerative disorders and its symptoms mostly appear in aged people. Catechol-o-methyltransferase (COMT) is one of the known target enzymes responsible for AD. With the use of 23 known inhibitors of COMT, a query has been generated and validated by screening against the database of 1500 decoys to obtain the GH score and enrichment value. The crucial features of the known inhibitors were evaluated by the online ZINC Pharmer to identify new leads from a ZINC database. Five hundred hits were retrieved from ZINC Pharmer and by ADMET (absorption, distribution, metabolism, excretion, and toxicity) filtering by using FAF-Drug-3 and 36 molecules were considered for molecular docking. From the COMT inhibitors, opicapone, fenoldopam, and quercetin were selected, while ZINC63625100_413 ZINC39411941_412, ZINC63234426_254, ZINC63637968_451, and ZINC64019452_303 were chosen for the molecular dynamics simulation analysis having high binding affinity and structural recognition. This study identified the potential COMT inhibitors through pharmacophore-based inhibitor screening leading to a more complete understanding of molecular-level interactions.

Shuffling Active Site Substate Populations Affects Catalytic Activity: The Case of Glucose Oxidase

Glucose oxidase has wide applications in the pharmaceutical, chemical, and food industries. Many recent studies have enhanced key properties of this enzyme using directed evolution, yet without being able to reveal why these mutations are actually beneficial. This work presents a synergistic combination of experimental and computational methods, indicating how mutations, even when distant from the active site, positively affect glucose oxidase catalysis. We have determined the crystal structures of glucose oxidase mutants containing molecular oxygen in the active site. The catalytically important His516 residue has been previously shown to be flexible in the wild-type enzyme. The molecular dynamics simulations performed in this work allow us to quantify this floppiness, revealing that His516 exists in two states: catalytic and noncatalytic. The relative populations of these two substates are almost identical in the wild-type enzyme, with His516 readily shuffling between them. In the glucose oxidase mutants, on the other hand, the mutations enrich the catalytic His516 conformation and reduce the flexibility of this residue, leading to an enhancement in their catalytic efficiency. This study stresses the benefit of active site preorganization with respect to enzyme conversion rates by reducing molecular reorientation needs. We further suggest that the computational approach based on Hamiltonian replica exchange molecular dynamics, used in this study, may be a general approach to screening in silico for improved enzyme variants involving flexible catalytic residues.

Mapping the Accessible Conformational Landscape of an Insect Carboxylesterase Using Conformational Ensemble Analysis and Kinetic Crystallography

The proper function of enzymes often depends upon their efficient interconversion between particular conformational sub-states on a free-energy landscape. Experimentally characterizing these sub-states is challenging, which has limited our understanding of the role of protein dynamics in many enzymes. Here, we have used a combination of kinetic crystallography and detailed analysis of crystallographic protein ensembles to map the accessible conformational landscape of an insect carboxylesterase (LcαE7) as it traverses all steps in its catalytic cycle. LcαE7 is of special interest because of its evolving role in organophosphate insecticide resistance. Our results reveal that a dynamically coupled network of residues extends from the substrate-binding site to a surface loop. Interestingly, the coupling of this network that is apparent in the apoenzyme appears to be reduced in the phosphorylated enzyme intermediate. Altogether, the results of this work highlight the importance of protein dynamics to enzyme function and the evolution of new activity.