Molecular dynamics simulation studies of a protein-RNA complex with a selectively modified binding interface

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.     

Simulation of the structure and dynamics of nonhelical RNA motifs

Computer simulation methods are increasingly being used to study possible conformations and dynamics of structural motifs in RNA. Recent results of molecular dynamics simulations and continum solvent studies of RNA structures and RNA-ligand complexes show promising agreement with experimental data. Combined with the ongoing progress in the experimental characterization of RNA structure and thermodynamics, these computational approaches can help to better understand the mechanism of RNA structure formation and the binding of ligands.     

Structure and thermodynamics of RNA-protein binding: using molecular dynamics and free energy analyses to calculate the free energies of binding and conformational change

An adaptive binding mechanism, requiring large conformational rearrangements, occurs commonly with many RNA-protein associations. To explore this process of reorganization, we have investigated the conformational change upon spliceosomal U1A-RNA binding with molecular dynamics (MD) simulations and free energy analyses. We computed the energetic cost of conformational change in U1A-hairpin and U1A-internal loop binding using a hybrid of molecular mechanics and continuum solvent methods. Encouragingly, in all four free energy comparisons (two slightly different proteins, two different RNAs), the free macromolecule was more stable than the bound form by the physically reasonable value of approximately 10 kcal/mol. We calculated the absolute binding free energies for both complexes to be in the same range as that found experimentally.     

Altering the RNA-binding mode of the U1A RBD1 protein

The N-terminal RNA-binding domain (RBD1) of the human U1A protein is evolutionarily designed to bind its RNA targets with great affinity and specificity. The physical mechanisms that modulate the coupling (local cooperativity) among amino acid residues on the extensive binding surface of RBD1 are investigated here, using mutants that replace a highly conserved glycine residue. This glycine residue, at the strand/loop junction of beta3/loop3, is found in U1A RBD1, and in most RBD domains, suggesting it has a specific role in modulation of RNA binding. Here, two RBD1 proteins are constructed in which that residue (Gly53) is replaced by either alanine or valine. These new proteins are shown by NMR methods and molecular dynamics simulations to be very similar to the wild-type RBD1, both in structure and in their backbone dynamics. However, RNA-binding assays show that affinity for the U1 snRNA stem-loop II RNA target is reduced by nearly 200-fold for the RBD1-G53A protein, and by 1.6 x 10(4)-fold for RBD1-G53V. The mode of RNA binding by RBD1-G53A is similar to that of RBD1-WT, displaying its characteristic non-additive free energies of base recognition and its salt-dependence. The binding mode of RBD1-G53V is altered, having lost its salt-dependence and displaying site-independence of base recognition. The molecular basis for this alteration in RNA-binding properties is proposed to result from the inability of the RNA to induce a change in the structure of the free protein to produce a high-affinity complex.     

Functional and structural insights revealed by molecular dynamics simulations of an essential RNA editing ligase in Trypanosoma brucei

RNA editing ligase 1 (TbREL1) is required for the survival of both the insect and bloodstream forms of Trypanosoma brucei, the parasite responsible for the devastating tropical disease African sleeping sickness. The type of RNA editing that TbREL1 is involved in is unique to the trypanosomes, and no close human homolog is known to exist. In addition, the high-resolution crystal structure revealed several unique features of the active site, making this enzyme a promising target for structure-based drug design. In this work, two 20 ns atomistic molecular dynamics (MD) simulations are employed to investigate the dynamics of TbREL1, both with and without the ATP substrate present. The flexibility of the active site, dynamics of conserved residues and crystallized water molecules, and the interactions between TbREL1 and the ATP substrate are investigated and discussed in the context of TbREL1's function. Differences in local and global motion upon ATP binding suggest that two peripheral loops, unique to the trypanosomes, may be involved in interdomain signaling events. Notably, a significant structural rearrangement of the enzyme's active site occurs during the apo simulations, opening an additional cavity adjacent to the ATP binding site that could be exploited in the development of effective inhibitors directed against this protozoan parasite. Finally, ensemble averaged electrostatics calculations over the MD simulations reveal a novel putative RNA binding site, a discovery that has previously eluded scientists. Ultimately, we use the insights gained through the MD simulations to make several predictions and recommendations, which we anticipate will help direct future experimental studies and structure-based drug discovery efforts against this vital enzyme.     

Understanding the origins of bacterial resistance to aminoglycosides through molecular dynamics mutational study of the ribosomal A-site

Paromomycin is an aminoglycosidic antibiotic that targets the RNA of the bacterial small ribosomal subunit. It binds in the A-site, which is one of the three tRNA binding sites, and affects translational fidelity by stabilizing two adenines (A1492 and A1493) in the flipped-out state. Experiments have shown that various mutations in the A-site result in bacterial resistance to aminoglycosides. In this study, we performed multiple molecular dynamics simulations of the mutated A-site RNA fragment in explicit solvent to analyze changes in the physicochemical features of the A-site that were introduced by substitutions of specific bases. The simulations were conducted for free RNA and in complex with paromomycin. We found that the specific mutations affect the shape and dynamics of the binding cleft as well as significantly alter its electrostatic properties. The most pronounced changes were observed in the U1406C∶U1495A mutant, where important hydrogen bonds between the RNA and paromomycin were disrupted. The present study aims to clarify the underlying physicochemical mechanisms of bacterial resistance to aminoglycosides due to target mutations.     

