Function of a viral genome packaging motor from bacteriophage T4 is insensitive to DNA sequence.

Many viruses employ ATP-powered motors during assembly to translocate DNA into procapsid shells. Previous reports raise the question if motor function is modulated by substrate DNA sequence: (i) the phage T4 motor exhibits large translocation rate fluctuations and pauses and slips; (ii) evidence suggests that the phage phi29 motor contacts DNA bases during translocation; and (iii) one theoretical model, the 'B-A scrunchworm', predicts that 'A-philic' sequences that transition more easily to A-form would alter motor function. Here, we use single-molecule optical tweezers measurements to compare translocation of phage, plasmid, and synthetic A-philic, GC rich sequences by the T4 motor. We observed no significant differences in motor velocities, even with A-philic sequences predicted to show higher translocation rate at high applied force. We also observed no significant changes in motor pausing and only modest changes in slipping. To more generally test for sequence dependence, we conducted correlation analyses across pairs of packaging events. No significant correlations in packaging rate, pausing or slipping versus sequence position were detected across repeated measurements with several different DNA sequences. These studies suggest that viral genome packaging is insensitive to DNA sequence and fluctuations in packaging motor velocity, pausing and slipping are primarily stochastic temporal events.

DNA Conformational Changes Play a Force-Generating Role during Bacteriophage Genome Packaging.

Motors that move DNA, or that move along DNA, play essential roles in DNA replication, transcription, recombination, and chromosome segregation. The mechanisms by which these DNA translocases operate remain largely unknown. Some double-stranded DNA (dsDNA) viruses use an ATP-dependent motor to drive DNA into preformed capsids. These include several human pathogens as well as dsDNA bacteriophages-viruses that infect bacteria. We previously proposed that DNA is not a passive substrate of bacteriophage packaging motors but is instead an active component of the machinery. We carried out computational studies on dsDNA in the channels of viral portal proteins, and they reveal DNA conformational changes consistent with that hypothesis. dsDNA becomes longer ("stretched") in regions of high negative electrostatic potential and shorter ("scrunched") in regions of high positive potential. These results suggest a mechanism that electrostatically couples the energy released by ATP hydrolysis to DNA translocation: The chemical cycle of ATP binding, hydrolysis, and product release drives a cycle of protein conformational changes. This produces changes in the electrostatic potential in the channel through the portal, and these drive cyclic changes in the length of dsDNA as the phosphate groups respond to the protein's electrostatic potential. The DNA motions are captured by a coordinated protein-DNA grip-and-release cycle to produce DNA translocation. In short, the ATPase, portal, and dsDNA work synergistically to promote genome packaging.

The entropic force generated by intrinsically disordered segments tunes protein function.

Protein structures are dynamic and can explore a large conformational landscape1,2. Only some of these structural substates are important for protein function (such as ligand binding, catalysis and regulation)3-5. How evolution shapes the structural ensemble to optimize a specific function is poorly understood3,4. One of the constraints on the evolution of proteins is the stability of the folded 'native' state. Despite this, 44% of the human proteome contains intrinsically disordered peptide segments greater than 30 residues in length6, the majority of which have no known function7-9. Here we show that the entropic force produced by an intrinsically disordered carboxy terminus (ID-tail) shifts the conformational ensemble of human UDP-α-D-glucose-6-dehydrogenase (UGDH) towards a substate with a high affinity for an allosteric inhibitor. The function of the ID-tail does not depend on its sequence or chemical composition. Instead, the affinity enhancement can be accurately predicted based on the length of the intrinsically disordered segment, and is consistent with the entropic force generated by an unstructured peptide attached to the protein surface10-13. Our data show that the unfolded state of the ID-tail rectifies the dynamics and structure of UGDH to favour inhibitor binding. Because this entropic rectifier does not have any sequence or structural constraints, it is an easily acquired adaptation. This model implies that evolution selects for disordered segments to tune the energy landscape of proteins, which may explain the persistence of intrinsic disorder in the proteome.

Transitions of Double-Stranded DNA Between the A- and B-Forms.

