Search

A working hypothesis for the pathogenesis of myotonic dystrophy type 1 (DM1) involves the aberrant sequestration of an alternative splicing regulator, MBNL1, by expanded CUG repeats, r(CUG)(exp). It has been suggested that a reversal of the myotonia and potentially other symptoms of the DM1 disease can be achieved by inhibiting the toxic MBNL1-r(CUG)(exp) interaction. Using rational design, we discovered an RNA-groove binding inhibitor (ligand 3) that contains two triaminotriazine units connected by a bisamidinium linker. Ligand 3 binds r(CUG)12 with a low micromolar affinity (K(d) = 8 ± 2 μM) and disrupts the MBNL1-r(CUG)12 interaction in vitro (K(i) = 8 ± 2 μM). In addition, ligand 3 is cell and nucleus permeable, exhibits negligible toxicity to mammalian cells, dissolves MBNL1-r(CUG)(exp) ribonuclear foci, and restores misregulated splicing of IR and cTNT in a DM1 cell culture model. Importantly, suppression of r(CUG)(exp) RNA-induced toxicity in a DM1 Drosophila model was observed after treatment with ligand 3. These results suggest ligand 3 as a lead for the treatment of DM1.

Myotonic dystrophy type 1 (DM1) is caused by an expanded CUG repeat (CUG(exp)) that sequesters muscleblind-like 1 protein (MBNL1), a protein that regulates alternative splicing. CUG(exp) RNA is a validated drug target for this currently untreatable disease. Herein, we develop a bioactive small molecule (1) that targets CUG(exp) RNA and is able to inhibit the CUG(exp)·MBNL1 interaction in cells that model DM1. The core of this small molecule is based on ligand 2, which was previously reported to be active in an in vitro assay. A polyamine-derivative side chain was conjugated to this core to make it aqueous-soluble and cell-penetrable. In a DM1 cell model this conjugate was found to disperse CUG(exp) ribonuclear foci, release MBNL1, and partially reverse the mis-splicing of the insulin receptor pre-mRNA. Direct evidence for ribonuclear foci dispersion by this ligand was obtained in a live DM1 cell model using time-lapse confocal microscopy.

Muscleblind-like proteins (MBNL) are RNA-binding proteins that bind to the poly(CUG) and poly(CCUG) sequences that are the causative agents of myotonic dystrophy. It has been suggested that as a result of binding to the repeating RNA sequences, MBNL1 is abnormally expressed and translocated, which leads to many of the misregulated events in myotonic dystrophy. In this work, steady-state fluorescence quenching experiments suggest that MBNL1 alters the structure of helical RNA targets upon binding, which may explain the selectivity of MBNL1 for less structured RNA sites. The removal of one pair of zinc fingers greatly impairs the binding affinity of MBNL1, which indicates that the two pairs of zinc fingers might possibly interact with RNA targets cooperatively. Alanine scanning mutagenesis results suggest that the binding energy may be distributed across the protein. Overall, the results presented here suggest that small molecules that stabilize the helical structure of poly(CUG) and poly(CCUG) RNAs will inhibit the formation of complexes with MBNL1., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Myotonic dystrophy type 2 (DM2) is an incurable neuromuscular disease caused by expanded CCUG repeats that may exhibit toxicity by sequestering the splicing regulator MBNL1. A series of triaminotriazine- and triaminopyrimidine-based small molecules (ligands 1-3) were designed, synthesized and tested as inhibitors of the MBNL1-CCUG interaction. Despite the structural similarities of the triaminotriazine and triaminopyrimidine units, the triaminopyrimidine-based ligands bind with low micromolar affinity to CCUG repeats (K(d) ∼ 0.1-3.6 µM) whereas the triaminotriazine ligands do not bind CCUG repeats. Importantly, these simple and small triaminopyrimidine ligands exhibit both strong inhibition (K(i) ∼ 2 µM) of the MBNL1-CCUG interaction and high selectivity for CCUG repeats over other RNA targets. These experiments suggest these compounds are potential lead agents for the treatment of DM2.

The effects of modifying the electronic characteristics of nonpolar base analogues substituted at positions involved in stacking interactions between SL2 RNA and the U1A protein are described. A surprisingly large difference in the stability between complexes formed with base analogues that differ only in the position of substitution of a single fluorine atom is observed. The results of high-level ab initio calculations of the interactions between the nonpolar base analogue and the amino acid side chain correlate with the experimentally observed trends in complex stability, which suggests that changes in stacking interactions that result from varying the position and degree of fluorine substitution contribute to the effects of fluorine substitution on the stability of the U1A-SL2 RNA complex.

