Protein engineering of nitrile hydratase activity of papain: molecular dynamics study of a mutant and wild-type enzyme

The mechanism of hydrolysis of the nitrile (N-acetyl-phenylalanyl-2-amino-propionitrile, I) catalyzed by Gln19Glu mutant of papain has been studied by nanosecond molecular dynamics (MD) simulations. MD simulations of the complex of mutant enzyme with I and of mutant enzyme covalently attached to both neutral (II) and protonated (III) thioimidate intermediates were performed. An MD simulation with the wild-type enzyme.I complex was undertaken as a reference. The ion pair between protonated His159 and thiolate of Cys25 is coplanar, and the hydrogen bonding interaction S(-)(25).HD1-ND1(159) is observed throughout MD simulation of the mutant enzyme.I complex. Such a sustained hydrogen bond is absent in nitrile-bound wild-type papain due to the flexibility of the imidazole ring of His159. The nature of the residue at position 19 plays a critical role in the hydrolysis of the covalent thioimidate intermediate. When position 19 represents Glu, the imidazolium ion of His159-ND1(+).Cys25-S(-) ion pair is distant, on average, from the nitrile nitrogen of substrate I. Near attack conformers (NACs) have been identified in which His159-ImH(+) is positioned to initiate a general acid-catalyzed addition of Cys-S(-) to nitrile. Though Glu19-CO(2)H is distant from nitrile nitrogen in the mutant.I structure, MD simulations of the mutant.II covalent adduct finds Glu19-CO(2)H hydrogen bonded to the thioimide nitrogen of II. This hydrogen bonded species is much less stable than the hydrogen bonded Glu19-CO(2)(-) with mutant-bound protonated thioimidate (III). This observation supports Glu19-CO(2)H general acid catalysis of the formation of mutant.III. This is the commitment step in the Gln19Glu mutant catalysis of nitrile hydrolysis.     

In silico studies of the mechanism of methanol oxidation by quinoprotein methanol dehydrogenase

The mechanism of bacterial methanol dehydrogenase involves hydride equivalent transfer from substrate to the ortho-quinone PQQ to provide a C5-reduced intermediate that subsequently rearranges to the hydroquinone PQQH(2). We have studied the PQQ reduction by molecular dynamic (MD) simulations in aqueous solution. Among the five simulated structures, either Asp297 or Glu171 or both are ionized. Reasonable structures are obtained only when both carboxyl groups are ionized. This is not unexpected since the kinetic pH optimum is 9.0. In the structure of the enzyme.PQQ.HOCH(3) complex, the hydrogen bonded Glu171-CO(2)(-).H-OCH(3) is in a position to act as a general base catalyst for hydride equivalent transfer to C5 of PQQ. We thus suggest that Glu171 plays the role of general base catalyst in PQQ reduction rather than Asp297 as previously suggested. The reduction is assisted by Arg324, which hydrogen bonds to the ortho-quinone moiety of PQQ. The rearrangement of the C5-reduced intermediate to provide hydroquinone PQQH(2) is also assisted by proton abstraction by Glu171-CO(2)(-) and the continuous hydrogen bonding of Arg324 throughout the entire reaction. These features as well as the mapping of the channel for substrate and water into the active site entrance are the observations of major importance.     

Mechanism of glucose oxidation by quinoprotein soluble glucose dehydrogenase: insights from molecular dynamics studies

We have generated 3 ns molecular dynamic (MD) simulations, in aqueous solution, of the bacterial soluble glucose dehydrogenase enzyme.PQQ.glucose complex and intermediates formed in PQQ reduction. In the MD structure of enzyme.PQQ.glucose complex the imidazole of His144 is hydrogen bonded to the hydroxyl hydrogen of H[bond]OC1(H) of glucose. The tightly hydrogen-bonded triad Asp163-His144-glucose (2.70 and 2.91 A) is involved in proton abstraction from glucose concerted with the hydride transfer from the C1[bond]H of glucose to the >C5[double bond]O quinone carbon of PQQ. The reaction is assisted by Arg228 hydrogen bonding to the carbonyl oxygen of >C5[double bond]O. The rearrangement of [bond](H)C5(O-)[bond]C4([double bond]O)[bond] of II to [bond]C5(OH)[double bond]C4(OH)[bond] of PQQH(2) hydroquinone is assisted by general acid protonatation of the >C4[double bond]O oxygen by protonated His144 and hydrogen bonds of Arg228 to the oxyanion O5. The continuous hydrogen bonding of the amide side chain of Asn229 to >C4[double bond]O4 oxygen and that of the O5 oxygen of the cofactor to Wat89 is observed throughout the entire reaction.     

