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

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.

Determining the Predominant Conformations of Mortiamides A-D in Solution Using NMR Data and Molecular Modeling Tools

Macrocyclic peptidomimetics have been seriously contributing to our arsenal of drugs to combat diseases. The search for nature's discoveries led us to mortiamides A-D (found in a novel fungus from Northern Canada), which is a family of cyclic peptides that clearly have demonstrated impressive pharmaceutical potential. This prompted us to learn more about their solution-state properties as these are central for binding to target molecules. Here, we secured and isolated mortiamide D, and then acquired high-resolution nuclear magnetic resonance (NMR) data to learn more about its structure and dynamics attributes. Sets of two-dimensional NMR experiments provided atomic-level (through-bond and through-space) data to confirm the primary structure, and NMR-driven molecular dynamics (MD) simulations suggested that more than one predominant three-dimensional (3D) structure exist in solution. Further steps of MD simulations are consistent with the finding that the backbones of mortiamides A-C also have at least two prominent macrocyclic shapes, but the side-chain structures and dynamics differed significantly. Knowledge of these solution properties can be exploited for drug design and discovery.

Varying molecular interactions explain aspects of crowder-dependent enzyme function of a viral protease

Biochemical processes in cells, including enzyme-catalyzed reactions, occur in crowded conditions with various background macromolecules occupying up to 40% of cytoplasm's volume. Viral enzymes in the host cell also encounter such crowded conditions as they often function at the endoplasmic reticulum membranes. We focus on an enzyme encoded by the hepatitis C virus, the NS3/4A protease, which is crucial for viral replication. We have previously found experimentally that synthetic crowders, polyethylene glycol (PEG) and branched polysucrose (Ficoll), differently affect the kinetic parameters of peptide hydrolysis catalyzed by NS3/4A. To gain understanding of the reasons for such behavior, we perform atomistic molecular dynamics simulations of NS3/4A in the presence of either PEG or Ficoll crowders and with and without the peptide substrates. We find that both crowder types make nanosecond long contacts with the protease and slow down its diffusion. However, they also affect the enzyme structural dynamics; crowders induce functionally relevant helical structures in the disordered parts of the protease cofactor, NS4A, with the PEG effect being more pronounced. Overall, PEG interactions with NS3/4A are slightly stronger but Ficoll forms more hydrogen bonds with NS3. The crowders also interact with substrates; we find that the substrate diffusion is reduced much more in the presence of PEG than Ficoll. However, contrary to NS3, the substrate interacts more strongly with Ficoll than with PEG crowders, with the substrate diffusion being similar to crowder diffusion. Importantly, crowders also affect the substrate-enzyme interactions. We observe that both PEG and Ficoll enhance the presence of substrates near the active site, especially near catalytic H57 but Ficoll crowders increase substrate binding more than PEG molecules.

Evaluation of interactions between the hepatitis C virus NS3/4A and sulfonamidobenzamide based molecules using molecular docking, molecular dynamics simulations and binding free energy calculations

The Hepatitis C Virus (HCV) NS3/4A is an attractive target for the treatment of Hepatitis C infection. Herein, we present an investigation of HCV NS3/4A inhibitors based on a sulfonamidobenzamide scaffold. Inhibitor interactions with HCV NS3/4A were explored by molecular docking, molecular dynamics simulations, and MM/PBSA binding free energy calculations. All of the inhibitors adopt similar molecular docking poses in the catalytic site of the protease that are stabilized by hydrogen bond interactions with G137 and the catalytic S139, which are known to be important for potency and binding stability. The quantitative assessments of binding free energies from MM/PBSA correlate well with the experimental results, with a high coefficient of determination, R2 of 0.92. Binding free energy decomposition analyses elucidate the different contributions of Q41, F43, H57, R109, K136, G137, S138, S139, A156, M485, and Q526 in binding different inhibitors. The importance of these sidechain contributions was further confirmed by computational alanine scanning mutagenesis. In addition, the sidechains of K136 and S139 show crucial but distinct contributions to inhibitor binding with HCV NS3/4A. The structural basis of the potency has been elucidated, demonstrating the importance of the R155 sidechain conformation. This extensive exploration of binding energies and interactions between these compounds and HCV NS3/4A at the atomic level should benefit future antiviral drug design.