This work describes the rational design, synthesis, and study of a ligand that selectively complexes CUG repeats in RNA (and CTG repeats in DNA) with high nanomolar affinity. This sequence is considered a causative agent of myotonic dystrophy type 1 (DM1) because of its ability to sequester muscleblind-like (MBNL) proteins. Ligand 1 was synthesized in two steps from commercially available compounds, and its binding to CTG and CUG repeats in oligonucleotides studied. Isothermal titration calorimetry studies of 1 with various sequences showed a preference toward the T-T mismatch (K(d) of 390 +/- 80 nM) with a 13-, 169-, and 85-fold reduction in affinity toward single C-C, A-A, and G-G mismatches, respectively. Binding and Job analysis of 1 to multiple CTG step sequences revealed high affinity binding to every other T-T mismatch with negative cooperativity for proximal T-T mismatches. The affinity of 1 for a (CUG)(4) step provided a K(d) of 430 nM with a binding stoichiometry of 1:1. The preference for the U-U in RNA was maintained with a 6-, >143-, and >143-fold reduction in affinity toward single C-C, A-A, and G-G mismatches, respectively. Ligand 1 destabilized the complexes formed between MBNL1N and (CUG)(4) and (CUG)(12) with IC(50) values of 52 +/- 20 microM and 46 +/- 7 microM, respectively, and K(i) values of 6 +/- 2 microM and 7 +/- 1 microM, respectively. These values were only minimally altered by the addition of competitor tRNA. Ligand 1 does not destabilize the unrelated RNA-protein complexes the U1A-SL2 RNA complex and the Sex lethal-tra RNA complex. Thus, ligand 1 selectively destabilizes the MBNL1N-poly(CUG) complex.

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.

Background: The RNA recognition motif (RRM) is one of the largest families of RNA binding domains. The RRM is modulated so that individual proteins containing RRMs can specifically recognize RNA targets with diverse sequences and structures. Understanding the principles governing this specificity will be important for the rational modification and design of RRM-RNA complexes., Results: In this paper we have investigated the origins of specificity of the N terminal RRM of the U1A protein for stem loop 2 (SL2) of U1 snRNA by substituting modified bases for essential purines in SL2 RNA. In one series of modified bases, hydrogen bond donors and acceptors were replaced by aliphatic groups to probe the importance of these functional groups to binding. In a second series of modified bases, hydrogen bond donors and acceptors were incorrectly placed on the purine bases to analyze the origins of discrimination between cognate and non-cognate RNA. The results of these experiments show that three different approaches are used by the U1A protein to gain specificity for purines. Specificity for the first base in the loop, A1, is based primarily on discrimination against RNA containing the incorrect base, specificity for the fourth base in the loop, G4, is based largely on recognition of the donors and acceptors of G4, while specificity for the sixth base in the loop, A6, results from a combination of direct recognition of the base and discrimination against incorrectly placed functional groups., Conclusion: These investigations identify different roles that hydrogen bond donors and acceptors on bases in both cognate and non-cognate RNA play in the specific recognition of RNA by the U1A protein. Taken together with investigations of other RNA-RRM complexes, the results contribute to a general understanding of the origins of RNA-RRM specificity and highlight, in particular, the contribution of steric and electrostatic repulsion to binding specificity.

An MM-GBSA computational protocol was used to investigate wild-type U1A-RNA and F56 U1A mutant experimental binding free energies. The trend in mutant binding free energies compared to wild-type is well-reproduced. Following application of a linear-response-like equation to scale the various energy components, the binding free energies agree quantitatively with observed experimental values. Conformational adaptation contributes to the binding free energy for both the protein and the RNA in these systems. Small differences in DeltaGs are the result of different and sometimes quite large relative contributions from various energetic components. Residual free energy decomposition indicates differences not only at the site of mutation, but throughout the entire protein. MM-GBSA and ab initio calculations performed on model systems suggest that stacking interactions may nearly, but not completely, account for observed differences in mutant binding affinities. This study indicates that there may be different underlying causes of ostensibly similar experimentally observed binding affinities of different mutants, and thus recommends caution in the interpretation of binding affinities and specificities purely by inspection.

Cooperative interactions play an important role in recognition and binding in macromolecular systems. In this study, we find that cross-correlated atomic fluctuations can be used to identify cooperative networks in a protein-RNA system. The dynamics of the RRM-containing protein U1A-stem loop 2 RNA complex have been calculated theoretically from a 10 ns molecular dynamics (MD) simulation. The simulation was analyzed by calculating the covariance matrix of all atomic fluctuations. These matrix elements are then presented in the form of a two-dimensional grid, which displays fluctuations on a per residue basis. The results indicate the presence of strong, selective cross-correlated fluctuations throughout the RRM in U1A-RNA. The atomic fluctuations correspond well with previous biophysical studies in which a multiplicity of cooperative networks have been reported and indicate that the various networks identified in separate individual experiments are fluctuationally correlated into a hyper-network encompassing most of the RRM. The calculated results also correspond well with independent results from a statistical covariance analysis of 330 aligned RRM sequences. This method has significant implications as a predictive tool regarding cooperativity in the protein-nucleic acid recognition process.

