Design, quality and validation of the EU-OPENSCREEN fragment library poised to a high-throughput screening collection

The EU-OPENSCREEN (EU-OS) European Research Infrastructure Consortium (ERIC) is a multinational, not-for-profit initiative that integrates high-capacity screening platforms and chemistry groups across Europe to facilitate research in chemical biology and early drug discovery. Over the years, the EU-OS has assembled a high-throughput screening compound collection, the European Chemical Biology Library (ECBL), that contains approximately 100 000 commercially available small molecules and a growing number of thousands of academic compounds crowdsourced through our network of European and non-European chemists. As an extension of the ECBL, here we describe the computational design, quality control and use case screenings of the European Fragment Screening Library (EFSL) composed of 1056 mini and small chemical fragments selected from a substructure analysis of the ECBL. Access to the EFSL is open to researchers from both academia and industry. Using EFSL, eight fragment screening campaigns using different structural and biophysical methods have successfully identified fragment hits in the last two years. As one of the highlighted projects for antibiotics, we describe the screening by Bio-Layer Interferometry (BLI) of the EFSL, the identification of a 35 μM fragment hit targeting the beta-ketoacyl-ACP synthase 2 (FabF), its binding confirmation to the protein by X-ray crystallography (PDB 8PJ0), its subsequent rapid exploration of its surrounding chemical space through hit-picking of ECBL compounds that contain the fragment hit as a core substructure, and the final binding confirmation of two follow-up hits by X-ray crystallography (PDB 8R0I and 8R1V).

NMR characterization and ligand binding site of the stem loop 2 motif (s2m) from the Delta variant of SARS-CoV-2

The stem loop 2 motif (s2m) in SARS-CoV-2 (SCoV-2) is located in the 3'-UTR. Although s2m has been reported to display characteristics of a mobile genomic element that might lead to an evolutionary advantage, its function has remained unknown. The secondary structure of the original SCoV-2 RNA sequence (Wuhan-Hu-1) was determined by NMR in late 2020, delineating the base pairing pattern and revealing substantial differences in secondary structure compared to SARS-CoV-1 (SCoV-1). The existence of a single G29742-A29756 mismatch in the upper stem of s2m leads to its destabilization and impedes a complete NMR analysis. With Delta, a variant of concern has evolved with one mutation compared to the original sequence that replaces G29742 by U29742. We show here that this mutation results in a more defined structure at ambient temperature accompanied by a rise in melting temperature. Consequently, we were able to identify over 90 % of the relevant NMR resonances using a combination of selective RNA labeling and filtered 2D NOESY as well as 4D NMR experiments. We present a comprehensive NMR analysis of the secondary structure, (sub-) nanosecond dynamics and ribose conformation of s2m Delta based on heteronuclear 13C NOE and T1 measurements and ribose carbon chemical shift-derived canonical coordinates. We further show that the G29742U mutation in Delta has no influence on the druggability of s2m compared to the Wuhan-Hu-1 sequence. With the assignment at hand, we identify the flexible regions of s2m as primary site for small molecule binding.

Guanosine-based hydrogel as a supramolecular scaffold for template-assisted macrocyclization

The G-quartet-like supramolecular assembly present in guanosine hydrogel templates macrocyclization between bis-azide and bis-alkyne fragments. The resulting macrocycle enhances viscoelastic properties, and strengthens the hydrogel network. This approach holds potential for the in situ synthesis of drugs and their simultaneous delivery in a stimuli-responsive manner.

Shine-Dalgarno Accessibility Governs Ribosome Binding to the Adenine Riboswitch

Translational riboswitches located in the 5' UTR of the messenger RNA (mRNA) regulate translation through variation of the accessibility of the ribosome binding site (RBS). These are the result of conformational changes in the riboswitch RNA governed by ligand binding. Here, we use a combination of single-molecule colocalization techniques (Single-Molecule Kinetic Analysis of RNA Transient Structure (SiM-KARTS) and Single-Molecule Kinetic Analysis of Ribosome Binding (SiM-KARB)) and microscale thermophoresis (MST) to investigate the adenine-sensing riboswitch in Vibrio vulnificus, focusing on the changes of accessibility between the ligand-free and ligand-bound states. We show that both methods faithfully report on the accessibility of the RBS within the riboswitch and that both methods identify an increase in accessibility upon adenine binding. Expanding on the regulatory context, we show the impact of the ribosomal protein S1 on the unwinding of the RNA secondary structure, thereby favoring ribosome binding even for the apo state. The determined rate constants suggest that binding of the ribosome is faster than the time required to change from the ON state to the OFF state, a prerequisite for efficient regulation decision.

Biophysical Investigation of RNA ⋅ DNA : DNA Triple Helix and RNA : DNA Heteroduplex Formation by the lncRNAs MEG3 and Fendrr

Long non-coding RNAs (lncRNAs) are important regulators of gene expression and can associate with DNA as RNA : DNA heteroduplexes or RNA ⋅ DNA : DNA triple helix structures. Here, we review in vitro biochemical and biophysical experiments including electromobility shift assays (EMSA), circular dichroism (CD) spectroscopy, thermal melting analysis, microscale thermophoresis (MST), single-molecule Förster resonance energy transfer (smFRET) and nuclear magnetic resonance (NMR) spectroscopy to investigate RNA ⋅ DNA : DNA triple helix and RNA : DNA heteroduplex formation. We present the investigations of the antiparallel triplex-forming lncRNA MEG3 targeting the gene TGFB2 and the parallel triplex-forming lncRNA Fendrr with its target gene Emp2. The thermodynamic properties of these oligonucleotides lead to concentration-dependent heterogeneous mixtures, where a DNA duplex, an RNA : DNA heteroduplex and an RNA ⋅ DNA : DNA triplex coexist and their relative populations are modulated in a temperature-dependent manner. The in vitro data provide a reliable readout of triplex structures, as RNA ⋅ DNA : DNA triplexes show distinct features compared to DNA duplexes and RNA : DNA heteroduplexes. Our experimental results can be used to validate computationally predicted triple helix formation between novel disease-relevant lncRNAs and their DNA target genes.

