Catalyzed decomposition of urea. Molecular dynamics simulations of the binding of urea to urease

We present the results of molecular dynamics simulations on the urea/urease system. The starting structure was prepared from the 2.0 angstrom crystal structure of Benini et al. [(1999) Struct. Folding Des. 7, 205-216] of DAP-inhibited urease (PDB code 3UBP), and the trimeric structure (2479 residues) resulted in 180K atoms after solvation by water. The force field parameters were derived using the bonded model approach described by Hoops et al. [(1991) J. Am. Chem. Soc. 113, 8262-8270]. Three different systems were analyzed, each one modeling a different protonation pattern for the His320 and His219 residues. In each case, the three monomers of urease have been analyzed separately. The time-averaged structures observed in the three monomers suggest that urease could follow two different competitive mechanisms. A "protein-assisted proton transfer" mechanism points to Asp221 as crucial for catalysis. An "Asp-mediated proton transfer" involves the transfer of a proton from the bridging OH to an NH2 moiety of urea, assisted by Asp360 in the active site. The impact of the simulation results on our understanding of urease catalysis is discussed in detail.