Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: Quantum mechanical/molecular mechanical free energy simulations

Sixue Zhang, Abir Ganguly, Puja Goyal, Jamie L. Bingaman, Philip C. Bevilacqua, Sharon Hammes-Schiffer

Research output: Contribution to journalArticle

38 Citations (Scopus)

Abstract

The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2′), which would be followed by proton transfer from G40(N1) to A-1(O2′). After this initial deprotonation, A-1(O2′) starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2′) to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2′) to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg2+ ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

Original languageEnglish (US)
Pages (from-to)784-798
Number of pages15
JournalJournal of the American Chemical Society
Volume137
Issue number2
DOIs
StatePublished - Jan 21 2015

Fingerprint

Catalytic RNA
Proton transfer
Guanine
Free energy
Protons
Rate constants
Catalytic Domain
Phosphates
Hydroxyl Radical
Glucosamine
Oxygen
Deprotonation
Energy barriers
Catalysis
Hydrogen
Ions
Acids

All Science Journal Classification (ASJC) codes

  • Catalysis
  • Chemistry(all)
  • Biochemistry
  • Colloid and Surface Chemistry

Cite this

@article{f662d22799c24669bc1a83ac54f6bbf2,
title = "Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: Quantum mechanical/molecular mechanical free energy simulations",
abstract = "The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2′), which would be followed by proton transfer from G40(N1) to A-1(O2′). After this initial deprotonation, A-1(O2′) starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2′) to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2′) to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg2+ ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.",
author = "Sixue Zhang and Abir Ganguly and Puja Goyal and Bingaman, {Jamie L.} and Bevilacqua, {Philip C.} and Sharon Hammes-Schiffer",
year = "2015",
month = "1",
day = "21",
doi = "10.1021/ja510387y",
language = "English (US)",
volume = "137",
pages = "784--798",
journal = "Journal of the American Chemical Society",
issn = "0002-7863",
publisher = "American Chemical Society",
number = "2",

}

Role of the active site guanine in the glmS ribozyme self-cleavage mechanism : Quantum mechanical/molecular mechanical free energy simulations. / Zhang, Sixue; Ganguly, Abir; Goyal, Puja; Bingaman, Jamie L.; Bevilacqua, Philip C.; Hammes-Schiffer, Sharon.

In: Journal of the American Chemical Society, Vol. 137, No. 2, 21.01.2015, p. 784-798.

Research output: Contribution to journalArticle

TY - JOUR

T1 - Role of the active site guanine in the glmS ribozyme self-cleavage mechanism

T2 - Quantum mechanical/molecular mechanical free energy simulations

AU - Zhang, Sixue

AU - Ganguly, Abir

AU - Goyal, Puja

AU - Bingaman, Jamie L.

AU - Bevilacqua, Philip C.

AU - Hammes-Schiffer, Sharon

PY - 2015/1/21

Y1 - 2015/1/21

N2 - The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2′), which would be followed by proton transfer from G40(N1) to A-1(O2′). After this initial deprotonation, A-1(O2′) starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2′) to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2′) to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg2+ ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

AB - The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2′), which would be followed by proton transfer from G40(N1) to A-1(O2′). After this initial deprotonation, A-1(O2′) starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2′) to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2′) to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg2+ ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

UR - http://www.scopus.com/inward/record.url?scp=84921519751&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84921519751&partnerID=8YFLogxK

U2 - 10.1021/ja510387y

DO - 10.1021/ja510387y

M3 - Article

C2 - 25526516

AN - SCOPUS:84921519751

VL - 137

SP - 784

EP - 798

JO - Journal of the American Chemical Society

JF - Journal of the American Chemical Society

SN - 0002-7863

IS - 2

ER -