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

Structural Basis of the Immunological Cross-Reactivity between Kiwi and Birch Pollen

Allergies related to kiwi consumption have become a growing health concern, with their prevalence on the rise. Many of these allergic reactions are attributed to cross-reactivity, particularly with the major allergen found in birch pollen. This cross-reactivity is associated with proteins belonging to the pathogenesis-related class 10 (PR-10) protein family. In our study, we determined the three-dimensional structures of the two PR-10 proteins in gold and green kiwi fruits, Act c 8 and Act d 8, using nuclear magnetic resonance (NMR) spectroscopy. The structures of both kiwi proteins closely resemble the major birch pollen allergen, Bet v 1, providing a molecular explanation for the observed immunological cross-reactivity between kiwi and birch pollen. Compared to Act d 11, however, a kiwi allergen that shares the same architecture as PR-10 proteins, structural differences are apparent. Moreover, despite both Act c 8 and Act d 8 containing multiple cysteine residues, no disulfide bridges are present within their structures. Instead, all the cysteines are accessible on the protein's surface and exposed to the surrounding solvent, where they are available for reactions with components of the natural food matrix. This structural characteristic sets Act c 8 and Act d 8 apart from other kiwi proteins with a high cysteine content. Furthermore, we demonstrate that pyrogallol, the most abundant phenolic compound found in kiwi, binds into the internal cavities of these two proteins, albeit with low affinity. Our research offers a foundation for further studies aimed at understanding allergic reactions associated with this fruit and exploring how interactions with the natural food matrix might be employed to enhance food safety.

Mapping N- to C-terminal allosteric coupling through disruption of a putative CD74 activation site in D-dopachrome tautomerase

The macrophage migration inhibitory factor (MIF) protein family consists of MIF and D-dopachrome tautomerase (also known as MIF-2). These homologs share 34% sequence identity while maintaining nearly indistinguishable tertiary and quaternary structure, which is likely a major contributor to their overlapping functions, including the binding and activation of the cluster of differentiation 74 (CD74) receptor to mediate inflammation. Previously, we investigated a novel allosteric site, Tyr99, that modulated N-terminal catalytic activity in MIF through a "pathway" of dynamically coupled residues. In a comparative study, we revealed an analogous allosteric pathway in MIF-2 despite its unique primary sequence. Disruptions of the MIF and MIF-2 N termini also diminished CD74 activation at the C terminus, though the receptor activation site is not fully defined in MIF-2. In this study, we use site-directed mutagenesis, NMR spectroscopy, molecular simulations, in vitro and in vivo biochemistry to explore the putative CD74 activation region of MIF-2 based on homology to MIF. We also confirm its reciprocal structural coupling to the MIF-2 allosteric site and N-terminal enzymatic site. Thus, we provide further insight into the CD74 activation site of MIF-2 and its allosteric coupling for immunoregulation.

J, Petermann O, Scapozza L, Orts J, Kirchmair J, Rabenau HF, Rollinger JM. 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.

Residence Times from Molecular Dynamics Simulations

In this work, efficient methods are proposed for the calculation, from molecular dynamics trajectories, of residence times (RTs) and related quantities. One of these was designed to obtain RT distributions, from which mean residence times (MRTs), residual time distributions, and mean residual times can be computed. This method does not require the assumptions and approximations made by the most commonly used methods. Its link to the most popular method in the literature is identified. It is shown how the much faster new method can be used as a replacement for the latter and the advantages of doing so. Also, a prescription for estimating the uncertainty in the MRTs obtained though the proposed method is provided. Another even faster method for the calculation of the MRTs, their uncertainties, and the mean residual times is also proposed. It yields exactly the same results as the first one but does not allow to obtain the mentioned distributions. Another very popular method, based on autocorrelation functions, for computing MRTs is analyzed in terms of its assumptions and approximations. An alternative, also based on autocorrelation functions, which is faster and requires fewer assumptions is presented. A prescription for the calculation of the uncertainty of the MRTs obtained with the latter method is also provided. In the literature, there are a few methods to discard short transient escapes. Here, an algorithm is suggested that is much faster than the most used one and allows a more integral treatment of the process. Also, it is more widely applicable because it is a preprocessing step that can be used in conjunction with any of the proposed methods mentioned above. The main disadvantage of these two approaches to discard brief escapes is that the maximum duration allowed for an escape to be considered transient appears as a parameter to be determined for the particular system under study. As an alternative, a parameter-free method of a similar character is also proposed to estimate the mean residence time of particles that reached a constant probability of leaving the region of interest.