Molecular dynamics and binding specificity analysis of the bovine immunodeficiency virus BIV Tat-TAR complex

We have performed molecular dynamics (MD) simulations, with particle-mesh Ewald, explicit waters, and counterions, and binding specificity analyses using combined molecular mechanics and continuum solvent (MM-PBSA) on the bovine immunodeficiency virus (BIV) Tat peptide-TAR RNA complex. The solution structure for the complex was solved independently by Patel and co-workers and Puglisi and co-workers. We investigated the differences in both structures and trajectories, particularly in the formation of the U-A-U base triple, the dynamic flexibility of the Tat peptide, and the interactions at the binding interface. We observed a decrease in RMSD in comparing the final average RNA structures and initial RNA structures of both trajectories, which suggests the convergence of the RNA structures to a MD equilibrated RNA structure. We also calculated the relative binding of different Tat peptide mutants to TAR RNA and found qualitative agreement with experimental studies.     

Minimum sequence requirements for the binding of paromomycin to the rRNA decoding site A

We have recently introduced a computational methodology that combines molecular dynamics (MD) simulations, free-energy calculations, and in vitro binding assays to predict the minimum RNA structural requirements for selective, high-affinity RNA binding to small-molecule ligands. Here, we show that this methodology can be applied to the conformationally flexible aminoglycoside antibiotic paromomycin. A RNA consisting of an 11-mer:10-mer duplex that contains one 16S ribosome RNA decoding A-site bound to paromomycin was simulated for 4 ns. The methodology predicts that the 11-mer:10-mer duplex binds to paromomycin with high affinity, whereas smaller RNA duplexes lose complex stability and the ability to bind paromomycin. The predicted high-affinity binding to paromomycin of the 11-mer:10-mer duplex was confirmed experimentally (EC(50) = 0.28 microM), as well as the inability of smaller complexes to bind. Our simulations show good agreement with experiment for dynamic and structural properties of the isolated A-site, including hydrogen-bonding networks and RNA structural rearrangements upon ligand binding. The results suggest that MD simulations can supplement in vitro methods as a tool for predicting minimum RNA-binding motifs for both small, rigid ligands, and large, flexible ligands when structural information is available.     

Correlated motions in the U1 snRNA stem/loop 2:U1A RBD1 complex

The complex formed by U1A RBD1 and the U1 snRNA stem/loop II is noted for its high affinity and exquisite specificity. Here, that complex is investigated by 5 ns molecular dynamics simulations and analyzed by reorientational eigenmode dynamics to determine the dynamic properties of the RNA:protein interface that could contribute to the binding mechanism. The analysis shows that there is extensive correlation between motions of the RNA and protein, involving 7 of the 10 RNA loop nucleotides, the protein beta-sheet surface, two of its loops, and its C-terminal tripeptide sequence. Order parameters of these regions of the complex are uniformly high, indicating restricted motion. However, several regions of both RNA and protein retain local flexibility, notably three nucleotides of the RNA loop and one loop of RBD1 that does not contact RNA. The highly correlated motions involving both molecules reflect the intricate network of interactions that characterize this complex and could account in part for the thermodynamic coupling observed for complex formation.     

Inhibition of the U1A-RNA complex by an aminoacridine derivative

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.     

Simulations of the dynamics at an RNA-protein interface.

Molecular dynamics simulations of the RNA-binding domain of the U1A spliceosomal protein in complex with its cognate RNA hairpin, performed at low and high ionic strength in aqueous solution, suggest a pathway for complex dissociation. First, cations condense around the RNA and compete with the protein for binding sites. Then solvated ions specifically destabilize residues at the RNA-protein interface. For a discrete cluster of residues at the complex interface, the simulations reveal an increased deviation from the crystal structure at high salt concentrations while the remaining protein scaffold is stabilized under these conditions. The microscopic picture of salt influence on the complex suggests guidelines for rational design of interface inhibitors targeted at RNA-protein complexes.     

Molecular dynamics simulation studies of induced fit and conformational capture in U1A–RNA binding: Do molecular substates code for specificity?

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.     

Spliceosomal SL1 RNA binding to U1-70K: the role of the extended RRM

The RNA recognition motif (RRM) occurs widely in RNA-binding proteins, but does not always by itself support full binding. For example, it is known that binding of SL1 RNA to the protein U1-70K in the U1 spliceosomal particle is reduced when a region flanking the RRM is truncated. How the RRM flanking regions that together with the RRM make up an ‘extended RRM’ (eRRM) contribute to complex stability and structural organization is unknown. We study the U1-70K eRRM bound to SL1 RNA by thermal dissociation and laser temperature jump kinetics; long-time molecular dynamics simulations interpret the experiments with atomistic resolution. Truncation of the helix flanking the RRM on its N-terminal side, ‘N-helix,’ strongly reduces overall binding, which is further weakened under higher salt and temperature conditions. Truncating the disordered region flanking the RRM on the C-terminal side, ‘C-IDR’, affects the local binding site. Surprisingly, all-atom simulations show that protein truncation enhances base stacking interactions in the binding site and leaves the overall number of hydrogen bonds intact. Instead, the flanking regions of the eRRM act in a distributed fashion via collective interactions with the RNA when external stresses such as temperature or high salt mimicking osmotic imbalance are applied.     

