血管紧张素转换酶与抑制肽结合模式的分子动力学研究

应用分子模拟方法研究了血管紧张素转换酶(Angiotensin-converting enzyme,ACE)C端结构域(C-domain)与两种抑制肽(RIGLF/AHEPVK)的结合机制,预测了两个体系的结合模式,提出在C-domain-RIGLF中His353,Asp377,Asp453,Phe457,His513,Tyr523和Phe527为RIGLF主要结合残基,而在C-domainAHEPVK中Gln281,His353,Ser355,Glu384,Lys511,His513和Tyr523等残基起关键作用.应用结合自由能计算比较了两个体系的结合能力,结果表明,RIGLF和AHEPVK均与C-domain活性位点残基存在较强作用,且AHEPVK对C-domain的结合能力较强,与实验结果一致.     

MC1R Gene Polymorphism Affects Skin Color and Phenotypic Features Related to Sun Sensitivity in a Population of French Adult Women

The melanocortin-1 receptor ( MC1R ) gene is known to play a major role in skin and hair pigmentation and to be highly polymorphic in Caucasians. This study was performed to investigate the relationships between MC1R gene polymorphisms and skin color in a large sample of French middle-aged Caucasian women. The codons 60 to 265 and the codon 294 of the MC1R gene were sequenced in 488 women. The skin color was measured on the inner side of the forearm using a spectrophotometric instrument. Fifteen variants were identified: Arg151Cys, Arg160Trp, Arg142His, Asp294His, Ile155Thr, Asp84Glu, Val60Leu, Val92Met, Arg163Gln, Ser83Pro, Thr95Met, Pro256Ser, Val265Ile, Ala166Ala and Gln233Gln. Women carrying Arg151Cys, Asp294His, Arg160Trp and Asp84Glu variants had a significantly higher reflectance in the red region, which indicates a lower level of functional melanin. This association was the most pronounced for women carrying Asp84Glu. In contrast, no significant difference was observed for other variants. Moreover, associations between MC1R polymorphisms and the risks of experiencing sunburn and of having freckles were found independently of skin color. Our findings support the hypothesis that MC1R polymorphisms do not necessarily alter the skin color but should sensitize the skin to UV-induced DNA damage.     

Synthesis and screening of libraries of synthetic tripodal receptor molecules with three different amino acid or peptide arms: identification of iron binders

A novel selectively deprotectable triazacyclophane scaffold was used for the design and split-mix synthesis of two libraries of solid-phase bound tripodal synthetic receptors possessing three different amino acid or peptidic arms. In the synthesis of the first library, the two outer arms consisted of amino acid Ala, Arg, Asp, Gln, Gly, Lys, Phe, Ser, Tyr, or Val and the middle arm consisted of amino acid Asn, Glu, His, Leu, or Pro. The second library contained amino acid and/or (di)peptide arms. The arms were different in all library members. The first outer arm consisted of amino acid(s) Ala, Arg, Gln, Phe, or Ser, the second outer arm consisted of amino acid(s) Asp, Gly, Lys, Tyr, or Val, and the middle arm consisted of amino acid(s) Asn, Glu, His, Leu, or Pro, leading to a 27 000 member library of synthetic tripodal receptor molecules. In on-bead screening experiments, a remarkable selectivity of some library members for Fe(3+) was observed and decoding of their structures by Edman degradation revealed consensus sequences with structural resemblance to non-heme iron proteins.     

Theoretical Studies on the Binding Mode and Reaction Mechanism of TLP Hydrolase kpHIUH