The structure of double-stranded DNA (dsDNA) is sensitive to solvent conditions. In solution, B-DNA is the favored conformation under physiological conditions, while A-DNA is the form found under low water activity. The A-form is induced locally in some protein-DNA complexes, and repeated transitions between the B- and A-forms have been proposed to generate the forces used to drive dsDNA into viral capsids during genome packaging. Here, we report analyses on previous molecular dynamics (MD) simulations on B-DNA, along with new MD simulations on the transition from A-DNA to B-DNA in solution. We introduce the A-B Index (ABI), a new metric along the A-B continuum, to quantify our results. When A-DNA is placed in an equilibrated solution at physiological ionic strength, there is no energy barrier to the transition to the B-form, which begins within about 1 ns. The transition is essentially complete within 5 ns, although occasionally a stretch of a few base pairs will remain A-like for up to ∼10 ns. A comparison of four sequences with a range of predicted A-phobicities shows that more A-phobic sequences make the transition more rapidly than less A-phobic sequences. Simulations on dsDNA with a region of roughly one turn locked in the A-form allow us to characterize the A/B junction, which has an average bend angle of 20-30°. Fluctuations in this angle occur with characteristic times of about 10 ns.

DNA Scrunching in the Packaging of Viral Genomes.

The motors that drive double-stranded DNA (dsDNA) genomes into viral capsids are among the strongest of all biological motors for which forces have been measured, but it is not known how they generate force. We previously proposed that the DNA is not a passive substrate but that it plays an active role in force generation. This "scrunchworm hypothesis" holds that the motor proteins repeatedly dehydrate and rehydrate the DNA, which then undergoes cyclic shortening and lengthening motions. These are captured by a coupled protein-DNA grip-and-release cycle to rectify the motion and translocate the DNA into the capsid. In this study, we examined the interactions of dsDNA with the dodecameric connector protein of bacteriophage ϕ29, using molecular dynamics simulations on four different DNA sequences, starting from two different conformations (A-DNA and B-DNA). In all four simulations starting with the protein equilibrated with A-DNA in the channel, we observed transitions to a common, metastable, highly scrunched conformation, designated A*. This conformation is very similar to one recently reported by Kumar and Grubmüller in much longer MD simulations on B-DNA docked into the ϕ29 connector. These results are significant for four reasons. First, the scrunched conformations occur spontaneously, without requiring lever-like protein motions often believed to be necessary for DNA translocation. Second, the transition takes place within the connector, providing the location of the putative "dehydrator". Third, the protein has more contacts with one strand of the DNA than with the other; the former was identified in single-molecule laser tweezer experiments as the "load-bearing strand". Finally, the spontaneity of the DNA-protein interaction suggests that it may play a role in the initial docking of DNA in motors like that of T4 that can load and package any sequence.

History of the ribosome and the origin of translation.

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.

Visualizing the global secondary structure of a viral RNA genome with cryo-electron microscopy.

The lifecycle, and therefore the virulence, of single-stranded (ss)-RNA viruses is regulated not only by their particular protein gene products, but also by the secondary and tertiary structure of their genomes. The secondary structure of the entire genomic RNA of satellite tobacco mosaic virus (STMV) was recently determined by selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE). The SHAPE analysis suggested a single highly extended secondary structure with much less branching than occurs in the ensemble of structures predicted by purely thermodynamic algorithms. Here we examine the solution-equilibrated STMV genome by direct visualization with cryo-electron microscopy (cryo-EM), using an RNA of similar length transcribed from the yeast genome as a control. The cryo-EM data reveal an ensemble of branching patterns that are collectively consistent with the SHAPE-derived secondary structure model. Thus, our results both elucidate the statistical nature of the secondary structure of large ss-RNAs and give visual support for modern RNA structure determination methods. Additionally, this work introduces cryo-EM as a means to distinguish between competing secondary structure models if the models differ significantly in terms of the number and/or length of branches. Furthermore, with the latest advances in cryo-EM technology, we suggest the possibility of developing methods that incorporate restraints from cryo-EM into the next generation of algorithms for the determination of RNA secondary and tertiary structures.

The scrunchworm hypothesis: transitions between A-DNA and B-DNA provide the driving force for genome packaging in double-stranded DNA bacteriophages.

Double-stranded DNA bacteriophages have motors that drive the genome into preformed capsids, using the energy released by hydrolysis of ATP to overcome the forces opposing DNA packaging. Viral packaging motors are the strongest of all biological motors, but it is not known how they generate these forces. Several models for the process of mechanochemical force generation have been put forward, but there is no consensus on which, if any, of these is correct. All the existing models assume that protein-generated forces drive the DNA forward. The scrunchworm hypothesis proposes that the DNA molecule is the active force-generating core of the motor, not simply a substrate on which the motor operates. The protein components of the motor dehydrate a section of the DNA, converting it from the B form to the A form and shortening it by about 23%. The proteins then rehydrate the DNA, which converts back to the B form. Other regions of the motor grip and release the DNA to capture the shortening-lengthening motions of the B→A→B cycle ("scrunching"), so that DNA is pulled into the motor and pushed forward into the capsid. This DNA-centric mechanism provides a quantitative physical explanation for the magnitude of the forces generated by viral packaging motors. It also provides a simple explanation for the fact that each of the steps in the burst cycle advances the DNA by 2.5 base pairs. The scrunchworm hypothesis is consistent with a large body of published data, and it makes four experimentally testable predictions.

