Toward an accurate theoretical framework for describing ensembles for proteins under strongly denaturing conditions

Our focus is on an appropriate theoretical framework for describing highly denatured proteins. In high concentrations of denaturants, proteins behave like polymers in a good solvent and ensembles for denatured proteins can be modeled by ignoring all interactions except excluded volume (EV) effects. To assay conformational preferences of highly denatured proteins, we quantify a variety of properties for EV-limit ensembles of 23 two-state proteins. We find that modeled denatured proteins can be best described as follows. Average shapes are consistent with prolate ellipsoids. Ensembles are characterized by large correlated fluctuations. Sequence-specific conformational preferences are restricted to local length scales that span five to nine residues. Beyond local length scales, chain properties follow well-defined power laws that are expected for generic polymers in the EV limit. The average available volume is filled inefficiently, and cavities of all sizes are found within the interiors of denatured proteins. All properties characterized from simulated ensembles match predictions from rigorous field theories. We use our results to resolve between conflicting proposals for structure in ensembles for highly denatured states.     

The separation between the 5′-3′ ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.     

Unified approach to partition functions of RNA secondary structures

RNA secondary structure formation is a field of considerable biological interest as well as a model system for understanding generic properties of heteropolymer folding. This system is particularly attractive because the partition function and thus all thermodynamic properties of RNA secondary structure ensembles can be calculated numerically in polynomial time for arbitrary sequences and homopolymer models admit analytical solutions. Such solutions for many different aspects of the combinatorics of RNA secondary structure formation share the property that the final solution depends on differences of statistical weights rather than on the weights alone. Here, we present a unified approach to a large class of problems in the field of RNA secondary structure formation. We prove a generic theorem for the calculation of RNA folding partition functions. Then, we show that this approach can be applied to the study of the molten-native transition, denaturation of RNA molecules, as well as to studies of the glass phase of random RNA sequences.     

Structural dynamics of a single-stranded RNA–helix junction using NMR

Many regulatory RNAs contain long single strands (ssRNA) that adjoin secondary structural elements. Here, we use NMR spectroscopy to study the dynamic properties of a 12-nucleotide (nt) ssRNA tail derived from the prequeuosine riboswitch linked to the 3′ end of a 48-nt hairpin. Analysis of chemical shifts, NOE connectivity, ¹³C spin relaxation, and residual dipolar coupling data suggests that the first two residues (A25 and U26) in the ssRNA tail stack onto the adjacent helix and assume an ordered conformation. The following U26-A27 step marks the beginning of an A₆-tract and forms an acute pivot point for substantial motions within the tail, which increase toward the terminal end. Despite substantial internal motions, the ssRNA tail adopts, on average, an A-form helical conformation that is coaxial with the helix. Our results reveal a surprising degree of structural and dynamic complexity at the ssRNA–helix junction, which involves a fine balance between order and disorder that may facilitate efficient pseudoknot formation on ligand recognition.     

RNA–LIM: A novel procedure for analyzing protein/single-stranded RNA propensity data with concomitant estimation of interface structure

RNA–LIM is a procedure that can analyze various pseudo-potentials describing the affinity between single-stranded RNA (ssRNA) ribonucleotides and surface amino acids to produce a coarse-grained estimate of the structure of the ssRNA at the protein interface. The search algorithm works by evolving an ssRNA chain, of known sequence, as a series of walks between fixed sites on a protein surface. Optimal routes are found by application of a set of minimal “limiting” restraints derived jointly from (i) selective sampling of the ribonucleotide amino acid affinity pseudo-potential data, (ii) limited surface path exploration by prior determination of surface arc lengths, and (iii) RNA structural specification obtained from a statistical potential gathered from a library of experimentally determined ssRNA structures. We describe the general approach using a NAST (Nucleic Acid Simulation Tool)-like approximation of the ssRNA chain and a generalized pseudo-potential reflecting the location of nucleic acid binding residues. Minimum and maximum performance indicators of the methodology are established using both synthetic data, for which the pseudo-potential defining nucleic acid binding affinity is systematically degraded, and a representative real case, where the RNA binding sites are predicted by the amplified antisense RNA (aaRNA) method. Some potential uses and extensions of the routine are discussed. RNA–LIM analysis programs along with detailed instructions for their use are available on request from the authors.     

DDX21 RNA解旋酶不同结构域的功能特点解析

