Yaroslava G. Yingling

RNAJunction: a database of RNA junctions and kissing loops for three-dimensional structural analysis and nanodesign

We developed a database called RNAJunction that contains structure and sequence information for RNA structural elements such as helical junctions, internal loops, bulges and loop–loop interactions. Our database provides a user-friendly way of searching structural elements by PDB code, structural classification, sequence, keyword or inter-helix angles. In addition, the structural data was subjected to energy minimization. This database is useful for analyzing RNA structures as well as for designing novel RNA structures on a nanoscale. The database can be accessed at: http://rnajunction.abcc.ncifcrf.gov/

Tertiary model of a plant cellulose synthase

A 3D atomistic model of a plant cellulose synthase (CESA) has remained elusive despite over forty years of experimental effort. Here, we report a computationally predicted 3D structure of 506 amino acids of cotton CESA within the cytosolic region. Comparison of the predicted plant CESA structure with the solved structure of a bacterial cellulose-synthesizing protein validates the overall fold of the modeled glycosyltransferase (GT) domain. The coaligned plant and bacterial GT domains share a six-stranded β-sheet, five α-helices, and conserved motifs similar to those required for catalysis in other GT-2 glycosyltransferases. Extending beyond the cross-kingdom similarities related to cellulose polymerization, the predicted structure of cotton CESA reveals that plant-specific modules (plant-conserved region and class-specific region) fold into distinct subdomains on the periphery of the catalytic region. Computational results support the importance of the plant-conserved region and/or class-specific region in CESA oligomerization to form the multimeric cellulose–synthesis complexes that are characteristic of plants. Relatively high sequence conservation between plant CESAs allowed mapping of known mutations and two previously undescribed mutations that perturb cellulose synthesis in Arabidopsis thaliana to their analogous positions in the modeled structure. Most of these mutation sites are near the predicted catalytic region, and the confluence of other mutation sites supports the existence of previously undefined functional nodes within the catalytic core of CESA. Overall, the predicted tertiary structure provides a platform for the biochemical engineering of plant CESAs.

Molecular dynamics simulations of matrix-assisted laser desorption—connections to experiment

The molecular dynamics (MD) simulation technique has been applied to investigate fundamental aspects of matrix-assisted laser desorption. In this paper, we focus on direct comparisons of the results from the simulations with experimental data and on establishing links between the measured or calculated parameters and the basic mechanisms of molecular ejection. The results on the fluence dependence of the ablation/desorption yields and composition of the ejected plume are compared with mass spectrometry and trapping plate experiments. Implications of the prediction of a fluence threshold for ablation are discussed. The strongly forward-peaked velocity and angular distributions of matrix and analyte molecules, predicted in the simulations, are related to the experimental distributions. The shapes and amplitudes of the acoustic waves transmitted from the absorption region through the irradiated sample are compared to recent photoacoustic measurements and related to the ejection mechanisms. The conformational changes during plume evolution and the ejection velocities of analyte molecules are studied and the directions for future investigations are discussed. Finally, we demonstrate that the MD simulation technique can be used to model other processes relevant to mass spectrometry applications, such as laser disintegration of aerosol particles and laser ablation in the presence of photochemical reactions.

Structural Domains within the 3′ Untranslated Region of Turnip Crinkle Virus

ABSTRACT The genomes of positive-strand RNA viruses undergo conformational shifts that complicate efforts to equate structures with function. We have initiated a detailed analysis of secondary and tertiary elements within the 3′ end of Turnip crinkle virus (TCV) that are required for viral accumulation in vivo. MPGAfold, a massively parallel genetic algorithm, suggested the presence of five hairpins (H4a, H4b, and previously identified hairpins H4, H5, and Pr) and one H-type pseudoknot (Ψ3) within the 3′-terminal 194 nucleotides (nt). In vivo compensatory mutagenesis analyses confirmed the existence of H4a, H4b, Ψ3 and a second pseudoknot (Ψ2) previously identified in a TCV satellite RNA. In-line structure probing of the 194-nt fragment supported the coexistence of H4, H4a, H4b, Ψ3 and a pseudoknot that connects H5 and the 3′ end (Ψ1). Stepwise replacements of TCV elements with the comparable elements from Cardamine chlorotic fleck virus indicated that the complete 142-nt 3′ end, and subsets containing Ψ3, H4a, and H4b or Ψ3, H4a, H4b, H5, and Ψ2, form functional domains for virus accumulation in vivo. A new 3-D molecular modeling protocol (RNA2D3D) predicted that H4a, H4b, H5, Ψ3, and Ψ2 are capable of simultaneous existence and bears some resemblance to a tRNA. The related Japanese iris necrotic ring virus does not have comparable domains. These results provide a framework for determining how interconnected elements participate in processes that require 3′ untranslated region sequences such as translation and replication.

Cellulose synthases: new insights from crystallography and modeling

