Abstract: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 R 2 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. RNA polymerase II (Pol II) synthesizes RNA with the help of an active site domain called the trigger loop (TL). Mutations in the TL cause changes in the activity of Pol II that range from gain-of-function (GOF, viable but hyperactive) to loss-of-function (LOF, viable but hypoactive) or lethal. This study provides a systematic characterization of the structural and functional outcomes of the TL mutations using molecular dynamics (MD) simulations and machine learning (ML). We obtained functional phenotypes of mutants by ML using genetic fitness scores (measure of growth defect strength) as input. We revealed that mutant TL sequences could predict the functional outcomes at a relatively high correlation. Then, we performed MD simulations to relate structural information to the phenotypes. The analysis of the MD data suggested that there are two subsets of lethal and LOF mutants, where one subset had increased distances between the TL and the substrate, while the other subset showed increased distances between TL and another active site domain called the bridge helix (BH). On the other hand, some of the GOF mutants altered a key hydrophobic pocket formed by interactions between residues near the active site. Overall, this study enhances our understanding of the effects of TL mutations to the Pol II function.

Connecting Mutations of the RNA Polymerase II C-Terminal Domain to Complex Phenotypic Changes Using Combined Gene Expression and Network Analyses

Mutations in MTO2 Related to Trna Modification Impair Mitochondrial Gene Expression and Protein Synthesis in the Presence of a Paromomycin Resistance Mutation in Mitochondrial 15 S Rrna

Comparative Genomics and Evolution of Proteins Associated with RNA Polymerase II C-terminal Domain.

Characterization of RNA polymerase II trigger loop mutations using molecular dynamics simulations and machine learning

Alternative stable conformation capable of protein misinteraction links tRNA synthetase to peripheral neuropathy.

Coordinated Evolution Of Transcriptional And Post-Transcriptional Regulation For Mitochondrial Functions In Yeast Strains

Cohesins coordinate gene transcriptions of related function within Saccharomyces cerevisiae

Dissecting the chemical interactions and substrate structural signatures governing RNA polymerase II trigger loop closure by synthetic nucleic acid analogues.

Widespread epistasis shapes RNA Polymerase II active site function and evolution

A New Point Mutation in the Nuclear Gene of Yeast Mitochondrial RNA Polymerase, RPO41, Identifies a Functionally Important Amino-Acid Residue in a Protein Region Conserved among Mitochondrial Core Enzymes

Variation of C-terminal domain governs RNA polymerase II genomic locations and alternative splicing in eukaryotic transcription

Sequence and structural determinants of RNAPII CTD phase-separation and phosphorylation by CDK7

Disordered C-terminal domain drives spatiotemporal confinement of RNAPII to enhance search for chromatin targets

Combining Nanopore direct RNA sequencing with genetics and mass spectrometry for analysis of T-loop base modifications across 42 yeast tRNA isoacceptors

Clustering of RNA Polymerase II C-Terminal Domain Models upon Phosphorylation

Functional and transcriptional profiling of non-coding RNAs in yeast reveal context-dependent phenotypes and widespread in trans effects on the protein regulatory network

Transcriptional dysregulation by a nucleus-localized aminoacyl-tRNA synthetase associated with Charcot-Marie-Tooth neuropathy

Dissection of a complex phenotype by functional genomics reveals roles for the yeast cyclin-dependent protein kinase Pho85 in stress adaptation and cell integrity

Repression of pervasive antisense transcription is the primary role of fission yeast RNA polymerase II CTD serine 2 phosphorylation

Higher-order epistasis within Pol II trigger loop haplotypes