Anton S. Petrov

Packaging Double-Helical DNA into Viral Capsids: Structures, Forces, and Energetics

Small, icosahedral double-stranded DNA bacteriophage pack their genomes tightly into preformed protein capsids using an ATP-driven motor. Coarse-grain molecular-mechanics models provide a detailed picture of DNA packaging in bacteriophage, revealing how conformation depends on capsid size and shape, and the presence or absence of a protein core. The forces that oppose packaging have large contributions from both electrostatic repulsions and the entropic penalty of confining the DNA into the capsid, whereas elastic deformations make only a modest contribution. The elastic deformation energy is very sensitive to the final conformation, whereas the electrostatic and entropic penalties are not, so the packaged DNA favors conformations that minimize the bending energy.

Evolution of the ribosome at atomic resolution

Significance Ribosomes exist in every cell and are responsible for translation from mRNA to protein. The structure of the ribosomal common core is highly conserved in all living species, while the outer regions of the ribosome are variable. Ribosomal RNA of eukaryotes contains expansion segments accreted onto the surface of the core, which is nearly identical in structure to that in prokaryotic ribosomes. Comparing eukaryotic and prokaryotic ribosomes allows us to identify 3D insertion fingerprints of the expansion segments. Similar fingerprints allow us to analyze the common core and detect ancestral expansion segments within it. We construct a molecular model of ribosomal evolution starting from primordial biological systems near the dawn of life, culminating with relatively recent changes specific to metazoans. The origins and evolution of the ribosome, 3–4 billion years ago, remain imprinted in the biochemistry of extant life and in the structure of the ribosome. Processes of ribosomal RNA (rRNA) expansion can be “observed” by comparing 3D rRNA structures of bacteria (small), yeast (medium), and metazoans (large). rRNA size correlates well with species complexity. Differences in ribosomes across species reveal that rRNA expansion segments have been added to rRNAs without perturbing the preexisting core. Here we show that rRNA growth occurs by a limited number of processes that include inserting a branch helix onto a preexisting trunk helix and elongation of a helix. rRNA expansions can leave distinctive atomic resolution fingerprints, which we call “insertion fingerprints.” Observation of insertion fingerprints in the ribosomal common core allows identification of probable ancestral expansion segments. Conceptually reversing these expansions allows extrapolation backward in time to generate models of primordial ribosomes. The approach presented here provides insight to the structure of pre-last universal common ancestor rRNAs and the subsequent expansions that shaped the peptidyl transferase center and the conserved core. We infer distinct phases of ribosomal evolution through which ribosomal particles evolve, acquiring coding and translocation, and extending and elaborating the exit tunnel.

Structural and electrostatic characterization of Pariacoto virus: Implications for viral assembly

We present the first all‐atom model for the structure of a T = 3 virus, pariacoto virus (PaV), which is a nonenveloped, icosahedral RNA virus and a member of the Nodaviridae family. The model is an extension of the crystal structure, which reveals about 88% of the protein structure but only about 35% of the RNA structure. New modeling methods, combining coarse‐grained and all‐atom approaches, were required for developing the model. Evaluation of alternative models confirms our earlier observation that the polycationic N‐ and C‐terminal tails of the capsid proteins must penetrate deeply into the core of the virus, where they stabilize the structure by neutralizing a substantial fraction of the RNA charge. This leads us to propose a model for the assembly of small icosahedral RNA viruses: nonspecific binding of the protein tails to the RNA leads to a collapse of the complex, in a fashion reminiscent of DNA condensation. The globular protein domains are excluded from the condensed phase but are tethered to it, so they accumulate in a shell around the condensed phase, where their concentration is high enough to trigger oligomerization and formation of the mature virus.

Secondary structure and domain architecture of the 23S and 5S rRNAs

