Theoretical study toward understanding the catalytic mechanism of pyruvate decarboxylase

Density functional calculations are employed to explore the mechanisms of all elementary reaction steps involved in the catalytic cycle of pyruvate decarboxylase (PDC). Different models are constructed for mimicking the involvement of some key residues in a certain step. The effect of the protein framework on the potential energy profiles of active site models is approximately modeled by fixing some freedoms, based on the crystal structure of the PDC enzyme from Saccharomyces cerevisiae (ScPDC). Our calculations confirm that Glu51 is the most important residue in the formation of the ylide and the release of acetaldehyde via the proton relay between Glu51, N1', and the 4'-amino group of thiamine diphosphate. The presence of Glu477 and Asp28 residues makes the decarboxylation of lactylthiamin diphosphate (LThDP) an endothermic process with a significant free energy barrier. The protonation of the alpha-carbanion to form 2-(1-hydroxyethyl)-thiamin diphosphate is found to go through a concerted double proton transfer transition state involving both Asp28 and His115 residues. The final step, acetaldehyde release, is likely to proceed through a concerted transition state involving carbon-carbon bond-breaking and the deprotonation of the alpha-hydroxyl group. The decarboxylation of LThDP and the protonation of the alpha-carbanion are two rate-limiting steps, relative to the facile occurrence of the ylide formation and acetaldehyde release. The catalytic roles of residues Glu51, Glu477, Asp28, and Gly417 in the active site of ScPDC in individual steps elucidated from the present study are in good agreement with those derived from site-directed mutagenesis.     

Characterizing the protonation states of the catalytic residues in apo and substrate-bound human T-cell leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) protease is an attractive target when developing inhibitors to treat HTLV-1 associated diseases. To study the catalytic mechanism and design novel HTLV-1 protease inhibitors, the protonation states of the two catalytic aspartic acid residues must be determined. Free energy simulations have been conducted to study the proton transfer reaction between the catalytic residues of HTLV-1 protease using a combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulation. The free energy profiles for the reaction in the apo-enzyme and in an enzyme – substrate complex have been obtained. In the apo-enzyme, the two catalytic residues are chemically equivalent and are expected to be both unprotonated. Upon substrate binding, the catalytic residues of HTLV-1 protease evolve to a singly protonated state, in which the OD1 of Asp32 is protonated and forms a hydrogen bond with the OD1 of Asp32′, which is unprotonated. The HTLV-1 protease–substrate complex structure obtained from this simulation can serve as the Michaelis complex structure for further mechanistic studies of HTLV-1 protease while providing a receptor structure with the correct protonation states for the active site residues toward the design of novel HTLV-1 protease inhibitors through virtual screening.     

Computational determination of fundamental pathway and activation barriers for acetohydroxyacid synthase‐catalyzed condensation reactions of α‐keto acids

Acetohydroxyacid synthase (AHAS) is the first common enzyme in the biosynthetic pathway leading to the production of various branched‐chain amino acids. AHAS is recognized as a promising target for new antituberculosis drugs, antibacterial drugs, and herbicides. Extensive first‐principles quantum mechanical (QM) and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations have enabled us, in this study, to uncover the fundamental reaction pathway, determine the activation barriers, and obtain valuable insights concerning the specific roles of key amino acid residues for the common steps of AHAS‐catalyzed condensation reactions of α‐keto acids. The computational results reveal that the rate‐determining step of the AHAS‐catalyzed reactions is the second reaction step and that the most important amino acid residues involved in the catalysis include Glu144′, Gln207′, Gly121′, and Gly511 that form favorable hydrogen bonds with the reaction center (consisting of atoms from the substrate and cofactor) during the reaction process. In addition, Glu144′ also accepts a proton from cofactor thiamin diphosphate (ThDP) through hydrogen bonding during the catalytic reaction. The favorable interactions between the reaction center and protein environment remarkably stabilize the transition state and, thus, lower the activation barrier for the rate‐determining reaction step by ～20 kcal/mol. The activation barrier calculated for the rate‐determining step is in good agreement with the experimental activation barrier. The detailed structural and mechanistic insights should be valuable for rational design of novel, potent AHAS inhibitors that may be used as promising new anti‐tuberculosis drugs, antibacterial drugs, and/or herbicides to overcome drug resistance problem. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Why substituting the asparagine at position 35 in Bacillus circulans xylanase with an aspartic acid remarkably improves the enzymatic catalytic activity? A quantum chemistry-based calculation study