Cross-Linked Poly-4-Acrylomorpholine: A Flexible and Reversibly Compressible Aligning Gel for Anisotropic NMR Analysis of Peptides and Small Molecules in Water

Determination of the solution conformation of both small organic molecules and peptides in water remains a substantial hurdle in using NMR solution conformations to guide drug design due to the lack of easy to use alignment media. Herein we report the design of a flexible compressible chemically cross-linked poly-4-acrylomorpholine gel that can be used for the alignment of both small molecules and cyclic peptides in water. To test the new gel, residual dipolar couplings (RDCs) and J-coupling constants were used in the configurational analysis of strychnine hydrochloride, a molecule that has been studied extensively in organic solvents as well as a small cyclic peptide that is known to form an α-helix in water. The conformational ensembles for each molecule with the best fit to the data are reported. Identification of minor conformers in water that cannot easily be determined by conventional NOE measurements will facilitate the use of RDC experiments in structure-based drug design.

Crowding affects structural dynamics and contributes to membrane association of the NS3/4A complex

Using molecular dynamics simulations, we describe how crowded environments affect the internal dynamics and diffusion of the hepatitis C virus proteases NS3/4A. This protease plays a key role in viral replication and is successfully used as a target for antiviral treatment. The NS3 enzyme requires a peptide cofactor, called NS4A, with its central part interacting with the NS3 β-sheet, and flexible, protruding terminal tails that are unstructured in water solution. The simulations describe the enzyme and water molecules at atomistic resolution, whereas crowders are modeled via either all-atom or coarse-grained models to emphasize different aspects of crowding. Crowders reflect the polyethylene glycol (PEG) molecules used in the experiments to mimic the crowded surrounding. A bead-shell model of folded coarse-grained PEG molecules considers mainly the excluded volume effect, whereas all-atom PEG models afford more protein-like crowder interactions. Circular dichroism spectroscopy experiments of the NS4A N-terminal tail show that a helical structure is formed in the presence of PEG crowders. The simulations suggest that crowding may assist in the formation of an NS4A helical fragment, positioned exactly where a transmembrane helix would fold upon the NS4A contact with the membrane. In addition, partially interactive PEGs help the NS4A N-tail to detach from the protease surface, thus enabling the process of helix insertion and potentially helping the virus establish a replication machinery needed to produce new viruses. Results point to an active role of crowding in assisting structural changes in disordered protein fragments that are necessary for their biological function.

Nuclear magnetic resonance free ligand conformations and atomic resolution dynamics

Knowledge of free ligand conformational preferences (energy minima) and conformational dynamics (rotational energy barriers) of small molecules in solution can guide drug design hypotheses and help rank ideas to bias syntheses towards more active compounds. Visualization of conformational exchange dynamics around torsion angles, by replica exchange with solute tempering molecular dynamics (REST-MD), gives results in agreement with high-resolution 1 H nuclear magnetic resonance (NMR) spectra and complements free ligand conformational analyses. Rotational energy barriers around individual bonds are comparable between calculated and experimental values, making the in-silico method relevant to ranking prospective design ideas in drug discovery programs, particularly across a series of analogs. Prioritizing design ideas, based on calculations and analysis of measurements across a series, efficiently guides rational discovery towards the "right molecules" for effective medicines.

Drug Design Strategies to Avoid Resistance in Direct-Acting Antivirals and Beyond

Drug resistance is prevalent across many diseases, rendering therapies ineffective with severe financial and health consequences. Rather than accepting resistance after the fact, proactive strategies need to be incorporated into the drug design and development process to minimize the impact of drug resistance. These strategies can be derived from our experience with viral disease targets where multiple generations of drugs had to be developed to combat resistance and avoid antiviral failure. Significant efforts including experimental and computational structural biology, medicinal chemistry, and machine learning have focused on understanding the mechanisms and structural basis of resistance against direct-acting antiviral (DAA) drugs. Integrated methods show promise for being predictive of resistance and potency. In this review, we give an overview of this research for human immunodeficiency virus type 1, hepatitis C virus, and influenza virus and the lessons learned from resistance mechanisms of DAAs. These lessons translate into rational strategies to avoid resistance in drug design, which can be generalized and applied beyond viral targets. While resistance may not be completely avoidable, rational drug design can and should incorporate strategies at the outset of drug development to decrease the prevalence of drug resistance.

