Design of an adenosine analogue that selectively improves the affinity of a mutant U1A protein for RNA

Interpreting the Dynamics of Binding Interactions of snRNA and U1A Using a Coarse-Grained Model

The binding interactions of small nuclear RNAs (snRNA) and the associated protein factors are critical to the function of spliceosomes in alternatively splicing primary RNA transcripts. Although molecular dynamics simulations are a powerful tool to interpret the mechanism of biological processes, the atomic-level simulations are, however, too expensive and with limited accuracy for the large-size systems, such as snRNA-protein complexes. We extend the coarse-grained Gaussian network model, which models the RNA-protein complexes as a harmonic chain of Cα, P, and O4' atoms, to investigating the impact of the snRNA-binding interaction on the dynamic stability of the human U1A protein, which is a major component of the spliceosomal U1 small nuclear ribonucleoprotein particle. The results reveal that the first and third loops and the C-terminal helix regions of the U1A domain undergo a significant loss of flexibility upon the RNA binding due to the forming of mostly electrostatic and hydrogen bond interactions with RNA 5' stem and loop. By examining the residues whose mutations significantly change the binding free energy between U1A and snRNA, the Gaussian network model-based calculations show that not only the residues at the binding sites that are traditionally considered to play a major role in U1A-RNA association but also those residues that are far away from the RNA-binding interface can participate in the long-range allosteric signal transmission; these calculations are quantitatively consistent with the data observed in the recent snRNA binding experiments. The study demonstrates a useful avenue to utilize the simplified elastic network model to investigate the dynamics characteristics of the biologically important macromolecular interactions.

Palladium-catalyzed modification of unprotected nucleosides, nucleotides, and oligonucleotides

Synthetic modification of nucleoside structures provides access to molecules of interest as pharmaceuticals, biochemical probes, and models to study diseases. Covalent modification of the purine and pyrimidine bases is an important strategy for the synthesis of these adducts. Palladium-catalyzed cross-coupling is a powerful method to attach groups to the base heterocycles through the formation of new carbon-carbon and carbon-heteroatom bonds. In this review, approaches to palladium-catalyzed modification of unprotected nucleosides, nucleotides, and oligonucleotides are reviewed. Polar reaction media, such as water or polar aprotic solvents, allow reactions to be performed directly on the hydrophilic nucleosides and nucleotides without the need to use protecting groups. Homogeneous aqueous-phase coupling reactions catalyzed by palladium complexes of water-soluble ligands provide a general approach to the synthesis of modified nucleosides, nucleotides, and oligonucleotides.

Structural analysis of zinc-finger (TTK) + [Cu(BPA)]2+ /[Cu(IDB)]2+ + DNA complexes: an investigation by molecular dynamics simulation

In the present study, the molecular dynamics simulation technique is employed to investigate the hydrogen abstraction possibility from sugar of DNA in two designed complexes of copper-based chemical nuclease [Cu(BPA)](2+) bis(2-pyridylmethyl) amine (BPA) or [Cu(IDB)](2+) N,N-bis(2-benzimidazolylmethyl) amine (IDB) bound to the zinc finger protein Tramtrack (TTK). The simulated results show that each of the designed complexes can form a stable conformation within 30 ns of simulation time with the substrate OOH(-) and an 18-base pair (bp) DNA segment and is located in the major groove of the DNA segment. The active terminal O atom of the OOH(-) substrate is found in close proximity to the target C2'H, C3'H, C4'H or C5'H proton of the DNA in TTK + [Cu(BPA)OOH](+) + DNA or TTK + [Cu(IDB)OOH](+) + DNA complex, which is crucial to propose the hydrogen abstraction possibility that is responsible for the DNA cleavage. The positions of copper-based chemical nucleases bound to TTK may substantially influence the hydrogen abstraction possibility. The structures and sizes of ligands in copper-based nucleases are also found to have influence on the order of difficulty of the hydrogen abstraction from the sugars of DNA.

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.

Characterization of photophysical and base-mimicking properties of a novel fluorescent adenine analogue in DNA

To increase the diversity of fluorescent base analogues with improved properties, we here present the straightforward click-chemistry-based synthesis of a novel fluorescent adenine-analogue triazole adenine (A^T) and its photophysical characterization inside DNA. A^T shows promising properties compared to the widely used adenine analogue 2-aminopurine. Quantum yields reach >20% and >5% in single- and double-stranded DNA, respectively, and show dependence on neighbouring bases. Moreover, A^T shows only a minor destabilization of DNA duplexes, comparable to 2-aminopurine, and circular dichroism investigations suggest that A^T only causes minimal structural perturbations to normal B-DNA. Furthermore, we find that A^T shows favourable base-pairing properties with thymine and more surprisingly also with normal adenine. In conclusion, A^T shows strong potential as a new fluorescent adenine analogue for monitoring changes within its microenvironment in DNA.

Characterization of two adenosine analogs as fluorescence probes in RNA

The fluorescence properties of two adenosine analogs, 2-(3-phenylpropyl)adenosine [A-3CPh] and 2-(4-phenylbutyl)adenosine [A-4CPh], are reported. As monomers, the quantum yields and the mean lifetimes are 0.011 and 6.22 ns for A-3CPh and 0.007 and 7.13 ns for A-4CPh, respectively. Surprisingly, the quantum yields of the two probes are enhanced 11- to 82-fold upon incorporation into RNA, while the mean lifetimes decrease 23-40%. The data suggest that a subpopulation of molecules is responsible for the fluorescence characteristics and that the distribution of emitting and non-emitting structures is altered upon incorporation of the probes into RNA. Thus, although both adenosine analogs have low quantum yields as monomers, their fluorescence signals are significantly enhanced in RNA. Thermodenaturation experiments and CD spectroscopy indicate that incorporation of the adenosine analogs into three different RNAs does not alter their global structure or stability. Therefore, these probes should be useful for probing events occurring close to the site of modification.

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements

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.

Recognition of essential purines by the U1A protein

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.

Affinity and specificity of protein U1A-RNA complex formation based on an additive component free energy model

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.

A study of collective atomic fluctuations and cooperativity in the U1A-RNA complex based on molecular dynamics simulations

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.

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.