Transferred 13C T1 Relaxation at Natural Isotopic Abundance: A Practical Method for Determining Site-Specific Changes in Ligand Flexibility upon Binding to a Macromolecule

NMR spectroscopy can help accelerate antiviral drug discovery programs

Small molecule drugs have an important role to play in combating viral infections, and biophysics support has been central for contributing to the discovery and design of direct acting antivirals. Perhaps one of the most successful biophysical tools for this purpose is NMR spectroscopy when utilized strategically and pragmatically within team workflows and timelines. This report describes some clear examples of how NMR applications contributed to the design of antivirals when combined with medicinal chemistry, biochemistry, X-ray crystallography and computational chemistry. Overall, these multidisciplinary approaches allowed teams to reveal and expose compound physical properties from which design ideas were spawned and tested to achieve the desired successes. Examples are discussed for the discovery of antivirals that target HCV, HIV and SARS-CoV-2.

Structure-thermodynamics-relationships of hepatitis C viral NS3 protease inhibitors

Thermodynamic parameters were determined for structurally-related inhibitors of HCV NS3 protease to assess how binding entropies and enthalpies vary with incremental changes at the P2 and P3 inhibitor subsites. Changing the heterocyclic substituent at P2 from a pyridyl to a 7-methoxy-2-phenyl-4-quinolyl group leads to a 710-fold increase in affinity. Annelating a benzene ring onto a pyridine ring leads to quinoline-derived inhibitors having higher affinities, but the individual enthalpy and entropy contributions are markedly different for each ligand pair. Introducing a phenyl group at C2 of the heterocyclic ring at P2 uniformly leads to higher affinity analogs with more favorable binding entropies, while adding a methoxy group at C7 of the quinoline ring at P2 provides derivatives with more favorable binding enthalpies. Significant enthalpy/entropy compensation is observed for structural changes made to inhibitors lacking a 2-phenyl substituent, whereas favorable changes in both binding enthalpies and entropies accompany structural modifications when a 2-phenyl group is present. Overall, binding energetics of inhibitors having a 2-phenyl-4-quinolyl group at P2 are dominated by entropic effects, whereas binding of the corresponding norphenyl analogs are primarily enthalpy driven. Notably, the reversal from an entropy driven association to an enthalpy driven one for this set of inhibitors also correlates with alternate binding modes. When the steric bulk of the side chain at P3 is increased from a hydrogen atom to a tert-butyl group, there is a 770-fold improvement in affinity. The 30-fold increase resulting from the first methyl group is solely the consequence of a more favorable change in entropy, whereas subsequent additions of methyl groups leads to modest increases in affinity that arise primarily from incremental improvements in binding enthalpies accompanied with smaller favorable entropic contributions.

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

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

2015 Philip S. Portoghese Medicinal Chemistry Lectureship. Curing Hepatitis C Virus Infection with Direct-Acting Antiviral Agents: The Arc of a Medicinal Chemistry Triumph

The development of direct-acting antiviral agents that can cure a chronic hepatitis C virus (HCV) infection after 8-12 weeks of daily, well-tolerated therapy has revolutionized the treatment of this insidious disease. In this article, three of Bristol-Myers Squibb's HCV programs are summarized, each of which produced a clinical candidate: the NS3 protease inhibitor asunaprevir (64), marketed as Sunvepra, the NS5A replication complex inhibitor daclatasvir (117), marketed as Daklinza, and the allosteric NS5B polymerase inhibitor beclabuvir (142), which is in late stage clinical studies. A clinical study with 64 and 117 established for the first time that a chronic HCV infection could be cured by treatment with direct-acting antiviral agents alone in the absence of interferon. The development of small molecule HCV therapeutics, designed by medicinal chemists, has been hailed as "the arc of a medical triumph" but may equally well be described as "the arc of a medicinal chemistry triumph".

Quantitative analysis of protein-ligand interactions by NMR

