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.

A Benchmark Study of Protein-Fragment Complex Structure Calculations with NMR2

Protein-fragment complex structures are particularly sought after in medicinal chemistry to rationally design lead molecules. These structures are usually derived using X-ray crystallography, but the failure rate is non-neglectable. NMR is a possible alternative for the calculation of weakly interacting complexes. Nevertheless, the time-consuming protein signal assignment step remains a barrier to its routine application. NMR Molecular Replacement (NMR2) is a versatile and rapid method that enables the elucidation of a protein-ligand complex structure. It has been successfully applied to peptides, drug-like molecules, and more recently to fragments. Due to the small size of the fragments, ca < 300 Da, solving the structures of the protein-fragment complexes is particularly challenging. Here, we present the expected performances of NMR2 when applied to protein-fragment complexes. The NMR2 approach has been benchmarked with the SERAPhic fragment library to identify the technical challenges in protein-fragment NMR structure calculation. A straightforward strategy is proposed to increase the method's success rate further. The presented work confirms that NMR2 is an alternative method to X-ray crystallography for solving protein-fragment complex structures.

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.

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.

An improved, time-efficient approach to extract accurate distance restraints for NMR2 structure calculation

Exact nuclear Overhauser enhancement (eNOE) yields highly accurate, ensemble averaged 1 H-1 H distance restraints with an accuracy of up to 0.1 Å for the multi-state structure determination of proteins as well as for nuclear magnetic resonance molecular replacement (N MR2 ) to determine the structure of the protein-ligand interaction site in a time-efficient manner. However, in the latter application, the acquired eNOEs lack the obtainable precision of 0.1 Å because of the asymmetrical nature of the filtered nuclear Overhauser enhancement spectroscopy (NOESY) experiment used in N MR2 . This error is further propagated to the eNOE equations used to fit and extract the distance restraints. In this work, a new analysis method is proposed to obtain inter-molecular distance restraints from the filtered NOESY spectrum more accurately and intuitively by dividing the NOE cross peak by the corresponding diagonal peak of the ligand. The method termed diagonal-normalised eNOEs was tested on the data acquired by on the complex of PIN1 and a small, weak-binding phenylimidazole fragment. N MR2 calculations performed using the distances derived from diagonal-normalised eNOEs yielded the right orientation of the fragment in the binding pocket and produced a structure that more closely resembles the benchmark X-ray structure (2XP6) with an average heavy-atom root-mean-square deviation (RMSD) of 1.681 Å with respect to it, when compared to the one produced with traditional N MR2 with an average heavy atom RMSD of 3.628 Å. This is attributed to the higher precision of the evaluated distance restraints.

NMR Molecular Replacement Provides New Insights into Binding Modes to Bromodomains of BRD4 and TRIM24

Structure-based drug discovery (SBDD) largely relies on structural information from X-ray crystallography because traditional NMR structure calculation methods are too time consuming to be aligned with typical drug discovery timelines. The recently developed NMR molecular replacement (NMR²) method dramatically reduces the time needed to generate ligand-protein complex structures using published structures (apo or holo) of the target protein and treating all observed NOEs as ambiguous restraints, bypassing the laborious process of obtaining sequence-specific resonance assignments for the protein target. We apply this method to two therapeutic targets, the bromodomain of TRIM24 and the second bromodomain of BRD4. We show that the NMR² methodology can guide SBDD by rationalizing the observed SAR. We also demonstrate that new types of restraints and selective methyl labeling have the potential to dramatically reduce "time to structure" and extend the method to targets beyond the reach of traditional NMR structure elucidation.

Large-Scale Recombinant Production of the SARS-CoV-2 Proteome for High-Throughput and Structural Biology Applications

The highly infectious disease COVID-19 caused by the Betacoronavirus SARS-CoV-2 poses a severe threat to humanity and demands the redirection of scientific efforts and criteria to organized research projects. The international COVID19-NMR consortium seeks to provide such new approaches by gathering scientific expertise worldwide. In particular, making available viral proteins and RNAs will pave the way to understanding the SARS-CoV-2 molecular components in detail. The research in COVID19-NMR and the resources provided through the consortium are fully disclosed to accelerate access and exploitation. NMR investigations of the viral molecular components are designated to provide the essential basis for further work, including macromolecular interaction studies and high-throughput drug screening. Here, we present the extensive catalog of a holistic SARS-CoV-2 protein preparation approach based on the consortium's collective efforts. We provide protocols for the large-scale production of more than 80% of all SARS-CoV-2 proteins or essential parts of them. Several of the proteins were produced in more than one laboratory, demonstrating the high interoperability between NMR groups worldwide. For the majority of proteins, we can produce isotope-labeled samples of HSQC-grade. Together with several NMR chemical shift assignments made publicly available on covid19-nmr.com, we here provide highly valuable resources for the production of SARS-CoV-2 proteins in isotope-labeled form.

