A non-radioactive, improved PAR-CLIP and small RNA cDNA library preparation protocol.

Crosslinking and immunoprecipitation (CLIP) methods are powerful techniques to interrogate direct protein-RNA interactions and dissect posttranscriptional gene regulatory networks. One widely used CLIP variant is photoactivatable ribonucleoside enhanced CLIP (PAR-CLIP) that involves in vivo labeling of nascent RNAs with the photoreactive nucleosides 4-thiouridine (4SU) or 6-thioguanosine (6SG), which can efficiently crosslink to interacting proteins using UVA and UVB light. Crosslinking of 4SU or 6SG to interacting amino acids changes their base-pairing properties and results in characteristic mutations in cDNA libraries prepared for high-throughput sequencing, which can be computationally exploited to remove abundant background from non-crosslinked sequences and help pinpoint RNA binding protein binding sites at nucleotide resolution on a transcriptome-wide scale. Here we present a streamlined protocol for fluorescence-based PAR-CLIP (fPAR-CLIP) that eliminates the need to use radioactivity. It is based on direct ligation of a fluorescently labeled adapter to the 3'end of crosslinked RNA on immobilized ribonucleoproteins, followed by isolation of the adapter-ligated RNA and efficient conversion into cDNA without the previously needed size fractionation on denaturing polyacrylamide gels. These improvements cut the experimentation by half to 2 days and increases sensitivity by 10-100-fold.

A SARS-CoV-2 mutant from B.1.258 lineage with ∆H69/∆V70 deletion in the Spike protein circulating in Central Europe in the fall 2020.

SARS-CoV-2 mutants carrying the ∆H69/∆V70 deletion in the amino-terminal domain of the Spike protein emerged independently in at least six lineages of the virus (namely, B.1.1.7, B.1.1.298, B.1.160, B.1.177, B.1.258, B.1.375). We analyzed SARS-CoV-2 samples collected from various regions of Slovakia between November and December 2020 that were presumed to contain B.1.1.7 variant due to drop-out of the Spike gene target in an RT-qPCR test caused by this deletion. Sequencing of these samples revealed that although in some cases the samples were indeed confirmed as B.1.1.7, a substantial fraction of samples contained another ∆H69/∆V70 carrying mutant belonging to the lineage B.1.258, which has been circulating in Central Europe since August 2020, long before the import of B.1.1.7. Phylogenetic analysis shows that the early sublineage of B.1.258 acquired the N439K substitution in the receptor-binding domain (RBD) of the Spike protein and, later on, also the deletion ∆H69/∆V70 in the Spike N-terminal domain (NTD). This variant was particularly common in several European countries including the Czech Republic and Slovakia but has been quickly replaced by B.1.1.7 early in 2021.

Identification of the RNA recognition element of the RBPMS family of RNA-binding proteins and their transcriptome-wide mRNA targets.

Recent studies implicated the RNA-binding protein with multiple splicing (RBPMS) family of proteins in oocyte, retinal ganglion cell, heart, and gastrointestinal smooth muscle development. These RNA-binding proteins contain a single RNA recognition motif (RRM), and their targets and molecular function have not yet been identified. We defined transcriptome-wide RNA targets using photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) in HEK293 cells, revealing exonic mature and intronic pre-mRNA binding sites, in agreement with the nuclear and cytoplasmic localization of the proteins. Computational and biochemical approaches defined the RNA recognition element (RRE) as a tandem CAC trinucleotide motif separated by a variable spacer region. Similar to other mRNA-binding proteins, RBPMS family of proteins relocalized to cytoplasmic stress granules under oxidative stress conditions suggestive of a support function for mRNA localization in large and/or multinucleated cells where it is preferentially expressed.

Crystal structure of a DNA containing the planar, phenoxazine-derived bi-functional spectroscopic probe C.

