3T MEGA-PRESS study of N-acetyl aspartyl glutamate and N-acetyl aspartate in activated visual cortex.

OBJECTIVE: To measure N-acetyl aspartyl glutamate (NAAG) and N-acetyl aspartate (NAA) concentrations in visual cortex activated by a continuous stimulation in a 3 T field. METHODS: NAAG and NAA spectra were obtained with MEGA-PRESS pulse sequence (TE/TR = 140/2000 ms; δONNAAG/δOFFNAAG = 4.61/4.15 ppm; δONNAA/δOFFNAA = 4.84/4.38 ppm) in 14 healthy volunteers at rest and upon stimulation by a radial checkerboard flickering at a frequency of 8 Hz. Spectra of all subjects were frequency and phase aligned and then averaged. Additionally, to obtain the time-dependency data, spectra were divided into time sections of 64 s each. The intensities of NAA, NAAG and lactate + macromolecular (Lac + MM) signals were defined by integration of the real part of spectra. The heights of the central resonance of NAAG and NAA signals were measured. RESULTS: The NAAG and NAA concentrations, measured with 2.5% and 0.5% error, respectively, were unaffected by visual activation. A significant increase in the Lac + MM signal by ~ 12% is clearly observed. No stimulation-induced time dependency was found for NAAG or NAA, while the increase in Lac + MM was gradual. The concentration values in visual cortex are in good agreement with the 7 T MRS measurements: [NAAG] = 1.55 mM, [NAA] = 11.95 mM. CONCLUSION: The MEGA-PRESS pulse sequence together with the spectral preprocessing techniques allowed to demonstrate that the concentrations of NAAG and NAA in the visual cortex remain constant during continuous visual stimulation within the margin of error. An increase in the lactate signal intensity signifies the activation of the anaerobic glycolysis in activated visual cortex.

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 (kcat/Km > 106 M-1min-1) 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.

Understanding the non-catalytic behavior of human butyrylcholinesterase silent variants: Comparison of wild-type enzyme, catalytically active Ala328Cys mutant, and silent Ala328Asp variant.

Conformational dynamics of wild-type human butyrylcholinesterase (BChE), two mutants of residue Ala328, the catalytically active Ala328Cys, and the catalytically inactive (silent) Ala328Asp, and their interactions with butyrylcholine were studied. The aim was to understand the molecular mechanisms by which point mutations may lead to silent BChE variant or alter catalytic activity. Importance of BChE natural variants is due to medical consequences, i.e. prolonged apnea, following administration of the myorelaxant esters, succinylcholine and mivacurium. Comparison of molecular dynamics (MD) simulations for the three model systems showed that: 1) the active mutant Ala328Cys mutant has some changes in configuration of catalytic residues, which do not prevent binding of butyrylcholine to the active site; 2) in the naturally-occurring silent variant Ala328Asp, the Asp328 carboxylate may either form a salt bridge with Lys339 or a H-bond with His438. In the first case, the Ω-loop swings off the gorge, disrupting the π-cation binding site and the catalytic triad. In the second case, binding of cationic substrates in the catalytic center is also impaired. MD simulations carried out in 0.15 M NaCl, close to physiological ionic strength conditions, favored the second situation. It was seen that Asp328 forms a H-bond with the catalytic triad His438, which in turn disrupts the catalytic machinery. Therefore, we concluded that the Ala328Asp variant is not catalytically active because of that dramatic event. Computational results, consistent with in vitro biochemical data and clinical observations, validate our MD approach.

Modeling the Complete Catalytic Cycle of Aspartoacylase.

The complete catalytic cycle of aspartoacylase (ASPA), a zinc-dependent enzyme responsible for cleavage of N-acetyl-l-aspartate, is characterized by the methods of molecular modeling. The reaction energy profile connecting the enzyme-substrate (ES) and the enzyme-product (EP) complexes is constructed by the quantum mechanics/molecular mechanics (QM/MM) method assisted by the molecular dynamics (MD) simulations with the QM/MM potentials. Starting from the crystal structure of ASPA complexed with the intermediate analogue, the minimum-energy geometry configurations and the corresponding transition states are located. The stages of substrate binding to the enzyme active site and release of the products are modeled by MD calculations with the replica-exchange umbrella sampling technique. It is shown that the first reaction steps, nucleophilic attack of a zinc-bound nucleophilic water molecule at the carbonyl carbon and the amide bond cleavage, are consistent with the glutamate-assisted mechanism hypothesized for the zinc-dependent hydrolases. The stages of formation of the products, acetate and l-aspartate, and regeneration of the enzyme are characterized for the first time. The constructed free energy diagram from the reactants to the products suggests that the enzyme regeneration, but not the nucleophilic attack of the catalytic water molecule, corresponds to the rate-determining stage of the full catalytic cycle of ASPA.

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.

Research on cholinesterases in the Soviet Union and Russia: a historical perspective.