Xylanases from Bacillus circulans (BCX) are known as configuration-retaining glycoside hydrolases, which hydrolyze xylans with two glutamic acid residues (Glu78 and Glu172) serving as catalytic active residues according to a double displacement mechanism. Existing experimental researches show that mutating the asparagines (Asn) to aspartic acid (Asp) at position 35 next to Glu172 can obviously improve the catalytic activity of BCX. To better understand the inherent mechanism for the experimental finding, we performed quantum chemistry calculations on two model systems to mimic the catalyses of wild-type and mutant BCXs. Geometrical structures and relative energies of intermediates and transition states involved in the hydrolysis reactions are given in detail. It is found that in the wild-type model system Asn35 interacts with Glu172 via a loose hydrogen bond, while in the mutant model system Asp35 forms a very tight hydrogen bond with Glu172. The glycosidic bond cleavage is proposed to be the rate-determining step for the hydrolysis reaction, whose barrier varies from 98 to 65 kJ mol⁻¹ when Asn35 is replaced by Asp35, showing the presence of Asp35 remarkably reduces the energy demand for the hydrolysis reaction. The present result provides a theoretical elucidation for why a single amino acid substitution can importantly influences catalytic activity of BCX.     

Proton exit channels in bovine cytochrome c oxidase

Cytochrome c oxidase (CcO) is the terminal transmembrane enzyme of the respiratory electron transport chain in aerobic cells. It catalyzes the reduction of oxygen to water and utilizes the free energy of the reduction reaction for proton pumping, a process which results in a membrane electrochemical proton gradient. Although the structure of the enzyme has been solved for several organisms, the molecular mechanism of proton pumping and proton exit pathways remain unknown. In our previous work, the continuum electrostatic calculations were employed to evaluate the electrostatic potential, energies, and protonation state of bovine cytochrome c oxidase for different redox states of the enzyme. A possible mechanism of oxygen reduction and proton pumping via His291 was proposed. In this paper, using electrostatic calculations, we examine the proton exit pathways in the enzyme. By monitoring the changes of the protonation states, proton affinities, and energies of electrostatic interactions between the titratable groups in different redox states of CcO, we identified the clusters of strongly interacting residues. Using these data, we detected four possible proton exit points on the periplasmic side of the membrane (Lys171B/Asp173B, His24B/Asp25B, Asp51, and Asp300). We then were able to trace the proton exit pathways and to evaluate the energy profiles along the paths. On the basis of energetic considerations and the conservation of the residues in a protein sequence, the most likely exit pathway is one via the Lys171B/Asp173B site. The obtained results are fully consistent with our His291 model of proton pumping, and provide a rationale for the absence of proton leaking in CcO between the pumping strokes.     

Class I ribonucleotide reductase revisited: the effect of removing a proton on Glu441

The substrate mechanism of class I ribonucleotide reductase has been revisited using the hybrid density functional B3LYP method. The molecular model used is based on the X-ray structure and includes all the residues of the R1 subunit commonly considered in the RNR substrate conversion scheme: Cys439 initiating the reaction as a thiyl radical, the redox-active cysteines Cys225 and Cys462, and the catalytically important Glu441 and Asn437. In contrast to previous theoretical studies of the overall mechanism, Glu441 is added as an anion. All relevant transition states have been optimized, including one where an electron is transferred 8 A from the disulfide to the substrate simultaneously with a proton transfer from Glu441. The calculated barrier for this step is 19.1 kcal/mol, which can be compared to the rate-limiting barrier indicated by experiments of about 17 kcal/mol. Even though the calculated barrier is somewhat higher than the experimental limit, the discrepancy is within the normal error bounds of B3LYP. The suggestion from the present modeling study is thus that a protonated Glu441 does not need to be present at the active site from the beginning of the catalytic cycle. However, the previously suggested mechanism with an initial protonation of Glu441 cannot be ruled out, because even with the cost added for protonation of Glu441 with a typical pK(a) of 4, the barrier for that mechanism is lower than the one obtained for the present mechanism. The results are compared to experimental results and suggestions.     

Factors governing the protonation state of Zn-bound histidine in proteins: a DFT/CDM study