Previously, we developed the deoxycytosine analog Ç (C-spin) as a bi-functional spectroscopic probe for the study of nucleic acid structure and dynamics using electron paramagnetic resonance (EPR) and fluorescence spectroscopy. To understand the effect of Ç on nucleic acid structure, we undertook a detailed crystallographic analysis. A 1.7 Å resolution crystal structure of Ç within a decamer duplex A-form DNA confirmed that Ç forms a non-perturbing base pair with deoxyguanosine, as designed. In the context of double-stranded DNA Ç adopted a planar conformation. In contrast, a crystal structure of the free spin-labeled base ç displayed a ∼ 20° bend at the oxazine linkage. Density function theory calculations revealed that the bent and planar conformations are close in energy and exhibit the same frequency for bending. These results indicate a small degree of flexibility around the oxazine linkage, which may be a consequence of the antiaromaticity of a 16-π electron ring system. Within DNA, the amplitude of the bending motion is restricted, presumably due to base-stacking interactions. This structural analysis shows that the Ç forms a planar, structurally non-perturbing base pair with G indicating it can be used with high confidence in EPR- or fluorescence-based structural and dynamics studies.

Vivid COVID-19 LAMP is an ultrasensitive, quadruplexed test using LNA-modified primers and a zinc ion and 5-Br-PAPS colorimetric detection system.

Sensitive and rapid point-of-care assays have been crucial in the global response to SARS-CoV-2. Loop-mediated isothermal amplification (LAMP) has emerged as an important diagnostic tool given its simplicity and minimal equipment requirements, although limitations exist regarding sensitivity and the methods used to detect reaction products. We describe the development of Vivid COVID-19 LAMP, which leverages a metallochromic detection system utilizing zinc ions and a zinc sensor, 5-Br-PAPS, to circumvent the limitations of classic detection systems dependent on pH indicators or magnesium chelators. We make important strides in improving RT-LAMP sensitivity by establishing principles for using LNA-modified LAMP primers, multiplexing, and conducting extensive optimizations of reaction parameters. To enable point-of-care testing, we introduce a rapid sample inactivation procedure without RNA extraction that is compatible with self-collected, non-invasive gargle samples. Our quadruplexed assay (targeting E, N, ORF1a, and RdRP) reliably detects 1 RNA copy/µl of sample (=8 copies/reaction) from extracted RNA and 2 RNA copies/µl of sample (=16 copies/reaction) directly from gargle samples, making it one of the most sensitive RT-LAMP tests and even comparable to RT-qPCR. Additionally, we demonstrate a self-contained, mobile version of our assay in a variety of high-throughput field testing scenarios on nearly 9,000 crude gargle samples. Vivid COVID-19 LAMP can be an important asset for the endemic phase of COVID-19 as well as preparing for future pandemics.

Multiplexed RNA-FISH-guided Laser Capture Microdissection RNA Sequencing Improves Breast Cancer Molecular Subtyping, Prognostic Classification, and Predicts Response to Antibody Drug Conjugates.

On a retrospective cohort of 1,082 FFPE breast tumors, we demonstrated the analytical validity of a test using multiplexed RNA-FISH-guided laser capture microdissection (LCM) coupled with RNA-sequencing (mFISHseq), which showed 93% accuracy compared to immunohistochemistry. The combination of these technologies makes strides in i) precisely assessing tumor heterogeneity, ii) obtaining pure tumor samples using LCM to ensure accurate biomarker expression and multigene testing, and iii) providing thorough and granular data from whole transcriptome profiling. We also constructed a 293-gene intrinsic subtype classifier that performed equivalent to the research based PAM50 and AIMS classifiers. By combining three molecular classifiers for consensus subtyping, mFISHseq alleviated single sample discordance, provided near perfect concordance with other classifiers (κ > 0.85), and reclassified 30% of samples into different subtypes with prognostic implications. We also use a consensus approach to combine information from 4 multigene prognostic classifiers and clinical risk to characterize high, low, and ultra-low risk patients that relapse early (< 5 years), late (> 10 years), and rarely, respectively. Lastly, to identify potential patient subpopulations that may be responsive to treatments like antibody drug-conjugates (ADC), we curated a list of 92 genes and 110 gene signatures to interrogate their association with molecular subtype and overall survival. Many genes and gene signatures related to ADC processing (e.g., antigen/payload targets, endocytosis, and lysosome activity) were independent predictors of overall survival in multivariate Cox regression models, thus highlighting potential ADC treatment-responsive subgroups. To test this hypothesis, we constructed a unique 19-feature classifier using multivariate logistic regression with elastic net that predicted response to trastuzumab emtansine (T-DM1; AUC = 0.96) better than either ERBB2 mRNA or Her2 IHC alone in the T-DM1 arm of the I-SPY2 trial. This test was deployed in a research-use only format on 26 patients and revealed clinical insights into patient selection for novel therapies like ADCs and immunotherapies and de-escalation of adjuvant chemotherapy.