Employing complementary spectroscopies to study the conformations of an epimeric pair of side-chain stapled peptides in aqueous solution

Understanding the conformational preferences of free ligands in solution is often necessary to rationalize structure-activity relationships in drug discovery. Herein, we examine the conformational behavior of an epimeric pair of side-chain stapled peptides that inhibit the FAD dependent amine oxidase lysine specific demethylase 1 (LSD1). The peptides differ only at a single stereocenter, but display a major difference in binding affinity. Their Raman optical activity (ROA) spectra are most likely dominated by the C-terminus, obscuring the analysis of the epimeric macrocycle. By employing NMR spectroscopy, we show a difference in conformational behavior between the two compounds and that the LSD1 bound conformation of the most potent compound is present to a measurable extent in aqueous solution. In addition, we illustrate that Molecular Dynamics (MD) simulations produce ensembles that include the most important solution conformations, but that it remains problematic to identify relevant conformations with no a priori knowledge from the large conformational pool. Furthermore, this work highlights the importance of understanding the scope and limitations of the available techniques for conducting conformational analyses. It also emphasizes the importance of conformational selection of a flexible ligand in molecular recognition.

NMR spectroscopy: the swiss army knife of drug discovery

Nuclear magnetic resonance (NMR) spectroscopy has evolved into a powerful tool within drug discovery over the last two decades. While traditionally being used by medicinal chemists for small molecule structure elucidation, it can also be a valuable tool for the identification of small molecules that bind to drug targets, for the characterization of target-ligand interactions and for hit-to-lead optimization. Here, we describe how NMR spectroscopy is integrated into the Pfizer drug discovery pipeline and how we utilize this approach to identify and validate initial hits and generate leads.

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.

Avoiding Drug Resistance by Substrate Envelope-Guided Design: Toward Potent and Robust HCV NS3/4A Protease Inhibitors

Hepatitis C virus (HCV) infects millions of people worldwide, causing chronic liver disease that can lead to cirrhosis, hepatocellular carcinoma, and liver transplant. In the last several years, the advent of direct-acting antivirals, including NS3/4A protease inhibitors (PIs), has remarkably improved treatment outcomes of HCV-infected patients. However, selection of resistance-associated substitutions and polymorphisms among genotypes can lead to drug resistance and in some cases treatment failure. A proactive strategy to combat resistance is to constrain PIs within evolutionarily conserved regions in the protease active site. Designing PIs using the substrate envelope is a rational strategy to decrease the susceptibility to resistance by using the constraints of substrate recognition. We successfully designed two series of HCV NS3/4A PIs to leverage unexploited areas in the substrate envelope to improve potency, specifically against resistance-associated substitutions at D168. Our design strategy achieved better resistance profiles over both the FDA-approved NS3/4A PI grazoprevir and the parent compound against the clinically relevant D168A substitution. Crystallographic structural analysis and inhibition assays confirmed that optimally filling the substrate envelope is critical to improve inhibitor potency while avoiding resistance. Specifically, inhibitors that enhanced hydrophobic packing in the S4 pocket and avoided an energetically frustrated pocket performed the best. Thus, the HCV substrate envelope proved to be a powerful tool to design robust PIs, offering a strategy that can be translated to other targets for rational design of inhibitors with improved potency and resistance profiles.IMPORTANCE Despite significant progress, hepatitis C virus (HCV) continues to be a major health problem with millions of people infected worldwide and thousands dying annually due to resulting complications. Recent antiviral combinations can achieve >95% cure, but late diagnosis, low access to treatment, and treatment failure due to drug resistance continue to be roadblocks against eradication of the virus. We report the rational design of two series of HCV NS3/4A protease inhibitors with improved resistance profiles by exploiting evolutionarily constrained regions of the active site using the substrate envelope model. Optimally filling the S4 pocket is critical to avoid resistance and improve potency. Our results provide drug design strategies to avoid resistance that are applicable to other quickly evolving viral drug targets.

Towards discovery of novel scaffold with potent antiangiogenic activity; design, synthesis of pyridazine based compounds, impact of hinge interaction, and accessibility of their bioactive conformation on VEGFR-2 activities