Protein-ligand interactions have been commonly studied through static structures of the protein-ligand complex. Recently, however, there has been increasing interest in investigating the dynamics of protein-ligand interactions both for fundamental understanding of the underlying mechanisms and for drug development. NMR is a versatile and powerful tool, especially because it provides site-specific quantitative information. NMR has widely been used to determine the dissociation constant (KD), in particular, for relatively weak interactions. The simplest NMR method is a chemical-shift titration experiment, in which the chemical-shift changes of a protein in response to ligand titration are measured. There are other quantitative NMR methods, but they mostly apply only to interactions in the fast-exchange regime. These methods derive the dissociation constant from population-averaged NMR quantities of the free and bound states of a protein or ligand. In contrast, the recent advent of new relaxation-based experiments, including R2 relaxation dispersion and ZZ-exchange, has enabled us to obtain kinetic information on protein-ligand interactions in the intermediate- and slow-exchange regimes. Based on R2 dispersion or ZZ-exchange, methods that can determine the association rate, kon, dissociation rate, koff, and KD have been developed. In these approaches, R2 dispersion or ZZ-exchange curves are measured for multiple samples with different protein and/or ligand concentration ratios, and the relaxation data are fitted to theoretical kinetic models. It is critical to choose an appropriate kinetic model, such as the two- or three-state exchange model, to derive the correct kinetic information. The R2 dispersion and ZZ-exchange methods are suitable for the analysis of protein-ligand interactions with a micromolar or sub-micromolar dissociation constant but not for very weak interactions, which are typical in very fast exchange. This contrasts with the NMR methods that are used to analyze population-averaged NMR quantities. Essentially, to apply NMR successfully, both the type of experiment and equation to fit the data must be carefully and specifically chosen for the protein-ligand interaction under analysis. In this review, we first explain the exchange regimes and kinetic models of protein-ligand interactions, and then describe the NMR methods that quantitatively analyze these specific interactions.

Conformation-based restrictions and scaffold replacements in the design of hepatitis C virus polymerase inhibitors: discovery of deleobuvir (BI 207127)

Conformational restrictions of flexible torsion angles were used to guide the identification of new chemotypes of HCV NS5B inhibitors. Sites for rigidification were based on an acquired conformational understanding of compound binding requirements and the roles of substituents in the free and bound states. Chemical bioisosteres of amide bonds were explored to improve cell-based potency. Examples are shown, including the design concept that led to the discovery of the phase III clinical candidate deleobuvir (BI 207127). The structure-based strategies employed have general utility in drug design.

Ligand bioactive conformation plays a critical role in the design of drugs that target the hepatitis C virus NS3 protease

A ligand-focused strategy employed NMR, X-ray, modeling, and medicinal chemistry to expose the critical role that bioactive conformation played in the design of a variety of drugs that target the HCV protease. The bioactive conformation (bound states) were determined for key inhibitors identified along our drug discovery pathway from the hit to clinical compounds. All adopt similar bioactive conformations for the common core derived from the hit peptide DDIVPC. A carefully designed SAR analysis, based on the advanced inhibitor 1 in which the P1 to P3 side chains and the N-terminal Boc were sequentially truncated, revealed a correlation between affinity and the relative predominance of the bioactive conformation in the free state. Interestingly, synergistic conformation effects on potency were also noted. Comparisons with clinical and recently marketed drugs from the pharmaceutical industry showed that all have the same core and similar bioactive conformations. This suggested that the variety of appendages discovered for these compounds also properly satisfy the bioactive conformation requirements and allowed for a large variety of HCV protease drug candidates to be designed.

Elucidating slow binding kinetics of a protein without observable bound resonances by longitudinal relaxation NMR spectroscopy

We developed a new method to elucidate the binding kinetics k(on) and k(off), and the dissociation constant K(D) (=k(off)/k(on)), of protein-protein interactions without observable bound resonances of the protein of interest due to high molecular weight in a complex with a large target protein. In our method, k(on) and k(off) rates are calculated from the analysis of longitudinal relaxation rates of free resonances measured for multiple samples containing different concentration ratios of (15)N-labeled protein and substoichiometric amounts of the target protein. The method is applicable to interactions that cannot be analyzed by relaxation dispersion spectroscopy due to slow interactions on millisecond to second timescale and/or minimal conformational (chemical shift) change upon binding. We applied the method to binding of the B1 domain of protein G (GB1) to immunoglobulin G, and derived the binding kinetics despite the absence of observable bound GB1 resonances.