In this work, we have investigated the binding conformations of the substrate in the active site of 5-HIU hydrolase kpHIUH and its catalytic hydrolysis mechanism. Docking calculations revealed that the substrate adopts a conformation in the active site with its molecular plane laying parallel to the binding interface of the protein dimer of kpHIUH, in which His7 and His92 are located adjacent to the hydrolysis site C6 and have hydrogen bond interactions with the lytic water. Based on this binding conformation, density functional theory calculations indicated that the optimal catalytic mechanism consists of two stages: (1) the lytic water molecule is deprotonated by His92 and carries out nucleophilic attack on C6=O of 5-HIU, resulting in an oxyanion intermediate; (2) by accepting a proton transferred from His92, C6–N5 bond is cleaved to completes the catalytic cycle. The roles of His7, His92, Ser108 and Arg49 in the catalytic reaction were revealed and discussed in detail.     

Amino acids of ricin and its polypeptides

Ricin and its corresponding polypeptides (A & B chain) were purified from castor seed. The molecular weight of ricin subunits were 29,000 and 28,000 daltons. The amino acids in ricin determined were Asp45 The22 Ser40 Glu53 Cys4 Gly96 His5 Ile21 Leu33 Lys20 Met4 Phe13 Pro37 Tyr11 Ala45 Val23 Arg20 indicating that ricin contains approximately 516 amino acid residues. The amino acids of the two subunits of ricin A and B chains were Asp23 The12 Ser21 Glu29 Cys2 Gly48 His3 Ile12, Leu17 Lys10 Met2 Phe6 Pro17 Tyr7 Ala35 Val13 Arg13 while in B chain the amino acids were Asp22 The10 Ser19 Glu25 Cys2 Gly47 His1 Ile10, Leu15 Lys11 Met1 Phe7 Pro6 Tyr5 Ala32Val11 Arg10. The total helical content of ricin came around 53.6% which is a new observation.     

Point Mutations of Nicotinic Receptor α1 Subunit Reveal New Molecular Features of G153S Slow-Channel Myasthenia

Slow-channel congenital myasthenic syndromes (SCCMSs) are rare genetic diseases caused by mutations in muscle nicotinic acetylcholine receptor (nAChR) subunits. Most of the known SCCMS-associated mutations localize at the transmembrane region near the ion pore. Only two SCCMS point mutations are at the extracellular domains near the acetylcholine binding site, α1(G153S) being one of them. In this work, a combination of molecular dynamics, targeted mutagenesis, fluorescent Ca²⁺ imaging and patch-clamp electrophysiology has been applied to G153S mutant muscle nAChR to investigate the role of hydrogen bonds formed by Ser 153 with C-loop residues near the acetylcholine-binding site. Introduction of L199T mutation to the C-loop in the vicinity of Ser 153 changed hydrogen bonds distribution, decreased acetylcholine potency (EC₅₀ 2607 vs. 146 nM) of the double mutant and decay kinetics of acetylcholine-evoked cytoplasmic Ca²⁺ rise (τ 14.2 ± 0.3 vs. 34.0 ± 0.4 s). These results shed light on molecular mechanisms of nAChR activation-desensitization and on the involvement of such mechanisms in channelopathy genesis.     

Intricate role of water in proton transport through cytochrome c oxidase

Cytochrome c oxidase (CytcO), the final electron acceptor in the respiratory chain, catalyzes the reduction of O(2) to H(2)O while simultaneously pumping protons across the inner mitochondrial or bacterial membrane to maintain a transmembrane electrochemical gradient that drives, for example, ATP synthesis. In this work mutations that were predicted to alter proton translocation and enzyme activity in preliminary computational studies are characterized with extensive experimental and computational analysis. The mutations were introduced in the D pathway, one of two proton-uptake pathways, in CytcO from Rhodobacter sphaeroides . Serine residues 200 and 201, which are hydrogen-bonded to crystallographically resolved water molecules halfway up the D pathway, were replaced by more bulky hydrophobic residues (Ser200Ile, Ser200Val/Ser201Val, and Ser200Val/Ser201Tyr) to query the effects of changing the local structure on enzyme activity as well as proton uptake, release, and intermediate transitions. In addition, the effects of these mutations on internal proton transfer were investigated by blocking proton uptake at the pathway entrance (Asp132Asn replacement in addition to the above-mentioned mutations). Even though the overall activities of all mutant CytcO's were lowered, both the Ser200Ile and Ser200Val/Ser201Val variants maintained the ability to pump protons. The lowered activities were shown to be due to slowed oxidation kinetics during the P(R) → F and F → O transitions (P(R) is the "peroxy" intermediate formed at the catalytic site upon reaction of the four-electron-reduced CytcO with O(2), F is the oxoferryl intermediate, and O is the fully oxidized CytcO). Furthermore, the P(R) → F transition is shown to be essentially pH independent up to pH 12 (i.e., the apparent pK(a) of Glu286 is increased from 9.4 by at least 3 pK(a) units) in the Ser200Val/Ser201Val mutant. Explicit simulations of proton transport in the mutated enzymes revealed that the solvation dynamics can cause intriguing energetic consequences and hence provide mechanistic insights that would never be detected in static structures or simulations of the system with fixed protonation states (i.e., lacking explicit proton transport). The results are discussed in terms of the proton-pumping mechanism of CytcO.     

