Relaxation Dispersion NMR Spectroscopy Research Articles

Identification of Natural Products Inhibiting SARS-CoV-2 by Targeting Viral Proteases: A Combined in Silico and in Vitro Approach.

In this study, an integrated in silico-in vitro approach was employed to discover natural products (NPs) active against SARS-CoV-2. The two SARS-CoV-2 viral proteases, i.e., main protease (Mpro) and papain-like protease (PLpro), were selected as targets for the in silico study. Virtual hits were obtained by docking more than 140,000 NPs and NP derivatives available in-house and from commercial sources, and 38 virtual hits were experimentally validated in vitro using two enzyme-based assays. Five inhibited the enzyme activity of SARS-CoV-2 Mpro by more than 60% at a concentration of 20 μM, and four of them with high potency (IC50 < 10 μM). These hit compounds were further evaluated for their antiviral activity against SARS-CoV-2 in Calu-3 cells. The results from the cell-based assay revealed three mulberry Diels-Alder-type adducts (MDAAs) from Morus alba with pronounced anti-SARS-CoV-2 activities. Sanggenons C (12), O (13), and G (15) showed IC50 values of 4.6, 8.0, and 7.6 μM and selectivity index values of 5.1, 3.1 and 6.5, respectively. The docking poses of MDAAs in SARS-CoV-2 Mpro proposed a butterfly-shaped binding conformation, which was supported by the results of saturation transfer difference NMR experiments and competitive 1H relaxation dispersion NMR spectroscopy.

Mutate-and-chemical-shift-fingerprint (MCSF) to characterize excited states in RNA using NMR spectroscopy.

It is important to understand the dynamics and higher energy structures of RNA, called excited states, to achieve better understanding of RNA function. R1ρ relaxation dispersion NMR spectroscopy (RD) determines chemical shift differences between the most stable, ground state and the short-lived, low-populated excited states. We describe a procedure for deducing the excited state structure from these chemical shift differences using the mutate-and-chemical-shift-fingerprint (MCSF) method, which requires ~2-6 weeks and moderate understanding of NMR and RNA structure. We recently applied the MCSF methodology to elucidate the excited state of microRNA 34a targeting the SIRT1 mRNA and use this example to demonstrate the analysis. The protocol comprises the following steps: (i) determination of the secondary structure of the excited state from RD chemical shift data, (ii) design of trapped excited state RNA, (iii) validation of the excited state structure by NMR, and (iv) MCSF analysis comparing the chemical shifts of the trapped excited state with the RD-derived chemical shift differences. MCSF enables observation of the short-lived RNA structures, which can be functionally and structurally characterized by entrapment.

Probing transient excited states of the bacterial cell division regulator MinE by relaxation dispersion NMR spectroscopy

Bacterial MinD and MinE form a standing oscillatory wave which positions the cell division inhibitor MinC, that binds MinD, everywhere on the membrane except at the midpoint of the cell, ensuring midcell positioning of the cytokinetic septum. During this process MinE undergoes fold switching as it interacts with different partners. We explore the exchange dynamics between major and excited states of the MinE dimer in 3 forms using 15N relaxation dispersion NMR: the full-length protein (6-stranded β-sheet sandwiched between 4 helices) representing the resting state; a 10-residue N-terminal deletion (Δ10) mimicking the membrane-binding competent state where the N-terminal helix is detached to interact with membrane; and N-terminal deletions of either 30 (Δ30) or 10 residues with an I24N mutation (Δ10/I24N), in which the β1-strands at the dimer interface are extruded and available to bind MinD, leaving behind a 4-stranded β-sheet. Full-length MinE samples 2 "excited" states: The first is similar to a full-length/Δ10 heterodimer; the second, also sampled by Δ10, is either similar to or well along the pathway toward the 4-stranded β-sheet form. Both Δ30 and Δ10/I24N sample 2 excited species: The first may involve destabilization of the β3- and β3'-strands at the dimer interface; changes in the second are more extensive, involving further disruption of secondary structure, possibly representing an ensemble of states on the pathway toward restoration of the resting state. The quantitative information on MinE conformational dynamics involving these excited states is crucial for understanding the oscillation pattern self-organization by MinD-MinE interaction dynamics on the membrane.