DEAD-box家族是在生物体内普遍存在的一类高度保守的RNA解旋酶,在RNA的合成和加工、细胞发育和细胞代谢等过程中都发挥着重要作用。DDX21 RNA解旋酶是DEAD-box家族成员之一,而目前为止DDX21的酶学功能及结构特征尚未被完全了解。本研究运用生物化学与生物物理学前沿技术,系统地研究了DDX21各结构域在不同功能中发挥的作用。首先重组构建并纯化了人的DDX21 RNA解旋酶及不同的截短蛋白质,利用动态激光散射和凝胶层析技术分析各蛋白质的寡聚形态,发现N-端的非功能区（N-端181aa）与C-端的4个FRGQR重复结构域对其结构有较大的影响;利用荧光偏振技术比较分析了各蛋白质与单链RNA的结合反应,结果显示,仅保留DEADc和HELICc结构域的截短蛋白质与单链RNA完全没有亲和性,缺失N-端181aa的截短蛋白质对ssRNA的结合能力与全长蛋白质基本一致,而仅缺失C-端的4个重复FRGQR结构域的截短蛋白质与单链RNA的亲和能力将显著下降;利用快速停流检测技术分析各截短蛋白质的解旋及退火活性,发现DEADc、HELICc及GUCT_RHII三个结构域共同参与DDX21的解旋功能,另一方面,缺失C-端4个FRGQR重复结构域的截短蛋白质导致退火能力的丧失。本研究揭示了DDX21的GUCT_RHII结构域及C-端4个FRGQR重复结构域在其结构及功能中发挥的重要作用,为今后研究DDX21的结构及其细胞功能提供了重要的理论依据。     

The palm subdomain-based active site is internally permuted in viral RNA-dependent RNA polymerases of an ancient lineage

Template-dependent polynucleotide synthesis is catalyzed by enzymes whose core component includes a ubiquitous alphabeta palm subdomain comprising A, B and C sequence motifs crucial for catalysis. Due to its unique, universal conservation in all RNA viruses, the palm subdomain of RNA-dependent RNA polymerases (RdRps) is widely used for evolutionary and taxonomic inferences. We report here the results of elaborated computer-assisted analysis of newly sequenced replicases from Thosea asigna virus (TaV) and the closely related Euprosterna elaeasa virus (EeV), insect-specific ssRNA+ viruses, which revise a capsid-based classification of these viruses with tetraviruses, an Alphavirus-like family. The replicases of TaV and EeV do not have characteristic methyltransferase and helicase domains, and include a putative RdRp with a unique C-A-B motif arrangement in the palm subdomain that is also found in two dsRNA birnaviruses. This circular motif rearrangement is a result of migration of approximately 22 amino acid (aa) residues encompassing motif C between two internal positions, separated by approximately 110 aa, in a conserved region of approximately 550 aa. Protein modeling shows that the canonical palm subdomain architecture of poliovirus (ssRNA+) RdRp could accommodate the identified sequence permutation through changes in backbone connectivity of the major structural elements in three loop regions underlying the active site. This permutation transforms the ferredoxin-like beta1alphaAbeta2beta3alphaBbeta4 fold of the palm subdomain into the beta2beta3beta1alphaAalphaBbeta4 structure and brings beta-strands carrying two principal catalytic Asp residues into sequential proximity such that unique structural properties and, ultimately, unique functionality of the permuted RdRps may result. The permuted enzymes show unprecedented interclass sequence conservation between RdRps of true ssRNA+ and dsRNA viruses and form a minor, deeply separated cluster in the RdRp tree, implying that other, as yet unidentified, viruses may employ this type of RdRp. The structural diversification of the palm subdomain might be a major event in the evolution of template-dependent polynucleotide polymerases in the RNA-protein world.     

Simple rules for efficient assembly predict the layout of a packaged viral RNA

Single-stranded RNA (ssRNA) viruses, which include major human pathogens, package their genomes as they assemble their capsids. We show here that the organization of the viral genomes within the capsids provides intriguing insights into the highly cooperative nature of the assembly process. A recent cryo-electron microscopy structure of bacteriophage MS2, determined with only 5-fold symmetry averaging, has revealed the asymmetric distribution of its encapsidated genome. Here we show that this RNA distribution is consistent with an assembly mechanism that follows two simple rules derived from experiment: (1) the binding of the MS2 maturation protein to the RNA constrains its conformation into a loop, and (2) the capsid must be built in an energetically favorable way. These results provide a new level of insight into the factors that drive efficient assembly of ssRNA viruses in vivo.     

Prediction of protein loop geometries in solution

