Binding of the Bacteriophage P22 N-Peptide to the boxB RNA Motif Studied by Molecular Dynamics Simulations

Protein-RNA interactions are important for many cellular processes. The Nut-utilization site (N)-protein of bacteriophages contains an N-terminal arginine-rich motif that undergoes a folding transition upon binding to the boxB RNA hairpin loop target structure. Molecular dynamics simulations were used to investigate the dynamics of the P22 N-peptide-boxB complex and to elucidate the energetic contributions to binding. In addition, the free-energy changes of RNA and peptide conformational adaptation to the bound forms, as well as the role of strongly bound water molecules at the peptide-RNA interface, were studied. The influence of peptide amino acid substitutions and the salt dependence of interaction were investigated and showed good agreement with available experimental results. Several tightly bound water molecules were found at the RNA-binding interface in both the presence and absence of N-peptide. Explicit consideration of the waters resulted in shifts of calculated contributions during the energetic analysis, but overall similar binding energy contributions were found. Of interest, it was found that the electrostatic field of the RNA has a favorable influence on the coil-to-α-helix transition of the N-peptide already outside of the peptide-binding site. This result may have important implications for understanding peptide-RNA complex formation, which often involves coupled folding and association processes. It indicates that electrostatic interactions near RNA molecules can lead to a shift in the equilibrium toward the bound form of an interacting partner before it enters the binding pocket.     

The dynamic structural basis of differential enhancement of conformational stability by 5'- and 3'-dangling ends in RNA

Unpaired bases at the end of an RNA duplex (dangling ends) can stabilize the core duplex in a sequence-dependent manner and are important determinants of RNA folding, recognition, and functions. Using 2-aminopurine as a dangling end purine base, we have employed femtosecond time-resolved fluorescence spectroscopy, combined with UV optical melting, to quantitatively investigate the physical and structural nature of the stacking interactions between the dangling end bases and the terminal base pairs. A 3'-dangling purine base has a large subpopulation that stacks on the guanine base of the terminal GC or UG pair, either intrastrand or cross-strand depending on the orientation of the pair, thus providing stabilization of different magnitudes. On the contrary, a 5'-dangling purine base only has a marginal subpopulation that stacks on the purine of the same strand (intrastrand) but has little cross-strand stacking. Thus a 5'-dangling purine does not provide significant stabilization. These stacking structures are not static, and a dangling end base samples a range of stacked and unstacked conformations with respect to the terminal base pair. Femtosecond time-resolved anisotropy decay reveals certain hindered base conformational dynamics that occur on the picosecond to nanosecond time scales, which allow the dangling base to sample these substates. When the dangling purine is opposite to a U and is able to form a potential base pair at the end of the duplex, there is an interplay of base stacking and hydrogen-bonding interactions that depends on the orientation of the base pair relative to the adjacent GC pair. By resolving these populations that are dynamically exchanging on fast time scales, we elucidated the correlation between dynamic conformational distributions and thermodynamic stability.     

Time-resolved fluorescence emission and excitation spectroscopy of d(TA) and d(AT) using synchrotron radiation

The photophysics of the sequence isomers d(TA) and d(AT) has been investigated at room temperature in 5 x 10(-5) M neutral aqueous solution using pulsed ultraviolet excitation from the ACO synchrotron and detection by time correlation or gated single-photon counting. Decay profiles of the emissions at 350, 400 and 460 have been analyzed both independently and globally by reiterative non-linear least-squares fitting to models of two and three independently emitting species. No evidence has been observed for excited-state reaction. Time-windowed spectra, both emission and excitation, have been collected for three time windows and have been deconvoluted to give time-resolved spectra using the lifetimes resulting from the decay analyses. Spectra are separated into two classes, with picosecond and nanosecond lifetimes, respectively. The picosecond spectra have the emission and excitation spectral characteristics of mixed monomer (A and T) fluorescences and are assigned as originating from the unstacked fractions of d(TA) and d(AT). The nanosecond emission spectra from d(TA) and d(AT) are both two-component, with lambda max approximately 350 and approximately 425 nm and lifetimes of 2.3 and 6.1 ns, respectively. The time-resolved excitation spectra for the nanosecond emissions are quite different from the isotropic absorption spectra of d(TA) and d(AT) but correlate with the anisotropic absorption for out-of-plane transitions between stacked bases of co-crystals of 9-methyladenine and 1-methylthymine reported by Stewart and Davidson. The nanosecond spectra thus represent the direct excitation and emission of stacked pairs of bases. These results provide no evidence for energy transfer and are probably related to sequence-specific photo-adduct formation.     