Studying sparsely populated conformational states in RNA combining chemical synthesis and solution NMR spectroscopy.

Using chemical synthesis and solution NMR spectroscopy, RNA structural ensembles including a major ground state and minor populated excited states can be studied at atomic resolution. In this work, atom-specific 13C labeled RNA building blocks - a 5-13C-uridine and a 2,8-13C2-adenosine building block - are used to introduce isolated 13C-1H-spin topologies into a target RNA to probe such structural ensembles via NMR spectroscopy. First, the 5-13C-uridine 2'-O-TBDMS-phosphoramidite building block was introduced into a 21 nucleotide (nt) tP5c stem construct of the tP5abc subdomain of the Tetrahymena group I ribozyme. Then, the 2,8-13C2-adenosine 2'-O-TBDMS-phosphoramidite building block was incorporated into a 9 kDa and a 15 kD construct derived from the epsilon (ε) RNA element of the duck Hepatitis B virus. The 2,8-13C2-adenosine resonances of the 9 kDa 28 nt sequence could be mapped to the full-length 53 nt construct. The isolated NMR active nuclei pairs were used to probe for low populated excited states (<10%) via 13C-Carr-Purcell-Meiboom-Gill (CPMG)-relaxation dispersion NMR spectroscopy. The 13C-CPMG relaxation dispersion experiment recapitulated a secondary structure switching event in the P5c hairpin of the group I intron construct previously revealed by 15N relaxation dispersion experiments. In the ε-HBV RNA an unfolding event occurring on the millisecond time scale was found in the upper stem in-line with earlier observations. This unpaired conformational state is presumed to be important for the binding of the epsilon reverse transcriptase (RT) enzyme. Thus, a full description of an RNA's folding landscape helps to obtain a deeper understanding of its function, as these high energy conformational states often represent functionally important intermediates involved in (un)folding or ribozyme catalysis.

Overview of Relaxation Dispersion NMR Spectroscopy to Study Protein Dynamics and Protein-Ligand Interactions.

Proteins and nucleic acids are central to all biological processes. NMR spectroscopy has proven to be excellent for studying the dynamics of these macromolecules over various timescales. Relaxation rates and heteronuclear nuclear Overhauser-effect values can resolve motion on pico- to nanosecond timescales, residual dipolar couplings provide information on submicro- to millisecond timescales, and even slower dynamics over seconds to hours can be resolved by hydrogen-exchange experiments. Relaxation dispersion experiments are especially valuable because they resolve motion on micro- to millisecond timescales, encompassing biomolecular motions associated with ligand binding, enzymatic catalysis, and domain-domain opening. These experiments provide structural, kinetic, and thermodynamic information on "invisible" excited conformational states. Relaxation dispersion can be applied not only to single biomolecules but also to protein-ligand complexes to study the kinetics and thermodynamics of association and dissociation. We review recent developments in relaxation dispersion methodology, outline the R1ρ relaxation dispersion experiment, and discuss application to biomolecular interactions. © 2018 by John Wiley & Sons, Inc.

Measurement and Characterization of Hydrogen-Deuterium Exchange Chemistry Using Relaxation Dispersion NMR Spectroscopy.