Pyridazine scaffolds are considered privileged structures pertaining to its novelty, chemical stability, and synthetic feasibility. In our quest towards the development of novel scaffolds for effective vascular endothelial growth 2 (VEGFR-2) inhibition with antiangiogenic activity, four novel series of pyridazines were designed and synthesised. Five of the synthesised compounds; namely (8c, 8f, 15, 18b, and 18c) exhibited potent VEGFR-2 inhibitory potency (>80%); with IC₅₀ values ranging from low micromolar to nanomolar range; namely compounds 8c, 8f, 15, 18c with (1.8 µM, 1.3 µM, 1.4 µM, 107 nM), respectively. Moreover, 3-[4-{(6-oxo-1,6-dihydropyridazin-3-yl)oxy}phenyl]urea derivative (18b) exhibited nanomolar potency towards VEGFR-2 (60.7 nM). In cellular assay, the above compounds showed excellent inhibition of VEGF-stimulated proliferation of human umbilical vein endothelial cells at 10 μM concentration. Finally, an extensive molecular simulation study was performed to investigate the probable interaction with VEGFR-2.

Free Ligand 1D NMR Conformational Signatures To Enhance Structure Based Drug Design of a Mcl-1 Inhibitor (AZD5991) and Other Synthetic Macrocycles

The three-dimensional conformations adopted by a free ligand in solution impact bioactivity and physicochemical properties. Solution 1D NMR spectra inherently contain information on ligand conformational flexibility and three-dimensional shape, as well as the propensity of the free ligand to fully preorganize into the bioactive conformation. Herein we discuss some key learnings, distilled from our experience developing potent and selective synthetic macrocyclic inhibitors, including Mcl-1 clinical candidate AZD5991. Case studies have been selected from recent oncology research projects, demonstrating how 1D NMR conformational signatures can complement X-ray protein-ligand structural information to guide medicinal chemistry optimization. Learning to extract free ligand conformational information from routinely available 1D NMR signatures has proven to be fast enough to guide medicinal chemistry decisions within design cycles for compound optimization.

Exploring small molecules with pan-genotypic inhibitory activities against hepatitis C virus NS3/4A serine protease

Among the many Hepatitis C virus (HCV) genotypes and subtypes, genotypes 1b and 3a are most prevalent in United States and Asia, respectively. A total of 132 commercially available analogs of a previous lead compound were initially investigated against wild-type HCV genotype 1b NS3/4A protease. Ten compounds showed inhibitory activities (IC₅₀ values) below 10 µM with comparable direct binding affinities (K_D values) determined by surface plasmon resonance (SPR). To identify pan-genotypic inhibitors, these ten selected compounds were tested against four additional genotypes (1a, 2a, 3a, and 4) and three drug-resistant mutants (A156S, R155K, and V36M). Four new analogs have been identified with better activities against all five tested genotypes than the prior lead compound. Further, the original lead compound did not show activity against genotype 3a NS3/4A, whereas four newly identified compounds exhibited IC₅₀ values below 33 µM against genotype 3a NS3/4A. Encouragingly, the best new compound F1813-0710 possessed promising activity toward genotype 3a, which is a huge improvement over the previous lead compound that had no effect on genotype 3a. This intriguing observation was further analyzed by molecular docking and molecular dynamics (MD) simulations to understand their different binding interactions, which should benefit future pan-genotypic inhibitor design and drug discovery.

Hepatitis C Virus NS3/4A Protease Inhibitors Incorporating Flexible P2 Quinoxalines Target Drug Resistant Viral Variants

A substrate envelope-guided design strategy is reported for improving the resistance profile of HCV NS3/4A protease inhibitors. Analogues of 5172-mcP1P3 were designed by incorporating diverse quinoxalines at the P2 position that predominantly interact with the invariant catalytic triad of the protease. Exploration of structure-activity relationships showed that inhibitors with small hydrophobic substituents at the 3-position of P2 quinoxaline maintain better potency against drug resistant variants, likely due to reduced interactions with residues in the S2 subsite. In contrast, inhibitors with larger groups at this position were highly susceptible to mutations at Arg155, Ala156, and Asp168. Excitingly, several inhibitors exhibited exceptional potency profiles with EC₅₀ values ≤5 nM against major drug resistant HCV variants. These findings support that inhibitors designed to interact with evolutionarily constrained regions of the protease, while avoiding interactions with residues not essential for substrate recognition, are less likely to be susceptible to drug resistance.