Importance of ligand bioactive conformation in the discovery of potent indole-diamide inhibitors of the hepatitis C virus NS5B

Significant advances have led to receptor induced-fit and conformational selection models for describing bimolecular recognition, but a more comprehensive view must evolve to also include ligand shape and conformational changes. Here, we describe an example where a ligand's "structural hinge" influences potency by inducing an "L-shape" bioactive conformation, and due to its solvent exposure in the complex, reasonable conformation-activity-relationships can be qualitatively attributed. From a ligand design perspective, this feature was exploited by successful linker hopping to an alternate "structural hinge" that led to a new and promising chemical series which matched the ligand bioactive conformation and the pocket bioactive space. Using a combination of X-ray crystallography, NMR and modeling with support from binding-site resistance mutant studies and photoaffinity labeling experiments, we were able to derive inhibitor-polymerase complexes for various chemical series.

Study of molecular interactions with 13C DNP-NMR

NMR spectroscopy is an established, versatile technique for the detection of molecular interactions, even when these interactions are weak. Signal enhancement by several orders of magnitude through dynamic nuclear polarization alleviates several practical limitations of NMR-based interaction studies. This enhanced non-equilibrium polarization contributes sensitivity for the detection of molecular interactions in a single NMR transient. We show that direct (13)C NMR ligand binding studies at natural isotopic abundance of (13)C gets feasible in this way. Resultant screens are easy to interpret and can be performed at (13)C concentrations below muM. In addition to such ligand-detected studies of molecular interaction, ligand binding can be assessed and quantified with enzymatic assays that employ hyperpolarized substrates at varying enzyme inhibitor concentrations. The physical labeling of nuclear spins by hyperpolarization thus provides the opportunity to devise fast novel in vitro experiments with low material requirement and without the need for synthetic modifications of target or ligands.

Dynamics of ligand binding from 13C NMR relaxation dispersion at natural abundance

We show that Carr-Purcell-Meiboom-Gill (CPMG) 13Calpha NMR relaxation dispersion measurements are a viable means for profiling mus-ms ligand dynamics involved in receptor binding. Critically, the dispersion is at natural 13C abundance; this matches typical pharmaceutical research settings in which ligand isotope-labeling is often impractical. The dispersion reveals ligand 13Calpha nuclei that experience mus-ms modulation of their chemical shifts due to binding. 13Calpha shifts are dominated by local torsion angles , psi, chi1; hence, these experiments identify flexible torsion angles that may assist complex formation. Since the experiments detect the ligand, they are viable even in the absence of a receptor structure. The mus-ms dynamic information gained helps establish flexibility-activity relationships. We apply these experiments to study the binding of a phospho-peptide substrate ligand to the peptidyl-prolyl isomerase Pin1.

Enhanced signal dispersion in saturation transfer difference experiments by conversion to a 1D-STD-homodecoupled spectrum

The saturation transfer difference (STD) experiment is a rich source of information on topological aspects of ligand binding to a receptor. The epitope mapping is based on a magnetization transfer after signal saturation from the receptor to the ligand, where interproton distances permit this process. Signal overlap in the STD spectrum can cause difficulties to correctly assign and/or quantitate the measured enhancements. To address this issue we report here a modified version of the routine experiment and a processing scheme that provides a 1D-STD homodecoupled spectrum (i.e. an experiment in which all STD signals appear as singlets) with line widths similar to those in original STD spectrum. These refinements contribute to alleviate problems of signal overlap. The experiment is based on 2D-J-resolved spectroscopy, one of the fastest 2D experiments under conventional data sampling in the indirect dimension, and provides excellent sensitivity, a key factor for the difference experiments.

NMR characterization of paclitaxel/poly (styrene-isobutylene-styrene) formulations