The ability to determine the structure of a protein in solution is a critical tool for structural biology, as proteins in their native state are found in aqueous environments. Using a physical chemistry based prediction protocol, we demonstrate the ability to reproduce protein loop geometries in experimentally derived solution structures. Predictions were run on loops drawn from (1)NMR entries in the Protein Databank (PDB), and from (2) the RECOORD database in which NMR entries from the PDB have been standardized and re-refined in explicit solvent. The predicted structures are validated by comparison with experimental distance restraints, a test of structural quality as defined by the WHAT IF structure validation program, root mean square deviation (RMSD) of the predicted loops to the original structural models, and comparison of precision of the original and predicted ensembles. Results show that for the RECOORD ensembles, the predicted loops are consistent with an average of 95%, 91%, and 87% of experimental restraints for the short, medium and long loops respectively. Prediction accuracy is strongly affected by the quality of the original models, with increases in the percentage of experimental restraints violated of 2% for the short loops, and 9% for both the medium and long loops in the PDB derived ensembles. We anticipate the application of our protocol to theoretical modeling of protein structures, such as fold recognition methods; as well as to experimental determination of protein structures, or segments, for which only sparse NMR restraint data is available.     

Global conformation of the Escherichia coli replication factor DnaC protein in absence and presence of nucleotide cofactors

Global conformational and oligomeric states of the Escherichia coli replicative factor DnaC protein in the absence and presence of magnesium and nucleotide cofactors, ATP and ADP, and their fluorescent analogues, MANT-ATP and MANT-ADP, have been examined using analytical sedimentation velocity and time-dependent fluorescence anisotropy techniques. In solution, the DnaC protein exists exclusively as a monomer over a large protein concentration range. The value of s(degrees) (20, w)= 2.45 +/- 0.07 S indicates that the protein molecule has an elongated shape. When modeled as a prolate ellipsoid of revolution, the hydrated DnaC protein has an axial ratio of 4.0 +/- 0.6 with long axis a = 112 A and the short axis b = 28 A, respectively. The presence of magnesium or nucleotide cofactors, ATP or ADP, does not affect the global conformation of the protein and its monomeric state. These data indicate that recently found cooperative interactions between the DnaC molecules, in the complex with the DnaB helicase, are induced by the binding to the helicase, i.e., they are not the intrinsic property of the DnaC protein. Fluorescence anisotropy decays of the DnaC-MANT-ATP and DnaC-MANT-ADP complexes indicate that the protein has a rigid global structure on the nanosecond time scale, little affected by the nucleotide cofactors. Nevertheless, the complex with ATP has a more flexible structure, while the complex with ADP is more rigid, with the protein molecule assuming a more elongated shape. Magnesium exerts control only on the complex with the ATP analogue. In the absence of magnesium, the ATP analogue is firmly held in the binding site. In the presence of Mg(2+), this fixed location is released and the analogue is allowed to assume a flexible conformational state. The significance of the results for the functioning of the DnaC protein is discussed.     

HCV NS3 protein helicase domain assists RNA structure conversion

NS3H, the helicase domain of HCV NS3, possesses RNA-stimulated ATPase and ATP hydrolysis-dependent dsRNA unwinding activities. Here, the ability of NS3H to facilitate RNA structural rearrangement is studied using relatively long RNA strands as the model substrates. NS3H promotes intermolecular annealing, resolves three-stranded RNA duplexes, and assists dsRNA and ssRNA inter-conversions to establish a steady state among RNA structures. NS3H facilitates RNA structure conversions in a mode distinct from an ATP-independent RNA chaperone. These findings expand the known function of HCV NS3 helicase and reveal a role for viral helicase in assisting RNA structure conversions during virus life cycle.     

Macrophages recognize streptococci through bacterial single-stranded RNA

Group B streptococcus (GBS) is a leading cause of both neonatal sepsis and meningitis, two diseases that are characterized by inflammation. However, the manner in which GBS organisms are recognized by monocytes and macrophages is poorly understood. In this study, we report that the recognition of GBS and other Gram-positive bacteria by macrophages and monocytes relies on bacterial single-stranded RNA (ssRNA). ssRNA interacts with a signalling complex, which comprises the Toll-like receptor adaptors MyD88 and UNC-93B, but not the established MyD88-dependent ssRNA sensors. The role of ssRNA in the recognition of Gram-positive bacteria--leading to the induction of inflammatory cytokines--has potential implications for sepsis pathogenesis, diagnosis and treatment.     

Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.     

Nanosecond fluorometric investigation of hydrodynamic properties of adenosine triphosphatase from thermophilic bacterium PS3