RiboVision suite for visualization and analysis of ribosomes.

RiboVision is a visualization and analysis tool for the simultaneous display of multiple layers of diverse information on primary (1D), secondary (2D), and three-dimensional (3D) structures of ribosomes. The ribosome is a macromolecular complex containing ribosomal RNA and ribosomal proteins and is a key component of life responsible for the synthesis of proteins in all living organisms. RiboVision is intended for rapid retrieval, analysis, filtering, and display of a variety of ribosomal data. Preloaded information includes 1D, 2D, and 3D structures augmented by base-pairing, base-stacking, and other molecular interactions. RiboVision is preloaded with rRNA secondary structures, rRNA domains and helical structures, phylogeny, crystallographic thermal factors, etc. RiboVision contains structures of ribosomal proteins and a database of their molecular interactions with rRNA. RiboVision contains preloaded structures and data for two bacterial ribosomes (Thermus thermophilus and Escherichia coli), one archaeal ribosome (Haloarcula marismortui), and three eukaryotic ribosomes (Saccharomyces cerevisiae, Drosophila melanogaster, and Homo sapiens). RiboVision revealed several major discrepancies between the 2D and 3D structures of the rRNAs of the small and large subunits (SSU and LSU). Revised structures mapped with a variety of data are available in RiboVision as well as in a public gallery (). RiboVision is designed to allow users to distill complex data quickly and to easily generate publication-quality images of data mapped onto secondary structures. Users can readily import and analyze their own data in the context of other work. This package allows users to import and map data from CSV files directly onto 1D, 2D, and 3D levels of structure. RiboVision has features in rough analogy with web-based map services capable of seamlessly switching the type of data displayed and the resolution or magnification of the display. RiboVision is available at .

Evolution of the ribosome at atomic resolution.

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.

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).

Comment on the letter by A. Ben-Shaul: "entropy, energy, and bending of DNA in viral capsids".

The conformational entropic penalty associated with packaging double-stranded DNA into viral capsids remains an issue of contention. So far, models based on a continuum approximation for DNA have either left the question unexamined, or they have assumed that the entropic penalty is negligible, following an early analysis by Riemer and Bloomfield. In contrast, molecular-dynamics (MD) simulations using bead-and-spring models consistently show a large penalty. A recent letter from Ben-Shaul attempts to reconcile the differences. While the letter makes some valid points, the issue of how to include conformational entropy in the continuum models remains unresolved. In this Comment, I show that the free energy decomposition from continuum models could be brought into line with the decomposition from the MD simulations with two adjustments. First, the entropy from Flory-Huggins theory should be replaced by the estimate of the entropic penalty given in Ben-Shaul's letter, which corresponds closely to that from the MD simulations. Second, the DNA-DNA repulsions are well described by the empirical relationship given by the Cal Tech group, but the strength of these should be reduced by about half, using parameters based on the Rau-Parsegian experiments, rather than treating them as "fitting parameters (tuned) to fit the data from (single molecule pulling) experiments."

Modeling the conformational changes underlying channel opening in CFTR.

Mutations in the gene encoding the cystic fibrosis transmembrane conductance regulator protein (CFTR) cause cystic fibrosis (CF), the most common life-shortening genetic disease among Caucasians. Although general features of the structure of CFTR have been predicted from homology models, the conformational changes that result in channel opening and closing have yet to be resolved. We created new closed- and open-state homology models of CFTR, and performed targeted molecular dynamics simulations of the conformational transitions in a channel opening event. The simulations predict a conformational wave that starts at the nucleotide binding domains and ends with the formation of an open conduction pathway. Changes in side-chain interactions are observed in all major domains of the protein, and experimental confirmation was obtained for a novel intra-protein salt bridge that breaks near the end of the transition. The models and simulation add to our understanding of the mechanism of ATP-dependent gating in this disease-relevant ion channel.

The icosahedral RNA virus as a grotto: organizing the genome into stalagmites and stalactites.

