Tertiary Motifs Revealed in Analyses of Higher-Order RNA Junctions

RNA junctions are secondary-structure elements formed when three or more helices come together. They are present in diverse RNA molecules with various fundamental functions in the cell. To better understand the intricate architecture of three-dimensional (3D) RNAs, we analyze currently solved 3D RNA junctions in terms of base-pair interactions and 3D configurations. First, we study base-pair interaction diagrams for solved RNA junctions with 5 to 10 helices and discuss common features. Second, we compare these higher-order junctions to those containing 3 or 4 helices and identify global motif patterns such as coaxial stacking and parallel and perpendicular helical configurations. These analyses show that higher-order junctions organize their helical components in parallel and helical configurations similar to lower-order junctions. Their sub-junctions also resemble local helical configurations found in three- and four-way junctions and are stabilized by similar long-range interaction preferences such as A-minor interactions. Furthermore, loop regions within junctions are high in adenine but low in cytosine, and in agreement with previous studies, we suggest that coaxial stacking between helices likely forms when the common single-stranded loop is small in size; however, other factors such as stacking interactions involving noncanonical base pairs and proteins can greatly determine or disrupt coaxial stacking. Finally, we introduce the ribo-base interactions: when combined with the along-groove packing motif, these ribo-base interactions form novel motifs involved in perpendicular helix-helix interactions. Overall, these analyses suggest recurrent tertiary motifs that stabilize junction architecture, pack helices, and help form helical configurations that occur as sub-elements of larger junction networks. The frequent occurrence of similar helical motifs suggest nature's finite and perhaps limited repertoire of RNA helical conformation preferences. More generally, studies of RNA junctions and tertiary building blocks can ultimately help in the difficult task of RNA 3D structure prediction.     

The role of RNA sequence and structure in RNA--protein interactions

We investigate the sequence and structural properties of RNA-protein interaction sites in 211 RNA-protein chain pairs, the largest set of RNA-protein complexes analyzed to date. Statistical analysis confirms and extends earlier analyses made on smaller data sets. There are 24.6% of hydrogen bonds between RNA and protein that are nucleobase specific, indicating the importance of both nucleobase-specific and -nonspecific interactions. While there is no significant difference between RNA base frequencies in protein-binding and non-binding regions, distinct preferences for RNA bases, RNA structural states, protein residues, and protein secondary structure emerge when nucleobase-specific and -nonspecific interactions are considered separately. Guanine nucleobase and unpaired RNA structural states are significantly preferred in nucleobase-specific interactions; however, nonspecific interactions disfavor guanine, while still favoring unpaired RNA structural states. The opposite preferences of nucleobase-specific and -nonspecific interactions for guanine may explain discrepancies between earlier studies with regard to base preferences in RNA-protein interaction regions. Preferences for amino acid residues differ significantly between nucleobase-specific and -nonspecific interactions, with nonspecific interactions showing the expected bias towards positively charged residues. Irregular protein structures are strongly favored in interactions with the protein backbone, whereas there is little preference for specific protein secondary structure in either nucleobase-specific interaction or -nonspecific interaction. Overall, this study shows strong preferences for both RNA bases and RNA structural states in protein-RNA interactions, indicating their mutual importance in protein recognition.     

Determinants of homodimerization specificity in histidine kinases

Two-component signal transduction pathways consisting of a histidine kinase and a response regulator are used by prokaryotes to respond to diverse environmental and intracellular stimuli. Most species encode numerous paralogous histidine kinases that exhibit significant structural similarity. Yet in almost all known examples, histidine kinases are thought to function as homodimers. We investigated the molecular basis of dimerization specificity, focusing on the model histidine kinase EnvZ and RstB, its closest paralog in Escherichia coli. Direct binding studies showed that the cytoplasmic domains of these proteins each form specific homodimers in vitro. Using a series of chimeric proteins, we identified specificity determinants at the base of the four-helix bundle in the dimerization and histidine phosphotransfer domain. Guided by molecular coevolution predictions and EnvZ structural information, we identified sets of residues in this region that are sufficient to establish homospecificity. Mutating these residues in EnvZ to the corresponding residues in RstB produced a functional kinase that preferentially homodimerized over interacting with EnvZ. EnvZ and RstB likely diverged following gene duplication to yield two homodimers that cannot heterodimerize, and the mutants we identified represent possible evolutionary intermediates in this process.     

Systematic Mutational Analysis of Peptide Inhibition of the p53-MDM2/MDMX Interactions