We present a de novo re-determination of the secondary (2°) structure and domain architecture of the 23S and 5S rRNAs, using 3D structures, determined by X-ray diffraction, as input. In the traditional 2° structure, the center of the 23S rRNA is an extended single strand, which in 3D is seen to be compact and double helical. Accurately assigning nucleotides to helices compels a revision of the 23S rRNA 2° structure. Unlike the traditional 2° structure, the revised 2° structure of the 23S rRNA shows architectural similarity with the 16S rRNA. The revised 2° structure also reveals a clear relationship with the 3D structure and is generalizable to rRNAs of other species from all three domains of life. The 2° structure revision required us to reconsider the domain architecture. We partitioned the 23S rRNA into domains through analysis of molecular interactions, calculations of 2D folding propensities and compactness. The best domain model for the 23S rRNA contains seven domains, not six as previously ascribed. Domain 0 forms the core of the 23S rRNA, to which the other six domains are rooted. Editable 2° structures mapped with various data are provided (http://apollo.chemistry.gatech.edu/RibosomeGallery).

RNA Folding and Catalysis Mediated by Iron (II)

Mg2+ shares a distinctive relationship with RNA, playing important and specific roles in the folding and function of essentially all large RNAs. Here we use theory and experiment to evaluate Fe2+ in the absence of free oxygen as a replacement for Mg2+ in RNA folding and catalysis. We describe both quantum mechanical calculations and experiments that suggest that the roles of Mg2+ in RNA folding and function can indeed be served by Fe2+. The results of quantum mechanical calculations show that the geometry of coordination of Fe2+ by RNA phosphates is similar to that of Mg2+. Chemical footprinting experiments suggest that the conformation of the Tetrahymena thermophila Group I intron P4–P6 domain RNA is conserved between complexes with Fe2+ or Mg2+. The catalytic activities of both the L1 ribozyme ligase, obtained previously by in vitro selection in the presence of Mg2+, and the hammerhead ribozyme are enhanced in the presence of Fe2+ compared to Mg2+. All chemical footprinting and ribozyme assays in the presence of Fe2+ were performed under anaerobic conditions. The primary motivation of this work is to understand RNA in plausible early earth conditions. Life originated during the early Archean Eon, characterized by a non-oxidative atmosphere and abundant soluble Fe2+. The combined biochemical and paleogeological data are consistent with a role for Fe2+ in an RNA World. RNA and Fe2+ could, in principle, support an array of RNA structures and catalytic functions more diverse than RNA with Mg2+ alone.

History of the ribosome and the origin of translation

Significance The ribosome, in analogy with a tree, contains a record of its history, spanning 4 billion years of life on earth. The information contained within ribosomes connects us to the prehistory of biology. Details of ribosomal RNA variation, observed by comparing three-dimensional structures of ribosomes across the tree of life, form the basis of our molecular-level model of the origins and evolution of the translational system. We infer many steps in the evolution of translation, mapping out acquisition of structure and function, revealing much about how modern biology originated from ancestral chemical systems. We present a molecular-level model for the origin and evolution of the translation system, using a 3D comparative method. In this model, the ribosome evolved by accretion, recursively adding expansion segments, iteratively growing, subsuming, and freezing the rRNA. Functions of expansion segments in the ancestral ribosome are assigned by correspondence with their functions in the extant ribosome. The model explains the evolution of the large ribosomal subunit, the small ribosomal subunit, tRNA, and mRNA. Prokaryotic ribosomes evolved in six phases, sequentially acquiring capabilities for RNA folding, catalysis, subunit association, correlated evolution, decoding, energy-driven translocation, and surface proteinization. Two additional phases exclusive to eukaryotes led to tentacle-like rRNA expansions. In this model, ribosomal proteinization was a driving force for the broad adoption of proteins in other biological processes. The exit tunnel was clearly a central theme of all phases of ribosomal evolution and was continuously extended and rigidified. In the primitive noncoding ribosome, proto-mRNA and the small ribosomal subunit acted as cofactors, positioning the activated ends of tRNAs within the peptidyl transferase center. This association linked the evolution of the large and small ribosomal subunits, proto-mRNA, and tRNA.

Secondary Structures of rRNAs from All Three Domains of Life

Accurate secondary structures are important for understanding ribosomes, which are extremely large and highly complex. Using 3D structures of ribosomes as input, we have revised and corrected traditional secondary (2°) structures of rRNAs. We identify helices by specific geometric and molecular interaction criteria, not by co-variation. The structural approach allows us to incorporate non-canonical base pairs on parity with Watson-Crick base pairs. The resulting rRNA 2° structures are up-to-date and consistent with three-dimensional structures, and are information-rich. These 2° structures are relatively simple to understand and are amenable to reproduction and modification by end-users. The 2° structures made available here broadly sample the phylogenetic tree and are mapped with a variety of data related to molecular interactions and geometry, phylogeny and evolution. We have generated 2° structures for both large subunit (LSU) 23S/28S and small subunit (SSU) 16S/18S rRNAs of Escherichia coli, Thermus thermophilus, Haloarcula marismortui (LSU rRNA only), Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens. We provide high-resolution editable versions of the 2° structures in several file formats. For the SSU rRNA, the 2° structures use an intuitive representation of the central pseudoknot where base triples are presented as pairs of base pairs. Both LSU and SSU secondary maps are available (http://apollo.chemistry.gatech.edu/RibosomeGallery). Mapping of data onto 2° structures was performed on the RiboVision server (http://apollo.chemistry.gatech.edu/RiboVision).

Bidentate RNA-magnesium clamps: On the origin of the special role of magnesium in RNA folding

Magnesium plays a special role in RNA function and folding. Although water is magnesiums most common first-shell ligand, the oxyanions of RNA have significant affinity for magnesium. Here we provide a quantum mechanical description of first-shell RNA-magnesium and DNA-magnesium interactions, demonstrating the unique features that characterize the energetics and geometry of magnesium complexes within large folded RNAs. Our work focuses on bidentate chelation of magnesium by RNA or DNA, where multiple phosphate oxyanions enter the first coordination shell of magnesium. These bidentate RNA clamps of magnesium occur frequently in large RNAs. The results here suggest that magnesium, compared to calcium and sodium, has an enhanced ability to form bidentate clamps with RNA. Bidentate RNA-sodium clamps, in particular, are unstable and spontaneously open. Due to magnesiums size and charge density it binds more intimately than other cations to the oxyanions of RNA, so that magnesium clamps are stabilized not only by electrostatic interactions, but also by charge transfer, polarization, and exchange interactions. These nonelectrostatic components of the binding are quite substantial with the high charge and small interatomic distances within the magnesium complexes, but are less pronounced for calcium due to its larger size, and for sodium due to its smaller charge. Additionally, bidentate RNA clamps of magnesium are more stable than those with DNA. The source of the additional stability of RNA complexes is twofold: there is a slightly attenuated energetic penalty for ring closure in the formation of RNA bidentate chelation complexes and elevated electrostatic interactions between the RNA and cations. In sum, it can be seen why sodium and calcium cannot replicate the structures or energetics of RNA-magnesium complexes.

The Role of DNA Twist in the Packaging of Viral Genomes

We performed molecular dynamics simulations of the genome packaging of bacteriophage P4 using two coarse-grained models of DNA. The first model, 1DNA6 (one pseudo-atom per six DNA basepairs), represents DNA as a string of beads, for which DNA torsions are undefined. The second model, 3DNA6 (three pseudo-atoms per six DNA basepairs), represents DNA as a series of base planes with torsions defined by the angles between successive planes. Bacteriophage P4 was packaged with 1DNA6, 3DNA6 in a torsionally relaxed state, and 3DNA6 in a torsionally strained state. We observed good agreement between the packed conformation of 1DNA6 and the packed conformations of 3DNA6. The free energies of packaging were in agreement, as well. Our results suggest that DNA torsions can be omitted from coarse-grained bacteriophage packaging simulations without significantly altering the DNA conformations or free energies of packaging that the simulations predict.

Molecular paleontology: a biochemical model of the ancestral ribosome

Ancient components of the ribosome, inferred from a consensus of previous work, were constructed in silico, in vitro and in vivo. The resulting model of the ancestral ribosome presented here incorporates ∼20% of the extant 23S rRNA and fragments of five ribosomal proteins. We test hypotheses that ancestral rRNA can: (i) assume canonical 23S rRNA-like secondary structure, (ii) assume canonical tertiary structure and (iii) form native complexes with ribosomal protein fragments. Footprinting experiments support formation of predicted secondary and tertiary structure. Gel shift, spectroscopic and yeast three-hybrid assays show specific interactions between ancestral rRNA and ribosomal protein fragments, independent of other, more recent, components of the ribosome. This robustness suggests that the catalytic core of the ribosome is an ancient construct that has survived billions of years of evolution without major changes in structure. Collectively, the data here support a model in which ancestors of the large and small subunits originated and evolved independently of each other, with autonomous functionalities.