Structural Insights into Novel 15-Prostaglandin Dehydrogenase Inhibitors

We discovered SW033291 in a high throughput chemical screen aimed at identifying 15-prostaglandin dehydrogenase (15-PGDH) modulators. The compound exhibited inhibitory activity in in vitro biochemical and cell-based assays of 15-PGDH activity. We subsequently demonstrated that this compound, and several analogs thereof, are effective in in vivo mouse models of bone marrow transplant, colitis, and liver regeneration, where increased levels of PGE2 positively potentiate tissue regeneration. To better understand the binding of SW033291, we carried out docking studies for both the substrate, PGE2, and an inhibitor, SW033291, to 15-PGDH. Our models suggest similarities in the ways that PGE2 and SW033291 interact with key residues in the 15-PGDH-NAD+ complex. We carried out molecular dynamics simulations (MD) of SW033291 bound to this complex, in order to understand the dynamics of the binding interactions for this compound. The butyl side chain (including the sulfoxide) of SW033291 participates in crucial binding interactions that are similar to those observed for the C15-OH and the C16-C20 alkyl chain of PGE2. In addition, interactions with residues Ser138, Tyr151, and Gln148 play key roles in orienting and stabilizing SW033291 in the binding site and lead to enantioselectivity for the R-enantiomer. Finally, we compare the binding mode of (R)-S(O)-SW033291 with the binding interactions of published 15-PGDH inhibitors.     

Synthesis,Structural Characterization,and In Vitro and In Silico Antifungal Evaluation of Azo-Azomethine Pyrazoles (PhN2(PhOH)CHN(C3N2(CH3)3)PhR,R = H or NO2)

The azo-azomethine imines, R¹-N=N-R²-CH=N-R³, are a class of active pharmacological ligands that have been prominent antifungal, antibacterial, and antitumor agents. In this study, four new azo-azomethines, R¹ = Ph, R² = phenol, and R³ = pyrazol-Ph-R’ (R = H or NO₂), have been synthesized, structurally characterized using X-ray, IR, NMR and UV–Vis techniques, and their antifungal activity evaluated against certified strains of Candida albicans and Cryptococcus neoformans. The antifungal tests revealed a high to moderate inhibitory activity towards both strains, which is regulated as a function of both the presence and the location of the nitro group in the aromatic ring of the series. These biological assays were further complemented with molecular docking studies against three different molecular targets from each fungus strain. Molecular dynamics simulations and binding free energy calculations were performed on the two best molecular docking results for each fungus strain. Better affinity for active sites for nitro compounds at the “meta” and “para” positions was found, making them promising building blocks for the development of new Schiff bases with high antifungal activity.     

Homology modeling and molecular dynamic simulation of UDP-N-acetylmuramoyl-l-alanine-d-glutamate ligase (MurD) from Mycobacterium tuberculosis H37Rv using in silico approach