Research on cholinesterases and effects of their inhibition in the USSR and Russia since 1930-1940s till present is exposed in historical aspects. The first physiological and toxicological effects of cholinesterase inhibition were reported by Alexander Ginetsinsky during World War II, when academic institutions were evacuated from Leningrad to Kazan. The main scientific schools that initiated research on chemistry, enzymology and physiology of cholinesterases and their inhibitors were leaded by Alexandr and Boris Arbuzovs, Victor Rozengart, Viktor Yakovlev, Michael Michelson, Martin Kabachnik, Mikhail Voronkov, Ivan Knunyants, Alexandr Bretskin and others. They investigated the main physiological effects of cholinesterase inhibitors, and analyzed the catalytic mechanisms of cholinesterases and related enzymes. Their contributions are landmarks in the history of cholinesterase research. At the present time revival of research on cholinesterases in different universities and institutes is vivid, in particular at the Moscow State University, research institutes of Russian Academy of Sciences and Kazan Scientific Center.

On quantum mechanical--molecular mechanical (QM/MM) approaches to model hydrolysis of acetylcholine by acetylcholinesterase.

We re-visited the results of quantum mechanics--molecular mechanics (QM/MM) approaches aiming to construct the reaction energy profile for the acylation stage of acetylcholine hydrolysis by acetylcholinesterase. The main emphasis of this study was on the energy of the first tetrahedral intermediate (TI) relative to the level of the enzyme-substrate (ES) complex for which contradictory data from different works had been reported. A new series of stationary points on the potential energy surface was calculated by using electronically embedding QM/MM schemes when starting from the crystal structure mimicking features of the reaction intermediate (PDB ID: 2VJA). A thoughtful analysis allows us to conclude that the energy of TI should be lower than that of ES, and a proper treatment of contributions from the oxyanion hole residues accounts for their relative positions.

Correlation between the substrate structure and the rate of acetylcholinesterase hydrolysis modeled with the combined quantum mechanical/molecular mechanical studies.

The combined quantum mechanical-molecular mechanical (QM/MM) based computational scheme for modeling the structure-reaction rate correlations was elaborated for the hydrolysis of the set of neutral esters in the active site of acetylcholinesterase (AChE). The energy barriers of hydrolysis were estimated on the basis of the equilibrium geometry configurations of the enzyme-substrate (ES) complexes. The obtained correlation between the rate of hydrolysis and the hydrophobicity of the substrate leaving group is consistent with experimental data. The developed method can be used to predict the substrate reactivity and to interpret the specific nature of the enzyme catalysis.

Free Energy Barriers for the N-Terminal Asparagine to Succinimide Conversion: Quantum Molecular Dynamics Simulations for the Fully Solvated Model.

Deamidation of asparagine residues represents one of the main routes for the post-translational modification of protein sequences. We computed the estimates of the free energy barriers for three stages of the deamidation process, deprotonation, cyclization, and deamination, of the conversion of asparagine to the succinimide intermediate within the fully solvated model with explicit water molecules. The Born-Oppenheimer molecular dynamics in the Gaussian and Plane Wave (GPW) approximation as implemented in the CP2K quantum chemistry package was utilized to sample the configurational space. By applying the metadynamics technique, the estimates of the free energy barriers were obtained for three separated stages of the reaction. In agreement with the experimental kinetic measurements, the estimated activation barriers do not exceed 21 kcal/mol. We demonstrate that the use of fully solvated models is the critical issue in theoretical studies of these reactions. We also conclude that more extensive sampling is necessary to obtain full free energy profiles and accurate barriers for the reaction stages.

Characterization of a complete cycle of acetylcholinesterase catalysis by ab initio QM/MM modeling.

The reaction mechanism of acetylcholine hydrolysis by acetylcholinesterase, including both acylation and deacylation stages from the enzyme-substrate (ES) to the enzyme-product (EP) molecular complexes, is examined by using an ab initio type quantum mechanical - molecular mechanical (QM/MM) approach. The density functional theory PBE0/aug-6-31+G* method for a fairly large quantum part trapped inside the native protein environment, and the AMBER force field parameters in the molecular mechanical part are employed in computations. All reaction steps, including the formation of the first tetrahedral intermediate (TI1), the acylenzyme (EA) complex, the second tetrahedral intermediate (TI2), and the EP complex, are modeled at the same theoretical level. In agreement with the experimental rate constants, the estimated activation energy barrier of the deacylation stage is slightly higher than that for the acylation phase. The critical role of the non-triad Glu202 amino acid residue in orienting lytic water molecule and in stabilizing the second tetrahedral intermediate at the deacylation stage of the enzymatic process is demonstrated.

A novel approach to distinguish between enzyme mechanisms: quasi-steady-state kinetic analysis of the prostaglandin H synthase peroxidase reaction.

A method of analysis for steady-state kinetic data has been developed that allows relationships between key partial reactions in the catalytic cycle of a functioning enzyme to be determined. The novel approach is based on a concept of scalar and vector 'kinetic connectivities' between enzyme intermediates in an arbitrary enzyme mechanism. The criterion for the agreement between experimental data and a proposed kinetic model is formulated as the kinetic connectivity of intermediate forms of the enzyme. This concept has advantages over conventional approaches and is better able to describe the complex kinetic behaviour of prostaglandin H synthase (PGHS) when catalysing the oxidation of adrenaline by H(2)O(2). To interpret the experimental data for PGHS, a generalized model for multi-substrate enzyme reactions was developed with provision for irreversible enzyme inactivation. This model showed that two enzyme intermediates must undergo inactivation during the catalytic cycle. These forms are proposed to be PGHS compound I and a compound I-adrenaline complex.