Conformation and dynamics of nucleotides in bulges and symmetric internal loops in duplex DNA studied by EPR and fluorescence spectroscopies.

The dynamics and conformation of base bulges and internal loops in duplex DNA were studied using the bifunctional spectroscopic probe Ç, which becomes fluorescent (Ç(f)) upon reduction of the nitroxide functional group, along with EPR and fluorescence spectroscopies. A one-base bulge was in a conformational equilibrium between looped-out and stacked states, the former favored at higher temperature and the latter at lower temperature. Stacking of bulge bases was favored in two- and three-base bulges, independent of temperature, resulting in DNA bending as evidenced by increased fluorescence of Ç(f). EPR spectra of Ç-labeled three-, four- and five-base symmetrical interior DNA bulges at 20 °C showed low mobility, indicating that the spin-label was stacked within the loop. The spin-label mobility at 37 °C increased as the loops became larger. A considerable variation in fluorescence between different loops was observed, as well as a temperature-dependence within constructs. Fluorescence unexpectedly increased as the size of the loop decreased at 2 °C. Fluorescence of the smallest loops, where a single T·T mismatch was located between the stem region and the probe, was even larger than for the single strand, indicating a considerable local structural deformation of these loops from regular B-DNA. These results show the value of combining EPR and fluorescence spectroscopy to study non-helical regions of nucleic acids.

Circulating Cell-Free DNA-Based Methylation Pattern in Saliva for Early Diagnosis of Head and Neck Cancer.

Head and neck cancer (HNC) remains one of the leading causes of mortality worldwide due to tumor diagnosis at a late stage, loco-regional aggression, and distant metastases. A standardized diagnostic procedure for HNC is a tissue biopsy that cannot faithfully portray the in-depth tumor dynamics. Therefore, there is an urgent need to develop simple, accurate, and non-invasive methods for cancer detection and follow-up. A saliva-based liquid biopsy allows convenient, non-invasive, and painless collection of high volumes of this biofluid, with the possibility of repetitive sampling, all enabling real-time monitoring of the disease. No approved clinical test for HNC has yet been established. However, epigenetic changes in saliva circulating cell-free DNA (cfDNA) have the potential for a wide range of clinical applications. Therefore, the aim of this review is to present an overview of cfDNA-based methylation patterns in saliva for early detection of HNC, with particular attention to circulating tumor DNA (ctDNA). Due to advancements in isolation and detection technologies, as well as next- and third-generation sequencing, recent data suggest that salivary biomarkers may be successfully applied for early detection of HNC in the future, but large prospective clinical trials are still warranted.

Sequential development of several RT-qPCR tests using LNA nucleotides and dual probe technology to differentiate SARS-CoV-2 from influenza A and B.

Sensitive and accurate RT-qPCR tests are the primary diagnostic tools to identify SARS-CoV-2-infected patients. While many SARS-CoV-2 RT-qPCR tests are available, there are significant differences in test sensitivity, workflow (e.g. hands-on-time), gene targets and other functionalities that users must consider. Several publicly available protocols shared by reference labs and public health authorities provide useful tools for SARS-CoV-2 diagnosis, but many have shortcomings related to sensitivity and laborious workflows. Here, we describe a series of SARS-CoV-2 RT-qPCR tests that are originally based on the protocol targeting regions of the RNA-dependent RNA polymerase (RdRp) and envelope (E) coding genes developed by the Charité Berlin. We redesigned the primers/probes, utilized locked nucleic acid nucleotides, incorporated dual probe technology and conducted extensive optimizations of reaction conditions to enhance the sensitivity and specificity of these tests. By incorporating an RNase P internal control and developing multiplexed assays for distinguishing SARS-CoV-2 and influenza A and B, we streamlined the workflow to provide quicker results and reduced consumable costs. Some of these tests use modified enzymes enabling the formulation of a room temperature-stable master mix and lyophilized positive control, thus increasing the functionality of the test and eliminating cold chain shipping and storage. Moreover, a rapid, RNA extraction-free version enables high sensitivity detection of SARS-CoV-2 in about an hour using minimally invasive, self-collected gargle samples. These RT-qPCR assays can easily be implemented in any diagnostic laboratory and can provide a powerful tool to detect SARS-CoV-2 and the most common seasonal influenzas during the vaccination phase of the pandemic.

