Fixation of the two Tabun isomers in acetylcholinesterase: a QM/MM study

De Novo Computational Design of a Lipase with Hydrolysis Activity towards Middle-Chained Fatty Acid Esters

Innovations in biocatalysts provide great prospects for intolerant environments or novel reactions. Due to the limited catalytic capacity and the long-term and labor-intensive characteristics of mining enzymes with the desired functions, de novo enzyme design was developed to obtain industrial application candidates in a rapid and convenient way. Here, based on the catalytic mechanisms and the known structures of proteins, we proposed a computational protein design strategy combining de novo enzyme design and laboratory-directed evolution. Starting with the theozyme constructed using a quantum-mechanical approach, the theoretical enzyme-skeleton combinations were assembled and optimized via the Rosetta "inside-out" protocol. A small number of designed sequences were experimentally screened using SDS-PAGE, mass spectrometry and a qualitative activity assay in which the designed enzyme 1a8uD₁ exhibited a measurable hydrolysis activity of 24.25 ± 0.57 U/g towards p-nitrophenyl octanoate. To improve the activity of the designed enzyme, molecular dynamics simulations and the RosettaDesign application were utilized to further optimize the substrate binding mode and amino acid sequence, thus keeping the residues of theozyme intact. The redesigned lipase 1a8uD₁-M8 displayed enhanced hydrolysis activity towards p-nitrophenyl octanoate-3.34 times higher than that of 1a8uD₁. Meanwhile, the natural skeleton protein (PDB entry 1a8u) did not display any hydrolysis activity, confirming that the hydrolysis abilities of the designed 1a8uD₁ and the redesigned 1a8uD₁-M8 were devised from scratch. More importantly, the designed 1a8uD₁-M8 was also able to hydrolyze the natural middle-chained substrate (glycerol trioctanoate), for which the activity was 27.67 ± 0.69 U/g. This study indicates that the strategy employed here has great potential to generate novel enzymes exhibiting the desired reactions.

Quantum mechanics implementation in drug-design workflows: does it really help?

The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?

On the Influence of the Protonation States of Active Site Residues on AChE Reactivation: A QM/MM Approach

Acetylcholinesterase (AChE), an enzyme of the serine hydrolase superfamily, is a mediator of signal transmission at cholinergic synapses by catalyzing acetylcholine cleavage into acetate and choline. This enzyme is vulnerable to covalent inhibition by organophosphate compounds (like VX). Covalent inhibition of AChE does not revert spontaneously. Known reactivator compounds have limited action in restoring catalytic activity. QM/MM simulations of VX-inhibited AChE reactivation by pralidoxime (2-PAM), a classical reactivator, were performed. These afforded a broad view of the effect of protonation states of active-site residues, and provide evidence for the role of Glu202, which needs to be protonated for reactivation to occur. In situ deprotonation of 2-PAM for both protonation states of Glu202 showed that His447 is able to deprotonate 2-PAM with the assistance of Glu202. Because the active site of serine hydrolases is highly conserved, this work provides new insights on the interplay between the catalytic triad residues and this glutamate, newly identified as protonatable.

How is acetylcholinesterase phosphonylated by soman? An ab initio QM/MM molecular dynamics study

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nerve system that hydrolyzes acetylcholine (ACh) and terminates synaptic signals by reducing the effective concentration of ACh in the synaptic clefts. Organophosphate compounds irreversibly inhibit AChEs, leading to irreparable damage to nerve cells. By employing Born–Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the covalent inhibition mechanism between AChE and the nerve toxin soman and determined its free energy profile for the first time. Our results indicate that phosphonylation of the catalytic serine by soman employs an addition–elimination mechanism, which is highly associative and stepwise: in the initial addition step, which is also rate-limiting, His440 acts as a general base to facilitate the nucleophilic attack of Ser200 on the soman’s phosphorus atom to form a trigonal bipyrimidal pentacovalent intermediate; in the subsequent elimination step, Try121 of the catalytic gorge stabilizes the leaving fluorine atom prior to its dissociation from the active site. Together with our previous characterization of the aging mechanism of soman inhibited AChE, our simulations have revealed detailed molecular mechanistic insights into the damaging function of the nerve agent soman.

Design of activated serine-containing catalytic triads with atomic-level accuracy