Targeting the Main Protease (Mpro, nsp5) by Growth of Fragment Scaffolds Exploiting Structure-Based Methodologies

The main protease Mpro, nsp5, of SARS-CoV-2 (SCoV2) is one of its most attractive drug targets. Here, we report primary screening data using nuclear magnetic resonance spectroscopy (NMR) of four different libraries and detailed follow-up synthesis on the promising uracil-containing fragment Z604 derived from these libraries. Z604 shows time-dependent binding. Its inhibitory effect is sensitive to reducing conditions. Starting with Z604, we synthesized and characterized 13 compounds designed by fragment growth strategies. Each compound was characterized by NMR and/or activity assays to investigate their interaction with Mpro. These investigations resulted in the four-armed compound 35b that binds directly to Mpro. 35b could be cocrystallized with Mpro revealing its noncovalent binding mode, which fills all four active site subpockets. Herein, we describe the NMR-derived fragment-to-hit pipeline and its application for the development of promising starting points for inhibitors of the main protease of SCoV2.

NMR 1H,19F-based screening of the four stem-looped structure 5_SL1-SL4 located in the 5'-untranslated region of SARS-CoV 2 RNA

Development of new antiviral medication against the beta-coronavirus SARS-CoV-2 (SCoV2) is actively being pursued. Both NMR spectroscopy and crystallography as structural screening technologies have been utilised to screen the viral proteome for binding to fragment libraries. Here, we report on NMR screening of elements of the viral RNA genome with two different ligand libraries using 1H-NMR-screening experiments and 1H and 19F NMR-screening experiments for fluorinated compounds. We screened against the 5'-terminal 119 nucleotides located in the 5'-untranslated region of the RNA genome of SCoV2 and further dissected the four stem-loops into its constituent RNA elements to test specificity of binding of ligands to shorter and longer viral RNA stretches. The first library (DRTL-F library) is enriched in ligands binding to RNA motifs, while the second library (DSI-poised library) represents a fragment library originally designed for protein screening. Conducting screens with two different libraries allows us to compare different NMR screening methodologies, describe NMR screening workflows, validate the two different fragment libraries, and derive initial leads for further downstream medicinal chemistry optimisation.

RNA G-quadruplex folding is a multi-pathway process driven by conformational entropy

The kinetics of folding is crucial for the function of many regulatory RNAs including RNA G-quadruplexes (rG4s). Here, we characterize the folding pathways of a G-quadruplex from the telomeric repeat-containing RNA by combining all-atom molecular dynamics and coarse-grained simulations with circular dichroism experiments. The quadruplex fold is stabilized by cations and thus, the ion atmosphere forming a double layer surrounding the highly charged quadruplex guides the folding process. To capture the ionic double layer in implicit solvent coarse-grained simulations correctly, we develop a matching procedure based on all-atom simulations in explicit water. The procedure yields quantitative agreement between simulations and experiments as judged by the populations of folded and unfolded states at different salt concentrations and temperatures. Subsequently, we show that coarse-grained simulations with a resolution of three interaction sites per nucleotide are well suited to resolve the folding pathways and their intermediate states. The results reveal that the folding progresses from unpaired chain via hairpin, triplex and double-hairpin constellations to the final folded structure. The two- and three-strand intermediates are stabilized by transient Hoogsteen interactions. Each pathway passes through two on-pathway intermediates. We hypothesize that conformational entropy is a hallmark of rG4 folding. Conformational entropy leads to the observed branched multi-pathway folding process for TERRA25. We corroborate this hypothesis by presenting the free energy landscapes and folding pathways of four rG4 systems with varying loop length.

Targeting EPHA2 with Kinase Inhibitors in Colorectal Cancer

The ephrin type-A 2 receptor tyrosine kinase (EPHA2) is involved in the development and progression of various cancer types, including colorectal cancer (CRC). There is also evidence that EPHA2 plays a key role in the development of resistance to the endothelial growth factor receptor (EGFR) monoclonal antibody Cetuximab used clinically in CRC. Despite the promising pharmacological potential of EPHA2, only a handful of specific inhibitors are currently available. In this concept paper, general strategies for EPHA2 inhibition with molecules of low molecular weight (small molecules) are described. Furthermore, available examples of inhibiting EPHA2 in CRC using small molecules are summarized, highlighting the potential of this approach.

PDP1 is a key metabolic gatekeeper and modulator of drug resistance in FLT3-ITD-positive acute myeloid leukemia