The soluble portion (TF1) of proton-translocating ATPase from thermophilic bacterium PS3 was labeled with a fluorescent dye N-(1-pyrene)maleimide. The decay of fluorescence anisotropy of the adduct showed that TF1 in aqueous solution was characterized by a volume of equivalent sphere of 1,120 nm3. This value is 2.4 times the volume calculated from the molecular weight and partial specific volume, indicating a non-spherical shape and/or extensive hydration. A prolate ellipsoid with an axial ratio of 2 to 3 is suggested as a first approximation of the shape of hydrated TF1. The presence or absence of ATP, ADP, or Mg2+ did not alter the volume of the equivalent sphere appreciably; the probable conformational change of TF1 induced by these ligands does not lead to a gross alteration of its hydrodynamic properties.     

Binding Site Identification and Flexible Docking of Single Stranded RNA to Proteins Using a Fragment-Based Approach

Protein-RNA docking is hampered by the high flexibility of RNA, and particularly single-stranded RNA (ssRNA). Yet, ssRNA regions typically carry the specificity of protein recognition. The lack of methodology for modeling such regions limits the accuracy of current protein-RNA docking methods. We developed a fragment-based approach to model protein-bound ssRNA, based on the structure of the protein and the sequence of the RNA, without any prior knowledge of the RNA binding site or the RNA structure. The conformational diversity of each fragment is sampled by an exhaustive RNA fragment library that was created from all the existing experimental structures of protein-ssRNA complexes. A systematic and detailed analysis of fragment-based ssRNA docking was performed which constitutes a proof-of-principle for the fragment-based approach. The method was tested on two 8-homo-nucleotide ssRNA-protein complexes and was able to identify the binding site on the protein within 10 Å. Moreover, a structure of each bound ssRNA could be generated in close agreement with the crystal structure with a mean deviation of ~1.5 Å except for a terminal nucleotide. This is the first time a bound ssRNA could be modeled from sequence with high precision.     

RNA multi-structure landscapes

Statistical properties of RNA folding landscapes obtained by the partition function algorithm (McCaskill 1990) are investigated in detail. The pair correlation of free energies as a function of the Hamming distance is used as a measure for the ruggedness of the landscape. The calculation of the partition function contains information about the entire ensemble of secondary structures as a function of temperature and opens the door to all quantities of thermodynamic interest, in contrast with the conventional minimal free energy approach. A metric distance of structure ensembles is introduced and pair correlations at the level of the structures themselves are computed. Just as with landscapes based on most stable secondary structure prediction, the landscapes defined on the full biophysical GCAU alphabet are much smoother than the landscapes restricted to pure GC sequences and the correlation lengths are almost constant fractions of the chain lengths. Correlation functions for multi-structure landscapes exhibit an increased correlation length, especially near the melting temperature. However, the main effect on evolution is rather an effective increase in sampling for finite populations where each sequence explores multiple structures. Correspondence to: P. Schuster     

Cryoelectron microscopy structures of rotavirus NSP2-NSP5 and NSP2-RNA complexes: implications for genome replication

The replication and packaging of the rotavirus genome, comprising 11 segments of double-stranded RNA, take place in specialized compartments called viroplasms, which are formed during infection and involve a coordinated interplay of multiple components. Two rotavirus nonstructural proteins, NSP2 (with nucleoside triphosphatase, single-stranded RNA [ssRNA] binding and helix-destabilizing activities) and NSP5, are essential in these events. Previous structural analysis of NSP2 showed that it is an octamer in crystals, obeying 4-2-2 crystal symmetry, with a large 35-A central hole along the fourfold axis and deep grooves at one of the twofold axes. To ascertain that the solution structure of NSP2 is the same as that in the crystals and investigate how NSP2 interacts with NSP5 and RNA, we carried out single-particle cryoelectron microscopy (cryo-EM) analysis of NSP2 alone and in complexes with NSP5 and ssRNA at subnanometer resolution. Because full-length NSP5 caused severe aggregation upon mixing with NSP2, the deletion construct NSP566-188 was used in these studies. Our studies show that the solution structure of NSP2 is same as the crystallographic octamer and that both NSP566-188 and ssRNA bind to the grooves in the octamer, which are lined by positively charged residues. The fitting of the NSP2 crystal structure to cryo-EM reconstructions of the complexes indicates that, in contrast to the binding of NSP566-188, the binding of RNA induces noticeable conformational changes in the NSP2 octamer. Consistent with the observation that both NSP5 and RNA share the same binding site on the NSP2 octamer, filter binding assays showed that NSP5 competes with ssRNA binding, indicating that one of the functions of NSP5 is to regulate NSP2-RNA interactions during genome replication.