Inhibition of the interaction between the tumor suppressor protein p53 and its negative regulators MDM2 and MDMX is of great interest in cancer biology and drug design. We previously reported a potent duodecimal peptide inhibitor, termed PMI (TSFAEYWNLLSP), of the p53-MDM2 and -MDMX interactions. PMI competes with p53 for MDM2 and MDMX binding at an affinity roughly 2 orders of magnitude higher than that of ^17-28p53 (ETFSDLWKLLPE) of the same length; both peptides adopt nearly identical α-helical conformations in the complexes, where the three highlighted hydrophobic residues Phe, Trp, and Leu dominate PMI or ^17-28p53 binding to MDM2 and MDMX. To elucidate the molecular determinants for PMI activity and specificity, we performed a systematic Ala scanning mutational analysis of PMI and ^17-28p53. The binding affinities for MDM2 and MDMX of a total of 35 peptides including 10 truncation analogs were quantified, affording a complete dissection of energetic contributions of individual residues of PMI and ^17-28p53 to MDM2 and MDMX association. Importantly, the N8A mutation turned PMI into the most potent dual-specific antagonist of MDM2 and MDMX reported to date, registering respective K_d values of 490 pM and 2.4 nM. The co-crystal structure of N8A-PMI-^25-109MDM2 was determined at 1.95 Å, affirming that high-affinity peptide binding to MDM2/MDMX necessitates, in addition to optimized intermolecular interactions, enhanced helix stability or propensity contributed by non-contact residues. The powerful empirical binding data and crystal structures present a unique opportunity for computational studies of peptide inhibition of the p53-MDM2/MDMX interactions.     

A Structural Model for the DEAD Box Helicase YxiN in Solution: Localization of the RNA Binding Domain

DEAD box proteins consist of a common helicase core formed by two globular RecA domains that are separated by a cleft. The helicase core acts as a nucleotide-dependent switch that alternates between open and closed conformations during the catalytic cycle of duplex separation, thereby providing basic helicase activity. Flanking domains can direct the helicase core to a specific RNA substrate by mediating high-affinity or high-specificity RNA binding. In addition, they may position RNA for the helicase core or may directly contribute to unwinding. While structures of different helicase cores have been determined previously, little is known about the orientation of flanking domains relative to the helicase core.YxiN is a DEAD box protein that consists of a helicase core and a C-terminal RNA binding domain (RBD) that mediates specific binding to hairpin 92 in 23S rRNA. To provide a framework for understanding the functional cooperation of the YxiN helicase core and the RBD, we mapped the orientation of the RBD in single-molecule fluorescence resonance energy transfer experiments. We present a model for the global conformation of YxiN in which the RBD lies above a slightly concave patch that is formed by flexible loops on the surface of the C-terminal RecA domain. The orientation of the RBD is different from the orientations of flanking domains in the Thermus thermophilus DEAD box protein Hera and in Saccharomyces cerevisiae Mss116p, in line with the different functions of these DEAD box proteins and of their RBDs. Interestingly, the corresponding patch on the C-terminal RecA domain that is covered by the YxiN RBD is also part of the interface between the translation factors eIF4A and eIF4G. Possibly, this region constitutes an adaptable interface that generally allows for the interaction of the helicase core with additional domains or interacting factors.     

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active “docked” conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.     

Structural Analysis of Peroxide-Soaked MnSOD Crystals Reveals Side-On Binding of Peroxide to Active-Site Manganese

The superoxide dismutase (SOD) enzymes are important antioxidant agents that protect cells from reactive oxygen species. The SOD family is responsible for catalyzing the disproportionation of superoxide radical to oxygen and hydrogen peroxide. Manganese- and iron-containing SOD exhibit product inhibition whereas Cu/ZnSOD does not. Here, we report the crystal structure of Escherichia coli MnSOD with hydrogen peroxide cryotrapped in the active site. Crystallographic refinement to 1.55 Å and close inspection revealed electron density for hydrogen peroxide in three of the four active sites in the asymmetric unit. The hydrogen peroxide molecules are in the position opposite His26 that is normally assumed by water in the trigonal bipyramidal resting state of the enzyme. Hydrogen peroxide is present in active sites B, C, and D and is side-on coordinated to the active-site manganese. In chains B and D, the peroxide is oriented in the plane formed by manganese and ligands Asp167 and His26. In chain C, the peroxide is bound, making a 70° angle to the plane. Comparison of the peroxide-bound active site with the hydroxide-bound octahedral form shows a shifting of residue Tyr34 towards the active site when peroxide is bound. Comparison with peroxide-soaked Cu/ZnSOD indicates end-on binding of peroxide when the SOD does not exhibit inhibition by peroxide and side-on binding of peroxide in the product-inhibited state of MnSOD.     

