Structures of additional crystal forms of Satellite tobacco mosaic virus grown from a variety of salts.

The structures of new crystal forms of Satellite tobacco mosaic virus (STMV) are described. These belong to space groups I2, P21212 (a low-resolution form), R3 (H3) and P23. The R3 crystals are 50%/50% twinned, as are two instances of the P23 crystals. The I2 and P21212 crystals were grown from ammonium sulfate solutions, as was one crystal in space group P23, while the R3 and the other P23 crystals were grown from sodium chloride, sodium bromide and sodium nitrate. The monoclinic and orthorhombic crystals have half a virus particle as the asymmetric unit, while the rhombohedral and cubic crystals have one third of a virus particle. RNA segments organized about the icosahedral twofold axes were present in crystals grown from ammonium sulfate and sodium chloride, as in the canonical I222 crystals (PDB entry 4oq8), but were not observed in crystals grown from sodium bromide and sodium nitrate. Bromide and nitrate ions generally replaced the RNA phosphates present in the I222 crystals, including the phosphates seen on fivefold axes, and were also found at threefold vertices in both the rhombohedral and cubic forms. An additional anion was also found on the fivefold axis 5 Šfrom the first anion, and slightly outside the capsid in crystals grown from sodium chloride, sodium bromide and sodium nitrate, suggesting that the path along the symmetry axis might be an ion channel. The electron densities for RNA strands at individual icosahedral dyads, as well as at the amino-terminal peptides of protein subunits, exhibited a diversity of orientations, in particular the residues at the ends.

Structural 3D Domain Reconstruction of the RNA Genome from Viruses with Secondary Structure Models.

Three-dimensional RNA domain reconstruction is important for the assembly, disassembly and delivery functionalities of a packed proteinaceus capsid. However, to date, the self-association of RNA molecules is still an open problem. Recent chemical probing reports provide, with high reliability, the secondary structure of diverse RNA ensembles, such as those of viral genomes. Here, we present a method for reconstructing the complete 3D structure of RNA genomes, which combines a coarse-grained model with a subdomain composition scheme to obtain the entire genome inside proteinaceus capsids based on secondary structures from experimental techniques. Despite the amount of sampling involved in the folded and also unfolded RNA molecules, advanced microscope techniques can provide points of anchoring, which enhance our model to include interactions between capsid pentamers and RNA subdomains. To test our method, we tackle the satellite tobacco mosaic virus (STMV) genome, which has been widely studied by both experimental and computational communities. We provide not only a methodology to structurally analyze the tertiary conformations of the RNA genome inside capsids, but a flexible platform that allows the easy implementation of features/descriptors coming from both theoretical and experimental approaches.

Investigation into the binding of dyes within protein crystals.

It was found that the crystals of at least a dozen different proteins could be thoroughly stained to an intense color with a panel of dyes. Many, if not most, of the stained protein crystals retained the dyes almost indefinitely when placed in large volumes of dye-free mother liquor. Dialysis experiments showed that most of the dyes that were retained in crystals also bound to the protein when free in solution; less frequently, some dyes bound only in the crystal. The experiments indicated a strong association of the dyes with the proteins. Four protein crystals were investigated by X-ray diffraction to ascertain the mode of binding. These were crystals of lysozyme, thaumatin, trypsin inhibited with benzamidine and satellite tobacco mosaic virus. In 30 X-ray analyses of protein crystal-dye complexes, in only three difference Fourier maps was any difference electron density present that was consistent with the binding of dye molecules, and even in these three cases (thaumatin plus thioflavin T, xylene cyanol and m-cresol purple) the amount of dye observed was inadequate to explain the intense color of the crystals. It was concluded that the dye molecules, which are clearly inside the crystals, are disordered but are paradoxically tightly bound to the protein. It is speculated that the dyes, which exhibit large hydrophobic cores and peripheral charged groups, may interact with the crystalline proteins in the manner of conventional detergents.

The Renormalization Group and Its Applications to Generating Coarse-Grained Models of Large Biological Molecular Systems.