Conformational flexibility of DNA.

Pulsed Electron-Electron Double Resonance (PELDOR) on double-stranded DNA (ds-DNA) was used to investigate the conformational flexibility of helical DNA. Stretching, twisting, and bending flexibility of ds-DNA was determined by incorporation of two rigid nitroxide spin labels into a series of 20 base pair (bp) DNA duplexes. Orientation-selective PELDOR experiments performed at both X-band (9 GHz/0.3 T) and G-band (180 GHz/6.4 T) with spin label distances in the range of 2-4 nm allowed us to differentiate between different simple models of DNA dynamics existing in the literature. All of our experimental results are in full agreement with a dynamic model for ds-DNA molecules, where stretching of the molecule leads to a slightly reduced radius of the helix induced by a cooperative twist-stretch coupling.

Folding of the cocaine aptamer studied by EPR and fluorescence spectroscopies using the bifunctional spectroscopic probe Ç.

The cocaine aptamer is a DNA molecule that binds cocaine at the junction of three helices. The bifunctional spectroscopic probe Ç was incorporated independently into three different positions of the aptamer and changes in structure and dynamics upon addition of the cocaine ligand were studied. Nucleoside Ç contains a rigid nitroxide spin label and can be studied directly by electron paramagnetic resonance (EPR) spectroscopy and fluorescence spectroscopy after reduction of the nitroxide to yield the fluoroside Ç(f). Both the EPR and the fluorescence data for aptamer 2 indicate that helix III is formed before cocaine binding. Upon addition of cocaine, increased fluorescence of a fully base-paired Ç(f), placed at the three-way junction in helix III, was observed and is consistent with a helical tilt from a coaxial stack of helices II and III. EPR and fluorescence data clearly show that helix I is formed upon addition of cocaine, concomitant with the formation of the Y-shaped three-way helical junction. The EPR data indicate that nucleotides in helix I are more mobile than nucleotides in regular duplex regions and may reflect increased dynamics due to the short length of helix I.

Surveillance of SARS-CoV-2 lineage B.1.1.7 in Slovakia using a novel, multiplexed RT-qPCR assay.

The emergence of a novel SARS-CoV-2 B.1.1.7 variant sparked global alarm due to increased transmissibility, mortality, and uncertainty about vaccine efficacy, thus accelerating efforts to detect and track the variant. Current approaches to detect B.1.1.7 include sequencing and RT-qPCR tests containing a target assay that fails or results in reduced sensitivity towards the B.1.1.7 variant. Since many countries lack genomic surveillance programs and failed assays detect unrelated variants containing similar mutations as B.1.1.7, we used allele-specific PCR, and judicious placement of LNA-modified nucleotides to develop an RT-qPCR test that accurately and rapidly differentiates B.1.1.7 from other SARS-CoV-2 variants. We validated the test on 106 clinical samples with lineage status confirmed by sequencing and conducted a country-wide surveillance study of B.1.1.7 prevalence in Slovakia. Our multiplexed RT-qPCR test showed 97% clinical sensitivity and retesting 6,886 SARS-CoV-2 positive samples obtained during three campaigns performed within one month, revealed pervasive spread of B.1.1.7 with an average prevalence of 82%. Labs can easily implement this test to rapidly scale B.1.1.7 surveillance efforts and it is particularly useful in countries with high prevalence of variants possessing only the ΔH69/ΔV70 deletion because current strategies using target failure assays incorrectly identify these as putative B.1.1.7 variants.