We have performed systematic theoretical studies to elucidate the factors governing the His protonation/deprotonation state in Zn-binding sites, especially those containing the ubiquitous Zn-His-Asp/Glu triad. Specifically, we have addressed the following three questions: (1) How does the transfer of the Zn-bound His imidazole proton to the second-shell Asp/Glu carboxylate oxygen depend on the composition of the other first-shell ligands and the solvent accessibility of the metal-binding site? (2) Can any second-shell ligand with a proton acceptor group such as the backbone carbonyl oxygen also act as a proton acceptor? (3) What is the effect of the Asp/Glu in the Zn-His-Asp/Glu triad on the Zn-bound water protonation state? To address these questions, we used a combination of quantum mechanical and continuum dielectric methods to compute the free energies for deprotonating a Zn-bound imidazole/water in various Zn complexes. The calculations show that whether the Zn-bound His is protonated or deprotonated depends on (1) the solvent accessibility of the metal-binding site, and (2) the Lewis acid ability of Zn, which is indirectly determined by both the first- and the second-shell Zn ligands. The calculations also show that the effect of the Zn-His-Asp/Glu interaction on the nucleophilicity of the Zn-bound water depends on the solvent accessibility of the catalytic Zn site. Furthermore, they show that the Asp/Glu side chain in the Zn-His-Asp/Glu triad can increase the negative charge of its partner, His, and create an anionic hole that may stabilize a cation in buried cavities, provided that the Zn complex is cationic/neutral. The findings of this work are in accord with available experimental data.     

Conversion of pyruvate decarboxylase into an enantioselective carboligase with biosynthetic potential

Pyruvate decarboxylase (PDC) catalyzes the decarboxylation of pyruvate into acetaldehyde and CO(2) and requires the cofactors thiamin diphosphate and Mg(2+) for activity. Owing to its catalytic promiscuity and relaxed substrate specificity, PDC catalyzes carboligation side reactions and is exploited for the asymmetric synthesis of 2-hydroxy ketones such as (R)-phenylacetyl carbinol, the precursor of (-)-ephedrine. Although PDC variants with enhanced carboligation efficiency were generated in the past, the native reaction, i.e., formation of aldehydes, is heavily favored over carboligation side reactions in all these biocatalysts. We characterized an active site variant (Glu473Gln) in which partitioning between aldehyde release versus carboligation is inverted with an up to 100-fold preference for the latter pathway. Due to a defective protonation of the central carbanion/enamine intermediate, substrate turnover stalls at this catalytic stage and addition of external aldehydes leads to quantitative and enantioselective formation of 2-hydroxy ketones as shown for (R)-phenylacetyl carbinol, which is afforded with unmatched yields, rates, and purity. This protein variant thus constitutes an example for the rational design of biocatalysts with greatly enhanced accidental catalytic promiscuity by selective blockage of the native reaction and accumulation of reactive intermediates under steady-state turnover conditions.     

Electron paramagnetic resonance detection of intermediates in the enzymatic cycle of an extradiol dioxygenase

Extradiol catecholic dioxygenases catalyze the cleavage of the aromatic ring of the substrate with incorporation of both oxygen atoms from O2. These enzymes are important in nature for the recovery of large amounts of carbon from aromatic compounds. The catalytic site contains either Fe or Mn coordinated by a facial triad of two His and one Glu or Asp residues. Previous studies have shown that Fe(II) and Mn(II) can be interchanged in enzymes from different organisms to catalyze similar substrate reactions. In combination, quantitative electron paramagnetic resonance spectroscopy and rapid freeze-quench experiments allow us to follow the concentrations of four different Mn species, including key metal intermediates in the catalytic cycle, as the enzyme turns over its natural substrate. Two intermediates are observed: a Mn(III)-radical species which is either Mn-superoxide or Mn-substrate radical, and a unique Mn(II) species which is involved in the rate-limiting step of the cycle and may be Mn-alkylperoxo.     

Differential effects of the Zn-His-Bkb vs Zn-His-[Asp/Glu] triad on Zn-core stability and reactivity