TAXUStrade mark is a coronary drug-eluting stent system utilizing a formulation consisting of cellular-target drug paclitaxel and poly (styrene-isobutylene-styrene) (SIBS). The present study investigates the interaction and interfacial dynamics of paclitaxel incorporated in a nano-polymeric matrix system. Solution and solid-state CP/MAS NMR experiments were designed to characterize the microstructure of heterogeneous drug-polymer mixtures in terms of its composition, molecular mobility, molecular order, paclitaxel-SIBS molecular interactions, and molecular mobility of the drug in the polymer matrix. The NMR spectra demonstrated unchanged chemical shifts between the neat and incorporated paclitaxel, and suggested that the level of the interactions between paclitaxel and SIBS is limited to non-bonding interactions or physical interactions between paclitaxel and SIBS when mixed in solution under NMR detection. Carbon spin-lattice relaxation time and proton spin-lattice relaxation time in the rotating frame offer further confirmation that the mobility of paclitaxel is increased in the paclitaxel-SIBS mixture. The results also indicate that a change occurs from crystalline packing to amorphous packing in paclitaxel due to its intermolecular interaction with SIBS. Our studies were used in understanding the detailed structure, morphology, and molecular motion of paclitaxel in the paclitaxel-SIBS system and to probe chemical and physical heterogeneity down to the nanometer scale.

A Systematic Approach to the Optimization of Substrate-Based Inhibitors of the Hepatitis C Virus NS3 Protease: Discovery of Potent and Specific Tripeptide Inhibitors

The inadequate efficacy and tolerability of current therapies for the infectious liver disease caused by the hepatitis C virus have warranted significant efforts in the development of new therapeutics. We have previously reported competitive peptide inhibitors of the NS3 serine protease based on the N-terminal cleavage products of peptide substrates. A detailed study of the interactions of these substrate-based inhibitors with the different subsites of the serine protease active site led to the discovery of novel residues that increased the affinity of the inhibitors. In this paper, we report the combination of the best binding residues in a tetrapeptide series that resulted in extremely potent inhibitors that bind exquisitely well to this enzyme. A substantial increase in potency was obtained with the simultaneous introduction of a 7-methoxy-2-phenyl-4-quinolinoxy moiety at the gamma-position of the P2 proline and a tert-leucine as a P3 residue. The increase in potency allowed for the further truncation and led to the identification of tripeptide inhibitors. Structure activity relationship studies on this inhibitor series led to the identification of carbamate-containing tripeptides that are able to inhibit replication of subgenomic HCV RNA in cell culture with potencies below 1 microM. This inhibitor series has the potential of becoming antiviral agents for the treatment of HCV infections.

Potent Inhibitors of the Hepatitis C Virus NS3 Protease: Design and Synthesis of Macrocyclic Substrate-Based β-Strand Mimics

The virally encoded NS3 protease is essential to the life cycle of the hepatitis C virus (HCV), an important human pathogen causing chronic hepatitis, cirrhosis of the liver, and hepatocellular carcinoma. The design and synthesis of 15-membered ring beta-strand mimics which are capable of inhibiting the interactions between the HCV NS3 protease enzyme and its polyprotein substrate will be described. The binding interactions between a macrocyclic ligand and the enzyme were explored by NMR and molecular dynamics, and a model of the ligand/enzyme complex was developed.

Using NMR for ligand discovery and optimization

Several recent technology-driven advances in the area of NMR have rekindled an interest in the application of the technology to problems in drug discovery and development. A unique aspect of NMR is that it has applicability in broadly different areas of the drug discovery and optimization processes. NMR techniques for screening aimed at the discovery of novel ligands or low molecular weight structures for fragment-based build up procedures are being applied commonly in the industry. Application of NMR in structure-guided drug design and metabonomics are also becoming routine. We present an overview of some of the most recent NMR developments in these areas.

The design of a potent inhibitor of the hepatitis C virus NS3 protease:BILN 2061?From the NMR tube to the clinic