High metabolic flexibility is pivotal for the persistence and therapy resistance of acute myeloid leukemia (AML). In 20-30% of AML patients, activating mutations of FLT3, specifically FLT3-ITD, are key therapeutic targets. Here, we investigated the influence of FLT3-ITD on AML metabolism. Nuclear Magnetic Resonance (NMR) profiling showed enhanced reshuffling of pyruvate towards the tricarboxylic acid (TCA) cycle, suggesting an increased activity of the pyruvate dehydrogenase complex (PDC). Consistently, FLT3-ITD-positive cells expressed high levels of PDP1, an activator of the PDC. Combining endogenous tagging of PDP1 with genome-wide CRISPR screens revealed that FLT3-ITD induces PDP1 expression through the RAS signaling axis. PDP1 knockdown resulted in reduced cellular respiration thereby impairing the proliferation of only FLT3-ITD cells. These cells continued to depend on PDP1, even in hypoxic conditions, and unlike FLT3-ITD-negative cells, they exhibited a rapid, PDP1-dependent revival of their respiratory capacity during reoxygenation. Moreover, we show that PDP1 modifies the response to FLT3 inhibition. Upon incubation with the FLT3 tyrosine kinase inhibitor quizartinib (AC220), PDP1 persisted or was upregulated, resulting in a further shift of glucose/pyruvate metabolism towards the TCA cycle. Overexpression of PDP1 enhanced, while PDP1 depletion diminished AC220 resistance in cell lines and peripheral blasts from an AC220-resistant AML patient in vivo. In conclusion, FLT3-ITD assures the expression of PDP1, a pivotal metabolic regulator that enhances oxidative glucose metabolism and drug resistance. Hence, PDP1 emerges as a potentially targetable vulnerability in the management of AML.

High-resolution structure of stem-loop 4 from the 5'-UTR of SARS-CoV-2 solved by solution state NMR

We present the high-resolution structure of stem-loop 4 of the 5'-untranslated region (5_SL4) of the severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2) genome solved by solution state nuclear magnetic resonance spectroscopy. 5_SL4 adopts an extended rod-like structure with a single flexible looped-out nucleotide and two mixed tandem mismatches, each composed of a G•U wobble base pair and a pyrimidine•pyrimidine mismatch, which are incorporated into the stem-loop structure. Both the tandem mismatches and the looped-out residue destabilize the stem-loop structure locally. Their distribution along the 5_SL4 stem-loop suggests a role of these non-canonical elements in retaining functionally important structural plasticity in particular with regard to the accessibility of the start codon of an upstream open reading frame located in the RNA's apical loop. The apical loop-although mostly flexible-harbors residual structural features suggesting an additional role in molecular recognition processes. 5_SL4 is highly conserved among the different variants of SARS-CoV-2 and can be targeted by small molecule ligands, which it binds with intermediate affinity in the vicinity of the non-canonical elements within the stem-loop structure.

Cell-free transcription-translation system: a dual read-out assay to characterize riboswitch function

Cell-free protein synthesis assays have become a valuable tool to understand transcriptional and translational processes. Here, we established a fluorescence-based coupled in vitro transcription-translation assay as a read-out system to simultaneously quantify mRNA and protein levels. We utilized the well-established quantification of the expression of shifted green fluorescent protein (sGFP) as a read-out of protein levels. In addition, we determined mRNA quantities using a fluorogenic Mango-(IV) RNA aptamer that becomes fluorescent upon binding to the fluorophore thiazole orange (TO). We utilized a Mango-(IV) RNA aptamer system comprising four subsequent Mango-(IV) RNA aptamer elements with improved sensitivity by building Mango arrays. The design of this reporter assay resulted in a sensitive read-out with a high signal-to-noise ratio, allowing us to monitor transcription and translation time courses in cell-free assays with continuous monitoring of fluorescence changes as well as snapshots of the reaction. Furthermore, we applied this dual read-out assay to investigate the function of thiamine-sensing riboswitches thiM and thiC from Escherichia coli and the adenine-sensing riboswitch ASW from Vibrio vulnificus and pbuE from Bacillus subtilis, which represent transcriptional and translational on- and off-riboswitches, respectively. This approach enabled a microplate-based application, a valuable addition to the toolbox for high-throughput screening of riboswitch function.

Integrated NMR/Molecular Dynamics Determination of the Ensemble Conformation of a Thermodynamically Stable CUUG RNA Tetraloop

Both experimental and theoretical structure determinations of RNAs have remained challenging due to the intrinsic dynamics of RNAs. We report here an integrated nuclear magnetic resonance/molecular dynamics (NMR/MD) structure determination approach to describe the dynamic structure of the CUUG tetraloop. We show that the tetraloop undergoes substantial dynamics, leading to averaging of the experimental data. These dynamics are particularly linked to the temperature-dependent presence of a hydrogen bond within the tetraloop. Interpreting the NMR data by a single structure represents the low-temperature structure well but fails to capture all conformational states occurring at a higher temperature. We integrate MD simulations, starting from structures of CUUG tetraloops within the Protein Data Bank, with an extensive set of NMR data, and provide a structural ensemble that describes the dynamic nature of the tetraloop and the experimental NMR data well. We thus show that one of the most stable and frequently found RNA tetraloops displays substantial dynamics, warranting such an integrated structural approach.

Fendrr synergizes with Wnt signalling to regulate fibrosis related genes during lung development via its RNA:dsDNA triplex element

Long non-coding RNAs are a very versatile class of molecules that can have important roles in regulating a cells function, including regulating other genes on the transcriptional level. One of these mechanisms is that RNA can directly interact with DNA thereby recruiting additional components such as proteins to these sites via an RNA:dsDNA triplex formation. We genetically deleted the triplex forming sequence (FendrrBox) from the lncRNA Fendrr in mice and found that this FendrrBox is partially required for Fendrr function in vivo. We found that the loss of the triplex forming site in developing lungs causes a dysregulation of gene programs associated with lung fibrosis. A set of these genes contain a triplex site directly at their promoter and are expressed in lung fibroblasts. We biophysically confirmed the formation of an RNA:dsDNA triplex with target promoters in vitro. We found that Fendrr with the Wnt signalling pathway regulates these genes, implicating that Fendrr synergizes with Wnt signalling in lung fibrosis.

NMR resonance assignment of a fibroblast growth factor 8 splicing isoform b

