Probing RNA dynamics via longitudinal exchange and CPMG relaxation dispersion NMR spectroscopy using a sensitive 13C-methyl label.

The refolding kinetics of bistable RNA sequences were studied in unperturbed equilibrium via (13)C exchange NMR spectroscopy. For this purpose a straightforward labeling technique was elaborated using a 2'-(13)C-methoxy uridine modification, which was prepared by a two-step synthesis and introduced into RNA using standard protocols. Using (13)C longitudinal exchange NMR spectroscopy the refolding kinetics of a 20 nt bistable RNA were characterized at temperatures between 298 and 310K, yielding the enthalpy and entropy differences between the conformers at equilibrium and the activation energy of the refolding process. The kinetics of a more stable 32 nt bistable RNA could be analyzed by the same approach at elevated temperatures, i.e. at 314 and 316 K. Finally, the dynamics of a multi-stable RNA able to fold into two hairpin- and a pseudo-knotted conformation was studied by (13)C relaxation dispersion NMR spectroscopy.

Structural insights into the dynamics and function of the C-terminus of the E. coli RNA chaperone Hfq.

The hexameric Escherichia coli RNA chaperone Hfq (Hfq(Ec)) is involved in riboregulation of target mRNAs by small trans-encoded RNAs. Hfq proteins of different bacteria comprise an evolutionarily conserved core, whereas the C-terminus is variable in length. Although the structure of the conserved core has been elucidated for several Hfq proteins, no structural information has yet been obtained for the C-terminus. Using bioinformatics, nuclear magnetic resonance spectroscopy, synchrotron radiation circular dichroism (SRCD) spectroscopy and small angle X-ray scattering we provide for the first time insights into the conformation and dynamic properties of the C-terminal extension of Hfq(Ec). These studies indicate that the C-termini are flexible and extend laterally away from the hexameric core, displaying in this way features typical of intrinsically disordered proteins that facilitate intermolecular interactions. We identified a minimal, intrinsically disordered region of the C-terminus supporting the interactions with longer RNA fragments. This minimal region together with rest of the C-terminal extension provides a flexible moiety capable of tethering long and structurally diverse RNA molecules. Furthermore, SRCD spectroscopy supported the hypothesis that RNA fragments exceeding a certain length interact with the C-termini of Hfq(Ec).

Mathematical treatment of adiabatic fast passage pulses for the computation of nuclear spin relaxation rates in proteins with conformational exchange.

Although originally designed for broadband inversion and decoupling in NMR spectroscopy, recent methodological developments have introduced adiabatic fast passage (AFP) pulses into the field of protein dynamics. AFP pulses employ a frequency sweep, and have not only superior inversion properties with respect to offset effects, but they are also easily implemented into a pulse sequence. As magnetization is dragged from the +z to the -z direction, Larmor precession is impeded since magnetization becomes spin-locked, which is a potentially useful feature for the investigation of microsecond to millisecond dynamics. A major drawback of these pulses as theoretical prediction is concerned, however, results from their time-dependent offset: simulations of spin density matrices under the influence of a time-dependent Hamiltonian with non-commuting elements are costly in terms of computational time, rendering data analysis impracticable. In this paper we suggest several ways to reduce the computational time without compromising accuracy with respect to effects such as cross-correlated relaxation and modulation of the chemical shift.

Longitudinal exchange: an alternative strategy towards quantification of dynamics parameters in ZZ exchange spectroscopy.

Longitudinal exchange experiments facilitate the quantification of the rates of interconversion between the exchanging species, along with their longitudinal relaxation rates, by analyzing the time-dependence of direct correlation and exchange cross peaks. Here we present a simple and robust alternative to this strategy, which is based on the combination of two complementary experiments, one with and one without resolving exchange cross peaks. We show that by combining the two data sets systematic errors that are caused by differential line-broadening of the exchanging species are avoided and reliable quantification of kinetic and relaxation parameters in the presence of additional conformational exchange on the ms-µs time scale is possible. The strategy is applied to a bistable DNA oligomer that displays different line-broadening in the two exchanging species.

Bias-invariant RNA-sequencing metadata annotation.

