Zentrum für Biomolekulare Magnetische Resonanz (BMRZ)
Refine
Year of publication
Document Type
- Article (111)
- Preprint (2)
- Part of a Book (1)
Has Fulltext
- yes (114)
Is part of the Bibliography
- no (114)
Keywords
- RNA (9)
- SARS-CoV-2 (8)
- Solution-state NMR (8)
- NMR spectroscopy (6)
- structural biology (6)
- rna (5)
- Biophysical chemistry (4)
- COVID19-NMR (4)
- NMR (4)
- rhodesain (4)
Institute
- Zentrum für Biomolekulare Magnetische Resonanz (BMRZ) (114)
- Biochemie und Chemie (63)
- Exzellenzcluster Makromolekulare Komplexe (39)
- Biochemie, Chemie und Pharmazie (30)
- Biowissenschaften (28)
- Sonderforschungsbereiche / Forschungskollegs (25)
- Medizin (4)
- Frankfurt Institute for Advanced Studies (FIAS) (3)
- Georg-Speyer-Haus (3)
- Pharmazie (3)
Mechanistic understanding of dynamic membrane proteins such as transporters, receptors, and channels requires accurate depictions of conformational ensembles, and the manner in which they interchange as a function of environmental factors including substrates, lipids, and inhibitors. Spectroscopic techniques such as electron spin resonance (ESR) pulsed electron–electron double resonance (PELDOR), also known as double electron–electron resonance (DEER), provide a complement to atomistic structures obtained from x-ray crystallography or cryo-EM, since spectroscopic data reflect an ensemble and can be measured in more native solvents, unperturbed by a crystal lattice. However, attempts to interpret DEER data are frequently stymied by discrepancies with the structural data, which may arise due to differences in conditions, the dynamics of the protein, or the flexibility of the attached paramagnetic spin labels. Recently, molecular simulation techniques such as EBMetaD have been developed that create a conformational ensemble matching an experimental distance distribution while applying the minimal possible bias. Moreover, it has been proposed that the work required during an EBMetaD simulation to match an experimentally determined distribution could be used as a metric with which to assign conformational states to a given measurement. Here, we demonstrate the application of this concept for a sodium-coupled transport protein, BetP. Because the probe, protein, and lipid bilayer are all represented in atomic detail, the different contributions to the work, such as the extent of protein backbone movements, can be separated. This work therefore illustrates how ranking simulations based on EBMetaD can help to bridge the gap between structural and biophysical data and thereby enhance our understanding of membrane protein conformational mechanisms.
Resonance assignments are challenging for membrane proteins due to the size of the lipid/detergent-protein complex and the presence of line-broadening from conformational exchange. As a consequence, many correlations are missing in the triple-resonance NMR experiments typically used for assignments. Herein, we present an approach in which correlations from these solution-state NMR experiments are supplemented by data from 13C unlabeling, single-amino acid type labeling, 4D NOESY data and proximity of moieties to lipids or water in combination with a structure of the protein. These additional data are used to edit the expected peaklists for the automated assignment protocol FLYA, a module of the program package CYANA. We demonstrate application of the protocol to the 262-residue proton pump from archaeal bacteriorhodopsin (bR) in lipid nanodiscs. The lipid-protein assembly is characterized by an overall correlation time of 44 ns. The protocol yielded assignments for 62% of all backbone (H, N, Cα, Cβ, C′) resonances of bR, corresponding to 74% of all observed backbone spin systems, and 60% of the Ala, Met, Ile (δ1), Leu and Val methyl groups, thus enabling to assign a large fraction of the protein without mutagenesis data. Most missing resonances stem from the extracellular half, likely due intermediate exchange line-broadening. Further analysis revealed that missing information of the amino acid type of the preceding residue is the largest problem, and that 4D NOESY experiments are particularly helpful to compensate for that information loss.
The facile synthesis and detailed investigation of a class of highly potent protease inhibitors based on 1,4-naphthoquinones with a dipeptidic recognition motif (HN-l-Phe-l-Leu-OR) in the 2-position and an electron-withdrawing group (EWG) in the 3-position is presented. One of the compound representatives, namely the acid with EWG = CN and with R = H proved to be a highly potent rhodesain inhibitor with nanomolar affinity. The respective benzyl ester (R = Bn) was found to be hydrolyzed by the target enzyme itself yielding the free acid. Detailed kinetic and mass spectrometry studies revealed a reversible covalent binding mode. Theoretical calculations with different density functionals (DFT) as well as wavefunction-based approaches were performed to elucidate the mode of action.
Current metabolomics approaches utilize cellular metabolite extracts, are destructive, and require high cell numbers. We introduce here an approach that enables the monitoring of cellular metabolism at lower cell numbers by observing the consumption/production of different metabolites over several kinetic data points of up to 48 hours. Our approach does not influence cellular viability, as we optimized the cellular matrix in comparison to other materials used in a variety of in‐cell NMR spectroscopy experiments. We are able to monitor real‐time metabolism of primary patient cells, which are extremely sensitive to external stress. Measurements are set up in an interleaved manner with short acquisition times (approximately 7 minutes per sample), which allows the monitoring of up to 15 patient samples simultaneously. Further, we implemented our approach for performing tracer‐based assays. Our approach will be important not only in the metabolomics fields, but also in individualized diagnostics.
Current metabolomics approaches utilize cellular metabolite extracts, are destructive, and require high cell numbers. We introduce here an approach that enables the monitoring of cellular metabolism at lower cell numbers by observing the consumption/production of different metabolites over several kinetic data points of up to 48 hours. Our approach does not influence cellular viability, as we optimized the cellular matrix in comparison to other materials used in a variety of in‐cell NMR spectroscopy experiments. We are able to monitor real‐time metabolism of primary patient cells, which are extremely sensitive to external stress. Measurements are set up in an interleaved manner with short acquisition times (approximately 7 minutes per sample), which allows the monitoring of up to 15 patient samples simultaneously. Further, we implemented our approach for performing tracer‐based assays. Our approach will be important not only in the metabolomics fields, but also in individualized diagnostics.
Amorphous formulation technologies to improve oral absorption of poorly soluble active pharmaceutical ingredients (APIs) have become increasingly prevalent. Currently, polymer-based amorphous formulations manufactured by spray drying, hot melt extrusion (HME), or co-precipitation are most common. However, these technologies have challenges in terms of the successful stabilization of poor glass former compounds in the amorphous form. An alternative approach is mesoporous silica, which stabilizes APIs in non-crystalline form via molecular adsorption inside nano-scale pores. In line with these considerations, two poor glass formers, haloperidol and carbamazepine, were formulated as polymer-based solid dispersion via HME and with mesoporous silica, and their stability was compared under accelerated conditions. Changes were monitored over three months with respect to solid-state form and dissolution. The results were supported by solid-state nuclear magnetic resonance spectroscopy (SS-NMR) and scanning electron microscopy (SEM). It was demonstrated that mesoporous silica was more successful than HME in the stabilization of the selected poor glass formers. While both drugs remained non-crystalline during the study using mesoporous silica, polymer-based HME formulations showed recrystallization after one week. Thus, mesoporous silica represents an attractive technology to extend the formulation toolbox to poorly soluble poor glass formers.
PELDOR (pulse electron-electron double resonance) is an established method to study intramolecular distances and can give evidence for conformational changes and flexibilities. However, it can also be used to study intermolecular interactions as for example oligerimization. Here, we used PELDOR to study the ‘end-to-end’ stacking of small double stranded (ds)RNAs. For this study, the dsRNA molecules were only singly labelled with the spin label TPA to avoid multi-spin effects and to measure only the intermolecular stacking interactions. It can be shown that small dsRNAs tend to assemble to rod-like structures due to π-π-interactions between the base pairs at the end of the strands. On the one hand, these interactions can influence or complicate measurements aimed at the determining of the structure and dynamics of the dsRNA molecule itself. On the other hand, it can be interesting to study such intermolecular stacking interactions in more detail, as for example their dependence on ion concentration. We quantitatively determined the stacking probability as a function of the monovalent NaCl salt and the dsRNA concentration. From this data the dissociation constant Kd was deduced and found to depend on the ratio between the NaCl salt and dsRNA concentrations. Additionally, the distances and distance distributions obtained predict a model for the stacking geometry of dsRNAs. Introducing a nucleotide overhangs at one end of the dsRNA molecule restricts the stacking to the other end, leading only to dimer formations. Introducing such an overhang at both ends of the dsRNA molecule fully suppresses stacking, as we could demonstrate by PELDOR experiments quantitatively.
In every established species, protein-protein interactions have evolved such that they are fit for purpose. However, the molecular details of the evolution of new protein-protein interactions are poorly understood. We have used nuclear magnetic resonance spectroscopy to investigate the changes in structure and dynamics during the evolution of a protein-protein interaction involving the intrinsically disordered CREBBP (CREB-binding protein) interaction domain (CID) and nuclear coactivator binding domain (NCBD) from the transcriptional coregulators NCOA (nuclear receptor coactivator) and CREBBP/p300, respectively. The most ancient low-affinity “Cambrian-like” [540 to 600 million years (Ma) ago] CID/NCBD complex contained less secondary structure and was more dynamic than the complexes from an evolutionarily younger “Ordovician-Silurian” fish ancestor (ca. 440 Ma ago) and extant human. The most ancient Cambrian-like CID/NCBD complex lacked one helix and several interdomain interactions, resulting in a larger solvent-accessible surface area. Furthermore, the most ancient complex had a high degree of millisecond-to-microsecond dynamics distributed along the entire sequences of both CID and NCBD. These motions were reduced in the Ordovician-Silurian CID/NCBD complex and further redistributed in the extant human CID/NCBD complex. Isothermal calorimetry experiments show that complex formation is enthalpically favorable and that affinity is modulated by a largely unfavorable entropic contribution to binding. Our data demonstrate how changes in structure and motion conspire to shape affinity during the evolution of a protein-protein complex and provide direct evidence for the role of structural, dynamic, and frustrational plasticity in the evolution of interactions between intrinsically disordered proteins.
1H-detected solid-state NMR experiments feasible at fast magic-angle spinning (MAS) frequencies allow accessing 1H chemical shifts of proteins in solids, which enables their interpretation in terms of secondary structure. Here we present 1H and 13C-detected NMR spectra of the RNA polymerase subunit Rpo7 in complex with unlabeled Rpo4 and use the 13C, 15N, and 1H chemical-shift values deduced from them to study the secondary structure of the protein in comparison to a known crystal structure. We applied the automated resonance assignment approach FLYA including 1H-detected solid-state NMR spectra and show its success in comparison to manual spectral assignment. Our results show that reasonably reliable secondary-structure information can be obtained from 1H secondary chemical shifts (SCS) alone by using the sum of 1Hα and 1HN SCS rather than by TALOS. The confidence, especially at the boundaries of the observed secondary structure elements, is found to increase when evaluating 13C chemical shifts, here either by using TALOS or in terms of 13C SCS.
The tetracycline-binding RNA aptamer (TC-aptamer) is a synthetic riboswitch that binds the antibiotic tetracycline (TC) with exceptionally high affinity. Although a crystal structure exists of the TC-bound state, little is known about the conformational dynamics and changes upon ligand binding. In this study, pulsed electron paramagnetic resonance techniques for measuring distances (PELDOR) in combination with rigid nitroxide spin labels (Çm spin label) were used to investigate the conformational flexibility of the TC-aptamer in the presence and absence of TC at different Mg2+ concentrations. TC was found to be the essential factor for stabilizing the tertiary structure at intermediate Mg2+ concentrations. At higher Mg2+ concentrations, Mg2+ alone is sufficient to stabilize the tertiary structure. In addition, the orientation of the two spin-labeled RNA helices with respect to each other was analyzed with orientation-selective PELDOR and compared to the crystal structure. These results demonstrate for the first time the unique value of the Çm spin label in combination with PELDOR to provide information about conformational flexibilities and orientations of secondary structure elements of biologically relevant RNAs.
Web spiders connect silk proteins, so-called spidroins, into fibers of extraordinary toughness. The spidroin N-terminal domain (NTD) plays a pivotal role in this process: it polymerizes spidroins through a complex mechanism of dimerization. Here we analyze sequences of spidroin NTDs and find an unusually high content of the amino acid methionine. We simultaneously mutate all methionines present in the hydrophobic core of a spidroin NTD from a nursery web spider’s dragline silk to leucine. The mutated NTD is strongly stabilized and folds at the theoretical speed limit. The structure of the mutant is preserved, yet its ability to dimerize is substantially impaired. We find that side chains of core methionines serve to mobilize the fold, which can thereby access various conformations and adapt the association interface for tight binding. Methionine in a hydrophobic core equips a protein with the capacity to dynamically change shape and thus to optimize its function.
As adapter molecules to convert the nucleic acid information into the amino acid sequence, tRNAs play a central role in protein synthesis. To fulfill this function in a reliable way, tRNAs exhibit highly conserved structural features common in all organisms and in all cellular compartments active in translation. However, in mitochondria of metazoans, certain dramatic deviations from the consensus tRNA structure are described, where some tRNAs lack the D- or T-arm without losing their function. In Enoplea, this miniaturization comes to an extreme, and functional mitochondrial tRNAs can lack both arms, leading to a considerable size reduction. Here, we investigate the secondary and tertiary structure of two such armless tRNAs from Romanomermis culicivorax. Despite their high AU content, the transcripts fold into a single and surprisingly stable hairpin structure, deviating from standard tRNAs. The three-dimensional form is boomerang-like and diverges from the standard L-shape. These results indicate that such unconventional miniaturized tRNAs can still fold into a tRNA-like shape, although their length and secondary structure are very unusual. They highlight the remarkable flexibility of the protein synthesis apparatus and suggest that the translational machinery of Enoplea mitochondria may show compensatory adaptations to accommodate these armless tRNAs for efficient translation.
The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop—a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.
LILBID and nESI : different native mass spectrometry techniques as tools in structural biology
(2018)
Native mass spectrometry is applied for the investigation of proteins and protein complexes worldwide. The challenge in native mass spectrometry is maintaining the features of the proteins of interest, such as oligomeric state, bound ligands, or the conformation of the protein complex, during transfer from solution to gas phase. This is an essential prerequisite to allow conclusions about the solution state protein complex, based on the gas phase measurements. Therefore, soft ionization techniques are required. Widely used for the analysis of protein complexes are nanoelectro spray ionization (nESI) mass spectrometers. A newer ionization method is laser induced liquid bead ion desorption (LILBID), which is based on the release of protein complexes from solution phase via infrared (IR) laser desorption. We use both methods in our lab, depending on the requirements of the biological system we are interested in. Here we benchmark the performance of our LILBID mass spectrometer in comparison to a nESI instrument, regarding sample conditions, buffer and additive tolerances, dissociation mechanism and applicability towards soluble and membrane protein complexes.
Global response of diacylglycerol kinase towards substrate binding observed by 2D and 3D MAS NMR
(2019)
Escherichia coli diacylglycerol kinase (DGK) is an integral membrane protein, which catalyses the ATP-dependent phosphorylation of diacylglycerol (DAG) to phosphatic acid (PA). It is a unique trimeric enzyme, which does not share sequence homology with typical kinases. It exhibits a notable complexity in structure and function despite of its small size. Here, chemical shift assignment of wild-type DGK within lipid bilayers was carried out based on 3D MAS NMR, utilizing manual and automatic analysis protocols. Upon nucleotide binding, extensive chemical shift perturbations could be observed. These data provide evidence for a symmetric DGK trimer with all of its three active sites concurrently occupied. Additionally, we could detect that the nucleotide substrate induces a substantial conformational change, most likely directing DGK into its catalytic active form. Furthermore, functionally relevant interprotomer interactions are identified by DNP-enhanced MAS NMR in combination with site-directed mutagenesis and functional assays.
HUWE1 E3 ligase promotes PINK1/PARKIN-independent mitophagy by regulating AMBRA1 activation via IKKα
(2018)
The selective removal of undesired or damaged mitochondria by autophagy, known as mitophagy, is crucial for cellular homoeostasis, and prevents tumour diffusion, neurodegeneration and ageing. The pro-autophagic molecule AMBRA1 (autophagy/beclin-1 regulator-1) has been defined as a novel regulator of mitophagy in both PINK1/PARKIN-dependent and -independent systems. Here, we identified the E3 ubiquitin ligase HUWE1 as a key inducing factor in AMBRA1-mediated mitophagy, a process that takes place independently of the main mitophagy receptors. Furthermore, we show that mitophagy function of AMBRA1 is post-translationally controlled, upon HUWE1 activity, by a positive phosphorylation on its serine 1014. This modification is mediated by the IKKα kinase and induces structural changes in AMBRA1, thus promoting its interaction with LC3/GABARAP (mATG8) proteins and its mitophagic activity. Altogether, these results demonstrate that AMBRA1 regulates mitophagy through a novel pathway, in which HUWE1 and IKKα are key factors, shedding new lights on the regulation of mitochondrial quality control and homoeostasis in mammalian cells.
Several peptides in clinical use are derived from non-ribosomal peptide synthetases (NRPS). In these systems multiple NRPS subunits interact with each other in a specific linear order mediated by specific docking domains (DDs), whose structures are not known yet, to synthesize well-defined peptide products. In contrast to classical NRPSs, single-module NRPS subunits responsible for the generation of rhabdopeptide/xenortide-like peptides (RXPs) can act in different order depending on subunit stoichiometry thereby producing peptide libraries. To define the basis for their unusual interaction patterns, we determine the structures of all N-terminal DDs (NDDs) as well as of an NDD-CDD complex and characterize all putative DD interactions thermodynamically for such a system. Key amino acid residues for DD interactions are identified that upon their exchange change the DD affinity and result in predictable changes in peptide production. Recognition rules for DD interactions are identified that also operate in other megasynthase complexes.
NMR and chromatography methods combined with mass spectrometry are the most important analytical techniques employed for plant metabolomics screening. Metabolomic analysis integrated to transcriptome screening add an important extra dimension to the information flow from DNA to RNA to protein. The most useful NMR experiment in metabolomics analysis is the proton spectra due the high receptivity of 1H and important structural information, through proton–proton scalar coupling. Routinely, databases have been used in identification of primary metabolites, however, there is currently no comparable data for identification of secondary metabolites, mainly, due to signal overlap in normal 1H NMR spectra and natural variation of plant. Related to spectra overlap, alternatively, better resolution can be find using 1H pure shift and 2D NMR pulse sequence in complex samples due to spreading the resonances in a second dimension. Thus, in data brief we provide a catalogue of metabolites and expression levels of genes identified in soy leaves and roots under flooding stress.
The ATP-binding cassette transporter TAPL translocates polypeptides from the cytosol into the lysosomal lumen. TAPL can be divided into two functional units: coreTAPL, active in ATP-dependent peptide translocation, and the N-terminal membrane spanning domain, TMD0, responsible for cellular localization and interaction with the lysosomal associated membrane proteins LAMP-1 and LAMP-2. Although the structure and function of ABC transporters were intensively analyzed in the past, the knowledge about accessory membrane embedded domains is limited. Therefore, we expressed the TMD0 of TAPL via a cell-free expression system and confirmed its correct folding by NMR and interaction studies. In cell as well as cell-free expressed TMD0 forms oligomers, which were assigned as dimers by PELDOR spectroscopy and static light scattering. By NMR spectroscopy of uniformly and selectively isotope labeled TMD0 we performed a complete backbone and partial side chain assignment. Accordingly, TMD0 has a four transmembrane helix topology with a short helical segment in a lysosomal loop. The topology of TMD0 was confirmed by paramagnetic relaxation enhancement with paramagnetic stearic acid as well as by nuclear Overhauser effects with c6-DHPC and cross-peaks with water.
The mitophagy receptor Nix interacts with LC3/GABARAP proteins, targeting mitochondria into autophagosomes for degradation. Here we present evidence for phosphorylation-driven regulation of the Nix:LC3B interaction. Isothermal titration calorimetry and NMR indicate a ~100 fold enhanced affinity of the serine 34/35-phosphorylated Nix LC3-interacting region (LIR) to LC3B and formation of a very rigid complex compared to the non-phosphorylated sequence. Moreover, the crystal structure of LC3B in complex with the Nix LIR peptide containing glutamic acids as phosphomimetic residues and NMR experiments revealed that LIR phosphorylation stabilizes the Nix:LC3B complex via formation of two additional hydrogen bonds between phosphorylated serines of Nix LIR and Arg11, Lys49 and Lys51 in LC3B. Substitution of Lys51 to Ala in LC3B abrogates binding of a phosphomimetic Nix mutant. Functionally, serine 34/35 phosphorylation enhances autophagosome recruitment to mitochondria in HeLa cells. Together, this study provides cellular, biochemical and biophysical evidence that phosphorylation of the LIR domain of Nix enhances mitophagy receptor engagement.
The field of dynamic nuclear polarization has undergone tremendous developments and diversification since its inception more than 6 decades ago. In this review we provide an in-depth overview of the relevant topics involved in DNP-enhanced MAS NMR spectroscopy. This includes the theoretical description of DNP mechanisms as well as of the polarization transfer pathways that can lead to a uniform or selective spreading of polarization between nuclear spins. Furthermore, we cover historical and state-of-the art aspects of dedicated instrumentation, polarizing agents, and optimization techniques for efficient MAS DNP. Finally, we present an extensive overview on applications in the fields of structural biology and materials science, which underlines that MAS DNP has moved far beyond the proof-of-concept stage and has become an important tool for research in these fields.
A key function of reversible protein phosphorylation is to regulate protein–protein interactions, many of which involve short linear motifs (3–12 amino acids). Motif‐based interactions are difficult to capture because of their often low‐to‐moderate affinities. Here, we describe phosphomimetic proteomic peptide‐phage display, a powerful method for simultaneously finding motif‐based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C‐terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD‐95/Dlg/ZO‐1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR. We uncover site‐specific phospho‐regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho‐regulation of motif‐based interactions on a large scale.
The p53 family of transcription factors (p53, p63 and p73) covers a wide range of functions critical for development, homeostasis and health of mammals across their lifespan. Beside the well-established tumor suppressor role, recent evidence has highlighted novel non-oncogenic functions exerted by p73. In particular, p73 is required for multiciliated cell (MCC) differentiation; MCCs have critical roles in brain and airways to move fluids across epithelial surfaces and to transport germ cells in the reproductive tract. This novel function of p73 provides a unifying cellular mechanism for the disparate inflammatory and immunological phenotypes of p73-deficient mice. Indeed, mice with Trp73 deficiency suffer from hydrocephalus, sterility and chronic respiratory tract infections due to profound defects in ciliogenesis and complete loss of mucociliary clearance since MCCs are essential for cleaning airways from inhaled pollutants, pathogens and allergens. Cross-species genomic analyses and functional rescue experiments identify TAp73 as the master transcriptional integrator of ciliogenesis, upstream of previously known central nodes. In addition, TAp73 shows a significant ability to regulate cellular metabolism and energy production through direct transcriptional regulation of several metabolic enzymes, such as glutaminase-2 and glucose-6 phosphate dehydrogenase. This recently uncovered role of TAp73 in the regulation of cellular metabolism strongly affects oxidative balance, thus potentially influencing all the biological aspects associated with p73 function, including development, homeostasis and cancer. Although through different mechanisms, p63 isoforms also contribute to regulation of cellular metabolism, thus indicating a common route used by all family members to control cell fate. At the structural level, the complexity of p73's function is further enhanced by its ability to form heterotetramers with some p63 isoforms, thus indicating the existence of an intrafamily crosstalk that determines the global outcome of p53 family function. In this review, we have tried to summarize all the recent evidence that have emerged on the novel non-oncogenic roles of p73, in an attempt to provide a unified view of the complex function of this gene within its family.
The bile acid activated transcription factor farnesoid X receptor (FXR) regulates numerous metabolic processes and is a rising target for the treatment of hepatic and metabolic disorders. FXR agonists have revealed efficacy in treating non-alcoholic steatohepatitis (NASH), diabetes and dyslipidemia. Here we characterize imatinib as first-in-class allosteric FXR modulator and report the development of an optimized descendant that markedly promotes agonist induced FXR activation in a reporter gene assay and FXR target gene expression in HepG2 cells. Differential effects of imatinib on agonist-induced bile salt export protein and small heterodimer partner expression suggest that allosteric FXR modulation could open a new avenue to gene-selective FXR modulators.
The identification of inhibitors of eukaryotic protein biosynthesis, which are targeting single translation factors, is highly demanded. Here we report on a small molecule inhibitor, gephyronic acid, isolated from the myxobacterium Archangium gephyra that inhibits growth of transformed mammalian cell lines in the nM range. In direct comparison, primary human fibroblasts were shown to be less sensitive to toxic effects of gephyronic acid than cancer-derived cells. Gephyronic acid is targeting the protein translation system. Experiments with IRES dual luciferase reporter assays identified it as an inhibitor of the translation initiation. DARTs approaches, co-localization studies and pull-down assays indicate that the binding partner could be the eukaryotic initiation factor 2 subunit alpha (eIF2α). Gephyronic acid seems to have a different mode of action than the structurally related polyketides tedanolide, myriaporone, and pederin and is a valuable tool for investigating the eukaryotic translation system. Because cancer derived cells were found to be especially sensitive, gephyronic acid could potentially find use as a drug candidate.
TEMPO spin labels protected with 2-nitrobenzyloxymethyl groups were attached to the amino residues of three different nucleosides: deoxycytidine, deoxyadenosine, and adenosine. The corresponding phosphoramidites could be incorporated by unmodified standard procedures into four different self-complementary DNA and two RNA oligonucleotides. After photochemical removal of the protective group, elimination of formic aldehyde and spontaneous air oxidation, the nitroxide radicals were regenerated in high yield. The resulting spin-labeled palindromic duplexes could be directly investigated by PELDOR spectroscopy without further purification steps. Spin–spin distances measured by PELDOR correspond well to the values obtained from molecular models.
Up to now, very small protein-coding genes have remained unrecognized in sequenced genomes. We identified an mRNA of 165 nucleotides (nt), which is conserved in Bradyrhizobiaceae and encodes a polypeptide with 14 amino acid residues (aa). The small mRNA harboring a unique Shine-Dalgarno sequence (SD) with a length of 17 nt was localized predominantly in the ribosome-containing P100 fraction of Bradyrhizobium japonicum USDA 110. Strong interaction between the mRNA and 30S ribosomal subunits was demonstrated by their co-sedimentation in sucrose density gradient. Using translational fusions with egfp, we detected weak translation and found that it is impeded by both the extended SD and the GTG start codon (instead of ATG). Biophysical characterization (CD- and NMR-spectroscopy) showed that synthesized polypeptide remained unstructured in physiological puffer. Replacement of the start codon by a stop codon increased the stability of the transcript, strongly suggesting additional posttranscriptional regulation at the ribosome. Therefore, the small gene was named rreB (ribosome-regulated expression in Bradyrhizobiaceae). Assuming that the unique ribosome binding site (RBS) is a hallmark of rreB homologs or similarly regulated genes, we looked for similar putative RBS in bacterial genomes and detected regions with at least 16 nt complementarity to the 3′-end of 16S rRNA upstream of sORFs in Caulobacterales, Rhizobiales, Rhodobacterales and Rhodospirillales. In the Rhodobacter/Roseobacter lineage of α-proteobacteria the corresponding gene (rreR) is conserved and encodes an 18 aa protein. This shows how specific RBS features can be used to identify new genes with presumably similar control of expression at the RNA level.
Investigating three-dimensional (3D) structures of proteins in living cells by in-cell nuclear magnetic resonance (NMR) spectroscopy opens an avenue towards understanding the structural basis of their functions and physical properties under physiological conditions inside cells. In-cell NMR provides data at atomic resolution non-invasively, and has been used to detect protein-protein interactions, thermodynamics of protein stability, the behavior of intrinsically disordered proteins, etc. in cells. However, so far only a single de novo 3D protein structure could be determined based on data derived only from in-cell NMR. Here we introduce methods that enable in-cell NMR protein structure determination for a larger number of proteins at concentrations that approach physiological ones. The new methods comprise (1) advances in the processing of non-uniformly sampled NMR data, which reduces the measurement time for the intrinsically short-lived in-cell NMR samples, (2) automatic chemical shift assignment for obtaining an optimal resonance assignment, and (3) structure refinement with Bayesian inference, which makes it possible to calculate accurate 3D protein structures from sparse data sets of conformational restraints. As an example application we determined the structure of the B1 domain of protein G at about 250 μM concentration in living E. coli cells.
Modification of SMN2 exon 7 (E7) splicing is a validated therapeutic strategy against spinal muscular atrophy (SMA). However, a target-based approach to identify small-molecule E7 splicing modifiers has not been attempted, which could reveal novel therapies with improved mechanistic insight. Here, we chose as a target the stem-loop RNA structure TSL2, which overlaps with the 5′ splicing site of E7. A small-molecule TSL2-binding compound, homocarbonyltopsentin (PK4C9), was identified that increases E7 splicing to therapeutic levels and rescues downstream molecular alterations in SMA cells. High-resolution NMR combined with molecular modelling revealed that PK4C9 binds to pentaloop conformations of TSL2 and promotes a shift to triloop conformations that display enhanced E7 splicing. Collectively, our study validates TSL2 as a target for small-molecule drug discovery in SMA, identifies a novel mechanism of action for an E7 splicing modifier, and sets a precedent for other splicing-mediated diseases where RNA structure could be similarly targeted.
Protein aggregation of the p63 transcription factor underlies severe skin fragility in AEC syndrome
(2018)
The p63 gene encodes a master regulator of epidermal commitment, development, and differentiation. Heterozygous mutations in the C-terminal domain of the p63 gene can cause ankyloblepharon-ectodermal defects-cleft lip/palate (AEC) syndrome, a life-threatening disorder characterized by skin fragility and severe, long-lasting skin erosions. Despite deep knowledge of p63 functions, little is known about mechanisms underlying disease pathology and possible treatments. Here, we show that multiple AEC-associated p63 mutations, but not those causative of other diseases, lead to thermodynamic protein destabilization, misfolding, and aggregation, similar to the known p53 gain-of-function mutants found in cancer. AEC mutant proteins exhibit impaired DNA binding and transcriptional activity, leading to dominant negative effects due to coaggregation with wild-type p63 and p73. Importantly, p63 aggregation occurs also in a conditional knock-in mouse model for the disorder, in which the misfolded p63 mutant protein leads to severe epidermal defects. Variants of p63 that abolish aggregation of the mutant proteins are able to rescue p63’s transcriptional function in reporter assays as well as in a human fibroblast-to-keratinocyte conversion assay. Our studies reveal that AEC syndrome is a protein aggregation disorder and opens avenues for therapeutic intervention.
Membrane proteins frequently assemble into higher order homo- or hetero-oligomers within their natural lipid environment. This complex formation can modulate their folding, activity as well as substrate selectivity. Non-disruptive methods avoiding critical steps, such as membrane disintegration, transfer into artificial environments or chemical modifications are therefore essential to analyze molecular mechanisms of native membrane protein assemblies. The combination of cell-free synthetic biology, nanodisc-technology and non-covalent mass spectrometry provides excellent synergies for the analysis of membrane protein oligomerization within defined membranes. We exemplify our strategy by oligomeric state characterization of various membrane proteins including ion channels, transporters and membrane-integrated enzymes assembling up to hexameric complexes. We further indicate a lipid-dependent dimer formation of MraY translocase correlating with the enzymatic activity. The detergent-free synthesis of membrane protein/nanodisc samples and the analysis by LILBID mass spectrometry provide a versatile platform for the analysis of membrane proteins in a native environment.
NMR spectroscopy is a powerful technique to study ribonucleic acids (RNAs) which are key players in a plethora of cellular processes. Although the NMR toolbox for structural studies of RNAs expanded during the last decades, they often remain challenging. Here, we show that solvent paramagnetic relaxation enhancements (sPRE) induced by the soluble, paramagnetic compound Gd(DTPA-BMA) provide a quantitative measure for RNA solvent accessibility and encode distance-to-surface information that correlates well with RNA structure and improves accuracy and convergence of RNA structure determination. Moreover, we show that sPRE data can be easily obtained for RNAs with any isotope labeling scheme and is advantageous regarding sample preparation, stability and recovery. sPRE data show a large dynamic range and reflect the global fold of the RNA suggesting that they are well suited to identify interaction surfaces, to score structural models and as restraints in RNA structure determination.
Structural and functional dissection of the DH and PH domains of oncogenic Bcr-Abl tyrosine kinase
(2017)
The two isoforms of the Bcr-Abl tyrosine kinase, p210 and p190, are associated with different leukemias and have a dramatically different signaling network, despite similar kinase activity. To provide a molecular rationale for these observations, we study the Dbl-homology (DH) and Pleckstrin-homology (PH) domains of Bcr-Abl p210, which constitute the only structural differences to p190. Here we report high-resolution structures of the DH and PH domains and characterize conformations of the DH–PH unit in solution. Our structural and functional analyses show no evidence that the DH domain acts as a guanine nucleotide exchange factor, whereas the PH domain binds to various phosphatidylinositol-phosphates. PH-domain mutants alter subcellular localization and result in decreased interactions with p210-selective interaction partners. Hence, the PH domain, but not the DH domain, plays an important role in the formation of the differential p210 and p190 Bcr-Abl signaling networks.
"Ästhetisch ist, was hilft"
(2017)
The full-length translation-regulating add adenine riboswitch (Asw) from Vibrio vulnificus has a more complex conformational space than its isolated aptamer domain. In addition to the predicted apo (apoA) and holo conformation that feature the conserved three-way junctional purine riboswitch aptamer, it adopts a second apo (apoB) conformation with a fundamentally different secondary structure. Here, we characterized the ligand-dependent conformational dynamics of the full-length add Asw by NMR and by single-molecule FRET (smFRET) spectroscopy. Both methods revealed an adenine-induced secondary structure switch from the apoB-form to the apoA-form that involves no tertiary structural interactions between aptamer and expression platform. This strongly suggests that the add Asw triggers translation by capturing the apoA-form secondary structure in the holo state. Intriguingly, NMR indicated a homogenous, docked aptamer kissing loop fold for apoA and holo, while smFRET showed persistent aptamer kissing loop docking dynamics between comparably stable, undocked and docked substates of the apoA and the holo conformation. Unraveling the folding of large junctional riboswitches thus requires the integration of complementary solution structural techniques such as NMR and smFRET.
ATP-binding cassette (ABC) transporters, a superfamily of integral membrane proteins, catalyse the translocation of substrates across the cellular membrane by ATP hydrolysis. Here we demonstrate by nucleotide turnover and binding studies based on 31P solid-state NMR spectroscopy that the ABC exporter and lipid A flippase MsbA can couple ATP hydrolysis to an adenylate kinase activity, where ADP is converted into AMP and ATP. Single-point mutations reveal that both ATPase and adenylate kinase mechanisms are associated with the same conserved motifs of the nucleotide-binding domain. Based on these results, we propose a model for the coupled ATPase-adenylate kinase mechanism, involving the canonical and an additional nucleotide-binding site. We extend these findings to other prokaryotic ABC exporters, namely LmrA and TmrAB, suggesting that the coupled activities are a general feature of ABC exporters.
Proton-pumping complex I of the mitochondrial respiratory chain is among the largest and most complex membrane protein complexes. The enzyme contributes substantially to oxidative energy-conversion in eukaryotic cells. Its malfunctions are implicated in many hereditary and degenerative disorders. Here, we report the X-ray structure of mitochondrial complex I at 3.6- 3.9 Å resolution describing in detail the central subunits that execute the bioenergetic function. A continuous axis of basic and acidic residues running centrally through the membrane arm connects the ubiquinone reduction site in the hydrophilic arm to four putative proton-pumping units. The binding position for a substrate analogous inhibitor and blockage of the predicted ubiquinone binding site provide a model for the ‘deactive’ form of the enzyme. The proposed transition into the active form is based on a concerted structural rearrangement at the ubiquinone reduction site rendering support for a two-state stabilization-change mechanism of protonpumping.
Mammalian oocytes are arrested in the dictyate stage of meiotic prophase I for long periods of time, during which the high concentration of the p53 family member TAp63α sensitizes them to DNA damage-induced apoptosis. TAp63α is kept in an inactive and exclusively dimeric state but undergoes rapid phosphorylation-induced tetramerization and concomitant activation upon detection of DNA damage. Here we show that the TAp63α dimer is a kinetically trapped state. Activation follows a spring-loaded mechanism not requiring further translation of other cellular factors in oocytes and is associated with unfolding of the inhibitory structure that blocks the tetramerization interface. Using a combination of biophysical methods as well as cell and ovary culture experiments we explain how TAp63α is kept inactive in the absence of DNA damage but causes rapid oocyte elimination in response to a few DNA double strand breaks thereby acting as the key quality control factor in maternal reproduction.
In bacteria, the regulation of gene expression by cis-acting transcriptional riboswitches located in the 5'-untranslated regions of messenger RNA requires the temporal synchronization of RNA synthesis and ligand binding-dependent conformational refolding. Ligand binding to the aptamer domain of the riboswitch induces premature termination of the mRNA synthesis of ligand-associated genes due to the coupled formation of 3'-structural elements acting as terminators. To date, there has been no high resolution structural description of the concerted process of synthesis and ligand-induced restructuring of the regulatory RNA element. Here, we show that for the guanine-sensing xpt-pbuX riboswitch from Bacillus subtilis, the conformation of the full-length transcripts is static: it exclusively populates the functional off-state but cannot switch to the on-state, regardless of the presence or absence of ligand. We show that only the combined matching of transcription rates and ligand binding enables transcription intermediates to undergo ligand-dependent conformational refolding.
Many naturally occurring or artificially created RNAs are capable of binding to guanine or guanine derivatives with high affinity and selectivity. They bind their ligands using very different recognition modes involving a diverse set of hydrogen bonding and stacking interactions. Apparently, the potential structural diversity for guanine, guanosine, and guanine nucleotide binding motifs is far from being fully explored. Szostak and coworkers have derived a large set of different GTP-binding aptamer families differing widely in sequence, secondary structure, and ligand specificity. The so-called class V–GTP aptamer from this set binds GTP with very high affinity and has a complex secondary structure. Here we use solution NMR spectroscopy to demonstrate that the class V aptamer binds GTP through the formation of an intermolecular two-layered G-quadruplex structure that directly incorporates the ligand and folds only upon ligand addition. Ligand binding and G-quadruplex formation depend strongly on the identity of monovalent cations present with a clear preference for potassium ions. GTP binding through direct insertion into an intermolecular G-quadruplex is a previously unobserved structural variation for ligand-binding RNA motifs and rationalizes the previously observed specificity pattern of the class V aptamer for GTP analogs.
The spliceosomal protein SF3b49, a component of the splicing factor 3b (SF3b) protein complex in the U2 small nuclear ribonucleoprotein, contains two RNA recognition motif (RRM) domains. In yeast, the first RRM domain (RRM1) of Hsh49 protein (yeast orthologue of human SF3b49) reportedly interacts with another component, Cus1 protein (orthologue of human SF3b145). Here, we solved the solution structure of the RRM1 of human SF3b49 and examined its mode of interaction with a fragment of human SF3b145 using NMR methods. Chemical shift mapping showed that the SF3b145 fragment spanning residues 598-631 interacts with SF3b49 RRM1, which adopts a canonical RRM fold with a topology of β1-α1-β2-β3-α2-β4. Furthermore, a docking model based on NOESY measurements suggests that residues 607-616 of the SF3b145 fragment adopt a helical structure that binds to RRM1 predominantly via α1, consequently exhibiting a helix-helix interaction in almost antiparallel. This mode of interaction was confirmed by a mutational analysis using GST pull-down assays. Comparison with structures of all RRM domains when complexed with a peptide found that this helix-helix interaction is unique to SF3b49 RRM1. Additionally, all amino acid residues involved in the interaction are well conserved among eukaryotes, suggesting evolutionary conservation of this interaction mode between SF3b49 RRM1 and SF3b145.
Mistakes in translation of messenger RNA into protein are clearly a detriment to the recombinant production of pure proteins for biophysical study or the biopharmaceutical market. However, they may also provide insight into mechanistic details of the translation process. Mistakes often involve the substitution of an amino acid having an abundant codon for one having a rare codon, differing by substitution of a G base by an A base, as in the case of substitution of a lysine (AAA) for arginine (AGA). In these cases one expects the substitution frequency to depend on the relative abundances of the respective tRNAs, and thus, one might expect frequencies to be similar for all sites having the same rare codon. Here we demonstrate that, for the ADP-ribosylation factor from yeast expressed in E. coli, lysine for arginine substitutions frequencies are not the same at the 9 sites containing a rare arginine codon; mis-incorporation frequencies instead vary from less than 1 to 16%. We suggest that the context in which the codons occur (clustering of rare sites) may be responsible for the variation. The method employed to determine the frequency of mis-incorporation involves a novel mass spectrometric analysis of the products from the parallel expression of wild type and codon-optimized genes in 15N and 14N enriched media, respectively. The high sensitivity and low material requirements of the method make this a promising technology for the collection of data relevant to other mis-incorporations. The additional data could be of value in refining models for the ribosomal translation elongation process.
Plant-released flavonoids induce the transcription of symbiotic genes in rhizobia and one of the first bacterial responses is the synthesis of so called Nod factors. They are responsible for the initial root hair curling during onset of root nodule development. This signal exchange is believed to be essential for initiating the plant symbiosis with rhizobia affiliated with the Alphaproteobacteria. Here, we provide evidence that in the broad host range strain Sinorhizobium fredii NGR234 the complete lack of quorum sensing molecules results in an elevated copy number of its symbiotic plasmid (pNGR234a). This in turn triggers the expression of symbiotic genes and the production of Nod factors in the absence of plant signals. Therefore, increasing the copy number of specific plasmids could be a widespread mechanism of specialized bacterial populations to bridge gaps in signaling cascades.
Relative orientation of POTRA domains from cyanobacterial Omp85 studied by pulsed EPR spectroscopy
(2016)
Many proteins of the outer membrane of Gram-negative bacteria and of the outer envelope of the endosymbiotically derived organelles mitochondria and plastids have a β-barrel fold. Their insertion is assisted by membrane proteins of the Omp85-TpsB superfamily. These proteins are composed of a C-terminal β-barrel and a different number of N-terminal POTRA domains, three in the case of cyanobacterial Omp85. Based on structural studies of Omp85 proteins, including the five POTRA-domain-containing BamA protein of Escherichia coli, it is predicted that anaP2 and anaP3 bear a fixed orientation, whereas anaP1 and anaP2 are connected via a flexible hinge. We challenged this proposal by investigating the conformational space of the N-terminal POTRA domains of Omp85 from the cyanobacterium Anabaena sp. PCC 7120 using pulsed electron-electron double resonance (PELDOR, or DEER) spectroscopy. The pronounced dipolar oscillations observed for most of the double spin-labeled positions indicate a rather rigid orientation of the POTRA domains in frozen liquid solution. Based on the PELDOR distance data, structure refinement of the POTRA domains was performed taking two different approaches: 1) treating the individual POTRA domains as rigid bodies; and 2) using an all-atom refinement of the structure. Both refinement approaches yielded ensembles of model structures that are more restricted compared to the conformational ensemble obtained by molecular dynamics simulations, with only a slightly different orientation of N-terminal POTRA domains anaP1 and anaP2 compared with the x-ray structure. The results are discussed in the context of the native environment of the POTRA domains in the periplasm.
We investigate complexes of two paramagnetic metal ions Gd3+ and Mn2+ to serve as polarizing agents for solid-state dynamic nuclear polarization (DNP) of 1H, 13C, and 15N at magnetic fields of 5, 9.4, and 14.1 T. Both ions are half-integer high-spin systems with a zero-field splitting and therefore exhibit a broadening of the mS = −1/2 ↔ +1/2 central transition which scales inversely with the external field strength. We investigate experimentally the influence of the chelator molecule, strong hyperfine coupling to the metal nucleus, and deuteration of the bulk matrix on DNP properties. At small Gd-DOTA concentrations the narrow central transition allows us to polarize nuclei with small gyromagnetic ratio such as 13C and even 15N via the solid effect. We demonstrate that enhancements observed are limited by the available microwave power and that large enhancement factors of >100 (for 1H) and on the order of 1000 (for 13C) can be achieved in the saturation limit even at 80 K. At larger Gd(III) concentrations (≥10 mM) where dipolar couplings between two neighboring Gd3+ complexes become substantial a transition towards cross effect as dominating DNP mechanism is observed. Furthermore, the slow spin-diffusion between 13C and 15N, respectively, allows for temporally resolved observation of enhanced polarization spreading from nuclei close to the paramagnetic ion towards nuclei further removed. Subsequently, we present preliminary DNP experiments on ubiquitin by site-directed spin-labeling with Gd3+ chelator tags. The results hold promise towards applications of such paramagnetically labeled proteins for DNP applications in biophysical chemistry and/or structural biology.
Structural biology and life sciences in general, and NMR in particular, have always been associated with advanced computing. The current challenges in the post-genomic era call for virtual research platforms that provide the worldwide research community with both user-friendly tools, platforms for data analysis and exchange, and an underlying e-Infrastructure. WeNMR, a three-year European Commission co-funded project started in November 2010, groups different research teams into a worldwide virtual research community. It builds on the established eNMR e-Infrastructure and its steadily growing virtual organisation, which is currently the second largest VO in the area of life sciences. WeNMR provides an e-Infrastructure platform and Science Gateway for structural biology. It involves researchers from around the world and will build bridges to other areas of structural biology.
Structured RNA regions are important gene control elements in prokaryotes and eukaryotes. Here, we show that the mRNA of a cyanobacterial heat shock gene contains a built-in thermosensor critical for photosynthetic activity under stress conditions. The exceptionally short 5´-untranslated region is comprised of a single hairpin with an internal asymmetric loop. It inhibits translation of the Synechocystis hsp17 transcript at normal growth conditions, permits translation initiation under stress conditions and shuts down Hsp17 production in the recovery phase. Point mutations that stabilized or destabilized the RNA structure deregulated reporter gene expression in vivo and ribosome binding in vitro. Introduction of such point mutations into the Synechocystis genome produced severe phenotypic defects. Reversible formation of the open and closed structure was beneficial for viability, integrity of the photosystem and oxygen evolution. Continuous production of Hsp17 was detrimental when the stress declined indicating that shutting-off heat shock protein production is an important, previously unrecognized function of RNA thermometers. We discovered a simple biosensor that strictly adjusts the cellular level of a molecular chaperone to the physiological need.
The mfl-riboswitch regulates expression of ribonucleotide reductase subunit in Mesoplasma florum by binding to 2´-deoxyguanosine and thereby promoting transcription termination. We characterized the structure of the ligand-bound aptamer domain by NMR spectroscopy and compared the mfl-aptamer to the aptamer domain of the closely related purine-sensing riboswitches. We show that the mfl-aptamer accommodates the extra 2´-deoxyribose unit of the ligand by forming a more relaxed binding pocket than these found in the purine-sensing riboswitches. Tertiary structures of the xpt-aptamer bound to guanine and of the mfl-aptamer bound to 2´-deoxyguanosine exhibit very similar features, although the sequence of the mfl-aptamer contains several alterations compared to the purine-aptamer consensus sequence. These alterations include the truncation of a hairpin loop which is crucial for complex formation in all purine-sensing riboswitches characterized to date. We further defined structural features and ligand binding requirements of the free mfl-aptamer and found that the presence of Mg2+ is not essential for complex formation, but facilitates ligand binding by promoting pre-organization of key structural motifs in the free aptamer.
Background: The automation of objectively selecting amino acid residue ranges for structure superpositions is important for meaningful and consistent protein structure analyses. So far there is no widely-used standard for choosing these residue ranges for experimentally determined protein structures, where the manual selection of residue ranges or the use of suboptimal criteria remain commonplace. Results: We present an automated and objective method for finding amino acid residue ranges for the superposition and analysis of protein structures, in particular for structure bundles resulting from NMR structure calculations. The method is implemented in an algorithm, CYRANGE, that yields, without protein-specific parameter adjustment, appropriate residue ranges in most commonly occurring situations, including low-precision structure bundles, multi-domain proteins, symmetric multimers, and protein complexes. Residue ranges are chosen to comprise as many residues of a protein domain that increasing their number would lead to a steep rise in the RMSD value. Residue ranges are determined by first clustering residues into domains based on the distance variance matrix, and then refining for each domain the initial choice of residues by excluding residues one by one until the relative decrease of the RMSD value becomes insignificant. A penalty for the opening of gaps favours contiguous residue ranges in order to obtain a result that is as simple as possible, but not simpler. Results are given for a set of 37 proteins and compared with those of commonly used protein structure validation packages. We also provide residue ranges for 6351 NMR structures in the Protein Data Bank. Conclusions: The CYRANGE method is capable of automatically determining residue ranges for the superposition of protein structure bundles for a large variety of protein structures. The method correctly identifies ordered regions. Global structure superpositions based on the CYRANGE residue ranges allow a clear presentation of the structure, and unnecessary small gaps within the selected ranges are absent. In the majority of cases, the residue ranges from CYRANGE contain fewer gaps and cover considerably larger parts of the sequence than those from other methods without significantly increasing the RMSD values. CYRANGE thus provides an objective and automatic method for standardizing the choice of residue ranges for the superposition of protein structures. Additional files Additional file 1: Dependence of Q on the order parameter rank. The quantity Qi is plotted against the order parameter rank i for 9 different protein structure bundles. Additional file 2: Dependence of P on the clustering stage. The quantity Pi is plotted against the clustering stage i for 9 different protein structure bundles. Additional file 3: Dependence of CYRANGE results on the minimal cluster size parameter my. The sequence coverage (red) and RMSD (blue) of the residue ranges determined by CYRANGE were plotted as a function of my for 9 different protein structure bundles. The dotted vertical line indicates the default value, my = 8. Where CYRANGE found two domains, the RMSD values of the individual domains are shown in light and dark blue. Additional file 4: Dependence of CYRANGE results on the domain boundary extension parameter m. See Additional File 3 for details. Additional file 5: Dependence of CYRANGE results on the minimal gap width g. See Additional File 3 for details. Additional file 6: Dependence of CYRANGE results on the relative RMSD decrease parameter delta. See Additional File 3 for details. Additional file 7: Dependence of CYRANGE results on the absolute RMSD decrease parameter delta abs. See Additional File 3 for details. Additional file 8: Dependence of CYRANGE results on the gap penalty parameter gamma. See Additional File 3 for details. Additional file 9: Correlation between the sequence coverage from CYRANGE, FindCore and PSVS, and the GDT total score, GDT_TS. Each data point represents a protein shown in Figures 3 and 4. The coverage is the percentage of amino acid residues included in the residue ranges found by the different methods. The GDT_TS value is defined by GDT_TS = (P1 + P2 + P4 + P8)/4, where Pd is the fraction of residues that can be superimposed under a distance cutoff of d Å. Additional file 10: Correlation between the RMSD value for the residue ranges from CYRANGE, FindCore and PSVS, and the GDT total score, GDT_TS. Each data point represents one protein domain. See Additional File 9 for details.
G-quadruplex topologies of telomeric repeat sequences from vertebrates were investigated in the presence of molecular crowding (MC) mimetics, namely polyethylene glycol 200 (PEG), Ficoll 70 as well as Xenopus laevis egg extract by CD and NMR spectroscopy and native PAGE. Here, we show that the conformational behavior of the telomeric repeats in X. laevis egg extract or in Ficoll is notably different from that observed in the presence of PEG. While the behavior of the telomeric repeat in X. laevis egg extract or in Ficoll resembles results obtained under dilute conditions, PEG promotes the formation of high-order parallel topologies. Our data suggest that PEG should not be used as a MC mimetic.
Human Transformer2-beta (hTra2-beta) is an important member of the serine/arginine-rich protein family, and contains one RNA recognition motif (RRM). It controls the alternative splicing of several pre-mRNAs, including those of the calcitonin/calcitonin gene-related peptide (CGRP), the survival motor neuron 1 (SMN1) protein and the tau protein. Accordingly, the RRM of hTra2-beta specifically binds to two types of RNA sequences [the CAA and (GAA)2 sequences]. We determined the solution structure of the hTra2-beta RRM (spanning residues Asn110–Thr201), which not only has a canonical RRM fold, but also an unusual alignment of the aromatic amino acids on the beta-sheet surface. We then solved the complex structure of the hTra2-beta RRM with the (GAA)2 sequence, and found that the AGAA tetra-nucleotide was specifically recognized through hydrogen-bond formation with several amino acids on the N- and C-terminal extensions, as well as stacking interactions mediated by the unusually aligned aromatic rings on the beta-sheet surface. Further NMR experiments revealed that the hTra2-beta RRM recognizes the CAA sequence when it is integrated in the stem-loop structure. This study indicates that the hTra2-beta RRM recognizes two types of RNA sequences in different RNA binding modes.
We present here a set of 13C-direct detected NMR experiments to facilitate the resonance assignment of RNA oligonucleotides. Three experiments have been developed: (1) the (H)CC-TOCSY-experiment utilizing a virtual decoupling scheme to assign the intraresidual ribose 13C-spins, (2) the (H)CPC-experiment that correlates each phosphorus with the C40 nuclei of adjacent nucleotides via J(C,P) couplings and (3) the (H)CPC-CCH-TOCSY-experiment that correlates the phosphorus nuclei with the respective C10,H10 ribose signals. The experiments were applied to two RNA hairpin structures. The current set of 13C-direct detected experiments allows direct and unambiguous assignment of the majority of the hetero nuclei and the identification of the individual ribose moieties following their sequential assignment. Thus, 13C-direct detected NMR methods constitute useful complements to the conventional 1H-detected approach for the resonance assignment of oligonucleotides that is often hindered by the limited chemical shift dispersion. The developed methods can also be applied to large deuterated RNAs. Keywords: NMR spectroscopy , Direct carbon , detection , RNA
Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain
(2007)
Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg2+-ions. Furthermore, a role for Mg2+-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg2+, but becomes partially pre-organized for ligand binding in the presence of Mg2+. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg2+ is more pronounced.
High-resolution NMR structure of an RNA model system : the 14-mer cUUCGg tetraloop hairpin RNA
(2009)
We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 Å) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2'-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.
In prokaryotes, RNA thermometers regulate a number of heat shock and virulence genes. These temperature sensitive RNA elements are usually located in the 5'-untranslated regions of the regulated genes. They repress translation initiation by base pairing to the Shine–Dalgarno sequence at low temperatures. We investigated the thermodynamic stability of the temperature labile hairpin 2 of the Salmonella fourU RNA thermometer over a broad temperature range and determined free energy, enthalpy and entropy values for the base-pair opening of individual nucleobases by measuring the temperature dependence of the imino proton exchange rates via NMR spectroscopy. Exchange rates were analyzed for the wild-type (wt) RNA and the A8C mutant. The wt RNA was found to be stabilized by the extraordinarily stable G14–C25 base pair. The mismatch base pair in the wt RNA thermometer (A8–G31) is responsible for the smaller cooperativity of the unfolding transition in the wt RNA. Enthalpy and entropy values for the base-pair opening events exhibit linear correlation for both RNAs. The slopes of these correlations coincide with the melting points of the RNAs determined by CD spectroscopy. RNA unfolding occurs at a temperature where all nucleobases have equal thermodynamic stabilities. Our results are in agreement with a consecutive zipper-type unfolding mechanism in which the stacking interaction is responsible for the observed cooperativity. Furthermore, remote effects of the A8C mutation affecting the stability of nucleobase G14 could be identified. According to our analysis we deduce that this effect is most probably transduced via the hydration shell of the RNA.
Nep1 (Emg1) is a highly conserved nucleolar protein with an essential function in ribosome biogenesis. A mutation in the human Nep1 homolog causes Bowen–Conradi syndrome—a severe developmental disorder. Structures of Nep1 revealed a dimer with a fold similar to the SPOUT-class of RNA-methyltransferases suggesting that Nep1 acts as a methyltransferase in ribosome biogenesis. The target for this putative methyltransferase activity has not been identified yet. We characterized the RNA-binding specificity of Methanocaldococcus jannaschii Nep1 by fluorescence- and NMR-spectroscopy as well as by yeast three-hybrid screening. Nep1 binds with high affinity to short RNA oligonucleotides corresponding to nt 910–921 of M. jannaschii 16S rRNA through a highly conserved basic surface cleft along the dimer interface. Nep1 only methylates RNAs containing a pseudouridine at a position corresponding to a previously identified hypermodified N1-methyl-N3-(3-amino-3-carboxypropyl) pseudouridine (m1acp3-Psi) in eukaryotic 18S rRNAs. Analysis of the methylated nucleoside by MALDI-mass spectrometry, HPLC and NMR shows that the methyl group is transferred to the N1 of the pseudouridine. Thus, Nep1 is the first identified example of an N1-specific pseudouridine methyltransferase. This enzymatic activity is also conserved in human Nep1 suggesting that Nep1 is the methyltransferase in the biosynthesis of m1acp3-Psi in eukaryotic 18S rRNAs.
The transcription factor p63 is expressed as at least six different isoforms, of which two have been assigned critical biological roles within ectodermal development and skin stem cell biology on the one hand and supervision of the genetic stability of oocytes on the other hand. These two isoforms contain a C-terminal inhibitory domain that negatively regulates their transcriptional activity. This inhibitory domain contains two individual components: one that uses an internal binding mechanism to interact with and mask the transactivation domain and one that is based on sumoylation. We have carried out an extensive alanine scanning study to identify critical regions within the inhibitory domain. These experiments show that a stretch of ~13 amino acids is crucial for the binding function. Further, investigation of transcriptional activity and the intracellular level of mutants that cannot be sumoylated suggests that sumoylation reduces the concentration of p63. We therefore propose that the inhibitory function of the C-terminal domain is in part due to direct inhibition of the transcriptional activity of the protein and in part due to indirect inhibition by controlling the concentration of p63. Keywords: p63, transcriptional regulation, auto-inhibition, sumoylation
The CUG-binding protein 1 (CUG-BP1) is a member of the CUG-BP1 and ETR-like factors (CELF) family or the Bruno-like family and is involved in the control of splicing, translation and mRNA degradation. Several target RNA sequences of CUG-BP1 have been predicted, such as the CUG triplet repeat, the GU-rich sequences and the AU-rich element of nuclear pre-mRNAs and/or cytoplasmic mRNA. CUG-BP1 has three RNA-recognition motifs (RRMs), among which the third RRM (RRM3) can bind to the target RNAs on its own. In this study, we solved the solution structure of the CUG-BP1 RRM3 by hetero-nuclear NMR spectroscopy. The CUG-BP1 RRM3 exhibited a noncanonical RRM fold, with the four-stranded b-sheet surface tightly associated with the N-terminal extension. Furthermore, we determined the solution structure of the CUG-BP1 RRM3 in the complex with (UG)3 RNA, and discovered that the UGU trinucleotide is specifically recognized through extensive stacking interactions and hydrogen bonds within the pocket formed by the b-sheet surface and the N-terminal extension. This study revealed the unique mechanism that enables the CUG-BP1 RRM3 to discriminate the short RNA segment from other sequences, thus providing the molecular basis for the comprehension of the role of the RRM3s in the CELF/Bruno-like family.
Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA's secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent ({approx}1 µs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.
A mild synthetic method for N-formyl-Met-Leu-Phe-OH (1) is described. After Fmoc solid phase peptide synthesis, on-bead formylation and HPLC purification, more than 30 mg of the fully 13C/15N-labelled tripeptide 1 could be isolated in a typical batch. This peptide can be easily crystallised and is therefore well suited as a standard sample for setting up solid-state NMR experiments.
In a combined NMR/MD study, the temperature-dependent changes in the conformation of two members of the RNA YNMG-tetraloop motif (cUUCGg and uCACGg) have been investigated at temperatures of 298, 317 and 325 K. The two members have considerable different thermal stability and biological functions. In order to address these differences, the combined NMR/MD study was performed. The large temperature range represents a challenge for both, NMR relaxation analysis (consistent choice of effective bond length and CSA parameter) and all-atom MD simulation with explicit solvent (necessity to rescale the temperature). A convincing agreement of experiment and theory is found. Employing a principle component analysis of the MD trajectories, the conformational distribution of both hairpins at various temperatures is investigated. The ground state conformation and dynamics of the two tetraloops are indeed found to be very similar. Furthermore, both systems are initially destabilized by a loss of the stacking interactions between the first and the third nucleobase in the loop region. While the global fold is still preserved, this initiation of unfolding is already observed at 317 K for the uCACGg hairpin but at a significantly higher temperature for the cUUCGg hairpin.
The degradation of the poly(A) tail is crucial for posttranscriptional gene regulation and for quality control of mRNA. Poly(A)-specific ribonuclease (PARN) is one of the major mammalian 3’ specific exo-ribonucleases involved in the degradation of the mRNA poly(A) tail, and it is also involved in the regulation of translation in early embryonic development. The interaction between PARN and the m7GpppG cap of mRNA plays a key role in stimulating the rate of deadenylation. Here we report the solution structures of the cap-binding domain of mouse PARN with and without the m7GpppG cap analog. The structure of the cap-binding domain adopts the RNA recognition motif (RRM) with a characteristic a-helical extension at its C-terminus, which covers the b-sheet surface (hereafter referred to as PARN RRM). In the complex structure of PARN RRM with the cap analog, the base of the N7-methyl guanosine (m7G) of the cap analog stacks with the solvent-exposed aromatic side chain of the distinctive tryptophan residue 468, located at the C-terminal end of the second b-strand. These unique structural features in PARN RRM reveal a novel cap-binding mode, which is distinct from the nucleotide recognition mode of the canonical RRM domains.
Oral presentations Background: We selected peptide ligands for the HIV-1 packaging signal PSI by screening phage displayed peptide libraries. Peptide ligands were optimized by screening spot synthesis peptide membranes. The aim of this study is the functional characterization of these peptide ligands with respect to inhibition of HIV-1 replication. Methods: Phage displayed peptide libraries were screened with PSI-RNA structures. The Trp-rich peptide motifs were optimized for specific binding on spot synthesis peptide membranes. The best binding peptide was expressed intracellularly in fusion with RFP or linked to a protein transduction domain (PTD) for intracellular delivery. The effects on virion production were analyzed using pseudotyped lentiviral particles. Results: After positive and negative selection rounds, phages binding specifically to PSI-RNA were identified by ELISA. Peptide inserts contained conserved motifs of aromatic amino acids known to be implicated in binding of PSI-RNA by the natural Gag ligand. The filter assay identified HKWPWW as the best binding ligand for PSI-RNA, which is delivered into several cell lines by addition of a PTD. Compared to a control peptide, the HKWPWW peptide inhibited HIV-1 replication as deduced from reduced titers of culture supernatants. As HKWPWW also binds to the TAR-RNA like the natural nucleocapsid PSI-RNA ligand, the effect on Tat-TAR inhibition will also be analyzed. Currently T-cell lines are established which stably express HKWPWW as well as a control peptide, which will be infected with HIV-1 to monitor the ability of HKWPWW to inhibit wild type HIV-1 replication. Conclusion: The selection of a peptide ligand for PSI-RNA able to inhibit HIV-1 replication proves the suitability of the phage display technology for the selection of peptides binding to RNA-structures. This enables the indentification of peptides serving as leads to interfere with additional targets in the HIV-1 replication cycle.
Ribosomal proteins are assumed to stabilize specific RNA structures and promote compact folding of the large rRNA. The conformational dynamics of the protein between the bound and unbound state play an important role in the binding process. We have studied those dynamical changes in detail for the highly conserved complex between the ribosomal protein L11 and the GTPase region of 23S rRNA. The RNA domain is compactly folded into a well defined tertiary structure, which is further stabilized by the association with the C-terminal domain of the L11 protein (L11ctd). In addition, the N-terminal domain of L11 (L11ntd) is implicated in the binding of the natural thiazole antibiotic thiostrepton, which disrupts the elongation factor function. We have studied the conformation of the ribosomal protein and its dynamics by NMR in the unbound state, the RNA bound state and in the ternary complex with the RNA and thiostrepton. Our data reveal a rearrangement of the L11ntd, placing it closer to the RNA after binding of thiostrepton, which may prevent binding of elongation factors. We propose a model for the ternary L11–RNA–thiostrepton complex that is additionally based on interaction data and conformational information of the L11 protein. The model is consistent with earlier findings and provides an explanation for the role of L11ntd in elongation factor binding.