Structures of the class D carbapenemase OXA-24 from Acinetobacter baumannii in complex with doripenem

The emergence of class D β-lactamases with carbapenemase activity presents an enormous challenge to health practitioners, particularly with regard to the treatment of infections caused by Gram-negative pathogens such as Acinetobacter baumannii. Unfortunately, class D β-lactamases with carbapenemase activity are resistant to β-lactamase inhibitors. To better understand the details of the how these enzymes bind and hydrolyze carbapenems, we have determined the structures of two deacylation-deficient variants (K84D and V130D) of the class D carbapenemase OXA-24 with doripenem bound as a covalent acyl-enzyme intermediate. Doripenem adopts essentially the same configuration in both OXA-24 variant structures, but varies significantly when compared to the non-carbapenemase class D member OXA-1/doripenem complex. The alcohol of the 6α hydroxyethyl moiety is directed away from the general base carboxy-K84, with implications for activation of the deacylating water. The tunnel formed by the Y112/M223 bridge in the apo form of OXA-24 is largely unchanged by the binding of doripenem. The presence of this bridge, however, causes the distal pyrrolidine/sulfonamide group to bind in a drastically different conformation compared to doripenem bound to OXA-1. The resulting difference in the position of the side-chain bridge sulfur of doripenem is consistent with the hypothesis that the tautomeric state of the pyrroline ring contributes to the different carbapenem hydrolysis rates of OXA-1 and OXA-24. These findings represent a snapshot of a key step in the catalytic mechanism of an important class D enzyme, and might be useful for the design of novel inhibitors.     

The importance of G.A hydrogen bonding in the metal ion- and protein-induced folding of a kink turn RNA

The kink turn (K-turn) is a common motif in RNA structure, found in many RNA species important in translation, RNA modification and splicing, and the control of gene expression. In general the K-turn comprises a three nucleotide bulge followed by trans sugar-Hoogsteen G·A pairs. The RNA adopts a tightly kinked conformation, and is a common target for binding proteins, exemplified by the L7Ae family. We have measured the rates of association and dissociation for the binding of L7Ae to the Kt-7 kink turn, from which we calculate an affinity of K_D = 10 pM. This high affinity is consistent with the role of this binding as the first stage in the assembly of key functional nucleoproteins such as box C/D snoRNP. Kink-turn RNA undergoes a two-state transition between the kinked conformation, and a more extended structure, and folding into the kinked form is induced by divalent metal ions, or by binding of proteins of the L7Ae class. The K-turn provides an excellent, simple model for RNA folding, which can be dissected at the atomic level. We have analyzed the contributions of the hydrogen bonds that form the G·A pairs to the ion- and protein-induced folding of the K-turn. We find that all four hydrogen bonds are important to the stability of the kinked form of the RNA, and we can now define all the important hydrogen bonding interactions that stabilize the K-turn. The high affinity of L7Ae binding is coupled to the induced folding of the K-turn, allowing some sub-optimal variants to adopt the kinked geometry. However, in all such cases the affinity is lowered, and the results underline the importance of both G·A pairs to the stability of the K-turn.     

Refinement of Protein Structures into Low-Resolution Density Maps Using Rosetta

We describe a method based on Rosetta structure refinement for generating high-resolution, all-atom protein models from electron cryomicroscopy density maps. A local measure of the fit of a model to the density is used to directly guide structure refinement and to identify regions incompatible with the density that are then targeted for extensive rebuilding. Over a range of test cases using both simulated and experimentally generated data, the method consistently increases the accuracy of starting models generated either by comparative modeling or by hand-tracing the density. The method can achieve near-atomic resolution starting from density maps at 4-6 Å resolution.     

Active-Site Monovalent Cations Revealed in a 1.55-Å-Resolution Hammerhead Ribozyme Structure

We have obtained a 1.55-Å crystal structure of a hammerhead ribozyme derived from Schistosoma mansoni under conditions that permit detailed observations of Na⁺ ion binding in the ribozyme's active site. At least two such Na⁺ ions are observed. The first Na⁺ ion binds to the N7 of G10.1 and the adjacent A9 phosphate in a manner identical with that previously observed for divalent cations. A second Na⁺ ion binds to the Hoogsteen face of G12, the general base in the hammerhead cleavage reaction, thereby potentially dissipating the negative charge of the catalytically active enolate form of the nucleotide base. A potential but more ambiguous third site bridges the A9 and scissile phosphates in a manner consistent with that of previous predictions. Hammerhead ribozymes have been observed to be active in the presence of high concentrations of monovalent cations, including Na⁺, but the mechanism by which monovalent cations substitute for divalent cations in hammerhead catalysis remains unclear. Our results enable us to suggest that Na⁺ directly and specifically substitutes for divalent cations in the hammerhead active site. The detailed geometry of the pre-catalytic active-site complex is also revealed with a new level of precision, thanks to the quality of the electron density maps obtained from what is currently the highest-resolution ribozyme structure in the Protein Data Bank.     