A challenge in the computational design of enzymes is that multiple properties, including substrate binding, transition state stabilization and product release, must be simultaneously optimized, and this has limited the absolute activity of successful designs. Here, we focus on a single critical property of many enzymes: the nucleophilicity of an active site residue that initiates catalysis. We design proteins with idealized serine-containing catalytic triads and assess their nucleophilicity directly in native biological systems using activity-based organophosphate probes. Crystal structures of the most successful designs show unprecedented agreement with computational models, including extensive hydrogen bonding networks between the catalytic triad (or quartet) residues, and mutagenesis experiments demonstrate that these networks are critical for serine activation and organophosphate reactivity. Following optimization by yeast display, the designs react with organophosphate probes at rates comparable to natural serine hydrolases. Co-crystal structures with diisopropyl fluorophosphate bound to the serine nucleophile suggest that the designs could provide the basis for a new class of organophosphate capture agents.

Molecular modeling evidence for His438 flip in the mechanism of butyrylcholinesterase hysteretic behavior

Cholinesterases display a hysteretic behavior with certain substrates and irreversible inhibitors. For years, this behavior has remained puzzling. However, several lines of evidence indicated that it is caused by perturbation of the catalytic triad and its water environment. In the present study, using molecular dynamics simulations of Ala328Cys BuChE mutant and wild-type BuChE in the absence and presence of a co-solvent (sucrose, glycerol), we provide evidence that hysteresis originates in a flip of the catalytic triad histidine (His438). This event is controlled by water molecules that interact with active site residues. The physiological significance of this phenomenon is still an issue.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

Aging mechanism of soman inhibited acetylcholinesterase

Acetylcholinesterase (AChE) is a crucial enzyme in the cholinergic nervous system that hydrolyzes neurotransmitter acetylcholine (ACh) and terminates synaptic signals. The catalytic serine of AChE can be phosphonylated by soman, one of the most potent nerve agents, and subsequently undergo an aging reaction. This phosphonylation and aging process leads to irreversible AChE inhibition, results in accumulation of excess ACh at the synaptic clefts, and causes neuromuscular paralysis. By employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations with umbrella sampling, a state-of-the-art approach to simulate enzyme reactions, we have characterized the aging mechanism of soman phosphonylated AChE and determined its free energy profile. This aging reaction starts with the scission of the O2-Cα bond, which is followed by methyl migration, and results in a tertiary carbenium intermediate. At the transition state, the scissile O2-Cα bond is already cleaved with an average O-C distance of 3.2 ± 0.3 Å and the migrating methyl group is shared between Cα and Cβ carbons with C-C distances of 1.9 ± 0.1 and 1.8 ± 0.1 Å, respectively. The negatively charged phosphonate group is stabilized by a salt bridge with the imidazole ring of the catalytic histidine. A major product of aging, 2,3-dimethyl-2-butanol can be formed swiftly by the reaction of a water molecule. Our characterized mechanism and simulation results provide new detailed insights into this important biochemical process.

Theoretical study on the aging and reactivation mechanism of tabun-inhibited acetylcholinesterase by using the quantum mechanical / molecular mechanical method

The organophosphorous compound tabun is highly neurotoxic because of its irreversible inhibition on acetylcholinesterase (AChE). It is wildly used as a warfare agent in the military. In this work, the aging and reactivation mechanism of tabun-inhibited AChE were studied by using the quantum mechanical / molecular mechanical (QM/MM) method. Geometry optimization of the stationary points were performed at the B3LYP/6–31G(d) level. Single-point energies were computed at the B3LYP/6–311++G(d,p) level. On the basis of the QM/MM results, a conclusion that the C–O bond scission is caused by water attack on the ethoxy group in the aging mechanism can be drawn. The reactivation process initialed by the antidotes CH2NO– or HLÖ-7 consists of three elemental steps, the nucleophilic attack on the P atom by the antidote, the dephosphorylation process, and the decomposition of the antidote–tabun complex. The highest energy barriers of the aging reaction, CH2NO–-induced reactivation, and HLÖ-7-induced reactivation are 19.9, 20.0, and 14.8 kcal/mol (1 cal = 4.184 J), respectively. The relative lower overall energy barrier of HLÖ-7-induced reactivation compared with that of the aging reaction indicates that HLÖ-7 is able to reactivate tabun-inhibited AChE. In addition, whether a newly designed antidote is able to reactivate tabun-inhibited AChE can be examined by the inequation X < 19.9 kcal/mol,where X means the highest energy barrier of the reactivation reaction of the newly designed antidote.