The splicing isoform b of human fibroblast growth factor 8 (FGF8b) is an important regulator of brain embryonic development. Here, we report the almost complete NMR chemical shift assignment of the backbone and aliphatic side chains of FGF8b. Obtained chemical shifts are in good agreement with the previously reported X-ray data, excluding the N-terminal gN helix, which apparently forms only in complex with the receptor. The reported data provide an NMR starting point for the investigation of FGF8b interaction with its receptors and with potential drugs or inhibitors.

Optimization of the Lead Compound NVP-BHG712 as a Colorectal Cancer Inhibitor

The ephrin type-A receptor 2 (EPHA2) kinase belongs to the largest family of receptor tyrosine kinases. There are several indications of an involvement of EPHA2 in the development of infectious diseases and cancer. Despite pharmacological potential, EPHA2 is an under-examined target protein. In this study, we synthesized a series of derivatives of the inhibitor NVP-BHG712 and triazine-based compounds. These compounds were evaluated to determine their potential as kinase inhibitors of EPHA2, including elucidation of their binding mode (X-ray crystallography), affinity (microscale thermophoresis), and selectivity (Kinobeads assay). Eight inhibitors showed affinities in the low-nanomolar regime (KD <10 nM). Testing in up to seven colon cancer cell lines that express EPHA2 reveals that several derivatives feature promising effects for the control of human colon carcinoma. Thus, we have developed a set of powerful tool compounds for fundamental new research on the interplay of EPH receptors in a cellular context.

Chemiluminescence Biosensor for the Determination of Cardiac Troponin I (cTnI)

Cardiac vascular diseases, especially acute myocardial infarction (AMI), are one of the leading causes of death worldwide. Therefore cardio-specific biomarkers such as cardiac troponin I (cTnI) play an essential role in the field of diagnostics. In order to enable rapid and accurate measurement of cTnI with the potential of online measurements, a chemiluminescence-based immunosensor is presented as a proof of concept. A flow cell was designed and combined with a sensitive CMOS camera allowing sensitive optical readout. In addition, a microfluidic setup was established, which achieved selective and quasi-online cTnI determination within ten minutes. The sensor was tested with recombinant cTnI in phosphate buffer and demonstrated cTnI measurements in the concentration range of 2-25 µg/L. With the optimized system, a limit of detection (LoD) of 0.6 µg/L (23 pmol/L) was achieved. Furthermore, the selectivity of the immunosensor was investigated with other recombinant proteins, such as cTnT, and cTnC, at a level of 16 µg/L. No cross-reactivity could be observed. Measurements with diluted blood plasma and serum resulted in an LoD of 60 µg/L (2.4 nmol/L) and 70 µg/L (2.9 nmol/L), respectively.

The COVID19-NMR Consortium: A Public Report on the Impact of this New Global Collaboration

The outbreak of COVID-19 in December 2019 required the formation of international consortia for a coordinated scientific effort to understand and combat the virus. In this Viewpoint Article, we discuss how the NMR community has gathered to investigate the genome and proteome of SARS-CoV-2 and tested them for binding to low-molecular-weight binders. External factors including extended lockdowns due to the global pandemic character of the viral infection triggered the transition from locally focused collaborative research conducted within individual research groups to digital exchange formats for immediate discussion of unpublished results and data analysis, sample sharing, and coordinated research between more than 50 groups from 18 countries simultaneously. We discuss key lessons that might pertain after the end of the pandemic and challenges that we need to address.

Identifying and Overcoming Artifacts in 1 H-Based Saturation Transfer NOE NMR Experiments

Magnetization transfer experiments are versatile nuclear magnetic resonance (NMR) tools providing site-specific information. We have recently discussed how saturation magnetization transfer (SMT) experiments could leverage repeated repolarizations arising from exchanges between labile and water protons to enhance connectivities revealed via the nuclear Overhauser effect (NOE). Repeated experience with SMT has shown that a number of artifacts may arise in these experiments, which may confound the information being sought - particularly when seeking small NOEs among closely spaced resonances. One of these pertains to what we refer to as "spill-over" effects, originating from the use of long saturation pulses leading to changes in the signals of proximate peaks. A second, related but in fact different effect, derives from what we describe as NOE "oversaturation", a phenomenon whereby the use of overtly intense RF fields overwhelms the cross-relaxation signature. The origin and ways to avoid these two effects are described. A final source of potential artifact arises in applications where the labile ¹Hs of interest are bound to ¹⁵N-labeled heteronuclei. SMT's long ¹H saturation times will then be usually implemented while under ¹⁵N decoupling based on cyclic schemes leading to decoupling sidebands. Although these sidebands usually remain invisible in NMR, they may lead to a very efficient saturation of the main resonance when touched by SMT frequencies. All of these phenomena are herein experimentally demonstrated, and solutions to overcome them are proposed.

Modulation of Aβ42 Aggregation Kinetics and Pathway by Low-Molecular-Weight Inhibitors

The aggregation of amyloid-β 42 (Aβ42) is directly related to the pathogenesis of Alzheimer's disease. Here, we have investigated the early stages of the aggregation process, during which most of the cytotoxic species are formed. Aβ42 aggregation kinetics, characterized by the quantification of Aβ42 monomer consumption, were tracked by real-time solution NMR spectroscopy (RT-NMR) allowing the impact that low-molecular-weight (LMW) inhibitors and modulators exert on the aggregation process to be analysed. Distinct differences in the Aβ42 kinetic profiles were apparent and were further investigated kinetically and structurally by using thioflavin T (ThT) and transmission electron microscopy (TEM), respectively. LMW inhibitors were shown to have a differential impact on early-state aggregation. Insight provided here could direct future therapeutic design based on kinetic profiling of the process of fibril formation.

Site-Specific Labeling of RNAs with Modified and 19 F-Labeled Nucleotides by Chemo-Enzymatic Synthesis