The most common partner of the Zn-bound His is the Asp/Glu carboxylate side chain in catalytic Zn sites and the backbone (Bkb) carbonyl group in structural Zn sites. To elucidate the factors governing the selection of the second-shell partner of the Zn-bound His in structural/catalytic Zn sites, systematic studies using density functional theory and continuum dielectric calculations were performed to determine the relative contributions of the second-shell Bkb carbonyl and the Asp/Glu carboxylate to the Zn-core stability and reactivity. The results show that the contributions of the second-shell Bkb carbonyl and Asp/Glu carboxylate to the Zn-core stability depend mainly on the solvent accessibility of the Zn-site and the composition of the Zn-core. They reveal the advantage of a second-shell Bkb carbonyl in anionic Zn cavities: it stabilizes anionic, buried Zn-cores more than the corresponding negatively charged Asp/Glu carboxylate, thus explaining the absence of the Zn-His-Asp/Glu triad in structural [Zn(Cys)3(His)]- cores. They also reveal the advantage of a second-shell Asp/Glu carboxylate in catalytic Zn-cores: relative to a Bkb carbonyl group, it increases (i) the HOMO energy of the cationic/neutral zinc core, (ii) the reactivity of the attacking Zn-bound OH-, (iii) electron transfer to the substrate, and (iv) the stability of the metal complex upon electron transfer. Furthermore, a second-shell Asp/Glu carboxylate could facilitate product release in the common cationic catalytic cores, by acting as a proton acceptor of the Zn-bound His creating an Asp...His- dyad that stabilizes the zinc dication more than the respective Bkb...His0 dyad.     

Copper-zinc superoxide dismutase: theoretical insights into the catalytic mechanism

The mechanism for the toxic superoxide radical disproportionation to molecular oxygen and hydrogen peroxide by copper-zinc superoxide dismutase (CuZnSOD) has been studied using the B3LYP hybrid density functional. On the basis of the X-ray structure of the enzyme, the molecular system investigated includes the first-shell protein ligands of the two metal centers as well as the second-shell ligand Asp122. The substrates of the model reaction are two superoxide radical anions, approaching the copper center at the beginning of two half-reactions: the first part of the catalytic cycle involving Cu+ oxidation and the second part reducing Cu2+ back to its initial state. The quantitative free energy profile of the reaction is obtained and discussed in connection with the experimental data on the reduction potentials and CuZnSOD kinetics. The optimized structures are analyzed and compared to the experimental ones. The two transition states alternate the protonation state of His61 and correspond to histidine Cu-His61-Zn bridge rupture/reformation. Modifications applied to the initial model allow the importance of Asp122 for catalysis to be estimated.     

Theoretical investigation on the dissociation of (R)‐benzoin catalyzed by benzaldehyde lyase

Benzaldehyde lyase (BAL) is a versatile thiamin diphosphate (THDP)‐dependent enzyme with widespread synthetic applications in industry. Besides lyase activity, BAL also performs the functions as carboligase and decarboxylase. Unlike many other THDP‐dependent enzymes, the active center of BAL is devoid of any acid‐base amino acid residues except Glu50 and His29, and therefore, the catalytic mechanism of BAL is unusual. In this article, the dissociation mechanism of (R)‐benzoin to benzaldehyde catalyzed by BAL has been studied by using density functional theory method. The calculation results indicate that the whole reaction consists of four elementary steps, and at least two steps contribute to rate‐limiting. A big difference with other THDP‐dependent enzymes is that, in the first stage of the reaction, the ligation of substrate and THDP ylide is not companied by proton transfer, and in the subsequent transition states and intermediates, the carbonyl oxygen always exists in the form of anion. Gln113, His29, and 4′‐amino group of THDP are found to have the function to stabilize the transition states and intermediates. His29 acts as the proton acceptor in step 2 and proton donor in step 3 using one water molecule as mediator. © 2013 Wiley Periodicals, Inc.     

The catalytic mechanism of mouse renin studied with QM/MM calculations