The present study aimed to identify the prospective inhibitors of MurD, a cytoplasmic enzyme that catalyzes the addition of d-glutamate to the UDP-N-acetylmuramoyl-l-alanine nucleotide precursor in Mycobacterium tuberculosis (MTB), using virtual screening, docking studies, pharmacokinetic analysis, Molecular Dynamic (MD) simulation, and Molecular Mechanics Generalized Born and Surface Area (MM-GBSA) analyses. The three dimensional (3D) structure was determined based on the homology technique using a template from Streptococcus agalactiae. The modeled structure had three binding sites, namely; substrate binding site (Val18, Thr19, Asp39, Asp40, Gly75, Asn147, Gln171 and His192), the ATP binding site (Gly123, Lys124, Thr125, Thr126, Glu166, Asp283, and Arg314) and the glutamic acid binding site (Arg382, Ser463, and Tyr470). These residues mentioned above play a critical role in the catalytic activity of the enzyme, and their inhibition could serve as a stumbling block to the normal function of the enzyme. A total of 10,344 obtained from virtual screened of Zinc and PubChem databases. These compounds further screened for Lipinski rule of five, docking studies and pharmacokinetic analysis. Four compounds with good binding energies (ZINC11881196 = −10.33 kcal/mol, ZINC12247644 = −8.90 kcal/mol, ZINC14995379 =−8.42 kcal/mol, and PubChem6185 = −8.20 kcal/mol), better than the binding energies of the ATP (−2.31 kcal/mol) and the ligand with known IC₅₀, Aminothiazole (−7.11 kcal/mol) were selected for the MD simulation and MM-GBSA analyses. The result of the analyses showed that all the four ligands formed a stable complex and had the binding free energies better than the binding energy of ATP. Therefore, these ligands considered as suitable prospective inhibitors of the MurD after experimental validation.     

Effects of mutation on the amyloidogenic propensity of apolipoprotein C-II(60-70) peptide

Using experimental and computational methods we identified the effects of mutation on the structure and dynamics of the amyloidogenic peptide apoC-II(60-70), in monomeric and oligomeric states. Methionine (Met60) substitutions to hydrophilic Gln, hydrophobic Val, and methionine sulfoxide residues were investigated and the results compared with observations of fibril formation by the wild-type, Met60Gln, Met60Val, and oxidised Met60 (oxi-Met) apoC-II(60-70) peptides. ThT fluorescence measurements showed fibril formation by all peptides, however with different kinetics. The wild-type and Met60Val peptides formed fibrils fastest, while oxi-Met and Met60Gln peptides exhibited significantly longer lag phases. Molecular dynamics simulations showed that the mutated monomers exhibited structural features consistent with fibril-forming propensity, such as β-hairpin conformation and a hydrophobic core. However, important differences to the wild-type were also noted, such as increased structural flexibility (oxi-Met and Met60Gln systems) and a broader distribution of the aromatic angle orientation, which could contribute to the different fibrillation kinetics observed in these peptides. Our results also showed that the critical nucleus size for fibril formation by apoC-II(60-70) may not be very large, since tetrameric oligomers in anti-parallel configuration were very stable within the 100 ns of simulations. The single-point mutations Met60Val and Met60Gln had no significant effect on the structural stability of the tetramer. The rate of fibril formation by apoC-II(60-70) peptides was generally much faster compared to longer apoC-II(56-76) peptides. Also, the effects of amino acid modifications on the kinetics of peptide fibril formation differ from the effects observed for apoC-II(56-76) and full-length apoC-II, suggesting that additional mechanisms are involved in fibril formation by mature apoC-II.     

Structure-Activity relationship of MDSA and its derivatives against Staphylococcus aureus Ser/Thr phosphatase Stp1

