Ekaterina D. Kots

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.

Macrocyclic derivatives of 6-methyluracil as ligands of the peripheral anionic site of acetylcholinesterase

Novel pyrimidinophanes possessing two o-nitrobenzylethyldialkylammonium heads bridging with different spacers were prepared. Pyrimidinophanes 2a, 2b and 3 are reversible inhibitors of cholinesterases. They show a very good selectivity for human acetylcholinesterase (AChE), with an inhibitory power 100–200 times higher than for human butyrylcholinesterase (BChE). Docking simulations indicate specific binding of pyrimidinophanes 2a and 4 onto the peripheral anionic site of AChE. Other compounds bind to the active center of AChE as well as to the peripheral anionic site. These compounds are dual binding site inhibitors. Pyrimidinophane 2b and its acyclic counterpart 1 were tested in the animal model of myasthenia gravis and may be considered as valuable candidates for the treatment of pathological muscle weakness syndromes.

Reaction Mechanism of Guanosine Triphosphate Hydrolysis by the Vision-Related Protein Complex Arl3–RP2

Complexes of small GTPases with GTPase-activating proteins have been intensively studied with the main focus on the complex of H-Ras with p120GAP (Ras-GAP). The detailed mechanism of GTP hydrolysis is still unresolved. To clarify it, we calculated the energy profile of GTP hydrolysis in the active site of a recently characterized vision-related member of this family, the Arl3-RP2 complex. The mechanism suggested in this study retains the main features of GTP hydrolysis by the Ras-GAP complex, but the relative energies of the corresponding intermediates are different and an additional intermediate exists in the Arl3-RP2 complex compared with the Ras-GAP. These differences arise from small deviations in the catalytic arginine conformation of the active site. In the Arl3-RP2 complex, the first two intermediates, corresponding to the Pγ-Oβγ bond cleavage and the glutamine-assisted proton transfer, are almost isoenergetic with the ES complex. Numerical simulations of the kinetic curves demonstrate that the concentrations of these intermediates are comparable with that of ES during the reaction. The calculated IR spectra reveal specific vibrational bands, corresponding to these intermediates. These specific features of the Arl3-RP2 complex open the opportunity to identify spectroscopically two more reaction intermediates in GTP hydrolysis in addition to the ES and EP complexes.

Role of Protein Dimeric Interface in Allosteric Inhibition of N-Acetyl-Aspartate Hydrolysis by Human Aspartoacylase

The results of molecular modeling suggest a mechanism of allosteric inhibition upon hydrolysis of N-acetyl-aspartate (NAA), one of the most abundant amino acid derivatives in brain, by human aspartoacylase (hAsp). Details of this reaction are important to suggest the practical ways to control the enzyme activity. Search for allosteric sites using the Allosite web server and SiteMap analysis allowed us to identify substrate binding pockets located at the interface between the subunits of the hAsp dimer molecule. Molecular docking of NAA to the pointed areas at the dimer interface predicted a specific site, in which the substrate molecule interacts with the Gly237, Arg233, Glu290, and Lys292 residues. Analysis of multiple long-scaled molecular dynamics trajectories (the total simulation time exceeded 1.5 μs) showed that binding of NAA to the identified allosteric site induced significant rigidity to the protein loops with the amino acid side chains forming gates to the enzyme active site. Application of the protein dynamical network algorithms showed that substantial reorganization of the signal propagation pathways of intersubunit communication in the dimer occurred upon allosteric NAA binding to the remote site. The modeling approaches provide an explanation to the observed decrease of the reaction rate of NAA hydrolysis by hAsp at high substrate concentrations.

Three Faces of N-Acetylaspartate: Activator, Substrate and Inhibitor of Human Aspartoacylase

Hydrolysis of N-acetylaspartate (NAA), one of the most concentrated metabolites in brain, catalyzed by human aspartoacylase (hAsp) shows a remarkable dependence of the reaction rate on substrate concentration. At low NAA concentrations, sigmoidal shape of kinetic curve is observed, followed by typical rate growth of the enzyme-catalyzed reaction, whereas at high NAA concentrations self-inhibition takes place. We show that this rate dependence is consistent with a molecular model, in which N-acetylaspartate appears to have three faces in the enzyme reaction, acting as activator at low concentrations, substrate at moderate concentrations, and inhibitor at high concentrations. To support this conclusion we identify binding sites of NAA at the hAsp dimer including those on the protein surface (activating sites) and at the dimer interface (inhibiting site). Using the Markov state model approach we demonstrate that population of either activating or inhibiting site shifts the equilibrium between the hAsp dimer conformations with the open and closed gates leading to the enzyme active site buried inside the protein. These conclusions are in accord with the calculated values of binding constants of NAA at the hAsp dimer, indicating that the activating site with a higher affinity to NAA should be occupied first, whereas the inhibiting site with a lower affinity to NAA should be occupied later. Application of the dynamical network analysis shows that communication pathways between the regulatory sites (activating or inhibiting) and the gates to the active site do not interfere. These considerations allow us to develop a kinetic mechanism and to derive the equation for the reaction rate covering the entire NAA concentration range. Perfect agreement between theoretical and experimental kinetic data provides strong support to the proposed catalytic model.

Molecular polymorphism of human enzymes as the basis of individual sensitivity to drugs. Supercomputer-assisted modeling as a tool for analysis of structural changes and enzymatic activity of proteins

The nature of individual sensitivity to drugs associated with molecular polymorphism of human enzymes is discussed. The influence of molecular polymorphism on the activity of key human esterases, in particular, cholinesterases and carboxylesterase, responsible for hydrolytic metabolism of ester-containing drugs, is analyzed. A method was developed, which involves supercomputer-assisted modeling as a tool for assessment of molecular mechanism of the impact of point mutations on the catalytic activity of enzymes. This work is a part of a study aimed at elaboration of the concept and methods of personalized medicine.

Supercomputer technologies for structural-kinetic study of mechanisms of enzyme catalysis: A quantum-chemical description of aspartoacylase catalysis

The results of modeling of the complete catalytic cycle of aspartoacylase-catalyzed N-acetylaspartate hydrolysis by the combined quantum mechanics/molecular mechanics method and with the use of umbrella sampling replica-exchange molecular dynamics simulations are reported. It has been shown that the decrease in the high-energy barriers of rate-limiting stages is achieved through the preceding equilibrium stages, such as proton transfer and conformational changes. General features of the catalytic behavior of enzymes have been formulated.