A Compare-and-Contrast NMR Dynamics Study of Two Related RRMs: U1A and SNF

The U1A/U2B″/SNF family of small nuclear ribonucleoproteins uses a phylogenetically conserved RNA recognition motif (RRM1) to bind RNA stemloops in U1 and/or U2 small nuclear RNA (snRNA). RRMs are characterized by their α/β sandwich topology, and these RRMs use their β-sheet as the RNA binding surface. Unique to this RRM family is the tyrosine-glutamine-phenylalanine (YQF) triad of solvent-exposed residues that are displayed on the β-sheet surface; the aromatic residues form a platform for RNA nucleobases to stack. U1A, U2B″, and SNF have very different patterns of RNA binding affinity and specificity, however, so here we ask how YQF in Drosophila SNF RRM1 contributes to RNA binding, as well as to domain stability and dynamics. Thermodynamic double-mutant cycles using tyrosine and phenylalanine substitutions probe the communication between those two residues in the free and bound states of the RRM. NMR experiments follow corresponding changes in the glutamine side-chain amide in both U1A and SNF, providing a physical picture of the RRM1 β-sheet surface. NMR relaxation and dispersion experiments compare fast (picosecond to nanosecond) and intermediate (microsecond-to-millisecond) dynamics of U1A and SNF RRM1. We conclude that there is a network of amino acid interactions involving Tyr-Gln-Phe in both SNF and U1A RRM1, but whereas mutations of the Tyr-Gln-Phe triad result in small local responses in U1A, they produce extensive microsecond-to-millisecond global motions throughout SNF that alter the conformational states of the RRM.     

Resurrection of an Urbilaterian U1A/U2B″/SNF Protein

The U1A/U2B″/SNF family of proteins found in the U1 and U2 spliceosomal small nuclear ribonucleoproteins is highly conserved. In spite of the high degree of sequence and structural conservation, modern members of this protein family have unique RNA binding properties. These differences have necessarily resulted from evolutionary processes, and therefore, we reconstructed the protein phylogeny in order to understand how and when divergence occurred and how protein function has been modulated. Contrary to the conventional understanding of an ancient human U1A/U2B″ gene duplication, we show that the last common ancestor of bilaterians contained a single ancestral protein (URB). The gene for URB was synthesized, the protein was overexpressed and purified, and we assessed RNA binding to modern snRNA sequences. We find that URB binds human and Drosophila U1 snRNA SLII and U2 snRNA SLIV with higher affinity than do modern homologs, suggesting that both Drosophila SNF and human U1A/U2B″ have evolved into weaker binders of one RNA or both RNAs.     

AMPK and SNF1: Snuffing Out Stress

AMP-activated protein kinase (AMPK) has attracted much attention for its key role in energy homeostasis. Three new papers providing structural information on mammalian and yeast AMPK homologs give insights into the binding of the regulatory nucleotides AMP and ATP and how mutations are associated with cardiac glycogen storage disorders.     

NMR studies of U1 snRNA recognition by the N-terminal RNP domain of the human U1A protein.

The RNP domain is a very common motif found in hundreds of proteins, including many protein components of the RNA processing machinery. The 70-90 amino acid domain contains two highly conserved stretches of 6-8 amino acids (RNP-1 and RNP-2) in the central strands of a four-stranded antiparallel beta-sheet, packed against two alpha-helices by a conserved hydrophobic core. Using multidimensional heteronuclear NMR, we have mapped intermolecular contacts between the human U1A protein 102 amino acid N-terminal RNP domain and a 31-mer oligonucleotide derived from stem-loop II of U1 snRNA. Chemical shift changes induced on the protein by the RNA define the surface of the beta-sheet as the recognition interface. The reverse face of the protein, with the two alpha-helices, remains exposed to the solvent in the presence of the RNA, and is potentially available for protein-protein contacts in spliceosome assembly or splice site selection. Protein-RNA contacts occur at the single-stranded apical loop of the hairpin, but also in the major groove of the helical stem at neighbouring U.G and U.U non-Watson-Crick base pairs. Examination of a proposed model for the complex in the light of the present results reveals several features of RNA recognition by RNP proteins. The quality of the spectra for this complex of 22 kDa demonstrates the feasibility of NMR investigation of RNA-protein complexes.     