BACKGROUND: Recent technological advances have resulted in an unprecedented increase in publicly available biomedical data, yet the reuse of the data is often precluded by experimental bias and a lack of annotation depth and consistency. Missing annotations makes it impossible for researchers to find datasets specific to their needs. FINDINGS: Here, we investigate RNA-sequencing metadata prediction based on gene expression values. We present a deep-learning-based domain adaptation algorithm for the automatic annotation of RNA-sequencing metadata. We show, in multiple experiments, that our model is better at integrating heterogeneous training data compared with existing linear regression-based approaches, resulting in improved tissue type classification. By using a model architecture similar to Siamese networks, the algorithm can learn biases from datasets with few samples. CONCLUSION: Using our novel domain adaptation approach, we achieved metadata annotation accuracies up to 15.7% better than a previously published method. Using the best model, we provide a list of >10,000 novel tissue and sex label annotations for 8,495 unique SRA samples. Our approach has the potential to revive idle datasets by automated annotation making them more searchable.

Pharmacophore mapping via cross-relaxation during adiabatic fast passage.

A novel NMR method is demonstrated for the investigation of protein ligand interactions. In this approach an adiabatic fast passage pulse, i.e. a long, weak pulse with a linear frequency sweep, is used to probe (1)H-(1)H NOEs. During the adiabatic fast passage the effective rotating-frame NOE is a weighted average of transverse and longitudinal cross-relaxation contributions that can be tuned by pulse power and frequency sweep rate. It is demonstrated that the occurrence of spin diffusion processes leads to sizable deviations from the theoretical relationship between effective relaxation rate and effective tilt angle in the spin lock frame and can be used to probe protein-ligand binding. This methodology comprises high sensitivity and ease of implementation. The feasibility of this technique is demonstrated with two protein complexes, vanillic acid bound to the quail lipocalin Q83 and NAD(+) and AMP binding to alcohol dehydrogenase (ADH).

Deep learning-based cell composition analysis from tissue expression profiles.

We present Scaden, a deep neural network for cell deconvolution that uses gene expression information to infer the cellular composition of tissues. Scaden is trained on single-cell RNA sequencing (RNA-seq) data to engineer discriminative features that confer robustness to bias and noise, making complex data preprocessing and feature selection unnecessary. We demonstrate that Scaden outperforms existing deconvolution algorithms in both precision and robustness. A single trained network reliably deconvolves bulk RNA-seq and microarray, human and mouse tissue expression data and leverages the combined information of multiple datasets. Because of this stability and flexibility, we surmise that deep learning will become an algorithmic mainstay for cell deconvolution of various data types. Scaden's software package and web application are easy to use on new as well as diverse existing expression datasets available in public resources, deepening the molecular and cellular understanding of developmental and disease processes.

Direct observation of the dynamic process underlying allosteric signal transmission.

Allosteric regulation is an effective mechanism of control in biological processes. In allosteric proteins a signal originating at one site in the molecule is communicated through the protein structure to trigger a specific response at a remote site. Using NMR relaxation dispersion techniques we directly observe the dynamic process through which the KIX domain of CREB binding protein communicates allosteric information between binding sites. KIX mediates cooperativity between pairs of transcription factors through binding to two distinct interaction surfaces in an allosteric manner. We show that binding the activation domain of the mixed lineage leukemia (MLL) transcription factor to KIX induces a redistribution of the relative populations of KIX conformations toward a high-energy state in which the allosterically activated second binding site is already preformed, consistent with the Monod-Wyman-Changeux (WMC) model of allostery. The structural rearrangement process that links the two conformers and by which allosteric information is communicated occurs with a time constant of 3 ms at 27 degrees C. Our dynamic NMR data reveal that an evolutionarily conserved network of hydrophobic amino acids constitutes the pathway through which information is transmitted.

Autocorrelation analysis of NOESY data provides residue compactness for folded and unfolded proteins.

A novel spectral entropy interpretation for protein NOESY data is presented for the investigation of the spatial distribution of residues in protein structures without the requirement of NOE cross peak assignments. In this approach individual traces S(i)(omega) from a 3D (15)N NOESY-HSQC taken at frequency positions corresponding to different amide groups (residue position i) are subjected to a self-convolution procedure thus leading to the autocorrelation function C(i)(omega) of the NOESY-trace for a particular backbone residue position. The characteristic spatial surrounding of a particular residue position is reflected in the corresponding autocorrelation function and can be quantified by taking the (spectral) entropy S(nu) as an information measure. The feasibility of this novel approach is demonstrated with applications to the proteins Cyclophilin D and Osteopontin and the protein complex between the lipocalin Q83 and the bacterial siderophore Enterobactin. Typically, large entropy values were found for residues located in structurally loosely defined regions, whereas small entropy values were found for residues in hydrophobic core regions of the protein with tightly interacting side chains and distinct chemical shift patterns. The applications to the unfolded Osteopontin and the Q83/Enterobactin protein complex indicated that both local compaction of the polypeptide chain due to transiently formed structural elements and subtle changes in side-chain packing can be efficiently probed by this novel approach.