A complete understanding of gene expression relies on a comprehensive understanding of the protein-RNA recognition process. However, the study of protein-RNA recognition is complicated by many factors that contribute to both binding affinity and specificity, including structure, energetics, dynamical motions, and cooperative interactions. Several recent studies have suggested that energetic coupling between residues contributes to formation of the complex between the U1A protein and stem loop 2 of U1 snRNA as a consequence of a cooperative network of interactions. We have performed molecular dynamics simulations on the U1A-RNA complex, including explicit water and counterions, and analyzed the results based on the calculated positional cross-correlations of atomic fluctuations. The results indicate that cross-correlations calculated on a per residue basis agree well with the observed inter-residue cooperativity and predict that the networks identified to date may also be coupled into an extensive hyper-network that reflects the intrinsic rigidity of the RNA recognition motif. In addition, we report a comparison of the MD calculated correlations with the results of a positional covariance analysis based on the sequences of 330 RNA recognition motifs, including U1A. The calculated inter-residue cross-correlations agree very well with the results of the sites exhibiting positional covariance. Collectively, these results strongly support the hypothesis that collective fluctuations contribute to cooperativity and the corresponding observed thermodynamic coupling. Predictions of additional sites in U1A that may be involved in cooperative networks are advanced.

The RNA recognition motif (RRM) is one of the most common RNA binding domains. We have investigated the contribution of three highly conserved aromatic amino acids to RNA binding by the N-terminal RRM of the U1A protein. Recently, we synthesized a modified base (A-4CPh) in which a phenyl group is tethered to adenine using a linker of 4 methylene groups. The substitution of this base for adenine in the target RNA selectively stabilizes the complex formed with a U1A protein in which one of the conserved aromatic amino acids is replaced with Ala (Phe56Ala). In this article, we report molecular dynamics (MD) simulations that probe the structural consequences of the substitution of A-4CPh for adenine in the wild type and Phe56Ala U1A-RNA complexes and in the free RNA. The simulations suggest that A-4CPh stabilizes the complex formed with Phe56Ala by adopting a folded conformation in which the tethered phenyl group fills the site occupied by the phenyl group of Phe56 in the wild-type complex. In contrast, an extended conformation of A-4CPh is predicted to be most stable in the complex formed with the wild-type protein. The calculations indicate A-4CPh is in an extended conformation in the free RNA. Therefore, preorganizing the structure of the phenyl-tethered base for binding may improve both the affinity and specificity of the RNA containing A-4CPh for the Phe56Ala U1A protein. Taken together, the previous experimental work and the calculations reported here suggest a general design strategy for altering RRM-RNA complex stability., (Copyright 2005 Wiley Periodicals, Inc.)

The RNA recognition motif (RRM), one of the most common RNA binding domains, contains three highly conserved aromatic amino acids that participate in stacking interactions with RNA bases. We have investigated the contribution of these highly conserved aromatic amino acids to the affinity of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA. Previously, we found that substitution of one of these conserved aromatic amino acids, Phe56, with Ala resulted in a large destabilization of the complex. Here, we have modified A6, the base in stem loop 2 RNA that stacks with Phe56, to compensate for a portion of the destabilization caused by the Phe56Ala mutation. We have designed two modified adenosines, A-3CPh and A-4CPh, in which a phenyl group is linked to the adenosine such that it may replace the phenyl group that is eliminated by the Phe56Ala mutation in the complex. We have found that incorporation of A-3CPh into stem loop 2 RNA stabilizes the complex formed with Phe56Ala by 0.6 kcal/mol, while incorporation of A-4CPh into stem loop 2 RNA stabilizes this complex by 1.8 kcal/mol. Either base modification destabilizes the wild-type complex by 0.8-0.9 kcal/mol. Experiments with other U1A mutant proteins suggest that the stabilization of the complex between the Phe56Ala U1A protein and stem loop 2 RNA is due to a specific interaction between the Phe56Ala U1A protein and A6-4CPh stem loop 2 RNA.

