Skip to main content


  • Oral presentation
  • Open Access

Investigation on the reaction mechanism of ribonuclease enzymes

  • 1,
  • 2,
  • 3 and
  • 1
Chemistry Central Journal20093 (Suppl 1) :O23

  • Published:


  • Bipyramidal
  • Free Energy Calculation
  • Transition State Structure
  • Vanadate Complex
  • Cyclic Phosphorane
Pancreatic ribonuclease A (RNase A) is an efficient enzyme for the hydrolysis of phosphate-ester linkages in single-stranded RNA. There is strong evidence that the catalytic process occurs via a two-step general acid-base mechanism. The first step involves intramolecular transphosphorylation leading to a trigonal bipyramidal cyclic phosphorane intermediate and displacement of the O5'-nucleotide product. In the subsequent step, the cyclic phosphorane is hydrolyzed to a 3'-monophosphate nucleotide as the final product. Three amino acid residues within the active site, Lys41, His12 and His119, are known to play significant roles in the catalytic process (Figure 1). They are believed to act as a general base and general acid. Because of the availability of high resolution X-ray and neutron scattering structural data [4], RNase A attracts considerable attention as a target for theoretical studies [1][2][3][5][6].

Figure 1

The purpose of this research program was to investigate the reaction mechanism of this enzyme using high level quantum chemistry methods of NWChem [7] with the highest feasible accuracy (DFT/B3LYP, Ahlrichs-pVDZ basis set) to provide an accurate representation of the chemistry of the active site (QM region) and molecular mechanics (Amber 99) to describe the remaining protein and solvent (MM region).

For this enzyme reaction there is no data providing a detailed picture of any structures along the reaction path [5]. However, there are structural data for a vanadate complex [4], which has been accepted as a transition state analogue. After evaluation of our computational analysis we substituted the central vanadium by phosphorus to model the intermediate state.

Afterwards we elongated the endocyclic and exocyclic P-O bonds to find reactant and product states allowing for relaxation at each step. In addition a thermally averaged NEB method has been implemented to produce an unbiased view of the movement of atoms as the system evolves along the reaction path. Free energy calculations have been performed to account for thermal fluctuations of the protein environment. This way we could calculate the energy profile of the hydrolysis step and identify the transition state structures along the reaction coordinate.

Authors’ Affiliations

University of Paderborn, Warburgerstr. 100, 33098 Paderborn, Germany
Pacific Northwest National Laboratory, Richland, WA-99352, USA
University of California San Diego, 9500 Gilman Dr, La Jolla, CA-92093, USA


  1. Richards FM, Wyckoff HW: The Enzymes. 1971, Academic Press, New York, 4:Google Scholar
  2. Blackburn P, Moore S: The Enzymes. 1982, Ac. Press, New York, 15:Google Scholar
  3. Lim C, Tole P: J Phys Chem. 1992, 97: 6212-J Am Chem Soc, 1992, 114: 7245. View ArticleGoogle Scholar
  4. Wlodawer A, Miller M, Sjölin L: Proc Natl Acad Sci USA. 1983, 80: 3628-3631. 10.1073/pnas.80.12.3628.View ArticleGoogle Scholar
  5. Lindquist RN, Lynn LJ, Lienhard GE: J Am Chem Soc. 1973, 95 (26): 8762-8768. 10.1021/ja00807a043.View ArticleGoogle Scholar
  6. Lopez X, York DM, Dejaeger A, Karplus M: Int J of Quant Chem. 2002, 86: 10-26. 10.1002/qua.1601.View ArticleGoogle Scholar
  7. Aprà E, Bylaska EJ, Valiev M, et al: NWChem, A Computational Chemistry Package for Parallel Computers, Version 4.7". 2005, Pacific Northwest National Laboratory, Richland, USAGoogle Scholar


© Elsässer et al; licensee BioMed Central Ltd. 2009

This article is published under license to BioMed Central Ltd.