More than 170 post-transcriptional modifications of RNAs have currently been identified. Detailed biophysical investigations of these modifications have been limited since large RNAs containing these post-transcriptional modifications are difficult to produce. Further, adequate readout of spectroscopic fingerprints are important, necessitating additional labeling procedures beyond the naturally occurring RNA modifications. Here, we report the chemo-enzymatic synthesis of RNA modifications and several structurally similar fluorine-modified analogs further optimizing a recently developed methodology.^[1] This chemo-enzymatic method allows synthesis of also large RNAs. We were able to incorporate 16 modified nucleotides and 6 ¹⁹ F-labeled nucleotides. To showcase the applicability of such modified large RNAs, we incorporated a ¹⁹ F-labeled cytidine into the aptamer domain of the 2'dG sensing riboswitch (2'dG-sw) from Mesoplasma florum, enabling characterizing RNA fold, ligand binding and kinetics. Thanks to the large chemical shift dispersion of ¹⁹ F, we can detect conformational heterogeneity in the apo state of the riboswitch.

The MECP2-TRD domain interacts with the DNMT3A-ADD domain at the H3-tail binding site

The DNMT3A DNA methyltransferase and MECP2 methylation reader are highly expressed in neurons. Both proteins interact via their DNMT3A-ADD and MECP2-TRD domains, and the MECP2 interaction regulates the activity and subnuclear localization of DNMT3A. Here, we mapped the interface of both domains using peptide SPOT array binding, protein pull-down, equilibrium peptide binding assays, and structural analyses. The region D529-D531 on the surface of the ADD domain was identified as interaction point with the TRD domain. This includes important residues of the histone H3 N-terminal tail binding site to the ADD domain, explaining why TRD and H3 binding to the ADD domain is competitive. On the TRD domain, residues 214-228 containing K219 and K223 were found to be essential for the ADD interaction. This part represents a folded patch within the otherwise largely disordered TRD domain. A crystal structure analysis of ADD revealed that the identified H3/TDR lysine binding pocket is occupied by an arginine residue from a crystallographic neighbor in the ADD apoprotein structure. Finally, we show that mutations in the interface of ADD and TRD domains disrupt the cellular interaction of both proteins in NIH3T3 cells. In summary, our data show that the H3 peptide binding cleft of the ADD domain also mediates the interaction with the MECP2-TRD domain suggesting that this binding site may have a broader role also in the interaction of DNMT3A with other proteins leading to complex regulation options by competitive and PTM specific binding.

Mechanistic Impact of Different Ligand Scaffolds on FXR Modulation Suggests Avenues to Selective Modulators

The bile-acid sensing nuclear farnesoid X receptor (FXR) is an attractive target for the treatment of hepatic and metabolic diseases, but application of this chemotherapeutic concept remains limited due to adverse effects of FXR activation observed in clinical trials. To elucidate the mechanistic basis of FXR activation at the molecular level, we have systematically studied FXR co-regulator interactions and dimerization in response to seven chemically diverse FXR ligands. Different molecular effects on FXR activation mediated by different scaffolds were evident and aligned with characteristic structural changes within the ligand binding domain of FXR. A partial FXR agonist acted mainly through co-repressor displacement from FXR and caused an FXR-regulated gene expression pattern markedly differing from FXR agonist effects. These results suggest selective modulation of FXR dimerization and co-regulator interactions for different ligands, offering a potential avenue for the design of gene- or tissue-selective FXR modulators.

A universal model of RNA.DNA:DNA triplex formation accurately predicts genome-wide RNA-DNA interactions

RNA.DNA:DNA triple helix (triplex) formation is a form of RNA-DNA interaction which regulates gene expression but is difficult to study experimentally in vivo. This makes accurate computational prediction of such interactions highly important in the field of RNA research. Current predictive methods use canonical Hoogsteen base pairing rules, which whilst biophysically valid, may not reflect the plastic nature of cell biology. Here, we present the first optimization approach to learn a probabilistic model describing RNA-DNA interactions directly from motifs derived from triplex sequencing data. We find that there are several stable interaction codes, including Hoogsteen base pairing and novel RNA-DNA base pairings, which agree with in vitro measurements. We implemented these findings in TriplexAligner, a program that uses the determined interaction codes to predict triplex binding. TriplexAligner predicts RNA-DNA interactions identified in all-to-all sequencing data more accurately than all previously published tools in human and mouse and also predicts previously studied triplex interactions with known regulatory functions. We further validated a novel triplex interaction using biophysical experiments. Our work is an important step towards better understanding of triplex formation and allows genome-wide analyses of RNA-DNA interactions.

Comprehensive Fragment Screening of the SARS‐CoV‐2 Proteome Explores Novel Chemical Space for Drug Development

SARS‐CoV‐2 (SCoV2) and its variants of concern pose serious challenges to the public health. The variants increased challenges to vaccines, thus necessitating for development of new intervention strategies including anti‐virals. Within the international Covid19‐NMR consortium, we have identified binders targeting the RNA genome of SCoV2. We established protocols for the production and NMR characterization of more than 80% of all SCoV2 proteins. Here, we performed an NMR screening using a fragment library for binding to 25 SCoV2 proteins and identified hits also against previously unexplored SCoV2 proteins. Computational mapping was used to predict binding sites and identify functional moieties (chemotypes) of the ligands occupying these pockets. Striking consensus was observed between NMR‐detected binding sites of the main protease and the computational procedure. Our investigation provides novel structural and chemical space for structure‐based drug design against the SCoV2 proteome.

Binding Adaptation of GS-441524 Diversifies Macro Domains and Downregulates SARS-CoV-2 de-MARylation Capacity