Hypertension is a chronic condition that affects nearly 25% of adults worldwide. As the Renin-Angiotensin-Aldosterone System is implicated in the control of blood pressure and body fluid homeostasis, its combined blockage is an attractive therapeutic strategy currently in use for the treatment of several cardiovascular conditions. We have performed QM/MM calculations to study the mouse renin catalytic mechanism in atomistic detail, using the N-terminal His6-Asn14 segment of angiotensinogen as substrate. The enzymatic reaction (hydrolysis of the peptidic bond between residues in the 10th and 11th positions) occurs through a general acid/base mechanism and, surprisingly, it is characterized by three mechanistic steps: it begins with the creation of a first very stable tetrahedral gem-diol intermediate, followed by protonation of the peptidic bond nitrogen, giving rise to a second intermediate. In a final step the peptidic bond is completely cleaved and both gem-diol hydroxyl protons are transferred to the catalytic dyad (Asp32 and Asp215). The final reaction products are two separate peptides with carboxylic acid and amine extremities. The activation energy for the formation of the gem-diol intermediate was calculated as 23.68 kcal mol(-1), whereas for the other steps the values were 15.51 kcal mol(-1) and 14.40 kcal mol(-1), respectively. The rate limiting states were the reactants and the first transition state. The associated barrier (23.68 kcal mol(-1)) is close to the experimental values for the angiotensinogen substrate (19.6 kcal mol(-1)). We have also tested the influence of the density functional on the activation and reaction energies. All eight density functionals tested (B3LYP, B3LYP-D3, X3LYP, M06, B1B95, BMK, mPWB1K and B2PLYP) gave very similar results.     

Catalytic mechanism of guanine deaminase: an ONIOM and molecular dynamics study

The catalytic mechanism of Bacillus subtilis guanine deaminase (bGD), a Zn metalloenzyme, has been investigated by a combination of quantum mechanical calculations using the multilayered ONIOM method and molecular dynamics simulations. In contrast to a previously proposed catalytic mechanism, which requires the bound guanine to assume a rare tautomeric state, the ONIOM calculations showed that the active-site residues of the enzyme do not affect the tautomeric state of guanine, and consequently the bound guanine is a tautomer that is the most abundant in aqueous solution. Two residues, Glutamate 55 and Aspartate 114, were found to play important roles in proton shuttling in the reaction. The proposed reaction path is initiated by proton transfer from a Zn-bound water to protonate Asp114. This process may be quite complex and rather dynamic in nature, as revealed by the molecular dynamics (MD) simulations, whereby another water may bridge the Zn-bound water and Asp114, which then is eliminated by positioning of guanine in the active site. The binding of guanine stabilizes protonated Asp114 by hydrogen bond formation. Asp114 can then transfer its proton to the N3 of the bound guanine, facilitating the nucleophilic attack on C2 of the guanine by the Zn-bound hydroxide to form a tetrahedral intermediate. This occurs with a rather low barrier. Glu55 then transfers a proton from the Zn-hydroxide to the amino group of the reaction intermediate and, at this point, the C2-N2 bond has lengthened by 0.2 A compared to guanine, making C2-N2 bond cleavage more facile. The C2-N2 bond breaks forming ammonia, with an energy barrier of approximately 8.8 kcal/mol. Ammonia leaves the active site, and xanthine is freed by the cleavage of the Zn-O2 bond, with a barrier approximately 8.4 kcal/mol. Along this reaction path, the highest barrier comes from C2-N2 bond cleavage, while the barrier from the cleavage of the Zn-O2 bond is slightly smaller. The Zn-O2 bond can be broken without the assistance of water during the release of xanthine.     

Direct NMR observation and pKa determination of the Asp102 side chain in a serine protease

The pK(a) value of aspartic acid in the catalytic triad of serine proteases has been a pivotal element in essentially every mechanism proposed for these enzymes over the past 40 years, but has, until now, eluded direct determination. Here, we have used the multinuclear 3D-NMR pulse programs HCACO and HCCH-TOCSY to directly identify and study the side-chain resonances of the aspartate and glutamate residues in uniformly (13)C-labeled α-lytic protease. Resonances from four of the six residues were detected and assigned, including that of Asp(102), which is notably the weakest of the four. pH titrations have shown all of the carboxylate (13)C signals to have unusually low pK(a) values: 2.0, 3.2, and 1.7 for Glu(129), Glu(174), and Glu(229), respectively, and an upper limit of 1.5 for Asp(102). The multiple H-bonds to Asp(102), long known from X-ray crystal studies, probably account for its unusually low pK(a) value through preferential stabilization of its anionic form. These H-bonds probably also contribute to the weakness of the NMR resonances of Asp(102) by restricting its mobility. The Asp(102)(13)C(γ) atom responds to the ionization of His(57) in the resting enzyme and to the inhibitor-derived oxyanion in a chloromethyl ketone complex, observations that strongly support the assignment. The low pK(a) value of Asp(102) would appear to be incompatible with mechanisms involving strong Asp(102)-His(57) H-bonds or high pK(a) values, but is compatible with mechanisms involving normal Asp(102)-His(57) H-bonds and moving His(57) imidazole rings, such as the reaction-driven ring flip.     