Rigid spin-labeled nucleoside C: a nonperturbing EPR probe of nucleic acid conformation.

Rigid spin-labeled nucleoside C, an analog of deoxycytidine that base-pairs with deoxyguanosine, was incorporated into DNA oligomers by chemical synthesis. Thermal denaturation experiments and circular dichroism (CD) measurements showed that C has a negligible effect on DNA duplex stability and conformation. Nucleoside C was incorporated into several positions within single-stranded DNA oligomers that can adopt two hairpin conformations of similar energy, each of which contains a four-base loop. The relative mobility of nucleotides in the alternating C/G hairpin loops, 5'-d(GCGC) and 5'-d(CGCG), was determined by electron paramagnetic resonance (EPR) spectroscopy. The most mobile nucleotide in the loop is the second one from the 5'-end, followed by the third, first and fourth nucleotides, consistent with previous NMR studies of DNA hairpin loops of different sequences. The EPR hairpin data were also corroborated by fluorescence spectroscopy using oligomers containing reduced C (C(f)), which is fluorescent. Furthermore, EPR spectra of duplex DNAs that contained C at the end of the helix showed features that indicated dipolar coupling between two spins. These data are consistent with end-to-end duplex stacking in solution, which was only observed when G was paired to C, but not when C was paired with A, C or T.

Identification of single-base mismatches in duplex DNA by EPR spectroscopy.

The spin-labeled nucleoside (T)C, containing 2,2,6,6-tetramethylpiperidine-1-oxyl (TEMPO) conjugated to the exocyclic amino group of C, was used to detect single-base mismatches in duplex DNA for the first time by electron paramagnetic resonance (EPR) spectroscopy. Furthermore, the EPR spectra of the fully base-paired duplex ((T)C.G) and the mismatches ((T)C.A, (T)C.C, and (T)C.T) were significantly different, showing that the probe can identify its base-pairing partner in DNA. At lower pH, the mobility of (T)C.A, (T)C.C, and (T)C.T became higher, consistent with increased protonation of the mismatched pairs. Although the duplexes for each of the three flanking sequences tested gave distinguishable EPR signals, the best discrimination between base pairs was achieved for sequences containing a flanking A.T base pair, in particular 5'-d(G(T)CA) and 5'-d(T(T)CA). This study shows that minor structural variations in nucleic acids can be detected with carefully chosen spin labels in conjunction with EPR spectroscopy.

Conformational equilibria of bulged sites in duplex DNA studied by EPR spectroscopy.

Conformational flexibility in nucleic acids provides a basis for complex structures, binding, and signaling. One-base bulges directly neighboring single-base mismatches in nucleic acids can be present in a minimum of two distinct conformations, complicating the examination of the thermodynamics by calorimetry or UV-monitored melting techniques. To provide additional information about such structures, we demonstrate how electron paramagnetic resonance (EPR) active spin-labeled base analogues, base-specifically incorporated into the DNA, are monitors of the superposition of different bulge-mismatch conformations. EPR spectra provide information about the dynamic environments of the probe. This information is cast in terms of "dynamic signatures" that have an underlying basis in structural variations. By examining the changes in the equilibrium of the different states across a range of temperatures, the enthalpy and entropy of the interconversion among possible conformations can be determined. The DNA constructs with a single bulge neighboring a single-base mismatch ("bulge-mismatches") may be approximately modeled as an equilibrium between two possible conformations. This structural information provides insight into the local composition of the bulge-mismatch sequences. Experiments on the bulge-mismatches show that basepairing across the helix can be understood in terms of purine and pyrimidine interactions, rather than specific bases. Measurements of the enthalpy and entropy of formation for the bulge-mismatches by differential scanning calorimetry and UV-monitored melting confirm that the formation of bulge-mismatches is in fact more complicated than a simple two-state process, consistent with the base-specific spectral data that bulge-mismatches exist in multiple conformations in the premelting temperature region. We find that the calculations with the nearest-neighbor (NN) model for the two likely conformations do not correlate well with the populations of structures and thermodynamic parameters inferred from the base-specific EPR dynamics probe. We report that the base-specific spin probes are able to identify a bistable, temperature dependent, switching between conformations for a particular complex bulged construct.

