Molecular dynamics study of active-site interactions with tetracoordinate transients in acetylcholinesterase and its mutants

2001 ◽  
Vol 353 (3) ◽  
pp. 645-653 ◽  
Author(s):  
Istvan J. ENYEDY ◽  
Ildiko M. KOVACH ◽  
Akos BENCSURA

The role of active-site residues in the dealkylation reaction in the PSCS diastereomer of 2-(3,3-dimethylbutyl)methylphosphonofluoridate (soman)-inhibited Torpedo californicaacetylcholinesterase (AChE) was investigated by full-scale molecular dynamics simulations using CHARMM: > 400ps equilibration was followed by 150–200ps production runs with the fully solvated tetracoordinate phosphonate adduct of the wild-type, Trp84Ala and Gly199Gln mutants of AChE. Parallel simulations were carried out with the tetrahedral intermediate formed between serine-200 Oγ of AChE and acetylcholine. We found that the NεH in histidine H+-440 is positioned to protonate the oxygen in choline and thus promote its departure. In contrast, NεH in histidine H+-440 is not aligned for a favourable proton transfer to the pinacolyl O to promote dealkylation, but electrostatic stabilization by histidine H+-440 of the developing anion on the phosphonate monoester occurs. Destabilizing interactions between residues and the alkyl fragment of the inhibitor enforce methyl migration from Cβ to Cα concerted with C—O bond breaking in soman-inhibited AChE. Tryptophan-84, phenyalanine-331 and glutamic acid-199 are within 3.7–3.9 Å (1 Å=10-10 m) from a methyl group in Cβ, 4.5–5.1 Å from Cβ and 4.8–5.8 Å from Cα, and can better stabilize the developing carbenium ion on Cβ than on Cα. The Trp84Ala mutation eliminates interactions between the incipient carbenium ion and the indole ring, but also reduces its interactions with phenylalanine-331 and aspartic acid-72. Tyrosine-130 promotes dealkylation by interacting with the indole ring of tryptophan-84. Glutamic acid-443 can influence the orientation of active-site residues through tyrosine-421, tyrosine-442 and histidine-440 in soman-inhibited AChE, and thus facilitate dealkylation.

2013 ◽  
Vol 12 (08) ◽  
pp. 1341002 ◽  
Author(s):  
XIN ZHANG ◽  
MING LEI

The deamination process of isoxanthopterin catalyzed by isoxanthopterin deaminase was determined using the combined QM(PM3)/MM molecular dynamics simulations. In this paper, the updated PM3 parameters were employed for zinc ions and the initial model was built up based on the crystal structure. Proton transfer and following steps have been investigated in two paths: Asp336 and His285 serve as the proton shuttle, respectively. Our simulations showed that His285 is more effective than Aap336 in proton transfer for deamination of isoxanthopterin. As hydrogen bonds between the substrate and surrounding residues play a key role in nucleophilic attack, we suggested mutating Thr195 to glutamic acid, which could enhance the hydrogen bonds and help isoxanthopterin get close to the active site. The simulations which change the substrate to pterin 6-carboxylate also performed for comparison. Our results provide reference for understanding of the mechanism of deaminase and for enhancing the deamination rate of isoxanthopterin deaminase.


2018 ◽  
Vol 115 (52) ◽  
pp. E12192-E12200 ◽  
Author(s):  
Haoran Yu ◽  
Paul A. Dalby

The directed evolution of enzymes for improved activity or substrate specificity commonly leads to a trade-off in stability. We have identified an activity–stability trade-off and a loss in unfolding cooperativity for a variant (3M) of Escherichia coli transketolase (TK) engineered to accept aromatic substrates. Molecular dynamics simulations of 3M revealed increased flexibility in several interconnected active-site regions that also form part of the dimer interface. Mutating the newly flexible active-site residues to regain stability risked losing the new activity. We hypothesized that stabilizing mutations could be targeted to residues outside of the active site, whose dynamics were correlated with the newly flexible active-site residues. We previously stabilized WT TK by targeting mutations to highly flexible regions. These regions were much less flexible in 3M and would not have been selected a priori as targets using the same strategy based on flexibility alone. However, their dynamics were highly correlated with the newly flexible active-site regions of 3M. Introducing the previous mutations into 3M reestablished the WT level of stability and unfolding cooperativity, giving a 10.8-fold improved half-life at 55 °C, and increased midpoint and aggregation onset temperatures by 3 °C and 4.3 °C, respectively. Even the activity toward aromatic aldehydes increased up to threefold. Molecular dynamics simulations confirmed that the mutations rigidified the active-site via the correlated network. This work provides insights into the impact of rigidifying mutations within highly correlated dynamic networks that could also be useful for developing improved computational protein engineering strategies.


2021 ◽  
Vol 14 ◽  
pp. 117864692110529
Author(s):  
Manon Mirgaux ◽  
Laurence Leherte ◽  
Johan Wouters

Protein dynamics governs most of the fundamental processes in the human body. Particularly, the dynamics of loops located near an active site can be involved in the positioning of the substrate and the reaction mechanism. The understanding of the functioning of dynamic loops is therefore a challenge, and often requires the use of a multi-disciplinary approach mixing, for example, crystallographic experiments and molecular dynamics simulations. In the present work, the dynamic behavior of the JK-loop of the human indoleamine 2,3-dioxygenase 1 hemoprotein, a target for immunotherapy, is investigated. To overcome the lack of knowledge on this dynamism, the study reported here is based on 3 crystal structures presenting different conformations of the loop, completed with molecular dynamics trajectories and MM-GBSA analyses, in order to trace the reaction pathway of the enzyme. In addition, the crystal structures identify an exo site in the small unit of the enzyme, that is populated redundantly by the substrate or the product of the reaction. The role of this newer reported exo site still needs to be investigated.


2017 ◽  
Vol 15 (43) ◽  
pp. 9095-9107 ◽  
Author(s):  
Fernanda Mendoza ◽  
José M. Lluch ◽  
Laura Masgrau

QM(DFT)/MM calculations and molecular dynamics simulations on wild-type retaining α1,4-N-acetylhexosaminyltransferase (EXTL2) and Arg293Ala, Asp246Ala, Arg293Ala/Asp246Ala and Asp246Glu mutants are used to understand the role of these two residues.


Soft Matter ◽  
2021 ◽  
Vol 17 (10) ◽  
pp. 2942-2956
Author(s):  
Rishabh D. Guha ◽  
Ogheneovo Idolor ◽  
Katherine Berkowitz ◽  
Melissa Pasquinelli ◽  
Landon R. Grace

We investigated the effect of temperature variation on the secondary bonding interactions between absorbed moisture and epoxies with different morphologies using molecular dynamics simulations.


Sign in / Sign up

Export Citation Format

Share Document