What is the active species of cytochrome P450 during camphor hydroxylation? QM/MM studies of different electronic states of compound I and of reduced and oxidized iron-oxo intermediates

We have investigated C-H hydroxylation of camphor by Compound I (Cpd I) of cytochrome P450cam in different electronic states and by its one-electron reduced and oxidized forms, using QM/MM calculations in the native protein/solvent environment. Cpd I species with five unpaired electrons (pentaradicaloids) are ca. 12 kcal/mol higher in energy than the ground state Cpd I species with three unpaired electrons (triradicaloids). The H-abstraction transition states of pentaradicaloids lie ca. 21 (9) kcal/mol above the triradicaloid (pentaradicaloid) reactants. Hydroxylation via pentaradicaloids is thus facile provided that they can react before relaxing to the ground-state triradicaloids. Excited states of Cpd I with an Fe(V)-oxo moiety lie more than 20 kcal/mol above the triradicaloid ground state in single-point gas-phase calculations, but these electronic configurations are not stable upon including the point-charge protein environment which causes SCF convergence to the triradicaloid ground state. One-electron reduced species (Cpd II) show sluggish reactivity compared with Cpd I in agreement with experimental model studies. One-electron oxidized species are more reactive than Cpd I but seem too high in energy to be accessible. The barriers to hydrogen abstraction for the various forms of Cpd I are generally not affected much by the chosen protonation states of the Asp297 and His355 residues near the propionate side chains of the heme or by the appearance of radical character at Asp297, His355, or the propionates.     

QM/MM study of mechanisms for compound I formation in the catalytic cycle of cytochrome P450cam

In the catalytic cycle of cytochrome P450cam, after molecular oxygen binds as a ligand to the heme iron atom to yield a ferrous dioxygen complex, there are fast proton transfers that lead to the formation of the active species, Compound I (Cpd I), which are not well understood because they occur so rapidly. In the present work, the conversion of the ferric hydroperoxo complex (Cpd 0) to Cpd I has been investigated by combined quantum-mechanical/molecular-mechanical (QM/MM) calculations. The residues Asp(251) and Glu(366) are considered as proton sources. In mechanism I, a proton is transported to the distal oxygen atom of the hydroperoxo group via a hydrogen bonding network to form protonated Cpd 0 (prot-Cpd0: FeOOH(2)), followed by heterolytic O-O bond cleavage that generates Cpd I and water. Although a local minimum is found for prot-Cpd0 in the Glu(366) channel, it is very high in energy (more than 20 kcal/mol above Cpd 0) and the barriers for its decay are only 3-4 kcal/mol (both toward Cpd 0 and Cpd I). In mechanism II, an initial O-O bond cleavage followed by a concomitant proton and electron transfer yields Cpd I and water. The rate-limiting step in mechanism II is O-O cleavage with a barrier of about 13-14 kcal/mol. According to the QM/MM calculations, the favored low-energy pathway to Cpd I is provided by mechanism II in the Asp(251) channel. Cpd 0 and Cpd I are of similar energies, with a slight preference for Cpd I.     

DFT/electrostatic calculations of pK(a) values in cytochrome c oxidase

Using classical electrostatic calculations, earlier we examined the dependence of the protonation state of bovine cytochrome c oxidase (CcO) on its redox state. Based on these calculations, we have proposed a model of CcO proton pumping that involves His291, one of the Cu(B) histidine ligands, which was found to respond to redox changes of the enzyme Fe(a)(3)-Cu(B) catalytic center. In this work, we employ combined density functional and continuum electrostatic calculations to evaluate the pK(a)() values of His291 and Glu242, two key residues of the model. The pK(a) values are calculated for different redox states of the enzyme, and the influence of different factors on the pK(a)'s is analyzed in detail. The calculated pK(a)() values of Glu242 are between 9.4 and 12.0, depending on the redox state of the protein, which is in excellent agreement with recent experimental measurements. Assuming the reduced state of heme a(3), His291 of the oxidized Cu(B) center possesses a pK(a)() between 2.1 and 4.0, while His291 of the reduced Cu(B) center has a pK(a) above 17. The obtained results support the proposal that the His291 ligand of the Cu(B) center in CcO is a proton pump element.     

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

应用分子模拟方法研究了血管紧张素转换酶(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的结合能力较强,与实验结果一致.