Detailed information about the structure and biochemical mechanisms of cellulose synthase (CelS) proteins remained elusive until a complex containing the catalytic subunit (BcsA) of CelS from Rhodobacter sphaeroides was crystalized. Additionally, a 3D structure of most of the cytosolic domain of a plant CelS (GhCESA1 from cotton, Gossypium hirsutum) was produced by computational modeling. This predicted structure contributes to our understanding of how plant CelS proteins may be similar and different as compared with BcsA. In this review, we highlight how these structures impact our understanding of the synthesis of cellulose and other extracellular polysaccharides. We show how the structures can be used to generate hypotheses for experiments testing mechanisms of glucan synthesis and translocation in plant CelS.

Weakly Charged Cationic Nanoparticles Induce DNA Bending and Strand Separation

Weakly charged cationic nanoparticles cause structural changes including local denaturing and compaction to DNA under mild conditions. The charged ligands bind to the phosphate backbone of DNA and the uncharged ligands penetrate the helix and disrupt base pairing. Mobility shifts in electrophoresis, molecular dynamics, and UV-vis spectrophotometry give clues to the details of the interactions.

Comparative Structural and Computational Analysis Supports Eighteen Cellulose Synthases in the Plant Cellulose Synthesis Complex.

A six-lobed membrane spanning cellulose synthesis complex (CSC) containing multiple cellulose synthase (CESA) glycosyltransferases mediates cellulose microfibril formation. The number of CESAs in the CSC has been debated for decades in light of changing estimates of the diameter of the smallest microfibril formed from the β-1,4 glucan chains synthesized by one CSC. We obtained more direct evidence through generating improved transmission electron microscopy (TEM) images and image averages of the rosette-type CSC, revealing the frequent triangularity and average cross-sectional area in the plasma membrane of its individual lobes. Trimeric oligomers of two alternative CESA computational models corresponded well with individual lobe geometry. A six-fold assembly of the trimeric computational oligomer had the lowest potential energy per monomer and was consistent with rosette CSC morphology. Negative stain TEM and image averaging showed the triangularity of a recombinant CESA cytosolic domain, consistent with previous modeling of its trimeric nature from small angle scattering (SAXS) data. Six trimeric SAXS models nearly filled the space below an average FF-TEM image of the rosette CSC. In summary, the multifaceted data support a rosette CSC with 18 CESAs that mediates the synthesis of a fundamental microfibril composed of 18 glucan chains.

Molecular Description of the LCST Behavior of an Elastin-Like Polypeptide

Elastin-like polypeptides (ELPs) with the repeat sequence of VPGVG are widely used as a model system for investigation of lower critical solution temperature (LCST) transition behavior. In this paper, the effect of temperature on the structure, dynamics and association of (VPGVG)18 in aqueous solution is investigated using atomistic molecular dynamics simulations. Our simulations show that as the temperature increases the ELP backbones undergo gradual conformational changes, which are attributed to the formation of more ordered secondary structures such as β-strands. In addition, increasing temperature changes the hydrophobicity of the ELP by exposure of hydrophobic valine-side chains to the solvent and hiding of proline residues. Based on our simulations, we conclude that the transition behavior of (VPGVG)18 can be attributed to a combination of thermal disruption of the water network that surrounds the polypeptide, reduction of solvent accessible surface area of the polypeptide, and increase in its hydrophobicity. Simulations of the association of two (VPGVG)18 molecules demonstrated that the observed gradual changes in the structural properties of the single polypeptide chain are enough to cause the aggregation of polypeptides above the LCST. These results lead us to propose that the LCST phase behavior of poly(VPGVG) is a collective phenomenon that originates from the correlated gradual changes in single polypeptide structure and the abrupt change in properties of hydration water around the peptide and is a result of a competition between peptide-peptide and peptide-water interactions. This is a computational study of an important intrinsically disordered peptide system that provides an atomic-level description of structural features and interactions that are relevant in the LCST phase behavior.

The role of the photochemical fragmentation in laser ablation: a molecular dynamics study

Despite numerous studies, the mechanistic understanding of the role of the photochemical processes and their coupling with the thermal processes in UV laser ablation is still far from being complete. In this work, the effects of the photochemical reactions on the laser ablation mechanism are delineated based on the results of molecular dynamics simulations of 248 nm laser irradiation of solid chlorobenzene. Photochemical reactions are found to release additional energy into the irradiated sample and decrease the average cohesive energy, therefore decreasing the value of the ablation threshold. The yield of emitted fragments becomes significant only above the ablation threshold. Below the ablation threshold only the most volatile photoproduct, HCl, is ejected in very small amounts, whereas the remainder of photoproducts are trapped inside the sample. Results of the simulations are in a good qualitative agreement with experimental data on the ejection of photoproducts in the laser ablation of chlorobenzene.

Protocols for the In Silico Design of RNA Nanostructures

Recent developments in the field of nanobiology have significantly expanded the possibilities for new modalities in the treatment of many diseases, including cancer. Ribonucleic acid (RNA) represents a relatively new molecular material for the development of these biologically oriented nanodevices. In addition, RNA nanobiology presents a relatively new approach for the development of RNA-based nanoparticles that can be used as crystallization substrates and scaffolds for RNA-based nanoarrays. Presented in this chapter are some methodological shaped-based protocols for the design of such RNA nanostructures. Included are descriptions and background materials describing protocols that use a database of three-dimensional RNA structure motifs; designed RNA secondary structure motifs; and a combination of the two approaches. An example is also given illustrating one of the protocols.