Disordered domains and high surface charge confer hubs with the ability to interact with multiple proteins in interaction networks

We investigate the structural properties of hubs that enable them to interact with several partners in protein-protein interaction networks. We find that hubs have more observed and predicted disordered residues with fewer loops/coils, and more charged residues on the surface as compared to non-hubs. Smaller hubs have fewer disordered residues and more charged residues on the surface than larger hubs. We conclude that the global flexibility provided by disordered domains, and high surface charge are complementary factors that play a significant role in the binding ability of hubs.     

A Structural Portrait of the PDZ Domain Family

PDZ (PSD-95/Discs-large/ZO1) domains are interaction modules that typically bind to specific C-terminal sequences of partner proteins and assemble signaling complexes in multicellular organisms. We have analyzed the existing database of PDZ domain structures in the context of a specificity tree based on binding specificities defined by peptide-phage binding selections. We have identified 16 structures of PDZ domains in complex with high-affinity ligands and have elucidated four additional structures to assemble a structural database that covers most of the branches of the PDZ specificity tree. A detailed comparison of the structures reveals features that are responsible for the diverse specificities across the PDZ domain family. Specificity differences can be explained by differences in PDZ residues that are in contact with the peptide ligands, but these contacts involve both side-chain and main-chain interactions. Most PDZ domains bind peptides in a canonical conformation in which the ligand main chain adopts an extended β-strand conformation by interacting in an antiparallel fashion with a PDZ β-strand. However, a subset of PDZ domains bind peptides with a bent main-chain conformation and the specificities of these non-canonical domains could not be explained based on canonical structures. Our analysis provides a structural portrait of the PDZ domain family, which serves as a guide in understanding the structural basis for the diverse specificities across the family.     

The capsid of the small RNA phage PRR1 is stabilized by metal ions

Many nonenveloped virus particles are stabilized by calcium ions bound in the interfaces between the protein subunits. These ions may have a role in the disassembly process. The small RNA phages of the Leviviridae family have T = 3 quasi-symmetry and are unique among simple viruses in that they have a coat protein with a translational repressor activity and a fold that has not been observed in other viruses. The crystal structure of phage PRR1 has been determined to 3.5 Å resolution. The structure shows a tentative binding site for a calcium ion close to the quasi-3-fold axis. The RNA-binding surface used for repressor activity is mostly conserved. The structure does not show any significant differences between quasi-equivalent subunits, which suggests that the assembly is not controlled by conformational switches as in many other simple viruses.     

On the Structure of the Proton-Binding Site in the Fo Rotor of Chloroplast ATP Synthases

The recently reported crystal structures of the membrane-embedded proton-dependent c-ring rotors of a cyanobacterial F₁F_o ATP synthase and a chloroplast F₁F_o ATP synthase have provided new insights into the mechanism of this essential enzyme. While the overall features of these c-rings are similar, a discrepancy in the structure and hydrogen-bonding interaction network of the H⁺ sites suggests two distinct binding modes, potentially reflecting a mechanistic differentiation. Importantly, the conformation of the key glutamate side chain to which the proton binds is also altered. To investigate the nature of these differences, we use molecular dynamics simulations of both c-rings embedded in a phospholipid membrane. We observe that the structure of the c₁₅ ring from Spirulina platensis is unequivocally stable within the simulation time. By contrast, the proposed structure of the H⁺ site in the chloroplast c₁₄ ring changes rapidly and consistently into that reported for the c₁₅ ring, indicating that the latter represents a common binding mode. To assess this hypothesis, we have remodeled the c₁₄ ring by molecular replacement using the published structure factors. The resulting structure provides clear evidence in support of a common binding site conformation and is also considerably improved statistically. These findings, taken together with a sequence analysis of c-subunits in the ATP synthase family, indicate that the so-called proton-locked conformation observed in the c₁₅ ring may be a common characteristic not only of light-driven systems such as chloroplasts and cyanobacteria but also of a selection of other bacterial species.     

Rotations of the 2B sub-domain of E. coli UvrD helicase/translocase coupled to nucleotide and DNA binding