Protein Allostery at Atomic Resolution

Protein allostery is a phenomenon involving the long range coupling between two distal sites in a protein. In order to elucidate allostery at atomic resoluion on the ligand-binding WW domain of the enzyme Pin1, multistate structures were calculated from exact nuclear Overhauser effect (eNOE). In its free form, the protein undergoes a microsecond exchange between two states, one of which is predisposed to interact with its parent catalytic domain. In presence of the positive allosteric ligand, the equilibrium between the two states is shifted towards domain-domain interaction, suggesting a population shift model. In contrast, the allostery-suppressing ligand decouples the side-chain arrangement at the inter-domain interface thereby reducing the inter-domain interaction. As such, this mechanism is an example of dynamic allostery. The presented distinct modes of action highlight the power of the interplay between dynamics and function in the biological activity of proteins.

Protein-ligand structure determination with the NMR molecular replacement tool, NMR2

We recently reported on a new method called NMR Molecular Replacement that efficiently derives the structure of a protein-ligand complex at the interaction site. The method was successfully applied to high and low affinity complexes covering ligands from peptides to small molecules. The algorithm used in the NMR Molecular Replacement program has until now not been described in detail. Here, we present a complete description of the NMR Molecular Replacement implementation as well as several new features that further reduce the time required for structure elucidation.

Protein-fragment complex structures derived by NMR molecular replacement

Recently we have established an NMR molecular replacement method, which is capable of solving the structure of the interaction site of protein-ligand complexes in a fully automated manner. While the method was successfully applied for ligands with strong and weak binding affinities, including small molecules and peptides, its applicability on ligand fragments remains to be shown. Structures of fragment-protein complexes are more challenging for the method since fragments contain only few protons. Here we show a successful application of the NMR molecular replacement method in solving structures of complexes between three derivatives of a ligand fragment and the protein receptor PIN1. We anticipate that this approach will find a broad application in fragment-based lead discovery.

Rational Structure-Based Design of Fluorescent Probes for Amyloid Folds

Amyloid fibrils are pathological hallmarks of various human diseases, including Parkinson's, Alzheimer's, amyotrophic lateral sclerosis (ALS or motor neurone disease), and prion diseases. Treatment of the amyloid diseases are hindered, among other factors, by timely detection and therefore, early detection of the amyloid fibrils would be beneficial for treatment against these disorders. Here, a small molecular fluorescent probe is reported that selectively recognize the fibrillar form of amyloid beta(1-42), α-synuclein, and HET-s(218-289) protein over their monomeric conformation. The rational design of the reporters relies on the well-known cross-β-sheet repetition motif, the key structural feature of amyloids.

Proteome-wide analysis of phospho-regulated PDZ domain interactions

A key function of reversible protein phosphorylation is to regulate protein-protein interactions, many of which involve short linear motifs (3-12 amino acids). Motif-based interactions are difficult to capture because of their often low-to-moderate affinities. Here, we describe phosphomimetic proteomic peptide-phage display, a powerful method for simultaneously finding motif-based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C-terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD-95/Dlg/ZO-1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR We uncover site-specific phospho-regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho-regulation of motif-based interactions on a large scale.

eNORA2 Exact NOE Analysis Program

We have recently developed an NMR protocol to extract exact distances between nuclei in proteins from an exact interpretation of NOESY buildup intensities (eNOEs). This enabled us to calculate multistate structural ensembles that exhibit realistic spatial sampling and long-range correlations. Our initial studies were laborious and required a deep understanding of the underlying spin dynamics. Here, we present a MatLab package that integrates all data processing steps required to convert intensities of assigned peaks in NOESY series into upper and lower distance limits for structure calculation. Those steps include organization of the data in object format, extraction of autorelaxation and cross-relaxation rate constants by fitting of diagonal peak decays and cross peak buildups, validation of the data, correction for spin diffusion, graphical display of the results, and generation of distance limits in CYANA compatible format. The analysis may be carried out using a full relaxation matrix or a simplified "divide and conquer" approach that allows for partial deuteration of protons. As the program does not require expertise beyond that of standard resonance assignment/structure calculation, it is suitable for experts and nonexperts alike.