Viral infection in cells triggers a cascade of molecular defense mechanisms to maintain host-cell homoeostasis. One of these mechanisms is ADP-ribosylation, a fundamental post-translational modification (PTM) characterized by the addition of ADP-ribose (ADPr) on substrates. Poly(ADP-ribose) polymerases (PARPs) are implicated in this process and they perform ADP-ribosylation on host and pathogen proteins. Some viral families contain structural motifs that can reverse this PTM. These motifs known as macro domains (MDs) are evolutionarily conserved protein domains found in all kingdoms of life. They are divided in different classes with the viral belonging to Macro-D-type class because of their properties to recognize and revert the ADP-ribosylation. Viral MDs are potential pharmaceutical targets, capable to counteract host immune response. Sequence and structural homology between viral and human MDs are an impediment for the development of new active compounds against their function. Remdesivir, is a drug administrated in viral infections inhibiting viral replication through RNA-dependent RNA polymerase (RdRp). Herein, GS-441524, the active metabolite of the remdesivir, is tested as a hydrolase inhibitor for several viral MDs and for its binding to human homologs found in PARPs. This study presents biochemical and biophysical studies, which indicate that GS-441524 selectively modifies SARS-CoV-2 MD de-MARylation activity, while it does not interact with hPARP14 MD2 and hPARP15 MD2. The structural investigation of MD•GS-441524 complexes, using solution NMR and X-ray crystallography, discloses the impact of certain amino acids in ADPr binding cavity suggesting that F360 and its adjacent residues tune the selective binding of the inhibitor to SARS-CoV-2 MD.

Synthesis of a Peptidoyl RNA Hairpin via a Combination of Solid-Phase and Template-Directed Chain Assembly

Peptidoyl RNAs are the products of ribosome‐free, single‐nucleotide translation. They contain a peptide in the backbone of the oligoribonucleotide and are interesting from a synthetic and a bioorganic point of view. A synthesis of a stabilized version of peptidoyl RNA, with an amide bond between the C‐terminus of a peptide and a 3′‐amino‐2′,3′‐dideoxynucleoside in the RNA chain was developed. The preferred synthetic route used an N‐Teoc‐protected aminonucleoside support and involved a solution‐phase coupling of the amino‐terminal oligonucleotide to a dipeptido dinucleotide. Exploratory UV‐melting and NMR analysis of the hairpin 5′‐UUGGCGAAAGCdC‐LeuLeu‐AA‐3′ indicated that the peptide‐linked RNA segments do not fold in a cooperative fashion. The synthetic access to doubly RNA‐linked peptides on a scale sufficient for structural biology opens the door to the exploration of their structural and biochemical properties.

The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

Oxidation of the Mycobacterium tuberculosis key virulence factor Protein Tyrosine Phosphatase A (MptpA) reduces its phosphatase activity

The Mycobacterium tuberculosis tyrosine-specific phosphatase MptpA and its cognate kinase PtkA are prospective targets for anti- tuberculosis drugs as they interact with the host defense response within the macrophages. Although both are structurally well characterized, the functional mechanism regulating their activity remains poorly understood. Here, we investigate the effect of post-translational oxidation in regulating the function of MptpA. Treatment of MptpA with H₂ O₂ /NaHCO₃ , mimicking cellular oxidative stress conditions, leads to oxidation of the catalytic cysteine (C11) and to a conformational rearrangement of the phosphorylation loop (D-loop) by repositioning the conserved tyrosine 128 (Y128) and generating a temporarily inactive pre-closed state of the phosphatase. Thus, the catalytic cysteine in the P-loop acts as a redox switch and regulates the phosphatase activity of MptpA.

1H, 13C and 15N chemical shift assignment of the stem-loops 5b + c from the 5'-UTR of SARS-CoV-2

The ongoing pandemic of the respiratory disease COVID-19 is caused by the SARS-CoV-2 (SCoV2) virus. SCoV2 is a member of the Betacoronavirus genus. The 30 kb positive sense, single stranded RNA genome of SCoV2 features 5'- and 3'-genomic ends that are highly conserved among Betacoronaviruses. These genomic ends contain structured cis-acting RNA elements, which are involved in the regulation of viral replication and translation. Structural information about these potential antiviral drug targets supports the development of novel classes of therapeutics against COVID-19. The highly conserved branched stem-loop 5 (SL5) found within the 5'-untranslated region (5'-UTR) consists of a basal stem and three stem-loops, namely SL5a, SL5b and SL5c. Both, SL5a and SL5b feature a 5'-UUUCGU-3' hexaloop that is also found among Alphacoronaviruses. Here, we report the extensive ¹H, ¹³C and ¹⁵N resonance assignment of the 37 nucleotides (nts) long sequence spanning SL5b and SL5c (SL5b + c), as basis for further in-depth structural studies by solution NMR spectroscopy.

The Extended Hadamard Transform: Sensitivity-Enhanced NMR Experiments Among Labile and Non-Labile 1 Hs of SARS-CoV-2-derived RNAs

Hadamard encoded saturation transfer can significantly improve the efficiency of NOE-based NMR correlations from labile protons in proteins, glycans and RNAs, increasing the sensitivity of cross-peaks by an order of magnitude and shortening experimental times by ≥100-fold. These schemes, however, fail when tackling correlations within a pool of labile protons - for instance imino-imino correlations in RNAs or amide-amide correlations in proteins. Here we analyze the origin of the artifacts appearing in these experiments and propose a way to obtain artifact-free correlations both within the labile pool as well as between labile and non-labile 1 Hs, while still enjoying the gains arising from Hadamard encoding and solvent repolarizations. The principles required for implementing what we define as the extended Hadamard scheme are derived, and its clean, artifact-free, sensitivity-enhancing performance is demonstrated on RNA fragments derived from the SARS-CoV-2 genome. Sensitivity gains per unit time approaching an order of magnitude are then achieved in both imino-imino and imino-amino/aromatic protons 2D correlations; similar artifact-free sensitivity gains can be observed when carrying out extended Hadamard encodings of 3D NOESY/HSQC-type experiments. The resulting spectra reveal significantly more correlations than their conventionally acquired counterparts, which can support the spectral assignment and secondary structure determination of structured RNA elements.