Comparative Population Dynamics of Two Closely Related Species Differing in Ploidy Level

Background

Many studies compare the population dynamics of single species within multiple habitat types, while much less is known about the differences in population dynamics in closely related species in the same habitat. Additionally, comparisons of the effect of habitat types and species are largely missing.

Methodology and Principal Findings

We estimated the importance of the habitat type and species for population dynamics of plants. Specifically, we compared the dynamics of two closely related species, the allotetraploid species Anthericum liliago and the diploid species Anthericum ramosum, occurring in the same habitat type. We also compared the dynamics of A. ramosum in two contrasting habitats. We examined three populations per species and habitat type. The results showed that single life history traits as well as the mean population dynamics of A. liliago and A. ramosum from the same habitat type were more similar than the population dynamics of A. ramosum from the two contrasting habitats.

Conclusions

Our findings suggest that when transferring knowledge regarding population dynamics between populations, we need to take habitat conditions into account, as these conditions appear to be more important than the species involved (ploidy level). However, the two species differ significantly in their overall population growth rates, indicating that the ploidy level has an effect on species performance. In contrast to what has been suggested by previous studies, we observed a higher population growth rate in the diploid species. This is in agreement with the wider range of habitats occupied by the diploid species.     

Functional analysis of SNF, the Drosophila U1A/U2B" homolog: identification of dispensable and indispensable motifs for both snRNP assembly and function in vivo

In Drosophila, the spliceosomal protein SNF fulfills the functions of two vertebrate proteins, U1 snRNP-UlA and U2 snRNP-U2B". The structure and sequence of SNF, U1A, and U2B" are nearly identical with two RNA recognition motifs (RRM) separated by a short linker region, yet they have different RNA-binding properties: U1A binds U1 snRNA, U2B" binds U2 snRNA, and SNF binds both snRNAs. Structure/function studies on the human proteins have identified motifs in the N-terminal RRM that are critical for RNA-binding specificity but have failed to identify a function for the C-terminal RRM. Interestingly, SNF is chimeric in these motifs, suggesting a basis for its dual specificity. Here, we test the importance of these motifs by introducing site-directed mutations in the snf coding region and examining the effects of these mutations on assembly into the snRNP and on snf function in vivo. We found that an N-terminal RRM mutant protein predicted to eliminate RNA binding still assembles into snRNPs and is capable of rescuing snf's lethal phenotype only if the normally dispensable C-terminal RRM is present. We also found that the mixed motif in the "RNA-specificity" domain is necessary for SNF's dual function whereas the mixed motif in the U2A'-protein-binding region is not. Finally, we demonstrate that animals carrying a snf mutation that converts SNF from a bifunctional protein to a U1 snRNP-specific protein are viable. This unexpected result suggests that SNF's presence within the U2 snRNP is not essential for splicing.     

Two functionally distinct steps mediate high affinity binding of U1A protein to U1 hairpin II RNA

Binding of the U1A protein to its RNA target U1 hairpin II has been extensively studied as a model for a high affinity RNA/protein interaction. However, the mechanism and kinetics by which this complex is formed remain largely unknown. Here we use real-time biomolecular interaction analysis to dissect the roles various protein and RNA structural elements play in the formation of the U1A.U1 hairpin II complex. We show that neutralization of positive charges on the protein or increasing the salt concentration slows the association rate, suggesting that electrostatic interactions play an important role in bringing RNA and protein together. In contrast, removal of hydrogen bonding or stacking interactions within the RNA/protein interface, or reducing the size of the RNA loop, dramatically destabilizes the complex, as seen by a strong increase in the dissociation rate. Our data support a binding mechanism consisting of a rapid initial association based on electrostatic interactions and a subsequent locking step based on close-range interactions that occur during the induced fit of RNA and protein. Remarkably, these two steps can be clearly distinguished using U1A mutants containing single amino acid substitutions. Our observations explain the extraordinary affinity of U1A for its target and may suggest a general mechanism for high affinity RNA/protein interactions.     

