E. coli RNA polymerase interacts homologously with two different promoters

We present and review experiments that identify points of close approach of the RNA polymerase to two promoters, lac UV5 and T7 A3. We identify the contacts to the phosphates along the DNA backbone, to the N7s of guanines in the major groove and the N3s of adenines in the minor groove, and to the methyl groups of thymines. These contacts to the two promoters are strikingly homologous in space, as shown on three-dimensional models, and identify major regions of interactions lying on one side of the DNA molecule (at -35 and -16), as well as further areas extending through the Pribnow box. Both promoters are unwound similarly by the polymerase, across a region of about twelve bases extending from the middle of the Pribnow box to just beyond the RNA start site. We discuss the areas of interaction in the context of promoter homologies and promoter mutations. The disposition of the contacts in space suggests a model for the pathway along which the RNA polymerase binds to promoters.     

A model for the LexA repressor DNA complex

A structural model for the interaction of the LexA repressor DNA binding domain (DBD) with operator DNA is derived by means of Monte Carlo docking. Protein–DNA complexes were generated by docking the LexA repressor DBD NMR solution structure onto both rigid and bent B-DNA structures while giving energy bonuses for contacts in agreement with experimental data. In the resulting complexes, helix III of the LexA repressor DBD is located in the major groove of the DNA and residues Asn-41, Glu-44, and Glu-45 form specific hydrogen bonds with bases of the CTGT DNA sequence. Ser-39, Ala-42, and Asn-41 are involved in a hydrophobic interaction with the methyl group of the first thymine base. Residues in the loop region connecting the two β-sheet strands are involved in nonspecific contacts near the dyad axis of the operator. The contacts observed in the docked complexes cover the entire consensus CTGT half-site DNA operator, thus explaining the specificity of the LexA repressor for such sequences. In addition, a large number of nonspecific interactions between protein and DNA is observed. The agreement between the derived model for the LexA repressor DBD/DNA complex and experimental biochemical results is discussed. © 1995 Wiley-Liss, Inc.     

Analysis of mixed DNA‐bisnaphthalimide interactions involving groove association and intercalation with surface‐based and solution methodologies

The bisnaphthalimide cytotoxic agent elinafide exhibits a mixed DNA binding mode including groove‐association and intercalation. We have compared the interaction of elinafide and two bisnaphthalimide analogues with various natural and modified DNA sequences using solution NMR and UV‐melting methods and surface plasmon resonance (SPR) experiments at different pH conditions. The combined data obtained with these techniques established a high‐affinity binding mode comprising intercalation and strong electrostatic contacts with guanine bases in the major groove, and a weaker interaction with A·T pairs likely involving groove association. However, the SPR binding constants and the NMR and UV‐melting binding parameters responded differently to variations in DNA bases and ligand intercalating moieties. The rates and equilibrium constants determined by SPR clearly responded to changes in pH and DNA groove composition, but were rather insensitive to alterations in drug rings and DNA bases affecting the intercalation process. Conversely, the intermolecular stacking interactions detected by NMR and the ligand‐induced thermal stabilizations measured by UV depended on both sets of factors and were controlled by the sequence‐dependent properties of the DNA helices, indicating that these data were modulated by naphthalimide stacking in addition to groove association. A two‐step binding process where a groove‐bound state is required prior to intercalation is proposed as an explanation for these observations. These findings may be useful for studying other classes of DNA‐ and RNA‐binding drugs, which frequently combine groove‐binding and stacking moieties. © 2012 Wiley Periodicals, Inc. Biopolymers 97:974–987, 2012.     

Hydrogen and hydration of DNA and RNA oligonucleotides

In this paper, hydrogen bonding interaction and hydration in crystal structures of both DNA and RNA oligonucleotides are discussed. Their roles in the formation and stabilization of oligonucleotides have been covered. Details of the Watson-Crick base pairs G.C and A.U in DNA and RNA are illustrated. The geometry of the wobble (mismatched) G.U base pairs and the cis and almost trans conformations of the mismatched U.U base pairs in RNA is described. The difference in hydration of the Watson-Crick base pairs G.C, A.U and the wobble G.U in different sequences of codon-anticodon interaction in double helical molecules are indicative of the effect of hydration. The hydration patterns of the phosphate, the 2'-hydroxyl groups, the water bridges linking the phosphate group, N7 (purine) and N4 of Cs or O4 of Us in the major groove, the water bridges between the 2'-hydroxyl group and N3 (purine) and O2 (pyrimidine) in the minor groove are discussed.     