Cell transformation by the v-myc oncogene abrogates c-Myc/Max-mediated suppression of a C/EBP beta-dependent lipocalin gene.

Using differential hybridization techniques, a cDNA clone (Q83) was isolated that corresponds to a highly abundant mRNA in quail embryo fibroblasts transformed by the v-myc oncogene. The deduced 178 amino acid protein product of Q83 contains an N-terminal signal sequence and a lipocalin sequence motif, the hallmark of a family of secretory proteins binding and transporting small hydrophobic molecules of diverse biological function, including retinoids and steroids. The quail Q83 protein displays 87% sequence identity with a developmentally regulated chicken protein, termed p20K or Ch21. Cell transformation specifically by v-myc, but not by other oncogenic agents, induces high-level expression of Q83 mRNA and of the Q83 protein. Nucleotide sequence analysis and transcriptional mapping revealed that the Q83 gene encompasses seven exons with the coding region confined to exons 1 through 6. The promoter region contains consensus binding sites for the transcriptional regulators Myc and C/EBP beta. Transcriptional activation of Q83 is principally dependent on C/EBP beta, but is blocked in normal cells by the endogenous c-Myc/Max/Mad transcription factor network. In v-myc-transformed cells, high-level expression of the v-Myc protein and formation of highly stable v-Myc/Max heterodimers leads to abrogation of Q83 gene suppression and activation by C/EBP beta. A 157 amino acid residue recombinant protein representing the secreted form of Q83 was used for structure determination by nuclear magnetic resonance spectroscopy. Q83 folds into a single globular domain of the lipocalin-type. The central part consists of an eight-stranded up-and-down beta-barrel core flanked by an N-terminal 3(10)-like helix and a C-terminal alpha-helix. The orientation of the C-terminal alpha-helix is partially determined by a disulfide bridge between Cys59 and Cys152. The three-dimensional structure determination of the Q83 protein will facilitate the identification of its authentic ligand and the assessment of its biological function, including the putative role in myc-induced cell transformation.

A novel paramagnetic relaxation enhancement tag for nucleic acids: a tool to study structure and dynamics of RNA.

In this work, we present a novel 2,2,6,6-tetramethylpiperidine 1-oxyl (TEMPO) radical phosphoramidite building block, which can be attached to the 5'-terminus of nucleic acids. To investigate the paramagnetic relaxation enhancement (PRE) emanating from this radical center, we incorporated the TEMPO label into various types of RNAs. We measured proton PREs for selectively (13)C-isotope labeled nucleotides to derive long-range distance restraints in a short 15 nucleotide stem-loop model system, underscoring the potential of the 5'-TEMPO tag to determine long-range distance restraints for solution structure determination. We subsequently applied the distance-dependent relaxation enhancement induced by the nitroxide radical to discern two folding states in a bistable RNA. Finally, we investigated the fast conformational sampling of the HIV-1 TAR RNA, a paradigm for structural flexibility in nucleic acids. With PRE NMR in combination with molecular dynamics simulations, the structural plasticity of this RNA was analyzed in the absence and presence of the ligand L-argininamide.

Direct methods and residue type specific isotope labeling in NMR structure determination and model-driven sequential assignment.

Direct methods in NMR based structure determination start from an unassigned ensemble of unconnected gaseous hydrogen atoms. Under favorable conditions they can produce low resolution structures of proteins. Usually a prohibitively large number of NOEs is required, to solve a protein structure ab-initio, but even with a much smaller set of distance restraints low resolution models can be obtained which resemble a protein fold. One problem is that at such low resolution and in the absence of a force field it is impossible to distinguish the correct protein fold from its mirror image. In a hybrid approach these ambiguous models have the potential to aid in the process of sequential backbone chemical shift assignment when (13)C(beta) and (13)C' shifts are not available for sensitivity reasons. Regardless of the overall fold they enhance the information content of the NOE spectra. These, combined with residue specific labeling and minimal triple-resonance data using (13)C(alpha) connectivity can provide almost complete sequential assignment. Strategies for residue type specific labeling with customized isotope labeling patterns are of great advantage in this context. Furthermore, this approach is to some extent error-tolerant with respect to data incompleteness, limited precision of the peak picking, and structural errors caused by misassignment of NOEs.

Generation and relaxation of high rank coherences in AX3 systems in a selectively methionine labelled SH2 domain.

