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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.