Fast NMR-Based Determination of the 3D Structure of the Binding Site of Protein-Ligand Complexes with Weak Affinity Binders

In early drug discovery approaches, screening hits are often weak affinity binders that are difficult to characterize in structural detail, particularly towards obtaining the 3D structure of protein-ligand complexes at atomic resolution. NMR is the outstanding technique to tackle such problems, yet suffers from a tedious structure calculation process. NMR² was recently developed to alleviate the laborious element of routine NMR structure calculation procedures and provides the structural information at protein-ligand interaction sites orders of magnitude faster than standard procedures. The NMR² method was extended to weak binders and applied to the oncoproteins HDM2 and MDMX. The structure of the MDMX-SJ212 complex is reported with a K_d of approximately 0.7 μm; the complex structure of HDM2 with the mm affinity ligand #845 exhibits a new scaffold.

Quenched hydrogen-deuterium exchange NMR of a disease-relevant Aβ(1-42) amyloid polymorph

Alzheimer’s disease is associated with the aggregation into amyloid fibrils of Aβ(1–42) and Aβ(1–40) peptides. Interestingly, these fibrils often do not obtain one single structure but rather show different morphologies, so-called polymorphs. Here, we compare quenched hydrogen-deuterium (H/D) exchange of a disease-relevant Aβ(1–42) fibril for which the 3D structure has been determined by solid-state NMR with H/D exchange previously determined on another structural polymorph. This comparison reveals secondary structural differences between the two polymorphs suggesting that the two polymorphisms can be classified as segmental polymorphs.

NMR-Based Determination of the 3D Structure of the Ligand-Protein Interaction Site without Protein Resonance Assignment

Molecular replacement in X-ray crystallography is the prime method for establishing structure-activity relationships of pharmaceutically relevant molecules. Such an approach is not available for NMR. Here, we establish a comparable method, called NMR molecular replacement (NMR(2)). The method requires experimentally measured ligand intramolecular NOEs and ligand-protein intermolecular NOEs as well as a previously known receptor structure or model. Our findings demonstrate that NMR(2) may open a new avenue for the fast and robust determination of the interaction site of ligand-protein complexes at atomic resolution.

The experimental accuracy of the uni-directional exact NOE

We have established protocols to calculate exact NOEs (eNOE) from NOE data. eNOEs lend unprecedented precision to the calculation of distance restraints used for structure calculation. Moreover, as eNOEs are averaged quantities over all conformations of a molecule, they may contain accessible information of the sampled conformational space. In practice, a prerequisite for an exact interpretation is the evaluation of both NOESY cross-peak buildups. For large molecular sizes, the fraction of NOEs which can only be obtained from one cross peak typically increases. Distance restraints derived from such NOEs must be used with a tolerance for errors associated with the broken symmetry of the individual magnetization transfer pathways. The correct choice of upper and lower limits is particularly important for multiple-state ensemble calculation, where too narrow tolerances may lead to incorrect spatial sampling. In order to dissect these pathways in heavy-atom resolved 3D NOESY experiments, we analyze 2D [(1)H, (1)H]-NOESY experiments, which are the fundamental building blocks of the former. In combination with an analysis of excitation and inversion profiles of pulses on heavy atoms and relaxation effects during HXQC elements, we derive a rule for the correct choice of upper and lower distance limits derived from such uni-directional NOEs. We show that normalization of the cross- to the diagonal-peak intensities of the spins of magnetization destination rather than origin leads to similar errors of the distance restraints. This opens up the prospect of extended collection of unidirectional eNOEs.

Solution NMR studies of recombinant Aβ(1-42): from the presence of a micellar entity to residual β-sheet structure in the soluble species