Relative orientation of rigid nitroxides by PELDOR: beyond distance measurements in nucleic acids.

Show me your angle: Incorporation of the rigid spin label C allows determination of both distance and orientation of two nitroxide spin labels in DNA by PELDOR experiments at common X-band frequencies. The orientational information is obtained by varying the position of the detection pulses over the nitroxide spectrum. Simulation of the set of time traces yields very precise distances and angles.

A nonafluoro nucleoside as a sensitive 19f NMR probe of nucleic acid conformation.

A nucleoside carrying a perfluorinated tert-butyl group ( 4) was prepared by a Sonogashira coupling of 5-iodo-2'-deoxyuridine with 4,4,4-trifluoro-3,3-bis(trifluoromethyl)-1-butyne in nearly quantitative yield and subsequently incorporated into DNA oligomers. Thermal denaturation studies showed that 4 had a negligible effect on duplex stability when compared to thymidine. Transition from single strand to duplex was monitored by (19)F NMR spectroscopy at micromolar concentrations of oligomers, demonstrating the sensitivity of 4 as an NMR reporter nucleoside.

Theory for spin-lattice relaxation of spin probes on weakly deformable DNA.

The weakly bending rod (WBR) model of double-stranded DNA (dsDNA) is adapted to analyze the internal dynamics of dsDNA as observed in electron paramagnetic resonance (EPR) measurements of the spin-lattice relaxation rate, R(1e), for spin probes rigidly attached to nucleic acid-bases. The WBR theory developed in this work models dsDNA base-pairs as diffusing rigid cylindrical discs connected by bending and twisting springs whose elastic force constants are kappa and alpha, respectively. Angular correlation functions for both rotational displacement and velocity are developed in detail so as to compute values for R(1e) due to four relaxation mechanisms: the chemical shift anisotropy (CSA), the electron-nuclear dipolar (END), the spin rotation (SR), and the generalized spin diffusion (GSD) relaxation processes. Measured spin-lattice relaxation rates in dsDNA under 50 bp in length are much faster than those calculated for the same DNAs modeled as rigid rods. The simplest way to account for this difference is by allowing for internal flexibility in models of DNA. Because of this discrepancy, we derive expressions for the spectral densities due to CSA, END, and SR mechanisms directly from a weakly bending rod model for DNA. Special emphasis in this development is given to the SR mechanism because of the lack of such detail in previous treatments. The theory developed in this paper provides a framework for computing relaxation rates from the WBR model to compare with magnetic resonance relaxation data and to ascertain the twisting and bending force constants that characterize DNA.

Single base interrogation by a fluorescent nucleotide: each of the four DNA bases identified by fluorescence spectroscopy.

Nucleoside C, which contains a rigid nitroxide spin label, is effectively reduced in DNA by sodium sulfide to the corresponding amine, yielding a fluorescent probe (Cf) that can report the identity of its base-pairing partner in duplex DNA.

Microscopic rate-constants for substrate binding and acylation in cold-adaptation of trypsin I from Atlantic cod.

Temperature imposes limits on where life can thrive and this is evident in the evolution of the basic structural properties of proteins. Cold-adaptation of enzymes is one example, where the catalytic rate constant (k(cat)) is increased compared with hot-acclimated homologous under identical assay conditions. Trypsin I from Atlantic cod (Gadus morhua) has catalytic efficiency (k(cat)/K(m)) for amide hydrolysis that is 17-fold larger than observed for bovine trypsin. Here, the individual rate-constants for association of substrate (k(1)), dissociation of substrate (k(-1)), and acylation of the enzyme (k(2)) have been determined using benzoyl-Arg-p-nitroanilide or benzyloxycarbonyl-Gly-Pro-Arg-p-nitroanilide as substrates. Rather unexpectedly, by far the largest difference (37-fold increase) was observed in k(1), the rate constant for binding of substrate. The cold-adaptation of the dissociation and catalytic steps were not as prominent (increased by 3.7-fold). The length of substrate did have an effect by increasing the reaction rate by 70-fold, and again, the step most affected was the initial binding-step.