Ser/thr phosphatase Stp1 is an important virulence factor for Staphylococcus aureus (S. aureus) and plays a key role in its infectivity, suggesting that it could serve as a potential target for treatment of S. aureus infection. Previous studies found that the activity of Stp1 was inhibited by MDSA and its derivatives. In this paper, we used molecular docking, molecular modeling, molecular dynamics simulations, binding free energy decomposition calculations, and hydrogen bond analyses to explore the structure-activity relationship. Energy decomposition indicated that MDSA, hydroxymethyl MDSA, carboxymethyl MDSA and methyl MDSA can bind to the catalytic pocket of Stp1. Furthermore, Met39, Ile163, Ile164, Val167, Gly195 and Asp233 were key residues in the Stp1-inhibitor complexes. Due to the lack of a double salicylate structure, salicylic acid cannot bind to the active site of Stpl, leading to loss of inhibitory activity. Based on these results, the structure-activity relationship at the atomic level was determined, which can promote the development of new and more effective anti-drug resistance inhibitors.     

AT1受体的中药活性成分筛选模型及其作用机理研究

传统中药对治疗心血管类疾病疗效显著,例如钩藤、黄芪、益母草等在临床应用广泛.现代药理研究表明钩藤碱可以降压;黄芪中毛蕊异黄酮能舒张血管平滑肌、保护心脑血管;益母草碱可扩张微血管,改善血液流变异常,但它们分子层面作用机制尚不明确.首先以牛视紫红质蛋白为模板,模建出心血管疾病主要靶点AT1受体的三维结构.然后将AT1受体拮抗剂和中药活性成分与受体模型结合的作用方式进行了对比研究,据此提出了中药活性成分治疗心血管疾病的作用机理,并建立了AT1受体的中药活性成分筛选模型.结果表明:黄芪毛蕊异黄酮等中药活性成分能与AT1受体活性口袋中的残基发生氢键作用,结合方式与AT1受体拮抗剂相似.每一种AT1受体拮抗剂均与His183,Lys199,His256,Gln257,Ser105,Ser109,Tyr113,Asn200中多个发生氢键作用;黄芪毛蕊异黄酮与Try113,Lys199,Gln257,Ser105发生氢键作用.本研究从分子层面上阐释了一些中药活性小分子的治病机理,为进一步挖掘中药资源,研究AT1受体相关的心脑血管类药物,合理设计和筛选AT1受体的拮抗剂提供重要依据.     

Investigation of the interactions between ginsenosides and amino acids by mass spectrometry and theoretical chemistry

In order to evaluate the essence of the interactions of ginsenosides and proteins which are composed by α-amino acids, electrospray ionization mass spectrometry was employed to study the noncovalent interactions between ginsenosides (Rb₂, Rb₃, Re, Rg₁ and Rh₁) and 18 kinds of α-amino acids (Asp, Glu, Asn, Phe, Gln, Thr, Ser, Met, Trp, Val, Gly, Ile, Ala, Leu, Pro, His, Lys and Arg). The 1:1 and 2:1 noncovalent complexes of ginsenosides and amino acids were observed in the mass spectra. The dissociation constants for the noncovalent complexes were directly calculated based on peak intensities of ginsenosides and the noncovalent complexes in the mass spectra. Based on the dissociation constants, it can be concluded that the acidic and the basic amino acids, Asp, Glu, Lys and Arg, bound to ginsenosides more strongly than other amino acids. The experimental results were verified by theoretical calculations of parameters of noncovalent interaction between ginsenoside Re and Arg which served as a representative example. Two kinds of binding forms, “head–tail” (“H–T”) and “head–head” (“H–H”), were proposed to explain the interaction between ginsenosides and amino acids. And the interaction in “H–T” form was stronger than that in “H–H” form.     

Isolation of Chalcomoracin as a Potential α-Glycosidase Inhibitor from Mulberry Leaves and Its Binding Mechanism