Conformational Studies of Two Bradykinin Antagonists by Using Two-dimensional NMR Techniques and Molecular Dynamics Simulations

Abstract

The solution conformations of two potent antagonists of bradykinin (Arg¹-Pro²-Pro³-Gly⁴- Phe⁵-Ser⁶-Pro⁷-Phe⁸-Arg⁹), [Aca-¹, DArg⁰, Hyp³, Thi⁵, DPhe⁷,(N-Bzl)Gly⁸]BK (1) and [Aaa- ¹, DArg⁰, Hyp³, Thi⁵,(2-DNal)⁷, Thi⁸]BK (2), were studied by using 2D NMR spectroscopy in DMSO-dg and molecular dynamics simulations. The NMR spectra of peptide 1 reveals the existence of at least two isomers arising from isomerization across the DPhe⁷-(N-Bzl)Gly⁸peptide bond. The more populated isomer possesses the cis peptide bond at this position. The ratio of cis/trans isomers amounted to 7:3. With both antagonists, the NMR data indicate a β-turn structure for the Hyp³-Gly⁴ residues. In addition, for peptide 2, position 2,3 is likely to be occupied by turn-like structures. The cis peptide bond between DPhe⁷ and (N- Bzl)Gly⁸ in analogue 1 suggests type VI β-turn at position 7,8. The molecular dynamics runs were performed on both peptides in DMSO solution. The results indicate that the structure of peptide 1 is characterized by type VIb β-turn comprising residues Ser-Arg⁹ and the βI or βII-turn involving the Pro²-Thi⁵ fragment, whereas peptide 2 shows the tendency towards the formation of type I β-turn at position 2,3. The structures of both antagonists are stabilized by a salt bridge between the guanidine moiety of Arg¹ and the carboxyl group of Arg⁹. Moreover, the side chain of DArg⁰ is apart of the rest of molecule and is not involved in structural elements except for a few calculated structures.     

Dynamics and metal ion binding in the U6 RNA intramolecular stem-loop as analyzed by NMR

The U6 RNA intramolecular stem-loop (ISL) is a conserved component of the spliceosome, and contains an essential metal ion binding site centered between a protonated adenine, A79, and U80. Correlated with protonation of A79, U80 undergoes a base-flipping conformational change accompanied by significant helical movement. We have investigated the dynamics of the U6 ISL by analyzing the power dependence of 13C NMR relaxation rates in the rotating frame. The data provide evidence that the conformational transition is centered around an exchange lifetime of 84 micros. The U80 nucleotide displays low internal mobility on the picosecond time-scale at pH 7.0 but high internal mobility at pH 6.0, in agreement with the global transition resulting in the base of U80 adopting a looped-out conformation with increased dynamic disorder. A kinetic analysis suggests that the conformational change, rather than adenine protonation, is the rate-limiting step in the pathway of the conformational transition. Two nucleotides, U70 and U80, were found from chemical shift perturbation mapping to interact with the magnesium ion, with apparent K(d) values in the micromolar to millimolar range. These nucleotides also displayed metal ion-induced elevation of R1 rates, which can be explained by a model that assumes dynamic metal ion coordination concomitant with an induced higher shielding anisotropy for the base 13C nuclei. Addition of Mg2+ shifts the conformational equilibrium toward the high-pH (base-stacked) structure, accompanied by a significant drop in the apparent pK(a) of A79.     

Isoforms of U1-70k Control Subunit Dynamics in the Human Spliceosomal U1 snRNP

