Publications by Author: Bercem Dutagaci

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg2+, and, more importantly, a putative third (site C) Mg2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Catalysis and fidelity of multisubunit RNA polymerases rely on a highly conserved active site domain called the trigger loop (TL), which achieves roles in transcription through conformational changes and interaction with NTP substrates. The mutations of TL residues cause distinct effects on catalysis including hypo- and hyperactivity and altered fidelity. We applied molecular dynamics simulation (MD) and machine learning (ML) techniques to characterize TL mutations in the Saccharomyces cerevisiae RNA Polymerase II (Pol II) system. We did so to determine relationships between individual mutations and phenotypes and to associate phenotypes with MD simulated structural alterations. Using fitness values of mutants under various stress conditions, we modeled phenotypes along a spectrum of continual values. We found that ML could predict the phenotypes with 0.68 R2 correlation from amino acid sequences alone. It was more difficult to incorporate MD data to improve predictions from machine learning, presumably because MD data is too noisy and possibly incomplete to directly infer functional phenotypes. However, a variational auto-encoder model based on the MD data allowed the clustering of mutants with different phenotypes based on structural details. Overall, we found that a subset of loss-of-function (LOF) and lethal mutations tended to increase distances of TL residues to the NTP substrate, while another subset of LOF and lethal substitutions tended to confer an increase in distances between TL and bridge helix (BH). In contrast, some of the gain-of-function (GOF) mutants appear to cause disruption of hydrophobic contacts among TL and nearby helices.