Molecular dynamics (MD) simulations on stem loop 2 of U1 small nuclear RNA and a construct of the U1A protein were carried out to obtain predictions of the structures for the unbound forms in solution and to elucidate dynamical aspects of induced fit upon binding. A crystal structure of the complex between the U1A protein and stem loop 2 RNA and an NMR structure for the uncomplexed form of the U1A protein are available from Oubridge et al. (Nature, 1994, Vol. 372, pp. 432-438) and Avis et al. (Journal of Molecular Biology, 1996, Vol. 257, pp. 398-411), respectively. As a consequence, U1A-RNA binding is a particularly attractive case for investigations of induced fit in protein-nucleic acid complexation. When combined with the available structural data, the results from simulations indicate that structural adaptation of U1A protein and RNA define distinct mechanisms for induced fit. For the protein, the calculations indicate that induced fit upon binding involves a non-native thermodynamic substate in which the structure is preorganized for binding. In contrast, induced fit of the RNA involves a distortion of the native structure in solution to an unstable form. However, the RNA solution structures predicted from simulation show evidence that structures in which groups of bases are favorably oriented for binding the U1A protein are thermally accessible. These results, which quantify with computational modeling recent proposals on induced fit and conformational capture by Leuillot and Varani (Biochemistry, 2001, Vol. 40, pp. 7947-7956) and by Williamson (Nature Structural Biology, 2000, Vol. 7, pp. 834-837) suggest an important role for intrinsic molecular architecture and substates other than the native form in the specificity of protein-RNA interactions., (Copyright 2002 Wiley Periodicals, Inc.)

The RNA recognition motif (RRM) binds to single-stranded RNA target sites of diverse sequences and structures. A conserved mode of base recognition by the RRM involves the simultaneous formation of a network of hydrogen bonds with the base functional groups and a stacking interaction between the base and a highly conserved aromatic amino acid. We have investigated the energetic contribution of the functional groups involved in the recognition of an essential adenine, A6, in stem-loop 2 of U1 snRNA by the N-terminal RRM of the U1A protein. Previously, we found that elimination of individual hydrogen bond donors and acceptors on A6 destabilized the complex by 0.8-1.9 kcal/mol, while mutation of the aromatic amino acid (Phe56) that stacks with A6 to Ala destabilized the complex by 5.5 kcal/mol. Here we continue to probe the contribution of A6 to complex stability through mutation of both the RNA and protein. We have removed two hydrogen-bonding functional groups by introducing a U1A mutation, Ser91Ala, and replacing A6 with tubercidin, purine, or 1-deazaadenine. We find that the complex is destabilized an additional 1.2-2.6 kcal/mol by the elimination of the second hydrogen bond donor or acceptor. Surprisingly, deletion of all of the functional groups involved in hydrogen bonds with the U1A protein by substituting adenine with 4-methylindole reduced the binding free energy by only 2.0 kcal/mol. Experiments with U1A proteins containing mutations of Phe56 suggested that improved stacking interactions due to the greater hydrophobicity of 4-methylindole than adenine may be partly responsible for the small destabilization of the complex upon substitution of 4-methylindole for A6. The data imply that hydrophobic interactions can compensate energetically for the disruption of the complex hydrogen-bonding network between nucleotide and protein.

The RNA recognition motif (RRM) is one of the most common RNA binding domains. There have been few investigations of small molecule inhibitors of RRM-RNA complexes, although these inhibitors could be valuable tools for probing biological processes involving RRM-RNA complexes and would have the potential to be effective drugs. In this paper, the inhibition by small molecules of the complex formed between the N-terminal RRM of the U1A protein and stem loop 2 of U1 snRNA has been investigated. An aminoacridine derivative has been found to promote dissociation of the U1A-stem loop 2 RNA complex with an IC(50) value of 1 microM. Fluorescence experiments indicate that two aminoacridine ligands bind to each RNA target site. RNase A footprinting suggests that one binding site may be near the base pair that closes the loop and the other may be in a more flexible region of the loop. The addition of the aminoacridine derivative to stem loop 2 RNA increases the susceptibility of other portions of the loop to digestion by RNase A, which implies that binding of the ligand changes the conformation or dynamics of the stem loop target site. Either direct binding to the RNA or indirect alteration of the structure or dynamics of the loop would be likely to inhibit binding of the U1A protein to this RNA.

Three highly conserved aromatic residues in RNA recognition motifs (RRM) participate in stacking interactions with RNA bases upon binding RNA. We have investigated the contribution of one of these aromatic residues, Phe56, to the complex formed between the N-terminal RRM of the spliceosomal protein U1A and stem-loop 2 of U1 snRNA. Previous work showed that the aromatic group is important for high affinity binding. Here we probe how mutation of Phe56 affects the kinetics of complex dissociation, the strength of the hydrogen bonds formed between U1A and the base that stacks with Phe56 (A6) and specific target site recognition. Substitution of Phe56 with Trp or Tyr increased the rate of dissociation of the complex, consistent with previously reported results. However, substitution of Phe56 with His decreased the rate of complex association, implying a change in the initial formation of the complex. Simultaneous modification of residue 56 and A6 revealed energetic coupling between the aromatic group and the functional groups of A6 that hydrogen bond to U1A. Finally, mutation of Phe56 to Leu reduced the ability of U1A to recognize stem-loop 2 correctly. Taken together, these experiments suggest that Phe56 contributes to binding affinity by stacking with A6 and participating in networks of energetically coupled interactions that enable this conserved aromatic amino acid to play a complex role in target site recognition.