Most human protein-encoding genes contain multiple exons that are spliced together, frequently in alternative arrangements, by the spliceosome. It is established that U1 snRNP is an essential component of the spliceosome, in human consisting of RNA and ten proteins, several of which are post-translationally modified and exist as multiple isoforms. Unresolved and challenging to investigate are the effects of these post translational modifications on the dynamics, interactions and stability of the particle. Using mass spectrometry we investigate the composition and dynamics of the native human U1 snRNP and compare native and recombinant complexes to isolate the effects of various subunits and isoforms on the overall stability. Our data reveal differential incorporation of four protein isoforms and dynamic interactions of subunits U1-A, U1-C and Sm-B/B''. Results also show that unstructured post-translationally modified C-terminal tails are responsible for the dynamics of Sm-B/B'' and U1-C and that their interactions with the Sm core are controlled by binding to different U1-70k isoforms and their phosphorylation status in vivo. These results therefore provide the important functional link between proteomics and structure as well as insight into the dynamic quaternary structure of the native U1 snRNP important for its function.     

Solution Conformation of [D-Pen2,D-Pen5]enkephalin in Water: A NMR and Molecular Dynamics Study

The solution conformation of [D -Pen²,D -Pen⁵] enkephalin (DPDPE), a highly potent δ-selective opioid agonist, was examined by means of NMR, molecular mechanics and molecular dynamics methods. The structural information in the solvent water was obtained employing one- and two-dimensional methods of ¹H and ¹³C-NMR spectroscopy. Based on the distance geometry technique using the ROE data as input, 400 conformers were obtained and considered in the structure analysis. Alternatively, about 2000 conformers were stochastically generated and related to the NMR data after energy minimization. The structure analysis provides one conformer in agreement with all NMR data, which belongs to the lowest energy conformation group. This structure may serve as a reference conformer for DPDPE analogues synthesized with the aim of activity increase.     

Physiological characterization of glucose repression in the strains with SNF1 and SNF4 genes deleted

We investigated the effect of Snf1 kinase and its regulatory subunit Snf4 on the regulation of glucose and galactose metabolism in the yeast Saccharomyces cerevisiae by physiologically characterizing Deltasnf1, Deltasnf4 and Deltasnf1Deltasnf4 in CEN.PK background in glucose and glucose-galactose-mixture batch cultivations. The main result of this study showed that delayed induction of galactose catabolism was SNF1 or SNF4 gene deletion specific. In comparison to the reference strain, growth delay on galactose was found to last 2.4 times (7 h), 3.1 times (10.5 h) and 9.6 times (43 h) longer for the Deltasnf4, Deltasnf1 and Deltasnf1Deltasnf4 strains, respectively. The maximum specific growth rates on galactose were determined to be two to three times lower for the recombinant strains compared to the reference strain (0.13 h(-1)) and were found to be 0.07, 0.08 and 0.04 h(-1) for the Deltasnf1, Deltasnf4 and Deltasnf1Deltasnf4 strains, respectively. The study showed that Snf1 kinase was not solely responsible for the derepression of galactose metabolism.     

Full Length Vpu from HIV-1: Combining Molecular Dynamics Simulations with NMR Spectroscopy

Abstract

Based on structures made available by solution NMR, molecular models of the protein Vpu from HIV-1 were built and refined by 6 ns MD simulations in a fully hydrated lipid bilayer. Vpu is an 81 amino acid type I integral membrane protein encoded by the human immunodeficiency virus type-1 (HIV-1) and closely related simian immunodeficiency viruses (SIVs). Its role is to amplify viral release. Upon phosphorylation, the cytoplasmic domain adopts a more compact shape with helices 2 and 3 becoming almost parallel to each other. A loss of helicity for several residues belonging to the helices adjacent to both ends of the loop region containing serines 53 and 57 is observed. A fourth helix, present in one of the NMR-based structures of the cytoplasmic domain and located near the C-terminus, is lost upon phosphorylation.     

Binding of Two Intrinsically Disordered Peptides to a Multi-Specific Protein: A Combined Monte Carlo and Molecular Dynamics Study