The Structural Flexibility of PR-10 Food Allergens

PR-10 proteins constitute a major cause of food allergic reactions. Birch-pollen-related food allergies are triggered by the immunologic cross-reactivity of IgE antibodies with structurally homologous PR-10 proteins that are present in birch pollen and various food sources. While the three-dimensional structures of PR-10 food allergens have been characterized in detail, only a few experimental studies have addressed the structural flexibility of these proteins. In this study, we analyze the millisecond-timescale structural flexibility of thirteen PR-10 proteins from prevalent plant food sources by NMR relaxation-dispersion spectroscopy, in a comparative manner. We show that all the allergens in this study have inherently flexible protein backbones in solution, yet the extent of the structural flexibility appears to be strikingly protein-specific (but not food-source-specific). Above-average flexibility is present in the two short helices, α1 and α2, which form a V-shaped support for the long C-terminal helix α3, and shape the internal ligand-binding cavity, which is characteristic for PR-10 proteins. An in-depth analysis of the NMR relaxation-dispersion data for the PR-10 allergen from peanut reveals the presence of at least two subglobal conformational transitions on the millisecond timescale, which may be related to the release of bound low-molecular-weight ligands from the internal cavity.

Application of Relaxation Dispersion of Hyperpolarized 13 C Spins to Protein-Ligand Binding

Nuclear spin relaxation dispersion parameters are proposed as indicators of the binding mode of a ligand to a protein. Hyperpolarization by dissolution dynamic nuclear polarization (D-DNP) provided a ¹³ C signal enhancement between 3000-6000 for the ligand 4-(trifluoromethyl) benzene-1-carboximidamide binding to trypsin. The measurement of ¹³ C R₂ relaxation dispersion was enabled without isotope enrichment, using a series of single-scan Carr-Purcell-Meiboom-Gill experiments with variable refocusing delays. The magnitude in dispersion for the spins of the ligand is correlated to the position with respect to the salt bridge between protein and the amidine group of the ligand, indicating the ligand binding orientation. Hyperpolarized relaxation dispersion is an alternative to chemical shift or NOE measurements for determining ligand binding modes.

NMR Reporter Assays for the Quantification of Weak-Affinity Receptor-Ligand Interactions

Biophysical methods are widely employed in academia and the pharmaceutical industry to detect and quantify weak molecular interactions. Such methods find broad application in fragment-based drug discovery (FBDD). In an FBDD campaign, a suitable affinity determination method is key to advancing a project beyond the initial screening phase. Protein-observed (PO) nuclear magnetic resonance (NMR) finds widespread use due to its ability to sensitively detect very weak interactions at residue-level resolution. When there are issues precluding the use of PO-NMR, ligand-observed (LO) NMR reporter assays can be a useful alternative. Such assays can measure affinities in a similar range to PO-NMR while offering some distinct advantages, especially with regard to protein consumption and compound throughput. In this paper, we take a closer look at setting up such assays for routine use, with the aim of getting high-quality, accurate data and good throughput. We assess some of the key characteristics of these assays in the mathematical framework established for fluorescence polarization assays with which the readers may be more familiar. We also provide guidance on setting up such assays and compare their performance with other affinity determination methods that are commonly used in drug discovery.

Structure and Zeatin Binding of the Peach Allergen Pru p 1