X-ray crystallographic snapshots of reaction intermediates in the G117H mutant of human butyrylcholinesterase, a nerve agent target engineered into a catalytic bioscavenger

OPs (organophosphylates) exert their acute toxicity through inhibition of acetylcholinesterase, by phosphylation of the catalytic serine residue. Engineering of human butyrylcholinesterase, by substitution of a histidine residue for the glycine residue at position 117, led to the creation of OP hydrolase activity. However, the lack of structural information and poor understanding of the hydrolytic mechanism of the G117H mutant has hampered further improvements in the catalytic activity. We have solved the crystallographic structure of the G117H mutant with a variety of ligands in its active site. A sulfate anion bound to the active site suggested the positioning for an OP prior to phosphylation. A fluoride anion was found in the active site when NaF was added to the crystallization buffer. In the fluoride complex, the imidazole ring from the His117 residue was substantially shifted, adopting a relaxed conformation probably close to that of the unliganded mutant enzyme. Additional X-ray structures were obtained from the transient covalent adducts formed upon reaction of the G117H mutant with the OPs echothiophate and VX [ethyl ({2-[bis(propan-2-yl)amino]ethyl}sulfanyl](methyl)phosphinate]. The position of the His117 residue shifted in response to the introduction of these adducts, overlaying the phosphylserine residue. These structural data suggest that the dephosphylation mechanism involves either a substantial conformational change of the His117 residue or an adjacent nucleophilic substitution by water.

Reaction profiles of the interaction between sarin and acetylcholinesterase and the S203C mutant: model nucleophiles and QM/MM potential energy surfaces

The phosphonylation mechanism of AChE and the S203C mutation by sarin (GB) is evaluated using two reaction schemes: a small model nucleophile (ethoxide, CH(3)CH(2)O(-)) and quantum mechanical/molecular mechanical (QM/MM) simulations. Calculations utilizing small model nucleophiles indicate that the reaction barrier for addition to GB is the rate-limiting step for both ethoxide and ethyl thiolate (CH(3)CH(2)S(-)); moreover, the activation barrier for addition to the phosphorus center of GB by ethyl thiolate is significantly larger (13.2 kcal/mol) than for ethoxide (8.3 kcal/mol). The decomposition transition state for both nucleophiles was determined to be approximately 1 kcal/mol. QM/MM simulations for AChE suggest a similar reaction mechanism for phosphonylation of the catalytic S203; however, the relative energetics are altered significantly compared to the isolated system. QM/MM results indicate that formation of the penta-coordinate intermediate is the rate-limiting step in the enzymatic system, with an activation barrier of 3.6 kcal/mol. Hydrogen-bonding interactions between the fluoride leaving group of GB with Y124 in AChE are observed throughout the reaction profile. The S203C mutation alters the relative energetics of the reaction, increasing the energy barrier for formation of the penta-coordinate intermediate to a value of 4.5 kcal/mol; moreover, the penta-coordinate intermediate (as product) is stabilized by an additional 6 kcal/mol when compared to wild-type AChE.

Aging mechanism of butyrylcholinesterase inhibited by an N-methyl analogue of tabun: implications of the trigonal-bipyramidal transition state rearrangement for the phosphylation or reactivation of cholinesterases

Cholinesterases are the main target of organophosphorus nerve agents (OPs). Their inhibition results in cholinergic syndrome and death. The enzymes are inhibited by phosphylation of the catalytic serine enzyme, but can be reactivated by oximes to some extent. However, phosphylated cholinesterases undergo a side reaction that progressively prevents their reactivatability. This unimolecular reaction, termed "aging", has been investigated for decades. It was shown that most OP-ChE conjugates aged by O-dealkylation of an alkoxy substituent of the phosphorus atom, a mechanism involving the stabilization of a transient carbocation. In this paper we present structural data supporting a substitution-based mechanism for aging of the huBChE conjugate of an N-mono-methyl analogue of tabun. This mechanism involves an adjacent nucleophilic attack followed by Berry pseudorotation. A similar adjacent attack and subsequent rearrangement of the transition state have been recently proposed for tabun phosphylation of AChE. We suggest that a similar mechanism is also possible for oxime reactivation of phosphylated cholinesterases. This opens new perspectives in terms of reactivator design.