One-dimensional heteronuclear relaxation dispersion NMR spectroscopy at 13C natural abundance successfully characterized the dynamics of the hydrogen-deuterium exchange reaction occurring at the Nε position in l-arginine by monitoring Cδ in varying amounts of D2O. A small equilibrium isotope effect was observed and quantified, corresponding to ΔG = -0.14 kcal mol-1. A bimolecular rate constant of kD = 5.1 × 109 s-1 M-1 was determined from the pH*-dependence of kex (where pH* is the direct electrode reading of pH in 10% D2O and kex is the nuclear spin exchange rate constant), consistent with diffusion-controlled kinetics. The measurement of ΔG serves to bridge the millisecond time scale lifetimes of the detectable positively charged arginine species with the nanosecond time scale lifetime of the nonobservable low-populated neutral arginine intermediate species, thus allowing for characterization of the equilibrium lifetimes of the various arginine species in solution as a function of fractional solvent deuterium content. Despite the system being in fast exchange on the chemical shift time scale, the magnitude of the secondary isotope shift due to the exchange reaction at Nε was accurately measured to be 0.12 ppm directly from curve-fitting D2O-dependent dispersion data collected at a single static field strength. These results indicate that relaxation dispersion NMR spectroscopy is a robust and general method for studying base-catalyzed hydrogen-deuterium exchange chemistry at equilibrium.

Effect of amino acid mutations on the conformational dynamics of amyloidogenic immunoglobulin light-chains: A combined NMR and in silico study

The conformational dynamics of a pathogenic κ4 human immunoglobulin light-chain variable domain, SMA, associated with AL amyloidosis, were investigated by 15N relaxation dispersion NMR spectroscopy. Compared to a homologous light-chain, LEN, which differs from SMA at eight positions but is non-amyloidogenic in vivo, we find that multiple residues in SMA clustered around the N-terminus and CDR loops experience considerable conformational exchange broadening caused by millisecond timescale protein motions, consistent with a destabilized dimer interface. To evaluate the contribution of each amino acid substitution to shaping the dynamic conformational landscape of SMA, NMR studies were performed for each SMA-like point mutant of LEN followed by in silico analysis for a subset of these proteins. These studies show that a combination of only three mutations located within or directly adjacent to CDR3 loop at the dimer interface, which remarkably include both destabilizing (Q89H and Y96Q) and stabilizing (T94H) mutations, largely accounts for the differences in conformational flexibility between LEN and SMA. Collectively, our studies indicate that a correct combination of stabilizing and destabilizing mutations is key for immunoglobulin light-chains populating unfolded intermediates that result in amyloid formation, and underscore the complex nature of correlations between light-chain conformational flexibility, thermodynamic stability and amyloidogenicity.

Structural Characterization of the Early Events in the Nucleation-Condensation Mechanism in a Protein Folding Process.

The nucleation-condensation mechanism represents a major paradigm to understand the folding process of many small globular proteins. Although substantial evidence has been acquired for this mechanism, it has remained very challenging to characterize the initial events leading to the formation of a folding nucleus. To achieve this goal, we used a combination of relaxation dispersion NMR spectroscopy and molecular dynamics simulations to determine ensembles of conformations corresponding to the denatured, transition, and native states in the folding of the activation domain of human procarboxypeptidase A2 (ADA2h). We found that the residues making up the folding nucleus tend to interact in the denatured state in a transient manner and not simultaneously, thereby forming incomplete and distorted versions of the folding nucleus. Only when all the contacts between these key residues are eventually formed can the protein reach the transition state and continue folding. Overall, our results elucidate the mechanism of formation of the folding nucleus of a protein and provide insights into how its folding rate can be modified during evolution by mutations that modulate the strength of the interactions between the residues forming the folding nucleus.

Measurement of Ligand-Target Residence Times by 1H Relaxation Dispersion NMR Spectroscopy.

A ligand-observed 1H NMR relaxation experiment is introduced for measuring the binding kinetics of low-molecular-weight compounds to their biomolecular targets. We show that this approach, which does not require any isotope labeling, is applicable to ligand–target systems involving proteins and nucleic acids of variable molecular size. The experiment is particularly useful for the systematic investigation of low affinity molecules with residence times in the micro- to millisecond time regime.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy.