The usefulness of selective isotope labelling patterns is demonstrated using the C-terminal SH2 domain of PLC-gamma1 selectively 13C labelled at methionine methyl groups. We demonstrate the generation and relaxation of coherences that are second rank in protons and first rank in carbons that derive from quadrupolar order in protons. The decay rates of second rank double quantum proton coherences are measured. These terms exhibit fewer channels for cross-correlated relaxation compared to single quantum coherences. Our results indicate the potential application of the measurement of high order proton coherences to the analysis of dynamics in methyl-bearing side chains.

An isolated helix persists in a sparsely populated form of KIX under native conditions.

NMR relaxation dispersion techniques were used to investigate conformational exchange of the three-helix bundle protein KIX under native conditions. These experiments provide site-resolved kinetic information about microsecond-to-millisecond time scale motions along with structural (chemical shift) information without requiring a perturbation of the equilibrium. All kinetic data are consistent with an apparent two-state transition between natively folded KIX and a partially unfolded high-energy state that is populated to 3.0 +/- 0.2% at 27 degrees C. By combining (13)C- and (15)N-based experiments that probe specific structural aspects, we show that the sparsely populated high-energy state displays a strong conformational preference. An isolated secondary structural element, C-terminal helix alpha3, is highly populated, while the hydrophobic core of the domain and the remainder of the protein backbone, including helices alpha1 and alpha2, are disordered and devoid of specific interactions. This high-energy state presumably represents the equilibrium analogue of a folding intermediate that is transiently populated in stopped-flow kinetic experiments [Horng, J. C., Tracz, S. M., Lumb, K. J., and Raleigh, D. P. (2002) Biochemistry 44, 627-634].

Automated NMR determination of protein backbone dihedral angles from cross-correlated spin relaxation.

The simultaneous interpretation of a suite of dipole-dipole and dipole-CSA cross-correlation rates involving the backbone nuclei 13Calpha, 1Halpha, 13CO, 15N and 1HN can be used to resolve the ambiguities associated with each individual cross-correlation rate. The method is based on the transformation of experimental cross-correlation rates via calculated values based on standard peptide plane geometry and solid-state 13CO CSA parameters into a dihedral angle probability surface. Triple resonance NMR experiments with improved sensitivity have been devised for the quantification of relaxation interference between 1Halpha(i)-13Calpha(i)/15N(i)-1HN(i) and 1Halpha(i-1)-13Calpha(i-1)/15N(i)-1HN(i) dipole-dipole mechanisms in 15N, 13C-labeled proteins. The approach is illustrated with an application to 13C, 15N-labeled ubiquitin.

Multiple-quantum relaxation dispersion NMR spectroscopy probing millisecond time-scale dynamics in proteins: theory and application.

New relaxation dispersion experiments are presented that probe millisecond time-scale dynamical processes in proteins. The experiments measure the relaxation of (1)H-(15)N multiple-quantum coherence as a function of the rate of application of either (1)H or (15)N refocusing pulses during a constant time relaxation interval. In contrast to the dispersion profiles generated from more conventional (15)N((1)H) single-quantum relaxation experiments that depend on changes in (15)N((1)H) chemical shifts between exchanging states, (1)H-(15)N multiple-quantum dispersions are sensitive to changes in the chemical environments of both (1)H and (15)N spins. The resulting multiple-quantum relaxation dispersion profiles can, therefore, be quite different from those generated by single-quantum experiments, so that an analysis of both single- and multiple-quantum profiles together provides a powerful approach for obtaining robust measures of exchange parameters. This is particularly the case in applications to protonated proteins where other methods for studying exchange involving amide proton spins are negatively influenced by contributions from neighboring protons. The methodology is demonstrated on protonated and perdeuterated samples of a G48M mutant of the Fyn SH3 domain that exchanges between folded and unfolded states in solution.

Probing slow dynamics in high molecular weight proteins by methyl-TROSY NMR spectroscopy: application to a 723-residue enzyme.

A new CPMG-based multiple quantum relaxation dispersion experiment is presented for measuring millisecond dynamic processes at side-chain methyl positions in high molecular weight proteins. The experiment benefits from a methyl-TROSY effect in which cancellation of intramethyl dipole fields occurs, leading to methyl (13)C-(1)H correlation spectra of high sensitivity and resolution (Tugarinov, V.; Hwang, P. M.; Ollerenshaw, J. E.; Kay, L. E. J. Am. Chem. Soc. 2003, 125, 10420-10428). The utility of the methodology is illustrated with an application to a highly deuterated, methyl-protonated sample of malate synthase G, an 82 kDa enzyme consisting of a single polypeptide chain. A comparison of the sensitivity obtained using the present approach relative to existing HSQC-type (13)C single quantum dispersion experiments shows a gain of a factor of 5.4 on average, significantly increasing the range of applications for this methodology.