The unique ability of intrinsically disordered proteins (IDPs) to fold upon binding to partner molecules makes them functionally well-suited for cellular communication networks. For example, the folding-binding of different IDP sequences onto the same surface of an ordered protein provides a mechanism for signaling in a many-to-one manner. Here, we study the molecular details of this signaling mechanism by applying both Molecular Dynamics and Monte Carlo methods to S100B, a calcium-modulated homodimeric protein, and two of its IDP targets, p53 and TRTK-12. Despite adopting somewhat different conformations in complex with S100B and showing no apparent sequence similarity, the two IDP targets associate in virtually the same manner. As free chains, both target sequences remain flexible and sample their respective bound, natively -helical states to a small extent. Association occurs through an intermediate state in the periphery of the S100B binding pocket, stabilized by nonnative interactions which are either hydrophobic or electrostatic in nature. Our results highlight the importance of overall physical properties of IDP segments, such as net charge or presence of strongly hydrophobic amino acids, for molecular recognition via coupled folding-binding.     

The Dynamics of Territory Acquisition: A Model of Two Coexisting Strategies

Two potential strategies for acquiring territories are (1) fighting and taking over a territory from its owner, or (2) waiting until a territory owner dies and then taking its place. In this paper we explore territory acquisition using these two strategies, using a population dynamical model. Factors expected to affect the predominance of each strategy are injury rate, rate of successful territory takeover, birth and death rates on the territory, and birth and death rates while non-territorial. We explore the effects of these parameters on numbers of territorial and nonterritorial fighters and waiters. Waiters predominate when injury rate, birth rate of nonterritorial individuals, and death rate of territory owners are high. Fighters predominate when rate of successful territory takeover, death rate of nonterritorial individuals, and birth rate of territory owners are high. We present supportive evidence for these preditions from the literature.     

Molecular analysis of the SNF4 gene of Saccharomyces cerevisiae: evidence for physical association of the SNF4 protein with the SNF1 protein kinase.

The SNF4 gene is required for expression of glucose-repressible genes in response to glucose deprivation in Saccharomyces cerevisiae. Previous evidence suggested that SNF4 is functionally related to SNF1, another essential gene in this global regulatory system that encodes a protein kinase. Increased SNF1 gene dosage partially compensates for a mutation in SNF4, and the SNF4 function is required for maximal SNF1 protein kinase activity in vitro. We have cloned SNF4 and identified its 1.2-kilobase RNA, which is not regulated by glucose repression. A 36-kilodalton SNF4 protein is predicted from the nucleotide sequence. Disruption of the chromosomal SNF4 locus revealed that the requirement for SNF4 function is less stringent at low temperature (23 degrees C). A bifunctional SNF4-lacZ gene fusion that includes almost the entire SNF4 coding sequence was constructed. The fusion protein was shown by immunofluorescence microscopy to be distributed throughout the cell, with partial localization to the nucleus. The SNF4-beta-galactosidase protein coimmunoprecipitated with the SNF1 protein kinase, thus providing evidence for the physical association of the two proteins.     

Drosophila SNF/D25 combines the functions of the two snRNP proteins U1A and U2B' that are encoded separately in human, potato, and yeast.

The plant and vertebrate snRP proteins U1A and U2B' are structurally closely related, but bind to different U snRNAs. Two additional related snRNP proteins, the yeast U2B' protein and Drosophila SNF/D25 protein, are analyzed here. We show that the previously described yeast open reading frame YIB9w encodes yeast U2B' as judged by the fact that the protein encoded by YIB9w bindsto stem-loop IV of yeast U2 snRNA in vitro and is part of the U2 snRNP in vivo. In contrast to the human U2B' protein, specific binding of yeast U2B' to RNA in vitro can occur in the absence of an accessory U2A' protein. The Drosophila SNF-D25 protein, unlike all other U1A/U2B' proteins studied to date, is shown to be a component of both U1 and U2 snRNPs. In vitro, SNF/D25 binds to U1 snRNA on itsown and to U2 snRNA in the presence of either the human U2A' protein or of Drosophila nuclear extract. Thus, its RNA-binding properties are the sum of those exhibited by human or potato U1A and U2B' proteins. Implications for the role of SNF/D25 in alternative splicing, and for the evolution of the U1A/U2B' protein family, are discussed.