Diabetes is a chronic metabolic disease, whereas α-glucosidases are key enzymes involved in the metabolism of starch and glycogen. There is a long history of the use of mulberry leaf (the leaf of Morus alba) as an antidiabetic herb in China, and we found that chalcomoracin, one of the specific Diels–Alder adducts in mulberry leaf, had prominent α-glucosidase inhibitory activity and has the potential to be a substitute for current hypoglycemic drugs such as acarbose, which have severe gastrointestinal side effects. In this study, chalcomoracin was effectively isolated from mulberry leaves, and its α-glucosidase inhibition was studied via enzymatic kinetics, isothermal titration (ITC) and molecular docking. The results showed that chalcomoracin inhibited α-glucosidase through both competitive and non-competitive manners, and its inhibitory activity was stronger than that of 1-doxymycin (1-DNJ) but slightly weaker than that of acarbose. ITC analysis revealed that the combination of chalcomoracin and α-glucosidase was an entropy-driven spontaneous reaction, and the molecular docking results also verified this conclusion. During the binding process, chalcomoracin went into the “pocket” of α-glucosidase via hydrophobic interactions, and it is linked with residues Val544, Asp95, Ala93, Gly119, Arg275 and Pro287 by hydrogen bonds. This study provided a potential compound for the prevention and treatment of diabetes and a theoretical basis for the discovery of novel candidates for α-glycosidase inhibitors.     

Anomalous structural dynamics of minimally frustrated residues in cardiac troponin C triggers hypertrophic cardiomyopathy

Cardiac TnC (cTnC) is highly conserved among mammals, and genetic variants can result in disease by perturbing Ca²⁺-regulation of myocardial contraction. Here, we report the molecular basis of a human mutation in cTnC''s αD-helix (TNNC1-p.C84Y) that impacts conformational dynamics of the D/E central-linker and sampling of discrete states in the N-domain, favoring the “primed” state associated with Ca²⁺ binding. We demonstrate cTnC''s αD-helix normally functions as a central hub that controls minimally frustrated interactions, maintaining evolutionarily conserved rigidity of the N-domain. αD-helix perturbation remotely alters conformational dynamics of the N-domain, compromising its structural rigidity. Transgenic mice carrying this cTnC mutation exhibit altered dynamics of sarcomere function and hypertrophic cardiomyopathy. Together, our data suggest that disruption of evolutionary conserved molecular frustration networks by a myofilament protein mutation may ultimately compromise contractile performance and trigger hypertrophic cardiomyopathy.

Cardiac TnC (cTnC) is highly conserved among mammals, and genetic variants can result in disease by perturbing Ca²⁺-regulation of myocardial contraction.     

谷氨酰胺结合蛋白的分子动力学模拟和自由能计算

谷氨酰胺结合蛋白(Glutamine-binding protein, GlnBp)是大肠杆菌透性酶系统中一个细胞外液底物专一性结合蛋白, 对于细胞外液中谷氨酰胺(Gln)的运输和传递至关重要. 本文运用分子动力学(Molecular dynamics, MD)模拟采样, 考察了GlnBp关键残基与底物Gln之间的相互作用和GlnBp两条铰链的功能差别; 并采用MM-PBSA方法计算了GlnBp与底物Gln的结合自由能. 结果表明: Ph13, Phe50, Thr118和Ile69与底物Gln的范德华相互作用和Arg75, Thr70, Asp157, Gly68, Lys115, Ala67, His156与底物Gln的静电相互作用是结合Gln的主要推动力; 复合物的铰链区85～89柔性大, 对构象开合提供了结构基础; 而铰链区181～185柔性小, 其作用更多是在功能上把底物Gln限制在口袋中; 自由能预测值与实验值吻合. 本研究很好地解释了GlnBp结构与功能的关系, 为进一步了解GlnBp的开合及转运Gln的机制提供了重要的结构信息.     

Insights into the Inhibitory Mechanism of Viniferifuran on Xanthine Oxidase by Multiple Spectroscopic Techniques and Molecular Docking