Amyloid-β (Aβ) peptide is the major component found in senile plaques of Alzheimer's disease patients. The 42-residue fragment Aβ(1-42) is proposed to be one of the most pathogenic species therein. Here, the soluble Aβ(1-42) species were analyzed by various liquid-state NMR methods. Transient formation of a micelle species was observed at the onset of the aggregation kinetics. This micelle is dissolved after approximately one day. Subsequent loss of this species and the formation of protofibrils are proposed to be the route of fibril formation. Consequently, the observed micelle species is suggested to be on an off-pathway mechanism. Furthermore, characterization of the NMR-observable soluble species shows that it is a random-coil-like entity with low propensities for four β-strands. These β-strands correlate with the β-strand segments observed in Aβ fibrils. This finding indicates that the 3D structure of the fibrils might already be predisposed in the soluble species.

Intermolecular detergent-membrane protein noes for the characterization of the dynamics of membrane protein-detergent complexes

The interaction between membrane proteins and lipids or lipid mimetics such as detergents is key for the three-dimensional structure and dynamics of membrane proteins. In NMR-based structural studies of membrane proteins, qualitative analysis of intermolecular nuclear Overhauser enhancements (NOEs) or paramagnetic resonance enhancement are used in general to identify the transmembrane segments of a membrane protein. Here, we employed a quantitative characterization of intermolecular NOEs between (1)H of the detergent and (1)H(N) of (2)H-perdeuterated, (15)N-labeled α-helical membrane protein-detergent complexes following the exact NOE (eNOE) approach. Structural considerations suggest that these intermolecular NOEs should show a helical-wheel-type behavior along a transmembrane helix or a membrane-attached helix within a membrane protein as experimentally demonstrated for the complete influenza hemagglutinin fusion domain HAfp23. The partial absence of such a NOE pattern along the amino acid sequence as shown for a truncated variant of HAfp23 and for the Escherichia coli inner membrane protein YidH indicates the presence of large tertiary structure fluctuations such as an opening between helices or the presence of large rotational dynamics of the helices. Detergent-protein NOEs thus appear to be a straightforward probe for a qualitative characterization of structural and dynamical properties of membrane proteins embedded in detergent micelles.

Towards a true protein movie: a perspective on the potential impact of the ensemble-based structure determination using exact NOEs

Confined by the Boltzmann distribution of the energies of the states, a multitude of structural states are inherent to biomolecules. For a detailed understanding of a protein's function, its entire structural landscape at atomic resolution and insight into the interconversion between all the structural states (i.e. dynamics) are required. Whereas dedicated trickery with NMR relaxation provides aspects of local dynamics, and 3D structure determination by NMR is well established, only recently have several attempts been made to formulate a more comprehensive description of the dynamics and the structural landscape of a protein. Here, a perspective is given on the use of exact NOEs (eNOEs) for the elucidation of structural ensembles of a protein describing the covered conformational space.

Stereospecific assignments in proteins using exact NOEs

Recently developed methods to measure distances in proteins with high accuracy by "exact" nuclear Overhauser effects (eNOEs) make it possible to determine stereospecific assignments, which are particularly important to fully exploit the accuracy of the eNOE distance measurements. Stereospecific assignments are determined by comparing the eNOE-derived distances to protein structure bundles calculated without stereospecific assignments, or an independently determined crystal structure. The absolute and relative CYANA target function difference upon swapping the stereospecific assignment of a diastereotopic group yields the respective stereospecific assignment. We applied the method to the eNOE data set that has recently been obtained for the third immunoglobulin-binding domain of protein G (GB3). The 884 eNOEs provide relevant data for 47 of the total of 75 diastereotopic groups. Stereospecific assignments could be established for 45 diastereotopic groups (96 %) using the X-ray structure, or for 27 diastereotopic groups (57 %) using structures calculated with the eNOE data set without stereospecific assignments, all of which are in agreement with those determined previously. The latter case is relevant for structure determinations based on eNOEs. The accuracy of the eNOE distance measurements is crucial for making stereospecific assignments because applying the same method to the traditional NOE data set for GB3 with imprecise upper distance bounds yields only 13 correct stereospecific assignments using the X-ray structure or 2 correct stereospecific assignments using NMR structures calculated without stereospecific assignments.

Discrete three-dimensional representation of macromolecular motion from eNOE-based ensemble calculation