Escherichia coli UvrD is a superfamily 1 DNA helicase and single-stranded DNA (ssDNA) translocase that functions in DNA repair and plasmid replication and as an anti-recombinase by removing RecA protein from ssDNA. UvrD couples ATP binding and hydrolysis to unwind double-stranded DNA and translocate along ssDNA with 3′-to-5′ directionality. Although a UvrD monomer is able to translocate along ssDNA rapidly and processively, DNA helicase activity in vitro requires a minimum of a UvrD dimer. Previous crystal structures of UvrD bound to a ssDNA/duplex DNA junction show that its 2B sub-domain exists in a “closed” state and interacts with the duplex DNA. Here, we report a crystal structure of an apo form of UvrD in which the 2B sub-domain is in an “open” state that differs by an ∼ 160° rotation of the 2B sub-domain. To study the rotational conformational states of the 2B sub-domain in various ligation states, we constructed a series of double-cysteine UvrD mutants and labeled them with fluorophores such that rotation of the 2B sub-domain results in changes in fluorescence resonance energy transfer. These studies show that the open and closed forms can interconvert in solution, with low salt favoring the closed conformation and high salt favoring the open conformation in the absence of DNA. Binding of UvrD to DNA and ATP binding and hydrolysis also affect the rotational conformational state of the 2B sub-domain, suggesting that 2B sub-domain rotation is coupled to the function of this nucleic acid motor enzyme.     

The conserved asparagine in the HNH motif serves an important structural role in metal finger endonucleases

The HNH motif is a small nucleic acid binding and cleavage module, widespread in metal finger endonucleases in all life kingdoms. Here we studied a non-specific endonuclease, the nuclease domain of ColE7 (N-ColE7), to decipher the role of the conserved asparagine and histidine residues in the HNH motif. We found, using fluorescence resonance energy transfer (FRET) assays, that the DNA hydrolysis activity of H545 N-ColE7 mutants was completely abolished while activities of N560 and H573 mutants varied from 6.9% to 83.2% of the wild-type activity. The crystal structures of three N-ColE7 mutants in complex with the inhibitor Im7, N560A-Im7, N560D-Im7 and H573A-Im7, were determined at a resolution of 1.9 A to 2.2 A. H573 is responsible for metal ion binding in the wild-type protein, as the zinc ion is still partially associated in the structure of H573A, suggesting that H573 plays a supportive role in metal binding. Both N560A and N560D contain a disordered loop in the HNH motif due to the disruption of the hydrogen bond network surrounding the side-chain of residue 560, and as a result, the imidazole ring of the general base residue H545 is tilted slightly and the scissile phosphate is shifted, leading to the large reductions in hydrolysis activities. These results suggest that the highly conserved asparagine in the HNH motif, in general, plays a structural role in constraining the loop in the metal finger structure and keeping the general base histidine and scissile phosphate in the correct position for DNA hydrolysis.     

The evolution of structural databases

Starting with the Protein Data Bank (PDB) as a common ancestor, the evolution of structural databases has been driven by the rapprochement of the structural world and the practical applications. The result is an impressive number of secondary structural databases that is welcomed by structural biologists and bioinformaticians but runs the risk of producing an embarrassment of riches among non-specialist users. Given that any profit depends on the number of customers, efficient interfaces between many structural data banks must be available to make their contents easily accessible. Increasing the information content of central structural repositories might be the best way to guide users through the many, sometimes overlapping databases.     

Mapping of protein-protein interaction sites by the 'absence of interference' approach

Protein-protein interactions are critical to most biological processes, and locating protein-protein interfaces on protein structures is an important task in molecular biology. We developed a new experimental strategy called the ‘absence of interference’ approach to determine surface residues involved in protein-protein interaction of established yeast two-hybrid pairs of interacting proteins. One of the proteins is subjected to high-level randomization by error-prone PCR. The resulting library is selected by yeast two-hybrid system for interacting clones that are isolated and sequenced. The interaction region can be identified by an absence or depletion of mutations. For data analysis and presentation, we developed a Web interface that analyzes the mutational spectrum and displays the mutational frequency on the surface of the structure (or a structural model) of the randomized protein†. Additionally, this interface might be of use for the display of mutational distributions determined by other types of random mutagenesis experiments. We applied the approach to map the interface of the catalytic domain of the DNA methyltransferase Dnmt3a with its regulatory factor Dnmt3L. Dnmt3a was randomized with high mutational load. A total of 76 interacting clones were isolated and sequenced, and 648 mutations were identified. The mutational pattern allowed to identify a unique interaction region on the surface of Dnmt3a, which comprises about 500-600 Å². The results were confirmed by site-directed mutagenesis and structural analysis. The absence-of-interference approach will allow high-throughput mapping of protein interaction sites suitable for functional studies and protein docking.