Discriminatory RNP remodeling by the DEAD-box protein DED1

DExH/D proteins catalyze NTP-driven rearrangements of RNA and RNA-protein complexes during most aspects of RNA metabolism. Although the vast majority of DExH/D proteins displays virtually no sequence-specificity when remodeling RNA complexes in vitro, the enzymes clearly distinguish between a large number of RNA and RNP complexes in a physiological context. It is unknown how this discrimination between potential substrates is achieved. Here we show one possible way by which a non-sequence specific DExH/D protein can discriminately remodel similar RNA complexes. We have measured in vitro the disassembly of model RNPs by two distinct DExH/D proteins, DED1 and NPH-II. Both enzymes displace the U1 snRNP from a tightly bound RNA in an active, ATP-dependent fashion. However, DED1 cannot actively displace the protein U1A from its binding site, whereas NPH-II can. The dissociation rate of U1A dictates the rate by which DED1 remodels RNA complexes with U1A bound. We further show that DED1 disassembles RNA complexes with slightly altered U1A binding sites at different rates, but only when U1A is bound to the RNA. These findings suggest that the "inability" to actively displace other proteins from RNA can provide non-sequence specific DExH/D proteins with the capacity to disassemble similar RNA complexes in a discriminatory fashion. In addition, our study illuminates possible mechanisms for protein displacement by DExH/D proteins.     

Multistep kinetics of the U1A-SL2 RNA complex dissociation

The U1A-SL2 RNA complex is a model system for studying interactions between RNA and the RNA recognition motif (RRM), which is one of the most common RNA binding domains. We report here kinetic studies of dissociation of the U1A-SL2 RNA complex, using laser temperature jump and stopped-flow fluorescence methods with U1A proteins labeled with the intrinsic chromophore tryptophan. An analysis of the kinetic data suggests three phases of dissociation with time scales of ∼ 100 μs, ∼ 50 ms, and ∼ 2 s. We propose that the first step of dissociation is a fast rearrangement of the complex to form a loosely bound complex. The intermediate step is assigned to be the dissociation of the U1A-SL2 RNA complex, and the final step is assigned to a reorganization of the U1A protein structure into the conformation of the free protein. These assignments are consistent with previous proposals based on thermodynamic, NMR, and surface plasmon resonance experiments and molecular dynamics simulations. Together, these results begin to build a comprehensive model of the complex dynamic processes involved in the formation and dissociation of an RRM-RNA complex.     

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Does 29-mer RNA hairpin of the HIV-1 TAR RNA sequence bind magnesium? 19F-NMR and modelling studies.

5-Fluorouridine residues have been introduced into functionally important bulge and loop regions of 29-mer HIV-1 TAR RNA hairpins I and II to study Mg2+ and Ca2+ binding using 19F-NMR spectroscopy. There was no substantial binding detected up to 20-molar excess in case of both cations, whereas association of argininamide, used as a reference ligand, could be detected at less than 1-molar excess. The deltadelta 19F value of 1.93 ppm observed for (F)U23 upon argininamide binding is in agreement with former NMR studies of TAR RNA/argininamide complex. However, obtained results do not confirm U38 x A27 x U23 base-triple formation. The unmodified HIV-1 TAR RNA hairpin resulted from 600 ps in aqua molecular dynamics simulation was subjected to a molecular mechanics modelling of Mg+ binding.     

U1A protein-stem loop 2 RNA recognition: prediction of structural differences from protein mutations

Molecular dynamics (MD) simulations were carried out to compare the free and bound structures of wild type U1A protein with several Phe56 mutant U1A proteins that bind the target stem loop 2 (SL2) RNA with a range of affinities. The simulations indicate the free U1A protein is more flexible than the U1A-RNA complex for both wild type and Phe56 mutant systems. A complete analysis of the hydrogen-bonding (HB) and non-bonded (VDW) interactions over the course of the MD simulations suggested that changes in the interactions in the free U1A protein caused by the Phe56Ala and Phe56Leu mutations may stabilize the helical character in loop 3, and contribute to the weak binding of these proteins to SL2 RNA. Compared with wild type, changes in HB and VDW interactions in Phe56 mutants of the free U1A protein are global, and include differences in β-sheet, loop 1 and loop 3 interactions. In the U1A-RNA complex, the Phe56Ala mutation leads to a series of differences in interactions that resonate through the complex, while the Phe56Leu and Phe56Trp mutations cause local differences around the site of mutation. The long-range networks of interactions identified in the simulations suggest that direct interactions and dynamic processes in both the free and bound forms contribute to complex stability.     

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.