Understanding the dynamics of biomolecules is the key to understanding their biological activities. Computational methods ranging from all-atom molecular dynamics simulations to coarse-grained normal-mode analyses based on simplified elastic networks provide a general framework to studying these dynamics. Despite recent successes in studying very large systems with up to a 100,000,000 atoms, those methods are currently limited to studying small- to medium-sized molecular systems due to computational limitations. One solution to circumvent these limitations is to reduce the size of the system under study. In this paper, we argue that coarse-graining, the standard approach to such size reduction, must define a hierarchy of models of decreasing sizes that are consistent with each other, i.e., that each model contains the information of the dynamics of its predecessor. We propose a new method, Decimate, for generating such a hierarchy within the context of elastic networks for normal-mode analysis. This method is based on the concept of the renormalization group developed in statistical physics. We highlight the details of its implementation, with a special focus on its scalability to large systems of up to millions of atoms. We illustrate its application on two large systems, the capsid of a virus and the ribosome translation complex. We show that highly decimated representations of those systems, containing down to 1% of their original number of atoms, still capture qualitatively and quantitatively their dynamics. Decimate is available as an OpenSource resource.

Symmetry-adapted digital modeling III. Coarse-grained icosahedral viruses.

Considered is the coarse-grained modeling of icosahedral viruses in terms of a three-dimensional lattice (the digital modeling lattice) selected among the projected points in space of a six-dimensional icosahedral lattice. Backbone atomic positions (Cα's for the residues of the capsid and phosphorus atoms P for the genome nucleotides) are then indexed by their nearest lattice point. This leads to a fine-grained lattice point characterization of the full viral chains in the backbone approximation (denoted as digital modeling). Coarse-grained models then follow by a proper selection of the indexed backbone positions, where for each chain one can choose the desired coarseness. This approach is applied to three viruses, the Satellite tobacco mosaic virus, the bacteriophage MS2 and the Pariacoto virus, on the basis of structural data from the Brookhaven Protein Data Bank. In each case the various stages of the procedure are illustrated for a given coarse-grained model and the corresponding indexed positions are listed. Alternative coarse-grained models have been derived and compared. Comments on related results and approaches, found among the very large set of publications in this field, conclude this article.

SPECTRUS: A Dimensionality Reduction Approach for Identifying Dynamical Domains in Protein Complexes from Limited Structural Datasets.

Identifying dynamical, quasi-rigid domains in proteins provides a powerful means for characterizing functionally oriented structural changes via a parsimonious set of degrees of freedom. In fact, the relative displacements of few dynamical domains usually suffice to rationalize the mechanics underpinning biological functionality in proteins and can even be exploited for structure determination or refinement purposes. Here we present SPECTRUS, a general scheme that, by solely using amino acid distance fluctuations, can pinpoint the innate quasi-rigid domains of single proteins or large complexes in a robust way. Consistent domains are usually obtained by using either a pair of representative structures or thousands of conformers. The functional insights offered by the approach are illustrated for biomolecular systems of very different size and complexity such as kinases, ion channels, and viral capsids. The decomposition tool is available as a software package and web server at spectrus.sissa.it.

Probing viral genomic structure: alternative viewpoints and alternative structures for satellite tobacco mosaic virus RNA.

Viral RNA structure prediction is a valuable tool for development of drugs against viral disease. This work discusses different approaches to predicting encapsidated viral RNA and highlights satellite tobacco mosaic virus (STMV) RNA as a model system with excellent crystallography data. Fundamentally important issues for debate include thermodynamic versus kinetic control of virus assembly and the possible consequences of quasi-species in the primary structure on RNA secondary structure prediction of a single structure or an ensemble of structures. Multiple computational tools and chemical reagents are now available for improved viral RNA structure prediction. Two different predicted structures for encapsidated STMV RNA result from differences in three main areas: a different approach and philosophy to studying encapsidated viral RNA, an emphasis on different RNA motifs, and technical differences in computational methods and chemical reagents. The experiments with traditional chemical probing and SHAPE reagents are compared in terms of chemistry, results, and interpretation for STMV RNA as well as other RNA protein assemblies, such as the 5'UTR of HIV and the ribosome. This discussion of the challenges of viral RNA structure prediction will lead to new experiments and improved future predictions for viral RNA.

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.

Space warping order parameters and symmetry: application to multiscale simulation of macromolecular assemblies.

Coarse-grained features of macromolecular assemblies are understood via a set of order parameters (OPs) constructed in terms of their all-atom configuration. OPs are shown to be slowly changing in time and capture the large-scale spatial features of macromolecular assemblies. The relationship of these variables to the classic notion of OPs based on symmetry breaking phase transitions is discussed. OPs based on space warping transformations are analyzed in detail as they naturally provide a connection between overall structure of an assembly and all-atom configuration. These OPs serve as the basis of a multiscale analysis that yields Langevin equations for OP dynamics. In this context, the characteristics of OPs and PCA modes are compared. The OPs enable efficient all-atom multiscale simulations of the dynamics of macromolecular assemblies in response to changes in microenvironmental conditions, as demonstrated on the structural transitions of cowpea chlorotic mottle virus capsid (CCMV) and RNA of the satellite tobacco mosaic virus (STMV).

Form, symmetry and packing of biomacromolecules. IV. Filled capsids of cowpea, tobacco, MS2 and pariacoto RNA viruses.

Four icosahedral RNA viruses are considered: the cowpea chlorotic mottle virus, the satellite tobacco mosaic virus, the pariacoto virus and the MS2 bacteriophage. The validity of the phenomenological rules derived in previous publications (crystallographic scaling, indexed forms enclosing axial-symmetric clusters, packing lattices of viral crystals) is confirmed and shown to apply equally well to the coat proteins as to the (ordered) RNA chains.

Form, symmetry and packing of biomacromolecules. V. Shells with boundaries at anti-nodes of resonant vibrations in icosahedral RNA viruses.

The RNA viruses cowpea chlorotic mottle, satellite tobacco mosaic, pariacoto and MS2, already considered in part IV of this series of papers [Janner, A. (2011a), Acta Cryst. A67, 517-520], are investigated further, with the aim to arrive at a possible physical basis for their structural properties. The shell structure of the filled capsid is analyzed in terms of successive spherical boundaries of the sets of icosahedral equivalent chains. By inversion in the sphere enclosing the capsid, the internal boundaries are transformed into external ones, which are more easily visualized. This graphical procedure reveals the presence of regularly spaced shells with boundaries fitting with anti-nodal surfaces of the virus considered as an elastic resonator. The centers of gravity of the various chains occur in the nodal regions of eigenvibrations with wavelength λ = R(0)/K(0), where R(0) is the radius of the virus and K(0) takes one of the values 12, 6, 4, 3, depending on the mode. The resonator model is consistent with practically all spherical shell boundaries, whereas deviations are observed for the icosahedral axial modes, which apparently play a secondary role with respect to the spherical ones. Both the spherical and the axial anti-nodal surfaces fit very well with the packed structure of the viruses in the crystal which, accordingly, is expected to have eigenfrequencies related to those of the virus. These results open the way to a better understanding of the possibility of breaking the capsid using resonant forced oscillations excited, for example, by an applied elastic shock or by irradiation with femtosecond laser pulses, as already realised by K.-T. Tsen and co-workers. An alternative `plywood' model connected to the extreme elastic properties of the capsid is also considered.

Ensemble of secondary structures for encapsidated satellite tobacco mosaic virus RNA consistent with chemical probing and crystallography constraints.

Viral genomic RNA adopts many conformations during its life cycle as the genome is replicated, translated, and encapsidated. The high-resolution crystallographic structure of the satellite tobacco mosaic virus (STMV) particle reveals 30 helices of well-ordered RNA. The crystallographic data provide global constraints on the possible secondary structures for the encapsidated RNA. Traditional free energy minimization methods of RNA secondary structure prediction do not generate structures consistent with the crystallographic data, and to date no complete STMV RNA basepaired secondary structure has been generated. RNA-protein interactions and tertiary interactions may contribute a significant degree of stability, and the kinetics of viral assembly may dominate the folding process. The computational tools, Helix Find & Combine, Crumple, and Sliding Windows and Assembly, evaluate and explore the possible secondary structures for encapsidated STMV RNA. All possible hairpins consistent with the experimental data and a cotranscriptional folding and assembly hypothesis were generated, and the combination of hairpins that was most consistent with experimental data is presented as the best representative structure of the ensemble. Multiple solutions to the genome packaging problem could be an evolutionary advantage for viruses. In such cases, an ensemble of structures that share favorable global features best represents the RNA fold.

Vibrational dynamics of icosahedrally symmetric biomolecular assemblies compared with predictions based on continuum elasticity.

Coarse-grained elastic network models elucidate the fluctuation dynamics of proteins around their native conformations. Low-frequency collective motions derived by simplified normal mode analysis are usually involved in biological function, and these motions often possess noteworthy symmetries related to the overall shape of the molecule. Here, insights into these motions and their frequencies are sought by considering continuum models with appropriate symmetry and boundary conditions to approximately represent the true atomistic molecular structure. We solve the elastic wave equations analytically for the case of spherical symmetry, yielding a symmetry-based classification of molecular motions together with explicit predictions for their vibrational frequencies. We address the case of icosahedral symmetry as a perturbation to the spherical case. Applications to lumazine synthase, satellite tobacco mosaic virus, and brome mosaic virus show that the spherical elastic model efficiently provides insights on collective motions that are otherwise obtained by detailed elastic network models. A major utility of the continuum models is the possibility of estimating macroscopic material properties such as the Young's modulus or Poisson's ratio for different types of viruses.

Elastic properties of viruses.

Viruses are compact biological nanoparticles whose elastic and dynamical properties are hardly known. Inelastic (Brillouin) light scattering was used to characterize these properties, from microcrystals of the Satellite Tobacco Mosaic Virus, a nearly spherical plant virus of 17-nm diameter. Longitudinal sound velocities in wet and dry Satellite Tobacco Mosaic Virus crystals were determined and compared to that of the well-known protein crystal, lysozyme. Localized vibrational modes of the viral particles (i.e., particle modes) were sought in the relevant frequency ranges, as derived assuming the viruses as full free nanospheres. Despite very favorable conditions, regarding virus concentration and expected low damping in dry microcrystals, no firm evidence of virus particle modes could be detected.

Molecular dynamics simulations of the complete satellite tobacco mosaic virus.

This work presents an all-atom molecular dynamics simulation of a complete virus, the satellite tobacco mosaic virus. Simulations with up to 1 million atoms for over 50 ns demonstrate the stability of the entire virion and of the RNA core alone, while the capsid without RNA exhibits a pronounced instability. Physical properties of the simulated virus particle including electrostatic potential, radial distribution of viral components, and patterns of correlated motion are analyzed, and the implications for the assembly and infection mechanism of the virus are discussed.

Molecular structures of viruses from Raman optical activity.

A vibrational Raman optical activity (ROA) study of a range of different structural types of virus exemplified by filamentous bacteriophage fd, tobacco mosaic virus, satellite tobacco mosaic virus, bacteriophage MS2 and cowpea mosaic virus has revealed that, on account of its sensitivity to chirality, ROA is an incisive probe of their aqueous solution structures at the molecular level. Protein ROA bands are especially prominent from which, as we have shown by comparison with the ROA spectra of proteins with known structures and by using a pattern recognition program, the folds of the major coat protein subunits may be deduced. Information about amino acid side-chain conformations, exemplified here by the determination of the sign and magnitude of the torsion angle chi(2,1) for tryptophan in fd, may also sometimes be obtained. By subtracting the ROA spectrum of the empty protein capsid (top component) of cowpea mosaic virus from those of the intact middle and bottom-upper components separated by means of a caesium chloride density gradient, the ROA spectrum of the viral RNA was obtained, which revealed that the RNA takes up an A-type single-stranded helical conformation and that the RNA conformations in the middle and bottom-upper components are very similar. This information is not available from the X-ray crystal structure of cowpea mosaic virus since no nucleic acid is visible.

Self-repair of biological fibers catalyzed by the surface of a virus crystal.

Helical fibers, presumably proteinaceous and of microbial origin, have been visualized by atomic force microscopy on the surfaces of crystals of satellite tobacco mosaic virus. If the crystals are growing, then the fibers are incorporated intact into the crystal lattice. If broken on the crystal surface, then within a few minutes, the fibers self-reassemble to reestablish continuity. This, we believe, is the first observation of such a crystal surface-catalyzed repair of a biological structure. The surfaces of virus crystals provide ideal workbenches for the visualization and manipulation of nanoscale objects, particularly extended structures such as these fibers.

Structural transitions of satellite tobacco mosaic virus particles.

Satellite tobacco mosaic virus (STMV) can undergo at least two physical transitions that significantly alter its mechanical and structural characteristics. At high pH the 17-nm STMV particles expand radially by about 5 A to yield particles having diameters of about 18 nm. This pH-induced transition is further promoted by aging of the virions and degradation of the RNA, so that swollen particles ultimately appear even at neutral pH. While the native 17-nm particles crystallize as orthorhombic or monoclinic crystals which diffract to high resolution (1.8 A), the enlarged 18-nm particles crystallize in a cubic form which diffracts to no better than 5 A. In the transition, not only do the capsid protein subunits move radially outward, but the helical RNA segments with which they interact do as well. This is noteworthy because it demonstrates that the RNA and the protein shell are capable of coordinated movement, and that neither structure is rigidly defined or independent of the other. Using atomic force microscopy, it can be shown that STMV particles, upon drying, lose their mechanical rigidity and undergo deformation. Virions initially 17 nm in diameter shrink to more uniform final sizes than do 18 nm, initially swollen particles. This transition appears to be irreversible, as the particles do not reassume their former size nor structural rigidity upon rehydration. Evidence is also presented that preparations of native virus and their crystals are naturally somewhat heterogeneous and contain a variety of particles of anomalous size.

Biophysical studies on the RNA cores of satellite tobacco mosaic virus.

Satellite tobacco mosaic virus (STMV) was probed using a variety of proteases. Consequences of the degradation were analyzed using gel electrophoresis, quasi-elastic light scattering (QELS), and atomic force microscopy (AFM). Proteolysis rates of 30 minutes for complete degradation of the protein capsid, up to many hours, were investigated. With each protease, degradation of virions 17 nm in diameter was shown by QELS to result in particles of 10 nm diameter, which is that of the RNA core observed in the virion by x-ray diffraction analysis. This was verified by direct visualization with atomic force microscopy. Using QELS, it was further shown that freshly prepared RNA cores remain as individual, stable, 10-nm condensed particles for 12 to 24 h. Clusters of particles then formed, followed by very large aggregates of 500 to 1000 nm diameter. AFM showed that the aggregates were composed of groups of the condensed RNA cores and were not due to unfolding of the nucleic acid. No unfolding of the core particles into extended conformation was seen by AFM until the samples were heated well beyond 90 degrees C. Mass spectrometry of RNA core particles revealed the presence of a major polypeptide whose amino acid sequence corresponded to residues 2 through 25 of the coat protein. Amino acids 13 through 25 were previously observed to be in direct contact with the RNA and are presumably protected from protease digestion. Low resolution difference Fourier analyses indicated the courses of the remainders of the amino terminal strands (amino acids 2-12) in intact virions. Any individual strand appears to have several choices of path, which accounts for the observed disorder at high resolution. These positively charged strands, serving as virtual polyamines, engage the helical segments of RNA. The intimate association of amino acid residues 2 through 25 with RNA likely contributes to the stability of the condensed conformation of the nucleic acid cores.

Satellite tobacco mosaic virus RNA: structure and implications for assembly.

The initial appearance of 45% of the single-stranded RNA of satellite tobacco mosaic virus in electron density maps suggested the entire RNA conformation could be delineated. Subsequent work has localized nearly 80% of the RNA as stem-loop elements. Connection of the stem-loops in the most efficient manner produces a persuasive model for the encapsidated RNA. This arrangement has significant implications for virus assembly and for the essential role of RNA.