In this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear 13C and 1H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site‐specific 2H and 13C labeling, spin topologies are introduced into DNA and RNA that make 1H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A‐site RNA was previously shown to undergo a two site exchange process in the micro‐ to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV‐1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for 1H RD experiments of nucleic acids.

Excited States of Nucleic Acids Probed by Proton Relaxation Dispersion NMR Spectroscopy

AbstractIn this work an improved stable isotope labeling protocol for nucleic acids is introduced. The novel building blocks eliminate/minimize homonuclear 13C and 1H scalar couplings thus allowing proton relaxation dispersion (RD) experiments to report accurately on the chemical exchange of nucleic acids. Using site‐specific 2H and 13C labeling, spin topologies are introduced into DNA and RNA that make 1H relaxation dispersion experiments applicable in a straightforward manner. The novel RNA/DNA building blocks were successfully incorporated into two nucleic acids. The A‐site RNA was previously shown to undergo a two site exchange process in the micro‐ to millisecond time regime. Using proton relaxation dispersion experiments the exchange parameters determined earlier could be recapitulated, thus validating the proposed approach. We further investigated the dynamics of the cTAR DNA, a DNA transcript that is involved in the viral replication cycle of HIV‐1. Again, an exchange process could be characterized and quantified. This shows the general applicablility of the novel labeling scheme for 1H RD experiments of nucleic acids.

PEGylated DOTA‐AHA‐Based GdIII Chelates: A Relaxometric Study

AbstractThree PEGylated derivatives of 6‐amino‐2‐[4,7,10‐tris(carboxymethyl)‐1,4,7,10‐tetrazacyclododec‐1‐yl]hexanoic acid (DOTA‐AHA) with different molecular weights were prepared and characterized. Their GdIII chelates were studied in aqueous solution through variable‐temperature 1H nuclear magnetic relaxation dispersion (NMRD) and 17O NMR spectroscopy to determine their relaxivities and the parameters that govern them. The relaxivity varied from 5.1 to 6.5 mM–1 s–1 (37 °C and 60 MHz) as the molecular weight of the polyethylene glycol (PEG) chain increased and was slightly higher than that of the parent chelate Gd(DOTA‐AHA), owing to a small contribution of a slow global rotation of the complexes. A variable‐temperature 1H NMR spectroscopy study of several LnIII chelates of DOTA‐A(PEG750)HA allowed the determination of the isomeric M/m ratio (M = square‐antiprismatic isomer and m = twisted square‐antiprismatic isomer; the latter presents a much faster water exchange), which was estimated to be ca. 1:0.2 for the GdIII chelate, very close to that of [Gd(DOTA)]–. This explains why the PEGylated GdIII chelate has a water rate exchange similar to that of [Gd(DOTA)]–. The predominance of the M isomer is a consequence of the bulky PEG moiety, which does not favor the stabilization of the m isomer in sterically crowded systems at the substituent site, in contrast to what happens with less packed asymmetrical DOTA‐type chelates with substitution at one of the four acetate C(α) atoms.

The Energetics of a Three‐State Protein Folding System Probed by High‐Pressure Relaxation Dispersion NMR Spectroscopy

AbstractThe energetic and volumetric properties of a three‐state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure‐dependent 15N‐relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol−1 and 40 μL mol−1 bar−1, respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on‐pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter.

The energetics of a three-state protein folding system probed by high-pressure relaxation dispersion NMR spectroscopy.

The energetic and volumetric properties of a three-state protein folding system, comprising a metastable triple mutant of the Fyn SH3 domain, have been investigated using pressure-dependent (15) N-relaxation dispersion NMR from 1 to 2500 bar. Changes in partial molar volumes (ΔV) and isothermal compressibilities (ΔκT ) between all the states along the folding pathway have been determined to reasonable accuracy. The partial volume and isothermal compressibility of the folded state are 100 mL mol(-1) and 40 μL mol(-1) bar(-1) , respectively, higher than those of the unfolded ensemble. Of particular interest are the findings related to the energetic and volumetric properties of the on-pathway folding intermediate. While the latter is energetically close to the unfolded state, its volumetric properties are similar to those of the folded protein. The compressibility of the intermediate is larger than that of the folded state reflecting the less rigid nature of the former relative to the latter.