ELECTROCHEMICAL LUMINESCENCE WITH N(5)-ETHYL-4a-HYDROXY-3-METHYL-4a, 5-DIHYDROLUMIFLAVIN. THE MECHANISM OF BACTERIAL LUCIFERASE

It has been proposed in the literature that the chemiluminescence of the flavoenzyme of bacterial luciferase is caused by a chemically initiated electron-exchange luminescence mechanism which provides an excited 4a-hydroxy-4a,5-dihydroflavin ([4a-FlHOH]*) as product of 1e- reduction of the radical 4a-FlHOH+.. Electrochemical/photon counting experiments were performed to assess the feasibility of this proposal. Potentials for step-wise oxidation of N(5)-ethyl-4a-hydroxy-4a,5-dihydroluminflavin (4a-FlEtOH) have been determined in dry N,N-dimethylformamide (DMF). Photon counting was carried out during the 1e- reduction of 4a-FlEtOH+.in both DMF and acetonitrile by use of an apparatus consisting of a photocell mounted below a Pt ring-disk electrode. By use of the ring-disk electrode a steady state concentration of [4a-FlEtOH]* could be maintained by continuous 1e- oxidation of 4a-FlEtOH----4a-FlEtOH+.and 1e- reduction of 4a-FlEtOH+.----4a-FlEtOH. A maximum of 14% collection (theoretical maximum is 18%) of FlEtOH.+ at the ring electrode was obtained below 5000 rotations per minute. Calibration of the apparatus using 9,10-diphenylanthracene allowed approximation of the quantum yield for 1e- reductive capture of 4a-FlEtOH+.as 10(-6) to 10(-4) in DMF and 10(-7) to 10(-5) in acetonitrile. No fluorescence for 4a-FlEtOH in DMF could be observed; if fluorescent, the efficiency of 4a-FlEtOH can be no greater than approximately 3 x 10(-5). No electrogenerated chemiluminescence is observed on the electrochemical recycling of FlEt+----FlEt2+ and FlEt2+----FlEt+.     

Computer simulations of trypanosomal nucleoside hydrolase: determination of the protonation state of the bound transition-state analogue

Inosine-uridine nucleoside hydrolase (IU-NH) catalyzes the hydrolysis of nucleosides into base and ribose moieties via a ribooxocarbenium ion transition state, which has been characterized using kinetic isotope effects. Protozoan parasites lack de novo purine and pyrimidine biosynthesis and depend on the purine salvage from the host. Vern Schramm and co-workers characterized p-aminophenyliminoribitol (pAPIR) to be a potent inhibitor of IU-NH from Crithidia fasciculata with K(d) of 30 nM. The cyclic amine function of the iminoribitol ring can be either protonated (pAPIRH(+)) or unprotonated (pAPIR). pAPIRH(+) resembles the charge and geometry of the ribooxocarbenium ion transition state and can be looked upon as a transition-state analogue inhibitor; however, it is known that the pAPIR species is initially bound to the enzyme. We have characterized the pAPIRH(+) species as resident of the active site using ab initio calculations and molecular dynamics simulations. This is a novel use of molecular dynamics to investigate the protonation state of the bound ligand to the active site. Nanosecond molecular dynamics simulations reveal a short hydrogen-bonding network between pAPIRH(+)-O2'-Asp14-His241 triad, which is not seen in the crystal structure. Other features discussed are: hydrogen bonding between pAPIRH(+) and Asn168, unusual geometry of the iminoribitol ring, and hydrophobic interactions.