Journal of the American Chemical Society, Vol.124, No.49, 14591-14600, 2002
Exploring nucleoside hydrolase catalysis in silico: Molecular dynamics study of enzyme-bound substrate and transition state
The mechanism of action of inosine-uridine nucleoside hydrolase has been investigated by long-term molecular dynamics (MD) simulation in TIP3P water using stochastic boundary conditions. Five MD studies have been performed with enzyme substrate complex (E.S), enzyme substrate complex with protonated His241 (EH.S), enzyme transition state complex (E.TS), enzyme transition state complex with protonated His241 (EH-TS), and His241 Ala transition state complex E(H241A).TS. Special attention has been given to the role of His241, which has been considered as the general acid catalyst to assist departure of the leaving nucleobase on the basis of its location in the active site in the X-ray crystal structure (2MAS). Yet on the basis of the location in the active site, Tyr229 is closer to the aniline ring of pAPIR as compared to His241. On initiation of MD simulations, His241 does not approach the nucleobase in the structures of EH.S, E.S, EH.TS, and E.TS. In the solvated enzyme, Tyr229, which is a member of the hydrogen bonding network inosine O2'-.Aspl4..His241..Tyr229..inosine N7, serves as a proton source to the leaving nucleobase. The loss of significant activity of His241 Ala mutant is shown to be related to the disruption of the above hydrogen bonded network and the distancing of Tyr229 from inosine N7. The structures of the enzyme complexes with substrate or TS are not visibly altered on protonation of His241, a most unusual outcome. The bell-shaped pH dependence upon pK(app)'s of 7.1 and 9.1 may be attributed to the necessity of the dissociation of Asp10 or Asp15 and the acid form of Tyr229, respectively. In TS, the residue IIe81 migrated closer, whereas Arg233 moved away from the nucleobase. The probability of ribooxocarbenium ion stabilization by Asn168 and Asp14 is discussed. The Asp14-CO2-is hydrogen bonded to the ribose 2'-OH for 96% of the MD simulation time. Nucleophilic addition of water138 to ribooxocarbenium ion is suggested to be assisted by the proton shuttle from water138-->Asp10-->Asp15--> water pool. An anticorrelation motion between Tyr229-OH and Asn168-OD1 in EH.S and E.S is observed. The relationship of this anticorrelated motion to mechanism, if any, deserves further exploration, perhaps the formation of a near attack conformation.