De novo computational RNA modeling into cryoEM maps of large ribonucleoprotein complexes

Abstract

RNA-protein assemblies carry out many critical biological functions including translation, RNA splicing, and telomere extension. Increasingly, cryo-electron microscopy (cryoEM) is used to determine the structures of these complexes, but nearly all maps determined with this method have regions in which the local resolution does not permit manual coordinate tracing. Because RNA coordinates typically cannot be determined by docking crystal structures of separate components and existing structure prediction algorithms cannot yet model RNA-protein complexes, RNA coordinates are frequently omitted from final models despite their biological importance. To address these omissions, we have developed a new framework for De novo Ribonucleoprotein modeling in Real-space through Assembly of Fragments Together with Electron density in Rosetta (DRRAFTER). We show that DRRAFTER recovers near-native models for a diverse benchmark set of small RNA-protein complexes, as well as for large RNA-protein machines, including the spliceosome, mitochondrial ribosome, and CRISPR-Cas9-sgRNA complexes where the availability of both high and low resolution maps enable rigorous tests. Blind tests on yeast U1 snRNP and spliceosomal P complex maps demonstrate that the method can successfully build RNA coordinates in real-world modeling scenarios. Additionally, to aid in final model interpretation, we present a method for reliable in situ estimation of DRRAFTER model accuracy. Finally, we apply this method to recently determined maps of telomerase, the HIV-1 reverse transcriptase initiation complex, and the packaged MS2 genome, demonstrating that DRRAFTER can be used to accelerate accurate model building in challenging cases.