Hybrid simulation approach incorporating microscopic interaction along with rigid body degrees of freedom for stacking between base pairs

Stacking interaction between the aromatic heterocyclic bases plays an important role in the double helical structures of nucleic acids. Considering the base as rigid body, there are total of 18 degrees of freedom of a dinucleotide step. Some of these parameters show sequence preferences, indicating that the detailed atomic interactions are important in the stacking. Large variants of non‐canonical base pairs have been seen in the crystallographic structures of RNA. However, their stacking preferences are not thoroughly deciphered yet from experimental results. The current theoretical approaches use either the rigid body degrees of freedom where the atomic information are lost or computationally expensive all atom simulations. We have used a hybrid simulation approach incorporating Monte‐Carlo Metropolis sampling in the hyperspace of 18 stacking parameters where the interaction energies using AMBER‐parm99bsc0 and CHARMM‐36 force‐fields were calculated from atomic positions. We have also performed stacking energy calculations for structures from Monte‐Carlo ensemble by Dispersion corrected density functional theory. The available experimental data with Watson–Crick base pairs are compared to establish the validity of the method. Stacking interaction involving A:U and G:C base pairs with non‐canonical G:U base pairs also were calculated and showed that these structures were also sequence dependent. This approach could be useful to generate multiscale modeling of nucleic acids in terms of coarse‐grained parameters where the atomic interactions are preserved. This method would also be useful to predict structure and dynamics of different base pair steps containing non Watson–Crick base pairs, as found often in the non‐coding RNA structures. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 212–226, 2016.     

Determination of the three-dimensional structure of the Mrf2-DNA complex using paramagnetic spin labeling

Understanding the mechanism of protein-DNA interactions at the molecular level is one of the main focuses in structural and molecular biological investigations. At present, NMR spectroscopy is the only approach that can provide atomic details of protein-DNA recognition in solution. However, determining the structures of protein-DNA complexes using NMR spectroscopy has been dependent on the observation of intermolecular nuclear Overhauser effects (NOE) and their assignments, which are difficult to obtain in many cases. In this study, we have shown that intermolecular distance constraints derived from a single spin-label in combination with docking calculations have defined many specific contacts of the complex between the AT-rich interaction domain (ARID) of Mrf2 and its target DNA. Mrf2 contacts DNA mainly using the two flexible loops, L1 and L2. While the L1 loop contacts the phosphate backbone, L2 and several residues in the adjacent helices interact with AT base pairs in the major groove of DNA. Despite the structural diversity in the ARID family of DNA-binding proteins, Mrf2 maintains contacts with DNA similar to those observed in the homologous Dri-DNA complex.     

Protein-RNA interactions: exploring binding patterns with a three-dimensional superposition analysis of high resolution structures

MOTIVATION: The recognition of specific RNA sequences and structures by proteins is critical to our understanding of RNA processing, gene expression and viral replication. The diversity of RNA structures suggests that RNA recognition is substantially different than that of DNA. RESULTS: The atomic coordinates of 41 protein-RNA complexes have been used to probe composite nucleoside binding pockets that form the structural and chemical underpinnings of base recognition. Composite nucleoside binding pockets were constructed using three-dimensional superpositions of each RNA nucleoside. Unlike protein-DNA interactions which are dominated by accessibility, RNA recognition frequently occurs in non-canonical and single-strand-like structures that allow interactions to occur from a much wider set of geometries and make fuller use of unique base shapes and hydrogen-bonding ability. By constructing composites that include all van der Waals, hydrogen-bonding, stacking and general non-polar interactions made to a particular nucleoside, the strategies employed are made readily visible. Protein-RNA interactions can result in the formation of a glove-like tight binding pocket around RNA bases, but the size, shape and non-polar binding patterns differ between specific RNA bases. We show that adenine can be distinguished from guanine based on the size and shape of the binding pocket and steric exclusion of the guanine N2 exocyclic amino group. The unique shape and hydrogen-bonding pattern for each RNA base allow proteins to make specific interactions through a very small number of contacts, as few as two in some cases. AVAILABILITY: The program ENTANGLE is available from http://www.bioc.rice.edu/~shamoo     

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.     