Structural analysis of temperature-dependent alternative splicing of HsfA2 pre-mRNA from tomato plants

Temperature-dependent alternative splicing was recently demonstrated for intron 2 of the gene coding for heat shock factor HsfA2 of the tomato plant Solanum lycopersicum, but the molecular mechanism regulating the abundance of such temperature-dependent splice variants is still unknown. We report here on regulatory pre-mRNA structures that could function as regulators by controlling the use of splice sites in a temperature-dependent manner. We investigate pre-mRNA structures at the splice sites of intron 2 of the gene coding for HsfA2 from S. lycopersicum using NMR- and CD-spectroscopy as well as in-line probing. The pre-mRNA undergoes conformational changes between two different secondary structures at the 3ʹ splice site of the intron in a temperature-dependent manner. Previously, it was shown that three single nucleotide polymorphisms (SNPs) in intron 2 of the HsfA2 pre-mRNA affect the splicing efficiency of its pre-mRNA and are linked to the thermotolerance in different tomato species. By comparing pre-mRNA fragments of the tomato species S. lycopersicum and S. peruvianum, we show that these SNPs result in substantial structural differences between the pre-mRNAs of the two species.

Mapping the conformational landscape of the neutral network of RNA sequences that connect two functional distinctly different ribozymes

During evolution of an RNA world, the development of enzymatic function was essential. Such enzymatic function was linked to RNA sequences capable of adopting specific RNA folds that possess catalytic pockets to promote catalysis. Within this primordial RNA world, initially evolved self-replicating ribozymes presumably mutated to ribozymes with new functions. In 2000, Schultes and Bartel investigated such conversion from one ribozyme to a new ribozyme with distinctly different catalytic functions. Within a neutral network that linked these two prototype ribozymes, a single RNA chain could be identified that exhibited both enzymatic functions. This system serves as a paradigm for an evolutionary system that allows neutral drifts by stepwise mutation from one ribozyme into a different ribozyme without loss of intermittent function. Here, we investigated this complex functional diversification by analyzing several RNA sequences within this neutral network between two ribozymes with class III ligase activity and with self-cleavage reactivity. We utilized rapid RNA sample preparation for NMR spectroscopic studies together with SHAPE analysis and in-line probing to characterize secondary structure changes within the neutral network. Our investigations allowed delineation of the 2 nd structure space and by comparison with the previously determined catalytic function allowed correlation of the structure-function relation of ribozyme function in this neutral network.

The Transcriptional Repressor Orphan Nuclear Receptor TLX Is Responsive to Xanthines

The orphan nuclear receptor tailless homologue (TLX) is expressed almost exclusively in neural stem cells acting as an essential factor for their survival and is hence considered as a promising drug target in neurodegeneration. However, few studies have characterized the roles of TLX due to the lack of ligands and limited functional understanding. Here, we identify xanthines including caffeine and istradefylline as TLX modulators that counteract the receptor's intrinsic repressor activity. Mutagenesis of residues lining a cavity within the TLX ligand binding domain altered the activity of these ligands, suggesting direct interactions with helix 5. Using xanthines as tool compounds, we observed a ligand-sensitive recruitment of the co-repressor silencing mediator for retinoid or thyroid-hormone receptors, TLX homodimerization, and heterodimerization with the retinoid X receptor. These protein-protein interactions evolve as factors that modulate the TLX function and suggest an unprecedented role of TLX in directly repressing other nuclear receptors.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS)

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine‐II riboswitch binds the small cation, guanidinium, and forms a kissing loop‐loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg²⁺ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small‐angle X‐ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg²⁺ and ligand binding.

Solid-Phase-Supported Chemoenzymatic Synthesis of a Light-Activatable tRNA Derivative

Herein, we present a multi‐cycle chemoenzymatic synthesis of modified RNA with simplified solid‐phase handling to overcome size limitations of RNA synthesis. It combines the advantages of classical chemical solid‐phase synthesis and enzymatic synthesis using magnetic streptavidin beads and biotinylated RNA. Successful introduction of light‐controllable RNA nucleotides into the tRNA^Met sequence was confirmed by gel electrophoresis and mass spectrometry. The methods tolerate modifications in the RNA phosphodiester backbone and allow introductions of photocaged and photoswitchable nucleotides as well as photocleavable strand breaks and fluorophores.

Repeating Aspartic Acid Residues Prefer Turn-like Conformations in the Unfolded State: Implications for Early Protein Folding

Protein folding can be described as a motion of the polypeptide chain in a potential energy funnel, where the conformational manifold is narrowed as the chain traverses from a completely unfolded state until it reaches the folded (native) state. The initial folding stages set the tone for this process by substantially narrowing the manifold of accessible conformations. In an ideally unfolded state with no long-range stabilizing forces, local conformations (i.e., residual structures) are likely to drive the folding process. While most amino acid residues tend to predominantly adopt extended structures in unfolded proteins and peptides, aspartic acid exhibits a relatively high intrinsic preference for turn-forming conformations. Regions in an unfolded polypeptide or protein that are rich in aspartic acid residues may therefore be crucial sites for protein folding steps. By combining NMR and vibrational spectroscopies, we observed that the conformational sampling of multiple sequentially neighbored aspartic acid residues in the model peptides GDDG and GDDDG even show an on average higher propensity for turn-forming structures than the intrinsic reference system D in GDG, which suggests that nearest neighbor interactions between adjacent aspartic acid residues stabilize local turn-forming structures. In the presence of the unlike neighbor phenylalanine, nearest neighbor interactions are of a totally different nature in that it they decrease the turn-forming propensities and mutually increase the sampling of polyproline II (pPII) conformations. We hypothesize the structural role of aspartic residues in intrinsically disordered proteins in general, and particularly in small linear motifs, that are very much determined by their respective neighbors.