Dinuclear DOTA‐Based GdIII Chelates – Revisiting a Straightforward Strategy for Relaxivity Improvement

AbstractThe need for magnetic resonance imaging (MRI) contrast agents with improved relaxivity maintains the development of new GdIII chelates as an intensive and demanding field of research. In this work, we introduce the new dimeric chelators bis{1,4,7,10‐tetraazacyclododecane‐1‐[(6‐amino)hexanoic]‐4,7,10‐triacetic acid}adipate [L2, bis(DOTA‐AHA)adipate] and bis{1,4,7,10‐tetraazacyclododecane‐1‐[(6‐amino)hexanoic]‐4,7,10‐triacetic acid}1,3‐phenyldiacetate [L3, bis(DOTA‐AHA)1,3‐phenyldiacetate], which are based on the bifunctional ligand 1,4,7,10‐tetraazacyclododecane‐1‐[(6‐amino)hexanoic]‐4,7,10‐triacetic acid (L1, DOTA‐AHA). Their GdIII chelates were studied by variable‐temperature 1H nuclear magnetic relaxation dispersion (NMRD) and 17O NMR spectroscopy to measure their relaxivities and the parameters that govern them. The rates of exchange of inner‐sphere water from the monomer GdL1 and from the two dinuclear chelates Gd2L2 and Gd2L3 are very similar (298kex ≈ 6.5 × 106 s–1) and slightly faster than that for [Gd(DOTA)H2O]– (298kex = 4.1 × 106 s–1). All three compounds form weakly bound aggregates with equilibrium constants 298K of 2.9, 15.6, and 14.6 for GdL1, Gd2L2, and Gd2L3, respectively. Even though the aggregates contain only 10 to 15 % of the total amount of GdIII ions, a marked increase in relaxivity between 30 and 100 MHz is observed. Furthermore, the distance between the two GdIII centers in the dinuclear compounds has been determined by double electron–electron resonance (DEER) spectroscopy experiments and by molecular modeling studies, which afforded comparable distances. The linkers between the chelating moieties allow GdIII–GdIII distances of ca. 3.0 nm for the completely stretched linker conformation and less than 1.9 nm for the conformation with the metal centers at a closer distance. These metal‐to‐metal distances by themselves cannot explain the considerably long tumbling times of the chelates in solution. Only a model consistent with some level of aggregation for the dinuclear chelates in aqueous solution could satisfactorily explain our results.

A minor conformation of a lanthanide tag on adenylate kinase characterized by paramagnetic relaxation dispersion NMR spectroscopy

NMR relaxation dispersion techniques provide a powerful method to study protein dynamics by characterizing lowly populated conformations that are in dynamic exchange with the major state. Paramagnetic NMR is a versatile tool for investigating the structures and dynamics of proteins. These two techniques were combined here to measure accurate and precise pseudocontact shifts of a lowly populated conformation. This method delivers valuable long-range structural restraints for higher energy conformations of macromolecules in solution. Another advantage of combining pseudocontact shifts with relaxation dispersion is the increase in the amplitude of dispersion profiles. Lowly populated states are often involved in functional processes, such as enzyme catalysis, signaling, and protein/protein interactions. The presented results also unveil a critical problem with the lanthanide tag used to generate paramagnetic relaxation dispersion effects in proteins, namely that the motions of the tag can interfere severely with the observation of protein dynamics. The two-point attached CLaNP-5 lanthanide tag was linked to adenylate kinase. From the paramagnetic relaxation dispersion only motion of the tag is observed. The data can be described accurately by a two-state model in which the protein-attached tag undergoes a 23° tilting motion on a timescale of milliseconds. The work demonstrates the large potential of paramagnetic relaxation dispersion and the challenge to improve current tags to minimize relaxation dispersion from tag movements.