1 A crystal structures of B-DNA reveal sequence-specific binding and groove-specific bending of DNA by magnesium and calcium

The 1 A resolution X-ray crystal structures of Mg(2+) and Ca(2+) salts of the B-DNA decamers CCAACGTTGG and CCAGCGCTGG reveal sequence-specific binding of Mg(2+) and Ca(2+) to the major and minor grooves of DNA, as well as non-specific binding to backbone phosphate oxygen atoms. Minor groove binding involves H-bond interactions between cross-strand DNA base atoms of adjacent base-pairs and the cations' water ligands. In the major groove the cations' water ligands can interact through H-bonds with O and N atoms from either one base or adjacent bases, and in addition the softer Ca(2+) can form polar covalent bonds bridging adjacent N7 and O6 atoms at GG bases. For reasons outlined earlier, localized monovalent cations are neither expected nor found.Ultra-high atomic resolution gives an unprecedented view of hydration in both grooves of DNA, permits an analysis of individual anisotropic displacement parameters, and reveals up to 22 divalent cations per DNA duplex. Each DNA helix is quite anisotropic, and alternate conformations, with motion in the direction of opening and closing the minor groove, are observed for the sugar-phosphate backbone. Taking into consideration the variability of experimental parameters and crystal packing environments among these four helices, and 24 other Mg(2+) and Ca(2+) bound B-DNA structures, we conclude that sequence-specific and strand-specific binding of Mg(2+) and Ca(2+) to the major groove causes DNA bending by base-roll compression towards the major groove, while sequence-specific binding of Mg(2+) and Ca(2+) in the minor groove has a negligible effect on helix curvature. The minor groove opens and closes to accommodate Mg(2+) and Ca(2+) without the necessity for significant bending of the overall helix.The program Shelxdna was written to facilitate refinement and analysis of X-ray crystal structures by Shelxl-97 and to plot and analyze one or more Curves and Freehelix output files.     

Structure of a photoactive rhodium complex intercalated into DNA

Intercalating complexes of rhodium(III) are strong photo-oxidants that promote DNA strand cleavage or electron transfer through the double helix. The 1.2 A resolution crystal structure of a sequence-specific rhodium intercalator bound to a DNA helix provides a rationale for the sequence specificity of rhodium intercalators. It also explains how intercalation in the center of an oligonucleotide modifies DNA conformation. The rhodium complex intercalates via the major groove where specific contacts are formed with the edges of the bases at the target site. The phi ligand is deeply inserted into the DNA base pair stack. The primary conformational change of the DNA is a doubling of the rise per residue, with no change in sugar pucker from B-form DNA. Based upon the five crystallographically independent views of an intercalated DNA helix observed in this structure, the intercalator may be considered as an additional base pair with specific functional groups positioned in the major groove.     

Mutations in 5S DNA and 5S RNA have different effects on the binding of Xenopus transcription factor IIIA

The effects on TFIIIA binding affinity of a series of substitution mutations in the Xenopus laevis oocyte 5S RNA gene were quantified. These data indicate that TFIIIA binds specifically to 5S DNA by forming sequence-specific contacts with three discrete sites located within the classical A and C boxes and the intermediate element of the internal control region. Substitution of the nucleotide sequence at any of the three sites significantly reduces TFIIIA binding affinity, with a 100-fold reduction observed for substitutions in the box C subregion. These results are consistent with a direct interaction of TFIIIA with specific base pairs within the major groove of the DNA. A comparison of the TFIIIA binding data for the same mutations expressed in 5S RNA indicates that the protein does not make any strong sequence-specific contacts with the RNA. Although the protein footprinting sites on the 5S DNA and 5S RNA are coincident, nucleotide substitutions in 5S RNA which moderately reduce TFIIIA binding affinity do not correspond at all to the three specific TFIIIA interaction sites within the gene. The implications of these results for models which attempt to reconcile the DNA and RNA binding activities of TFIIIA by proposing a common structural motif for the two nucleic acids are discussed.     