Folding dynamics of polymorphic G-quadruplex structures

G-quadruplexes (G4), found in numerous places within the human genome, are involved in essential processes of cell regulation. Chromosomal DNA G4s are involved for example, in replication and transcription as first steps of gene expression. Hence, they influence a plethora of downstream processes. G4s possess an intricate structure that differs from canonical B-form DNA. Identical DNA G4 sequences can adopt multiple long-lived conformations, a phenomenon known as G4 polymorphism. A detailed understanding of the molecular mechanisms that drive G4 folding is essential to understand their ambivalent regulatory roles. Disentangling the inherent dynamic and polymorphic nature of G4 structures thus is key to unravel their biological functions and make them amenable as molecular targets in novel therapeutic approaches. We here review recent experimental approaches to monitor G4 folding and discuss structural aspects for possible folding pathways. Substantial progress in the understanding of G4 folding within the recent years now allows drawing comprehensive models of the complex folding energy landscape of G4s that we herein evaluate based on computational and experimental evidence.

1H, 13C and 15N assignment of stem-loop SL1 from the 5'-UTR of SARS-CoV-2

The stem-loop (SL1) is the 5'-terminal structural element within the single-stranded SARS-CoV-2 RNA genome. It is formed by nucleotides 7-33 and consists of two short helical segments interrupted by an asymmetric internal loop. This architecture is conserved among Betacoronaviruses. SL1 is present in genomic SARS-CoV-2 RNA as well as in all subgenomic mRNA species produced by the virus during replication, thus representing a ubiquitous cis-regulatory RNA with potential functions at all stages of the viral life cycle. We present here the ¹H, ¹³C and ¹⁵N chemical shift assignment of the 29 nucleotides-RNA construct 5_SL1, which denotes the native 27mer SL1 stabilized by an additional terminal G-C base-pair.

Backbone chemical shift spectral assignments of SARS coronavirus-2 non-structural protein nsp9

As part of an International consortium aiming at the characterization by NMR of the proteins of the SARS-CoV-2 virus, we have obtained the virtually complete assignment of the backbone atoms of the non-structural protein nsp9. This small (12 kDa) protein is encoded by ORF1a, binds to RNA and seems to be essential for viral RNA synthesis. The crystal structures of the SARS-CoV-2 protein and other homologues suggest that the protein is dimeric as also confirmed by analytical ultracentrifugation and dynamic light scattering. Our data constitute the prerequisite for further NMR-based characterization, and provide the starting point for the identification of small molecule lead compounds that could interfere with RNA binding and prevent viral replication.

1H, 13C, and 15N backbone chemical-shift assignments of SARS-CoV-2 non-structural protein 1 (leader protein)

The current COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has become a worldwide health crisis, necessitating coordinated scientific research and urgent identification of new drug targets for treatment of COVID-19 lung disease. The covid19-nmr consortium seeks to support drug development by providing publicly accessible NMR data on the viral RNA elements and proteins. The SARS-CoV-2 genome comprises a single RNA of about 30 kb in length, in which 14 open reading frames (ORFs) have been annotated, and encodes approximately 30 proteins. The first two-thirds of the SARS-CoV-2 genome is made up of two large overlapping open-reading-frames (ORF1a and ORF1b) encoding a replicase polyprotein, which is subsequently cleaved to yield 16 so-called non-structural proteins. The non-structural protein 1 (Nsp1), which is considered to be a major virulence factor, suppresses host immune functions by associating with host ribosomal complexes at the very end of its C-terminus. Furthermore, Nsp1 facilitates initiation of viral RNA translation via an interaction of its N-terminal domain with the 5' untranslated region (UTR) of the viral RNA. Here, we report the near-complete backbone chemical-shift assignments of full-length SARS-CoV-2 Nsp1 (19.8 kDa), which reveal the domain organization, secondary structure and backbone dynamics of Nsp1, and which will be of value to further NMR-based investigations of both the biochemical and physiological functions of Nsp1.

1H, 13C, 15N and 31P chemical shift assignment for stem-loop 4 from the 5'-UTR of SARS-CoV-2

The SARS-CoV-2 virus is the cause of the respiratory disease COVID-19. As of today, therapeutic interventions in severe COVID-19 cases are still not available as no effective therapeutics have been developed so far. Despite the ongoing development of a number of effective vaccines, therapeutics to fight the disease once it has been contracted will still be required. Promising targets for the development of antiviral agents against SARS-CoV-2 can be found in the viral RNA genome. The 5'- and 3'-genomic ends of the 30 kb SCoV-2 genome are highly conserved among Betacoronaviruses and contain structured RNA elements involved in the translation and replication of the viral genome. The 40 nucleotides (nt) long highly conserved stem-loop 4 (5_SL4) is located within the 5'-untranslated region (5'-UTR) important for viral replication. 5_SL4 features an extended stem structure disrupted by several pyrimidine mismatches and is capped by a pentaloop. Here, we report extensive ¹H, ¹³C, ¹⁵N and ³¹P resonance assignments of 5_SL4 as the basis for in-depth structural and ligand screening studies by solution NMR spectroscopy.