Three-dimensional structural data and description of dynamics are fundamental to infer and understand protein function. Structure determination by NMR follows well-established protocols while NMR relaxation phenomena provide insights into local molecular dynamics. However, methods to detect concerted motion were not pursued until very recently. Here, we present an ensemble-based structure determination protocol using ensemble-averaged distance restraints obtained from exact NOE (eNOE) rate constants. An application of our protocol to the model protein GB3 established an ensemble of structures that reveals correlated motion across the β-sheet and concerted motion between the backbone and side chains localized in the core. Furthermore, the data repudiate concerted conformational exchange between the β-sheet and the α-helix.

The description of protein internal motions aids selection of ligand binding poses by the INPHARMA method

Protein internal motions influence observables of NMR experiments. The effect of internal motions occurring at the sub-nanosecond timescale can be described by NMR order parameters. Here, we report that the use of order parameters derived from Molecular Dynamics (MD) simulations of two holo-structures of Protein Kinase A increase the discrimination power of INPHARMA, an NMR based methodology that selects docked ligand orientations by maximizing the correlation of back-calculated to experimental data. By including internal motion in the back-calculation of the INPHARMA transfer, we obtain a more realistic description of the system, which better represents the experimental data. Furthermore, we propose a set of generic order parameters, derived from MD simulations of globular proteins, which can be used in the back-calculation of INPHARMA NOEs for any protein-ligand complex, thus by-passing the need of obtaining system-specific order parameters for new protein-ligand complexes.

Relaxation Matrix Analysis of Spin Diffusion for the NMR Structure Calculation with eNOEs

NMR structure determination is usually based on distance restraints extracted semiquantitatively from cross peak volumes or intensities in NOESY spectra. The recent introduction of exact NOEs (eNOE) by Vogeli et al. opens an avenue for the ensemble-based structure determination of proteins on the basis of eNOE-derived quantitative distance restraints. We present an approach to extract eNOE from build-up curve intensities. For the determination of eNOEs, spin diffusion is a major source of errors. A full relaxation matrix analysis is used to calculate the spin diffusion contribution to the NOESY cross peaks of each individual spin pair of interest. A software program is written, which requires as input the peak intensities from the various NOESY spectra as well as a 3D structure of the protein. This structure can be either an X-ray structure or an NMR structure determined with the conventional approach. The outputs of the program are the eNOE rates, the autorelaxation rates, as well as graphs and quality factors from the individual NOE build-up curves for semiautomated analysis of the derived rates. The protocol is straightforward, and the program integrates well into the current structure calculation workflow.

An NMR-based scoring function improves the accuracy of binding pose predictions by docking by two orders of magnitude

Low-affinity ligands can be efficiently optimized into high-affinity drug leads by structure based drug design when atomic-resolution structural information on the protein/ligand complexes is available. In this work we show that the use of a few, easily obtainable, experimental restraints improves the accuracy of the docking experiments by two orders of magnitude. The experimental data are measured in nuclear magnetic resonance spectra and consist of protein-mediated NOEs between two competitively binding ligands. The methodology can be widely applied as the data are readily obtained for low-affinity ligands in the presence of non-labelled receptor at low concentration. The experimental inter-ligand NOEs are efficiently used to filter and rank complex model structures that have been pre-selected by docking protocols. This approach dramatically reduces the degeneracy and inaccuracy of the chosen model in docking experiments, is robust with respect to inaccuracy of the structural model used to represent the free receptor and is suitable for high-throughput docking campaigns.

The INPHARMA technique for pharmacophore mapping: A theoretical guide to the method

During the process of drug discovery, INPHARMA can be used to derive the structure of receptor/lead compound complexes binding to each other with a K(d) in the microM to mM range. To be successful, the methodology needs adjustment of various parameters that depend on the physical constants of the binding event and on the receptor size. Here we present a thorough theoretical analysis of the INPHARMA interligand NOE effect in dependence of experimental parameters and physical constants. This analysis helps the experimentalist to choose the correct experimental parameters and consequentially to achieve optimal performance of the methodology.

Specimen collection time for enzyme analysis in urine

Catalytic activity of N-acetyl-beta-D-glucosaminidase, alanine aminopeptidase, lactate dehydrogenase, isoenzyme 1 of lactate dehydrogenase, lysozyme, gamma-glutamyl transferase and alkaline phosphatase in urine specimens collected between 6 a.m. and 9 a.m. were determined in 25 patients with acute renal failure. We found no statistical differences (Wilcoxon's t test) between specimens collected at 6 a.m. and 9 a.m. We conclude that, in renal patients, the first morning specimen (overnight urine) may be used for enzyme analysis.