Stacking-unstacking of the dinucleoside monophosphate guanylyl-3'',5''-uridine studied with molecular dynamics.

Molecular dynamics simulations were carried out on two conformations of the dinucleoside monophosphate guanylyl-3',5'-uridine (GpU) in aqueous solution with one sodium counterion. One stacked conformation and one with the C3'-O3'-P-O5' backbone torsion angle twisted 180 degrees to create an unstacked conformation. We observed a relatively stable behavior of the stacked conformation, which remained stacked throughout the simulation, whereas the unstacked conformation showed major changes in the backbone torsion and glycosidic angles. During the simulation the unstacked conformation transformed into a more stacked form and then back again to an unstacked one. The calculated correlation times for rotational diffusion from the molecular dynamics simulations are in agreement with fluorescence anisotropy and nuclear magnetic resonance data. As expected, the correlation times for rotational diffusion of the unstacked conformation were observed to be longer than for the stacked conformation. The 2'OH group may contribute in stabilizing the stacked conformation, where the O2'-H...O4' hydrogen bond occurred in 82.7% of the simulation.     

Base stacking and unstacking as determined from a DNA decamer containing a fluorescent base

Time-resolved fluorescence decay of a single-stranded DNA decamer d(CTGAAT5CAG), where d5 is the fluorescent base 1-(beta-D-2'-deoxyribosyl)-5-methyl-2-pyrimidinone, was measured and analyzed at several temperatures. The d5 base in the decamer is resolved into three states according to their fluorescence decay lifetime characteristics and temperature dependence of their associated amplitudes: fully extended and completely unstacked state, loosely associated state, and fully stacked state. These states are in slow exchange compared to their fluorescence decay rates. The population of the fully extended and completely unstacked state is small and decreases further with increasing temperature. The loosely associated state, whose fluorescence can still be efficiently quenched by other DNA bases, occupies a large portion of the conventionally defined unstacked state. Stacking enthalpy and entropy for the d5 base with thymine or cytosine bases in the DNA decamer are calculated to be -6.6 kcal/mol and -22 cal/mol.K, respectively. This work shows that fluorescent bases in DNA can be useful to the study of local conformations of bases.     

Intrinsic flexibility of snRNA hairpin loops facilitates protein binding

Stem–loop II of U1 snRNA and Stem–loop IV of U2 snRNA typically have 10 or 11 nucleotides in their loops. The fluorescent nucleobase 2-aminopurine was used as a substitute for the adenines in each loop to probe the local and global structures and dynamics of these unusually long loops. Using steady-state and time-resolved fluorescence, we find that, while the bases in the loops are stacked, they are able to undergo significant local motion on the picosecond/nanosecond timescale. In addition, the loops have a global conformational change at low temperatures that occurs on the microsecond timescale, as determined using laser T-jump experiments. Nucleobase and loop motions are present at temperatures far below the melting temperature of the hairpin stem, which may facilitate the conformational change required for specific protein binding to these RNA loops.     

RNA2D3D: A program for Generating,Viewing, and Comparing 3-Dimensional Models of RNA

Abstract

Using primary and secondary structure information of an RNA molecule, the program RNA2D3D automatically and rapidly produces a first-order approximation of a 3-dimensional conformation consistent with this information. Applicable to structures of arbitrary branching complexity and pseudoknot content, it features efficient interactive graphical editing for the removal of any overlaps introduced by the initial generating procedure and for making conformational changes favorable to targeted features and subsequent refinement. With emphasis on fast exploration of alternative 3D conformations, one may interactively add or delete base-pairs, adjacent stems can be coaxially stacked or unstacked, single strands can be shaped to accommodate special constraints, and arbitrary subsets can be defined and manipulated as rigid bodies. Compaction, whereby base stacking within stems is optimally extended into connecting single strands, is also available as a means of strategically making the structures more compact and revealing folding motifs. Subsequent refinement of the first-order approximation, of modifications, and for the imposing of tertiary constraints is assisted with standard energy refinement techniques. Previously determined coordinates for any part of the molecule are readily incorporated, and any part of the modeled structure can be output as a PDB or XYZ file. Illustrative applications in the areas of ribozymes, viral kissing loops, viral internal ribosome entry sites, and nanobiology are presented.     

Sequence analysis of an artificial family of RNA-binding peptides

Diverse peptide sequences recognizing the lambda boxB RNA hairpin were previously isolated from a library encoding the 22-residue lambda N peptide with random amino acids at positions 13-22 using mRNA display. We have statistically analyzed amino acid distributions in 65 unique sequences from rounds 11 and 12 of this selection and evaluated the resulting structural and functional predictions by alanine-scanning mutagenesis and circular dichroism spectrometry. This artificial sequence family has a consensus structure that continues the bent alpha helix of lambda N up to position 17 when bound to lambda boxB. A charge pair (E(14)R(15)) and hydrophobic patch (A(21)L(22) or V(21)L(22)) have important functional roles in this context. Notably, amino acid covariance reveals six specific pairs of random region positions with >95% significant linkage and strong overall helical (i+1, i+3, and i+4) couplings. The covariance analysis suggests that (1) the sequence context of every residue in each insert has been optimized, (2) selected sequences are local optima on a rugged fitness landscape, and (3) it is possible to detect more subtle structural features with artificial protein sequence families than natural homologs. Our results provide a framework for investigating the structures of in vitro selected proteins by functional minimization, reselection, and covariance analysis.     

Photosystem II solubilizes as a monomer by mild detergent treatment of unstacked thylakoid membranes*

We studied the aggregation state of Photosystem II in stacked and unstacked thylakoid membranes from spinach after a quick and mild solubilization with the non-ionic detergent n-dodecyl-α,D-maltoside, followed by analysis by diode-array-assisted gel filtration chromatography and electron microscopy. The results suggest that Photosystem II (PS II) isolates either as a paired, appressed membrane fragment or as a dimeric PS II-LHC II supercomplex upon mild solubilization of stacked thylakoid membranes or PS II grana membranes, but predominantly as a core monomer upon mild solubilization of unstacked thylakoid membranes. Analysis of paired grana membrane fragments reveals that the number of PS II dimers is strongly reduced in single membranes at the margins of the grana membrane fragments. We suggest that unstacking of thylakoid membranes results in a spontaneous disintegration of the PS II-LHC II supercomplexes into separated PS II core monomers and peripheral light-harvesting complexes. This revised version was published online in June 2006 with corrections to the Cover Date.     

Assignment and modeling of the Rev Response Element RNA bound to a Rev peptide using 13C-heteronuclear NMR

Summary The Rev Response Element (RRE) RNA-Rev protein interaction is important for regulation of gene expression in the human immunodeficiency virus. A model system for this interaction, which includes stem IIB of the RRE RNA and an arginine-rich peptide from the RNA-binding domain of Rev, was studied using multidimensional heteronuclear NMR. Assignment of the RNA when bound to the peptide was obtained from NMR experiments utilizing uniformly and specifically ¹³C-labeled RNA. Isotopic filtering experiments on the specifically labeled RNA enabled unambiguous assignment of unusual nonsequential NOE patterns present in the internal loop of the RRE. A three-dimensional model of the RNA in the complex was obtained using restrained molecular dynamics calculations. The internal loop contains two purine-purine base pairs, which are stacked to form one continuous helix flanked by two A-form regions. The formation of a G-G base pair in the internal loop requires an unusual structure of the phosphate backbone. This structural feature is consistent with mutational data as being important for the binding of Rev to the RRE. The G-G base pair may play an important role in opening the normally narrow major groove of A-form RNA to permit binding of the Rev basic domain.     

Nucleotide-mediated conformational changes of monomeric actin and Arp3 studied by molecular dynamics simulations

Members of the actin family of proteins exhibit different biochemical properties when ATP, ADP-P_i, ADP, or no nucleotide is bound. We used molecular dynamics simulations to study the effect of nucleotides on the behavior of actin and actin-related protein 3 (Arp3). In all of the actin simulations, the nucleotide cleft stayed closed, as in most crystal structures. ADP was much more mobile within the cleft than ATP, despite the fact that both nucleotides adopt identical conformations in actin crystal structures. The nucleotide cleft of Arp3 opened in most simulations with ATP, ADP, and no bound nucleotide. Deletion of a C-terminal region of Arp3 that extends beyond the conserved actin sequence reduced the tendency of the Arp3 cleft to open. When the Arp3 cleft opened, we observed multiple instances of partial release of the nucleotide. Cleft opening in Arp3 also allowed us to observe correlated movements of the phosphate clamp, cleft mouth, and barbed-end groove, providing a way for changes in the nucleotide state to be relayed to other parts of Arp3. The DNase binding loop of actin was highly flexible regardless of the nucleotide state. The conformation of Ser14/Thr14 in the P1 loop was sensitive to the presence of the γ-phosphate, but other changes observed in crystal structures were not correlated with the nucleotide state on nanosecond timescales. The divalent cation occupied three positions in the nucleotide cleft, one of which was not previously observed in actin or Arp2/3 complex structures. In sum, these simulations show that subtle differences in structures of actin family proteins have profound effects on their nucleotide-driven behavior.     

The extent to which the spatial separation between photosystems I and II associated with granal formation limits noncyclic electron flow in isolated lettuce chloroplasts

Uncoupled noncyclic electron flow in stacked (granal) chloroplasts with a lateral heterogeneity in the distribution of the two photosystems has been compared with that in unstacked (agranal) chloroplasts with a near-uniform distribution. Chloroplasts were maintained in either structural state in the same assay medium so as to equalize effects of ionic composition which may influence reaction rates. The assay medium, an ion-deficient solution, was capable of supporting high rates of electron flow from water to methyl viologen. At high irradiance, unstacked chloroplasts exhibited an uncoupled rate which was 30% (in chloroplasts isolated from lettuce grown in low light) or 55% (in chloroplasts isolated from lettuce grown in high light) higher than that of stacked chloroplasts; the percentage remained relatively constant in the temperature range 7 to 22 degrees C for both high-light and low-light chloroplasts. At low irradiance, stacked low-light chloroplasts, despite the spatial separation of the two photosystems, gave higher rates of electron flow than did unstacked low-light chloroplasts. The addition of MgCl2 to stacked chloroplasts increased the uncoupled rate of noncyclic electron flow, but only at relatively high irradiances. The differences observed for stacked and unstacked chloroplasts, and for high-light and low-light chloroplasts are discussed. The approach taken in this work should be useful in other comparisons of stacked and unstacked chloroplasts.     

The secondary structure of histone IV and its stability.

The secondary structure within histone IV and its fragments obtained by cyanogen bromide (CNBr) and cleavage at Met 84 has been examined by circular dichroism and spectophotometric pH titration measurements. These studies have confirmed the existence of stable secondary structure within the C-terminal fragment of histone IV (C-peptide which can be perturbed only by 6M urea at pH greater than 8 or 8 M guanidine-HCL. In contrast, the N-terminal fragment (N-peptide) appears to lack significant secondary structure at low ionic strengths but acquires approximately 15% betasheet conformation and 5% alpha-helix upon aggregation at ionic strengths larger than or equal to 0.4. The rates of nitration of the N- and C-peptides by tetranitromethane (TNM) have also been measured as a function of ionic strengths. Under comparable conditions, the rate constant for nitration of the N-peptide was found to be about six times greater than that for the C-peptide, further evidence in support of the presence of stable secondary structure within the C-terminal region of histone IV. After binding these histone IV fragments to DNA, however, the nitration reaction rate constants for the N- and C-peptide in the bound form are found to be 2% and 27% of the corresponding free peptides. Reconstituted nucleohistone IV is about 10% as reactive to TNM as histone IV at comparable ionic strength.     

Alternative conformations at the RNA-binding surface of the N-terminal U2AF(65) RNA recognition motif

The essential pre-mRNA splicing factor, U2 auxiliary factor 65KD (U2AF(65)) recognizes the polypyrimidine tract (Py-tract) consensus sequence of the pre-mRNA using two RNA recognition motifs (RRMs), the most prevalent class of eukaryotic RNA-binding domain. The Py-tracts of higher eukaryotic pre-mRNAs are often interrupted with purines, yet U2AF(65) must identify these degenerate Py-tracts for accurate pre-mRNA splicing. Previously, the structure of a U2AF(65) variant in complex with poly(U) RNA suggested that rearrangement of flexible side-chains or bound water molecules may contribute to degenerate Py-tract recognition by U2AF(65). Here, the X-ray structure of the N-terminal RRM domain of U2AF(65) (RRM1) is described at 1.47 A resolution in the absence of RNA. Notably, RNA-binding by U2AF(65) selectively stabilizes pre-existing alternative conformations of three side-chains located at the RNA interface (Arg150, Lys225, and Arg227). Additionally, a flexible loop connecting the beta2/beta3 strands undergoes a conformational change to interact with the RNA. These pre-existing alternative conformations may contribute to the ability of U2AF(65) to recognize a variety of Py-tract sequences. This rare, high-resolution view of an important member of the RRM class of RNA-binding domains highlights the role of alternative side-chain conformations in RNA recognition.     

A Novel Stable RNA Pentaloop that Interacts Specifically with A Motif Peptide of Lambda-N Protein

To achieve a novel specific peptide–nucleic acid binding model, we designed an in vitro selection procedure to decrease the energetic contribution of the electrostatic interaction in the total binding energy and to increase the contribution of hydrogen bonding and π–π stacking. After the selection of hairpin-loop RNAs that specifically bound to a model peptide of lambda N protein (N peptide), a new thermostable pentaloop RNA motif (N binding thermostable RNA hairpin: NTS RNA) was revealed. The obtained NTS RNA was able to bind to the N peptide with superior specificity to the boxB RNA, which is the naturally occurring partner of the lambda N protein.     

Structural basis for recognition of AU-rich element RNA by the HuD protein

Hu proteins bind to adenosine-uridine (AU)-rich elements (AREs) in the 3' untranslated regions of many short-lived mRNAs, thereby stabilizing them. Here we report the crystal structures of the first two RNA recognition motif (RRM) domains of the HuD protein in complex with an 11-nucleotide fragment of a class I ARE (the c-fos ARE; to 1.8 A), and with an 11-nucleotide fragment of a class II ARE (the tumor necrosis factor alpha ARE; to 2.3 A). These structures reveal a consensus RNA recognition sequence that suggests a preference for pyrimidine-rich sequences and a requirement for a central uracil residue in the clustered AUUUA repeats found in class II AREs. Comparison to structures of other RRM domain-nucleic acid complexes reveals two base recognition pockets in all the structures that interact with bases using residues in conserved ribonucleoprotein motifs and at the C-terminal ends of RRM domains. Different conformations of nucleic acid can be bound by RRM domains by using different combinations of base recognition pockets and multiple RRM domains.