Peach (Prunus persica) is among the fruits most frequently reported to cause food allergies. Allergic reactions commonly result from previous sensitization to the birch pollen allergen Bet v 1, followed by immunological cross-reactivity of IgE antibodies to structurally related proteins in peach. In this study, we present the three-dimensional NMR solution structure of the cross-reactive peach allergen Pru p 1 (isoform Pru p 1.0101). This 17.5 kDa protein adopts the canonical Bet v 1 fold, composed of a seven-stranded β-sheet and three α-helices enclosing an internal cavity. In Pru p 1, the inner surface of the cavity contains an array of hydroxyl-bearing amino acids surrounded by a hydrophobic patch, constituting a docking site for amphiphilic molecules. NMR-guided docking of the cytokinin molecule zeatin to the internal cavity of Pru p 1 provides a structure-based rationale for the effect that zeatin binding has on the protein's RNase activity.

Relaxation Dispersion NMR to Reveal Fast Dynamics in Brønsted Acid Catalysis: Influence of Sterics and H-Bond Strength on Conformations and Substrate Hopping

NMR provides both structural and dynamic information, which is key to connecting intermediates and to understanding reaction pathways. However, fast exchanging catalytic intermediates are often inaccessible by conventional NMR due its limited time resolution. Here, we show the combined application of the 1H off-resonance R1ρ NMR method and low temperature (185-175 K) to resolve intermediates exchanging on a μs time scale (ns at room temperature). The potential of the approach is demonstrated on chiral phosphoric acid (CPA) catalysts in their complexes with imines. The otherwise inaccessible exchange kinetics of the E-I ⇌ E-II imine conformations and thermodynamic E-I:E-II imine ratios inside the catalyst pocket are experimentally determined and corroborated by calculations. The E-I ⇌ E-II exchange rate constants (kex185 K) for different catalyst-substrate binary complexes varied between 2500 and 19 000 s-1 (τex = 500-50 μs). Theoretical analysis of these exchange rate constants revealed the involvement of an intermediary tilted conformation E-III, which structurally resembles the hydride transfer transition state. The main E-I and E-II exchange pathway is a hydrogen bond strength dependent tilting-switching-tilting mechanism via a bifurcated hydrogen bond as a transition state. The reduction in the sterics of the catalyst showed an accelerated switching process by at least an order of magnitude and enabled an additional rotational pathway. Hence, the exchange process is mainly a function of the intrinsic properties of the 3,3'-substituents of the catalyst. Overall, we believe that the present study opens a new dimension in catalysis via experimental access to structures, populations, and kinetics of catalyst-substrate complexes on the μs time scale by the 1H off-resonance R1ρ method.

Ligand-observed NMR techniques to probe RNA-small molecule interactions

A growing understanding of the structure and function of RNA has revealed it as a key regulator of gene expression and disease. A multitude of noncoding functions apart from the central roles of RNA in coding for and facilitating protein biogenesis has stimulated research into RNA as a pharmacological target. Despite many exciting advances, RNA remains an understudied target for small molecules, and techniques to investigate RNA-binding molecules are still emerging. A key stumbling block in this area has been validation of RNA-small molecule interactions. Our laboratory has recently used multiple ligand-observed NMR techniques in this regard, including CPMG and WaterLOGSY. This work describes methods to use these techniques in the context of studying RNA-ligand interactions.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

Ligand Proton Pseudocontact Shifts Determined from Paramagnetic Relaxation Dispersion in the Limit of NMR Intermediate Exchange

Delineation of protein-ligand interaction modes is key for rational drug discovery. The availability of complex crystal structures is often limited by the aqueous solubility of the compounds, while lead-like compounds with micromolar affinities normally fall into the NMR intermediate exchange regime, in which severe line broadening to beyond the detection of interfacial resonances limits NMR applications. Here, we developed a new method to retrieve low-populated bound-state ¹H pseudocontact shifts (PCSs) using paramagnetic relaxation dispersion (RD). We evaluated using a ¹H PCS-RD approach in a BRM bromodomain lead-like inhibitor to filter molecular docking poses using multiple intermolecular structural restraints. Considering the universal presence of proton atoms in druglike compounds, our work will have wide application in structure-guided drug discovery even under an extreme condition of NMR intermediate exchange and low aqueous solubility of ligands.