Searching for bioactive conformations of drug-like ligands with current force fields: how good are we?

Drug-like ligands obtained from protein-ligand complexes deposited in the Protein Databank were subjected to conformational searching using various force fields and solvation settings. For each ligand, the resulting conformer pool was examined for the presence of the bioactive (crystal pose-like) conformation. Similarity of conformers toward the crystal-pose was quantified as the best achievable root mean squared deviation (RMSD, heavy atoms only). Analyzing the conformer pools generated by various force fields revealed only small differences in the likelihood of finding a crystal pose-like conformation. However, employing different solvents in the conformational search was found to be very important for achieving RMSDs below 1.0 Å. The best statistical values of likelihood were observed with a recently released force field covering a large portion of dihedral angles occurring in drug-like compounds in combination with the water as solvent. In order to enable computational chemists and modelers to efficiently use available software tools, we have additionally performed several focused analyses on ligands, grouped according to descriptors most relevant for the rational drug design.

Small molecule inhibitors of mesotrypsin from a structure-based docking screen

PRSS3/mesotrypsin is an atypical isoform of trypsin, the upregulation of which has been implicated in promoting tumor progression. To date there are no mesotrypsin-selective pharmacological inhibitors which could serve as tools for deciphering the pathological role of this enzyme, and could potentially form the basis for novel therapeutic strategies targeting mesotrypsin. A virtual screen of the Natural Product Database (NPD) and Food and Drug Administration (FDA) approved Drug Database was conducted by high-throughput molecular docking utilizing crystal structures of mesotrypsin. Twelve high-scoring compounds were selected for testing based on lowest free energy docking scores, interaction with key mesotrypsin active site residues, and commercial availability. Diminazene (CID22956468), along with two similar compounds presenting the bis-benzamidine substructure, was validated as a competitive inhibitor of mesotrypsin and other human trypsin isoforms. Diminazene is the most potent small molecule inhibitor of mesotrypsin reported to date with an inhibitory constant (Ki) of 3.6±0.3 μM. Diminazene was subsequently co-crystalized with mesotrypsin and the crystal structure was solved and refined to 1.25 Å resolution. This high resolution crystal structure can now offer a foundation for structure-guided efforts to develop novel and potentially more selective mesotrypsin inhibitors based on similar molecular substructures.

Validation and correction of Zn-CysxHisy complexes

Many crystal structures in the Protein Data Bank contain zinc ions in a geometrically distorted tetrahedral complex with four Cys and/or His ligands. A method is presented to automatically validate and correct these zinc complexes. Analysis of the corrected zinc complexes shows that the average Zn-Cys distances and Cys-Zn-Cys angles are a function of the number of cysteines and histidines involved. The observed trends can be used to develop more context-sensitive targets for model validation and refinement.

Unusual Amino Acids in Medicinal Chemistry

Unusual amino acids are fundamental building blocks of modern medicinal chemistry. The combination of readily functionalized amine and carboxyl groups attached to a chiral central core along with one or two potentially diverse side chains provides a unique three-dimensional structure with a high degree of functionality. This makes them invaluable as starting materials for syntheses of complex molecules, highly diverse elements for SAR campaigns, integral components of peptidomimetic drugs, and potential drugs on their own. This Perspective highlights the diversity of unnatural amino acid structures found in hit-to-lead and lead optimization campaigns and clinical stage and approved drugs, reflecting their increasingly important role in medicinal chemistry.

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

PepComposer: computational design of peptides binding to a given protein surface

There is a wide interest in designing peptides able to bind to a specific region of a protein with the aim of interfering with a known interaction or as starting point for the design of inhibitors. Here we describe PepComposer, a new pipeline for the computational design of peptides binding to a given protein surface. PepComposer only requires the target protein structure and an approximate definition of the binding site as input. We first retrieve a set of peptide backbone scaffolds from monomeric proteins that harbor the same backbone arrangement as the binding site of the protein of interest. Next, we design optimal sequences for the identified peptide scaffolds. The method is fully automatic and available as a web server at http://biocomputing.it/pepcomposer/webserver.