Ligand‐Detected Relaxation Dispersion NMR Spectroscopy: Dynamics of preQ1–RNA Binding

AbstractAn NMR‐based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand‐detected Carr‐Purcell‐Meiboom‐Gill (CPMG) relaxation dispersion experiments is described. A 15N‐modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand–receptor interaction: the off rate and the population of the small molecule–receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch’s function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand‐detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.

Ligand-detected relaxation dispersion NMR spectroscopy: dynamics of preQ1-RNA binding.

An NMR-based approach to characterizing the binding kinetics of ligand molecules to biomolecules, like RNA or proteins, by ligand-detected Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion experiments is described. A (15)N-modified preQ1 ligand is used to acquire relaxation dispersion experiments in the presence of low amounts of the Fsu class I preQ1 aptamer RNA, and increasing ligand concentrations to probe the RNA small molecule interaction. Our experimental data strongly support the conformational selection mechanism postulated. The approach gives direct access to two parameters of a ligand-receptor interaction: the off rate and the population of the small molecule-receptor complex. A detailed description of the kinetics underlying the ligand binding process is of crucial importance to fully understanding a riboswitch's function and to evaluate potential new antibiotics candidates targeting the noncoding RNA species. Ligand-detected NMR relaxation dispersion experiments represent a valuable diagnostic tool for the characterization of binding mechanisms.

A Similar In Vitro and In Cell Lysate Folding Intermediate for the FF Domain

Understanding the mechanisms by which proteins fold into their three-dimensional structures, including a description of the intermediates that are formed during the folding process, remains a goal of protein science. Most studies are performed under carefully controlled conditions in which the folding reaction is monitored in a buffer solution that is far from the natural milieu of the cell. Here, we have used 13C and 1H relaxation dispersion NMR spectroscopy to study folding of the FF domain in both Escherichia coli and Saccharomyces cerevisiae cellular lysates. We find that a conformationally excited state is populated in both lysates, which is very similar in structure to a folding intermediate observed in previous studies in buffer, with the kinetics and thermodynamics of the interconversion between native and intermediate conformers somewhat changed. The results point to the importance of extending folding studies beyond the test tube yet emphasize that insights can be obtained through careful experiments recorded in controlled buffer solutions.

Uncovering pH-dependent transient states of proteins with buried ionizable residues.

The role of pH in regulating biological activity is ubiquitous, and understanding pH-mediated activity has traditionally relied on analyzing static biomolecular structures of highly populated ground states solved near physiological pH. However, recent advances have shown the increasing importance of transiently populated states, the characterization of which is extremely challenging but made plausible with the development of techniques such as relaxation dispersion NMR spectroscopy. To unlock the pH dependence of these transient states with atomistic-level details, we applied the recently developed explicit solvent constant pH molecular dynamics (CPHMDMSλD) framework to a series of staphylococcal nuclease (SNase) mutants with buried ionizable residues and probed their dynamics in different pH environments. Among our key findings is the existence of open states in all SNase mutants containing “buried” residues with highly shifted pKa’s, where local solvation around the protonation site was observed. The calculated pKa demonstrated good agreement with experimental pKa’s, with a low average unsigned error of 1.3 pKa units and correlation coefficient R2 = 0.78. Sampling both open and closed states in their respective pH range, where they are expected to be dominant, was necessary to reproduce experimental pKa’s, and in the most extreme examples of pKa shifts measured, it can be interpreted that the open-state structures are transient at physiological pH, contributing a small population of 1–2%. This suggests that buried ionizable residues can trigger conformational fluctuations that may be observed as transient-state structures at physiological pH. Furthermore, the coupled relationship of both open and closed states and their role in recapitulating macroscopic experimental observables suggest that structural analysis of buried residues may benefit from looking at structural pairs, as opposed to the conventional approach of looking at a single static ground-state conformation.