The virally encoded serine protease NS3/NS4A is essential to the life cycle of the hepatitis C virus (HCV), an important human pathogen causing chronic hepatitis, cirrhosis of the liver, and hepatocellular carcinoma. Until very recently, the design of inhibitors for the HCV NS3 protease was limited to large peptidomimetic compounds with poor pharmacokinetic properties, making drug discovery an extremely challenging endeavor. In our quest for the discovery of a small-molecule lead that could block replication of the hepatitis C virus by binding to the HCV NS3 protease, the critical protein-polypeptide interactions between the virally encoded NS3 serine protease and its polyprotein substrate were investigated. Lead optimization of a substrate-based hexapeptide, guided by structural data, led to the understanding of the molecular dynamics and electronic effects that modulate the affinity of peptidomimetic ligands for the active site of this enzyme. Macrocyclic beta-strand scaffolds were designed that allowed the discovery of potent, highly selective, and orally bioavailable compounds. These molecules were the first HCV NS3 protease inhibitors reported that inhibit replication of HCV subgenomic RNA in a cell-based replicon assay at low nanomolar concentrations. Optimization of their biopharmaceutical properties led to the discovery of the clinical candidate BILN 2061. Oral administration of BILN 2061 to patients infected with the hepatitis C genotype 1 virus resulted in an impressive reduction of viral RNA levels, establishing proof-of-concept for HCV NS3 protease inhibitors as therapeutic agents in humans.

New probes of ligand flexibility in drug design: transferred (13)C CSA-dipolar cross-correlated relaxation at natural abundance

Understanding the impact of molecular flexibility remains an important outstanding problem in rational drug design. Toward this end, we present new NMR relaxation methods that describe ligand flexibility at the atomic level. Specifically, we measure natural abundance (13)C cross-correlated relaxation parameters for ligands in rapid exchange between the free and receptor-bound states. The rapid exchange transfers the bound state relaxation parameters to the free state, such that a comparison of relaxation rates in the absence and presence of protein receptor yields site-specific information concerning the bound ligand flexibility. We perform these measurements for aromatic carbons, which are highly prevalent in drug-like molecules and demonstrate significant cross-correlated relaxation between the (13)C-(1)H dipole-dipole (DD) and (13)C chemical shift anisotropy (CSA) relaxation mechanisms. Our use of natural abundance measurements addresses the practical difficulties of obtaining isotope-labeled ligands in pharmaceutical research settings. We demonstrate our methods on a small ligand of the 42 kDa kinase domain of the p38 MAP kinase. We show that exchange-transferred cross-correlated relaxation measurements are not only sensitive probes of bound ligand flexibility but also offer complementary advantages over standard R(1) = 1/T(1) and R(2) = 1/T(2) measurements. The ligand flexibility profiles obtained from the relaxation data can help assess the influence of dynamics on ligand potency or pharmacokinetic properties or both, and thereby include inherent molecular flexibility in drug design.

The effect of relaxation on the epitope mapping by saturation transfer difference NMR

The effect of longitudinal relaxation of ligand protons on saturation transfer difference (STD) was investigated by using a known binding system, dihydrofolate reductase and trimethoprim. The results indicate that T1 relaxation of ligand protons has a severe interference on the epitope map derived from a STD measurement. When the T1s of individual ligand protons are distinctly different, STD experiments may not give an accurate epitope map for the ligand-target interactions. Measuring the relaxation times prior to mapping is strongly advised. A saturation time shorter than T1s is suggested for improving the potential epitope map. Reduction in temperature was seen to enhance the saturation efficiency in small to medium size targets.

Applications of NMR in drug discovery

NMR, already some 50 years old, has long been an invaluable analytical method in industry for verification of chemical synthesis and compound characterisation. The range of molecular information accessible through NMR, however, offers a far larger horizon of applications. Of these, ligand screening by NMR has emerged as a very promising new method in drug discovery. Its unmatched screening sensitivity, combined with the abundance of available information on the structure and nature of molecular binding, justifies the growing interest in this dynamically expanding NMR application.