Efficient Detection of Structure and Dynamics in Unlabeled RNAs: The SELOPE Approach

The knowledge of structure and dynamics is crucial to explain the function of RNAs. While nuclear magnetic resonance (NMR) is well suited to probe these for complex biomolecules, it requires expensive, isotopically labeled samples, and long measurement times. Here we present SELOPE, a new robust, proton-only NMR method that allows us to obtain site-specific overview of structure and dynamics in an entire RNA molecule using an unlabeled sample. SELOPE simplifies assignment and allows for cost-effective screening of the response of nucleic acids to physiological changes (e.g. ion concentration) or screening of drugs in a high throughput fashion. This single technique allows us to probe an unprecedented range of exchange time scales (the whole μs to ms motion range) with increased sensitivity, surpassing all current experiments to detect chemical exchange. For the first time we could describe an RNA excited state using an unlabeled RNA.

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

One-dimensional heteronuclear relaxation dispersion NMR spectroscopy at ¹³C 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 D₂O. A small equilibrium isotope effect was observed and quantified, corresponding to ΔG = -0.14 kcal mol^-1. A bimolecular rate constant of k_D = 5.1 × 10⁹ s^-1 M^-1 was determined from the pH*-dependence of k_ex (where pH* is the direct electrode reading of pH in 10% D₂O and k_ex 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 D₂O-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.

Simultaneous determination of fast and slow dynamics in molecules using extreme CPMG relaxation dispersion experiments

Molecular dynamics play a significant role in how molecules perform their function. A critical method that provides information on dynamics, at the atomic level, is NMR-based relaxation dispersion (RD) experiments. RD experiments have been utilized for understanding multiple biological processes occurring at micro-to-millisecond time, such as enzyme catalysis, molecular recognition, ligand binding and protein folding. Here, we applied the recently developed high-power RD concept to the Carr-Purcell-Meiboom-Gill sequence (extreme CPMG; E-CPMG) for the simultaneous detection of fast and slow dynamics. Using a fast folding protein, gpW, we have shown that previously inaccessible kinetics can be accessed with the improved precision and efficiency of the measurement by using this experiment.

19F multiple-quantum coherence NMR spectroscopy for probing protein–ligand interactions

A new ¹⁹F NMR method is presented which can be used to detect weak protein binding of small molecules with up to mM affinity. The method capitalizes on the synthetic availability of unique SF₅ containing compounds and the generation of five-quantum coherences (5QC). Given the high sensitivity of 5QC relaxation to exchange events (i.e. reversible protein binding) fragments which bind to the target with weak affinity can be identified. The utility of the method in early stage drug discovery programs is demonstrated with applications to two model proteins, the neurotoxic NGAL and the prominent tumor target β-catenin.

An NMR experiment is presented that allows identification of weak binders typically found in early stages of drug discovery programs.

Label-free NMR-based dissociation kinetics determination

Understanding the dissociation of molecules is the basis to modulate interactions of biomedical interest. Optimizing drugs for dissociation rates is found to be important for their efficacy, selectivity, and safety. Here, we show an application of the high-power relaxation dispersion (RD) method to the determination of the dissociation rates of weak binding ligands from receptors. The experiment probes proton RD on the ligand and, therefore, avoids the need for any isotopic labeling. The large ligand excess eases the detection significantly. Importantly, the use of large spin-lock fields allows the detection of faster dissociation rates than other relaxation approaches. Moreover, this experimental approach allows to access directly the off-rate of the binding process without the need for analyzing a series of samples with increasing ligand saturation. The validity of the method is shown with small molecule interactions using two macromolecules, bovine serum albumin and tubulin heterodimers.