There are two important problems in the assembly of small, icosahedral RNA viruses. First, how does the capsid protein select the viral RNA for packaging, when there are so many other candidate RNA molecules available? Second, what is the mechanism of assembly? With regard to the first question, there are a number of cases where a particular RNA sequence or structure--often one or more stem-loops--either promotes assembly or is required for assembly, but there are others where specific packaging signals are apparently not required. With regard to the assembly pathway, in those cases where stem-loops are involved, the first step is generally believed to be binding of the capsid proteins to these "fingers" of the RNA secondary structure. In the mature virus, the core of the RNA would then occupy the center of the viral particle, and the stem-loops would reach outward, towards the capsid, like stalagmites reaching up from the floor of a grotto towards the ceiling. Those viruses whose assembly does not depend on protein binding to stem-loops could have a different structure, with the core of the RNA lying just under the capsid, and the fingers reaching down into the interior of the virus, like stalactites. We review the literature on these alternative structures, focusing on RNA selectivity and the assembly mechanism, and we propose experiments aimed at determining, in a given virus, which of the two structures actually occurs.

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 with iron(II) as a cofactor catalyses electron transfer.

Mg(2+) is essential for RNA folding and catalysis. However, for the first 1.5 billion years of life on Earth RNA inhabited an anoxic Earth with abundant and benign Fe(2+). We hypothesize that Fe(2+) was an RNA cofactor when iron was abundant, and was substantially replaced by Mg(2+) during a period known as the 'great oxidation', brought on by photosynthesis. Here, we demonstrate that reversing this putative metal substitution in an anoxic environment, by removing Mg(2+) and replacing it with Fe(2+), expands the catalytic repertoire of RNA. Fe(2+) can confer on some RNAs a previously uncharacterized ability to catalyse single-electron transfer. We propose that RNA function, in analogy with protein function, can be understood fully only in the context of association with a range of possible metals. The catalysis of electron transfer, requisite for metabolic activity, may have been attenuated in RNA by photosynthesis and the rise of O2.

Effects of pulling forces, osmotic pressure, condensing agents and viscosity on the thermodynamics and kinetics of DNA ejection from bacteriophages to bacterial cells: a computational study.

In this work, we report on simulations of double-stranded DNA (dsDNA) ejection from bacteriophage φ29 into a bacterial cell. The ejection was studied with a coarse-grained model, in which viral dsDNA was represented by beads on a torsion-less string. The bacteriophage's capsid and the bacterial cell were defined by sets of spherical constraints. To account for the effects of the viscous medium inside the bacterial cell, the simulations were carried out using a Langevin dynamics protocol. Our simplest simulations (involving constant viscosity and no external biasing forces) produced results compatible with the push-pull model of DNA ejection, with an ejection rate significantly higher in the first part of ejection than in the latter parts. Additionally, we performed more complicated simulations, in which we included additional factors such as external forces, osmotic pressure, condensing agents and ejection-dependent viscosity. The effects of these factors (independently and in combination) on the thermodynamics and kinetics of DNA ejection were studied. We found that, in general, the dependence of ejection forces and ejection rates on the amount of DNA ejected becomes more complex if the ejection is modeled with a broader, more realistic set of parameters and influences (such as variation in the solvent's viscosity and the application of an external force). However, certain combinations of factors and numerical parameters led to the opposition of some ejection-driving and ejection-inhibiting influences, ultimately causing an apparent simplification of the ejection profiles.

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.

In vitro secondary structure of the genomic RNA of satellite tobacco mosaic virus.

Satellite tobacco mosaic virus (STMV) is a T = 1 icosahedral virus with a single-stranded RNA genome. It is widely accepted that the RNA genome plays an important structural role during assembly of the STMV virion. While the encapsidated form of the RNA has been extensively studied, less is known about the structure of the free RNA, aside from a purported tRNA-like structure at the 3' end. Here we use selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis to examine the secondary structure of in vitro transcribed STMV RNA. The predicted secondary structure is unusual in the sense that it is highly extended, which could be significant for protecting the RNA from degradation. The SHAPE data are also consistent with the previously predicted tRNA-like fold at the 3' end of the molecule, which is also known to hinder degradation. Our data are not consistent with the secondary structure proposed for the encapsidated RNA by Schroeder et al., suggesting that, if the Schroeder structure is correct, either the RNA is packaged as it emerges from the replication complex, or the RNA undergoes extensive refolding upon encapsidation. We also consider the alternative, i.e., that the structure of the encapsidated STMV RNA might be the same as the in vitro structure presented here, and we examine how this structure might be organized in the virus. This possibility is not rigorously ruled out by the available data, so it remains open to examination by experiment.