Viniferifuran was investigated for its potential to inhibit the activity of xanthine oxidase (XO), a key enzyme catalyzing xanthine to uric acid. An enzyme kinetics analysis showed that viniferifuran possessed a strong inhibition on XO in a typical anti-competitive manner with an IC₅₀ value of 12.32 μM (IC₅₀ for the first-line clinical drug allopurinol: 29.72 μM). FT-IR and CD data analyses showed that viniferifuran could induce a conformational change of XO with a decrease in the α-helix and increases in the β-sheet, β-turn, and random coil structures. A molecular docking analysis revealed that viniferifuran bound to the amino acid residues located within the activity cavity of XO by a strong hydrophobic interaction (for Ser1214, Val1011, Phe914, Phe1009, Leu1014, and Phe649) and hydrogen bonding (for Asn768, Ser876, and Tyr735). These findings suggested that viniferifuran might be a promising XO inhibitor with a favorable mechanism of action.     

Insights into the Allosteric Effect of SENP1 Q597A Mutation on the Hydrolytic Reaction of SUMO1 via an Integrated Computational Study

Small ubiquitin-related modifier (SUMO)-specific protease 1 (SENP1) is a cysteine protease that catalyzes the cleavage of the C-terminus of SUMO1 for the processing of SUMO precursors and deSUMOylation of target proteins. SENP1 is considered to be a promising target for the treatment of hepatocellular carcinoma (HCC) and prostate cancer. SENP1 Gln597 is located at the unstructured loop connecting the helices α4 to α5. The Q597A mutation of SENP1 allosterically disrupts the hydrolytic reaction of SUMO1 through an unknown mechanism. Here, extensive multiple replicates of microsecond molecular dynamics (MD) simulations, coupled with principal component analysis, dynamic cross-correlation analysis, community network analysis, and binding free energy calculations, were performed to elucidate the detailed mechanism. Our MD simulations showed that the Q597A mutation induced marked dynamic conformational changes in SENP1, especially in the unstructured loop connecting the helices α4 to α5 which the mutation site occupies. Moreover, the Q597A mutation caused conformational changes to catalytic Cys603 and His533 at the active site, which might impair the catalytic activity of SENP1 in processing SUMO1. Moreover, binding free energy calculations revealed that the Q597A mutation had a minor effect on the binding affinity of SUMO1 to SENP1. Together, these results may broaden our understanding of the allosteric modulation of the SENP1−SUMO1 complex.     

Two Sides of Quantum-Based Modeling of Enzyme-Catalyzed Reactions: Mechanistic and Electronic Structure Aspects of the Hydrolysis by Glutamate Carboxypeptidase

We report the results of a computational study of the hydrolysis reaction mechanism of N-acetyl-l-aspartyl-l-glutamate (NAAG) catalyzed by glutamate carboxypeptidase II. Analysis of both mechanistic and electronic structure aspects of this multistep reaction is in the focus of this work. In these simulations, model systems are constructed using the relevant crystal structure of the mutated inactive enzyme. After selection of reaction coordinates, the Gibbs energy profiles of elementary steps of the reaction are computed using molecular dynamics simulations with ab initio type QM/MM potentials (QM/MM MD). Energies and forces in the large QM subsystem are estimated in the DFT(PBE0-D3/6-31G**) approximation. The established mechanism includes four elementary steps with the activation energy barriers not exceeding 7 kcal/mol. The models explain the role of point mutations in the enzyme observed in the experimental kinetic studies; namely, the Tyr552Ile substitution disturbs the “oxyanion hole”, and the Glu424Gln replacement increases the distance of the nucleophilic attack. Both issues diminish the substrate activation in the enzyme active site. To quantify the substrate activation, we apply the QTAIM-based approaches and the NBO analysis of dynamic features of the corresponding enzyme-substrate complexes. Analysis of the 2D Laplacian of electron density maps allows one to define structures with the electron density deconcentration on the substrate carbon atom, i.e., at the electrophilic site of reactants. The similar electronic structure element in the NBO approach is a lone vacancy on the carbonyl carbon atom in the reactive species. The electronic structure patterns revealed in the NBO and QTAIM-based analyses consistently clarify the reactivity issues in this system.