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.

Methyl probes in proteins for determining ligand binding mode in weak protein-ligand complexes.

Structures of protein-ligand complexes provide critical information for drug design. Most protein-ligand complex structures are determined using X-ray crystallography, but where crystallography is not able to generate a structure for a complex, NMR is often the best alternative. However, the available tools to enable rapid and robust structure determination of protein-ligand complexes by NMR are currently limited. This leads to situations where projects are either discontinued or pursued without structural data, rendering the task more difficult. We previously reported the NMR Molecular Replacement (NMR2) approach that allows the structure of a protein-ligand complex to be determined without requiring the cumbersome task of protein resonance assignment. Herein, we describe the NMR2 approach to determine the binding pose of a small molecule in a weak protein-ligand complex by collecting sparse protein methyl-to-ligand NOEs from a selectively labeled protein sample and an unlabeled ligand. In the selective labeling scheme all methyl containing residues of the protein are protonated in an otherwise deuterated background. This allows measurement of intermolecular NOEs with greater sensitivity using standard NOESY pulse sequences instead of isotope-filtered NMR experiments. This labelling approach is well suited to the NMR2 approach and extends its utility to include larger protein-ligand complexes.

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 (NMR2) 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 NMR2 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.

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

Exact nuclear Overhauser enhancement (eNOE) yields highly accurate, ensemble averaged 1H-1H 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 (NMR2) 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 NMR2. 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. NMR2 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 NMR2 with an average heavy atom RMSD of 3.628 Å. This is attributed to the higher precision of the evaluated distance restraints.

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.

Structure determination of protein-ligand complexes by NMR in solution.

In this paper, we discuss methods for determining structures of protein-ligand complexes by NMR in solution. Our discussion is based on small ligands (<2 kDa) as for example drugs, metabolites or oligo-peptides, but most of the considerations also apply to more general cases. In NMR in solution, the kinetics of association and dissociation of the complex - the exchange rate - determines the optimal sample preparation and the NMR experimental approach. Additionally, depending on the part of the complex that will be studied (only the bound ligand, the protein, the protein-ligand interface or the entire protein-ligand complex structure), different types of NMR experiments are needed. Therefore, the choice of a combination of the appropriate experiment and a suitable sample preparation in terms of ligand to protein ratios are discussed in detail. Also, considerations for practically preparing samples of protein-ligand complexes and carrying out experiments including trouble shooting are described. For structure determination, the scope of this paper is limited to NOE-based methods and some of the most recent approaches will be covered.

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.

The Exact Nuclear Overhauser Enhancement: Recent Advances.

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

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. NMR2 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 NMR2 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 Kd 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.

A Structural Ensemble for the Enzyme Cyclophilin Reveals an Orchestrated Mode of Action at Atomic Resolution.

For enzyme activity, an exact structural and motional orchestration of the active site and its surroundings is believed to be key. In order to reveal such possible phenomena at atomic resolution on the basis of experimental evidence, an experimental restraint driven two-state ensemble of the prototypical enzyme cyclophilin was determined by using a recently introduced exact NOE approach. The ensemble description reveals the presence of an open and a closed state of cyclophilin, which is indicative of large-scale correlated motion. In the open state, the catalytic site is preorganized for catalysis, thus suggesting the mechanism of action to be conformational sampling, while the ligand-binding loop appears to act through an induced fit mechanism. This finding is supported by affinity measurements of a cyclophilin designed to be more open. Overall, more than 60-70 % of the side-chain conformations of cyclophilin appear to be correlated.

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.