The P1 phage replication protein RepA contacts an otherwise inaccessible thymine N3 proton by DNA distortion or base flipping

The RepA protein from bacteriophage P1 binds DNA to initiate replication. RepA covers one face of the DNA and the binding site has a completely conserved T that directly faces RepA from the minor groove at position +7. Although all four bases can be distinguished through contacts in the major groove of B-form DNA, contacts in the minor groove cannot easily distinguish between A and T bases. Therefore the 100% conservation at this position cannot be accounted for by direct contacts approaching into the minor groove of B-form DNA. RepA binding sites with modified base pairs at position +7 were used to investigate contacts with RepA. The data show that RepA contacts the N3 proton of T at position +7 and that the T=A hydrogen bonds are already broken in the DNA before RepA binds. To accommodate the N3 proton contact the T₊₇ /A_+7^′ base pair must be distorted. One possibility is that T₊₇ is flipped out of the helix. The energetics of the contact allows RepA to distinguish between all four bases, accounting for the observed high sequence conservation. After protein binding, base pair distortion or base flipping could initiate DNA melting as the second step in DNA replication.     

Flexibility of the DNA enhances promoter affinity of Escherichia coli RNA polymerase.

Two types of mechanisms are discussed for the formation of active protein-DNA complexes: contacts with specific bases and interaction via specific DNA structures within the cognate DNA. We have studied the effect of a single nucleoside deletion on the interaction of Escherichia coli RNA polymerase with a strong promoter. This study reveals three patterns of interaction which can be attributed to different sites of the promoter, (i) direct base contact with the template strand in the '-35 region' (the 'recognition domain'), (ii) a DNA structure dependent interaction in the '-10 region' (the 'melting domain'), and (iii) an interaction which is based on a defined spatial relationship between the two domains of a promoter, namely the 'recognition domain' and the 'melting domain'.     

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.     

Crystal structure of an RNA quadruplex containing inosine tetrad: implications for the roles of NH2 group in purine tetrads

Polyinosinic acid has been known to adopt the four-stranded helical structure but its basic unit, inosine tetrad (I tetrad), has not been determined at the atomic level. Here we report the crystal structure of an RNA quadruplex containing an I tetrad at 1.4 A resolution. The I tetrad has one cyclic hydrogen bond N1...O6 with the bond length of 2.7 A. A water bridge is observed in the minor groove side of the base tetrad. Even though it is sandwiched by guanine tetrads (G tetrads), the I tetrad is buckled towards the 3' side of the tetrad plane, which results from the different interaction strength with K ions on two sides of the tetrad plane. Comparison with both G tetrad and adenine tetrad indicates that lack of NH2 in the C2 position makes the I tetrad prone to buckle for interactions with ligands. Two U*(G-G-G-G) base pentads are observed at the junction of the 5' termini of two quadruplexes. The uridine residue in the base pentad is engaged in two hydrogen bonding interactions (N2(G)-H...O2(U) and O2'(G)-H...O4(U)) and a water-mediated interaction (N3(G) and N3(U)) with the G tetrad. We also discuss the roles of amino group in purine tetrads and the inter-quadruplex interactions in RNA molecules. These quadruplexes may interact with each other by stacking, groove binding and intercalation.     

Recent advances in RNA-protein interaction studies

Conclusions The RNA binding sites for several small proteins have been characterised. These sites include double helical regions with hairpins, bulged bases and internal loops. As seen in Flock House virus structure, some proteins may recognise phosphate backbone of the canonical A-form helix not in a sequence-specific manner. If sequence-specific base contacts are to be made, then the A-helic major groove must be widened. This can be accomplished by introducing bulges, internal loops and hairpin loops into double helical regions. In these cases proteins may recognise both distorted backbone conformations and read out base sequences in a widened major groove. Crystallographic studies on complexes of aminoacyl-tRNA synthetase and tRNA showed that even RNAs with stable tertiary fold undergo substantial structural changes upon binding to the synthetases. The structural variability of RNA as well as the ability of RNA to distort upon protein binding may be crucial in RNA-protein interactions.     

The role of lysine 55 in determining the specificity of the purine repressor for its operators through minor groove interactions.

The interaction of the dimeric Escherichia coli purine repressor (PurR) with its cognate sequences leads to a 45 degrees to 50 degrees kink at a central CpG base step towards the major groove, as dyad-related leucine side-chains interdigitate between these bases from the minor groove. The resulting broadening of the minor groove increases the accessibility of the six central base-pairs towards minor groove interactions with residues from PurR. It has been shown that lysine 55 of PurR makes a direct contact with the adenine base (Ade8) directly 5' to the central CpG base-pair step in the high-affinity purF operator sequence. We have investigated the importance of this interaction in the specificity and affinity of wild-type PurR (WT) for its operators and we have studied a mutant of PurR in which Lys55 is replaced with alanine (K55A). Complexes of WT and K55A with duplex DNA containing pur operator sequences varied at position 8 were investigated crystallographically, and binding studies were performed using fluorescence anisotropy. The structures of the protein-DNA complexes reveal a relatively unperturbed global conformation regardless of the identity of the base-pair at position 8 or residue 55. In all structures the combination of higher resolution and a palindromic purF operator site allowed several new PurR.DNA interactions to be observed, including contacts by Thr15, Thr16 and His20. The side-chain of Lys55 makes productive, though varying, interactions with the adenine, thymine or cytosine base at position 8 that result in equilibrium dissociation constants of 2.6 nM, 10 nM and 35 nM, respectively. However, the bulk of the lysine side-chain apparently blocks high-affinity binding of operators with guanine at position 8 (Kd620 nM). Also, the high-affinity binding conformation appears blocked, as crystals of WT bound to DNA with guanine at position 8 could not be grown. In complexes containing K55A, the alanine side-chain is too far removed to engage in van der Waals interactions with the operator, and, with the loss of the general electrostatic interaction between the phosphate backbone and the ammonium group of lysine, K55A binds each operator weakly. However, the mutation leads to a swap of specificity of PurR for the base at position 8, with K55A exhibiting a twofold preference for guanine over adenine. In addition to defining the role of Lys55 in PurR minor groove binding, these studies provide structural insight into the minor groove binding specificities of other LacI/GalR family members that have either alanine (e.g. LacI, GalR, CcpA) or a basic residue (e.g. RafR, ScrR, RbtR) at the comparable position.     

Structure-based analysis of protein-RNA interactions using the program ENTANGLE

Until recently, drawing general conclusions about RNA recognition by proteins has been hindered by the paucity of high-resolution structures. We have analyzed 45 PDB entries of protein-RNA complexes to explore the underlying chemical principles governing both specific and non-sequence specific binding. To facilitate the analysis, we have constructed a database of interactions using ENTANGLE, a JAVA-based program that uses available structural models in their PDB format and searches for appropriate hydrogen bonding, stacking, electrostatic, hydrophobic and van der Waals interactions. The resulting database of interactions reveals correlations that suggest the basis for the discrimination of RNA from DNA and for base-specific recognition. The data illustrate both major and minor interaction strategies employed by families of proteins such as tRNA synthetases, ribosomal proteins, or RNA recognition motifs with their RNA targets. Perhaps most surprisingly, specific RNA recognition appears to be mediated largely by interactions of amide and carbonyl groups in the protein backbone with the edge of the RNA base. In cases where a base accepts a proton, the dominant amino acid donor is arginine, whereas in cases where the base donates a proton, the predominant acceptor is the backbone carbonyl group, not a side-chain group. This is in marked contrast to DNA-protein interactions, which are governed predominantly by amino acid side-chain interactions with functional groups that are presented in the accessible major groove. RNA recognition often proceeds through loops, bulges, kinks and other irregular structures that permit use of all the RNA functional groups and this is seen throughout the protein-RNA interaction database.