A Comprehensive Computational Analysis for the Binding Modes of Hepatitis C Virus NS5A Inhibitors: The Question of Symmetry

Leveraging structural and 2D-QSAR to investigate the role of functional group substitutions, conserved surface residues and desolvation in triggering the small molecule-induced dimerization of hPD-L1

Small molecules are rising as a new generation of immune checkpoints’ inhibitors, with compounds targeting the human Programmed death-ligand 1 (hPD-L1) protein are pioneering this area of research. Promising examples include the recently disclosed compounds from Bristol-Myers-Squibb (BMS). These molecules bind specifically to hPD-L1 through a unique mode of action. They induce dimerization between two hPD-L1 monomers through the hPD-1 binding interface in each monomer, thereby inhibiting the PD-1/PD-L1 axis. While the recently reported crystal structures of such small molecules bound to hPD-L1 reveal valuable insights regarding their molecular interactions, there is still limited information about the dynamics driving this unusual complex formation. The current study provides an in-depth computational structural analysis to study the interactions of five small molecule compounds in complex with hPD-L1. By employing a combination of molecular dynamic simulations, binding energy calculations and computational solvent mapping techniques, our analyses quantified the dynamic roles of different hydrophilic and lipophilic residues at the surface of hPD-L1 in mediating these interactions. Furthermore, ligand-based analyses, including Free-Wilson 2D-QSAR was conducted to quantify the impact of R-group substitutions at different sites of the phenoxy-methyl biphenyl core. Our results emphasize the importance of a terminal phenyl ring that must be present in any hPD-L1 small molecule inhibitor. This phenyl moiety overlaps with a very unfavorable hydration site, which can explain the ability of such small molecules to trigger hPD-L1 dimerization.

Redesigning of the cap conformation and symmetry of the diphenylethyne core to yield highly potent pan-genotypic NS5A inhibitors with high potency and high resistance barrier

Herein, we report the discovery of several NS5A inhibitors with potency against HCV genotype 1b in the picomolar range. Compounds (15, 33) were of extremely high potency against HCV genotype 1b (EC₅₀ ≈ 1 pM), improved activity against genotype 3a (GT 3a) and good metabolic stability. We studied the impact of changing the cap conformation relative to the diphenylethyne core and/or compound symmetry on both potency and metabolic stability. The analogs obtained exhibited improved potency against HCV genotypes 1a, 1b, 3a and 4a compared to the clinically approved candidate daclatasvir with EC₅₀ values in the low picomolar range and SI₅₀s > 7 orders of magnitude. Compound 15, a symmetrically m-, m'-substituted diphenyl ethyne analog, was 150-fold more potent than daclatasvir against GT 3a, while compound 33, an asymmetrically m-, p-substituted diphenyl ethyne analog, was 35-fold more potent than daclatasvir against GT 3a. In addition, compound 15 exhibited a higher resistance barrier than daclatasvir against genotype 1b.

Serinol-Based Benzoic Acid Esters as New Scaffolds for the Development of Adenovirus Infection Inhibitors: Design, Synthesis, and In Vitro Biological Evaluation

Over the years, human adenovirus (HAdV) has progressively been recognized as a significant viral pathogen. Traditionally associated with self-limited respiratory, gastrointestinal, and conjunctival infections, mainly in immunocompromised patients, HAdV is currently considered to be a pathogen presenting significant morbidity and mortality in both immunosuppressed and otherwise healthy individuals. Currently available therapeutic options are limited because of their lack of effectivity and related side effects. In this context, there is an urgent need to develop effective anti-HAdV drugs with suitable therapeutic indexes. In this work, we identified new serinol-derived benzoic acid esters as novel scaffolds for the inhibition of HAdV infections. A set of 38 compounds were designed and synthesized, and their antiviral activity and cytotoxicity were evaluated. Four compounds (13, 14, 27, and 32) inhibited HAdV infection at low micromolar concentrations (2.82-5.35 μM). Their half maximal inhibitory concentration (IC₅₀) values were lower compared to that of cidofovir, the current drug of choice. All compounds significantly reduced the HAdV DNA replication process, while they did not block any step of the viral entry. Our results showed that compounds 13, 14, and 32 seem to be targeting the expression of the E1A early gene. Moreover, all four derivatives demonstrated a significant inhibition of human cytomegalovirus (HCMV) DNA replication. This new scaffold may represent a potential tool useful for the development of effective anti-HAdV drugs.

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.

Synthesis, biological evaluation and in silico modeling of novel pan-genotypic NS5A inhibitors

A series of novel small-molecule pan-genotypic hepatitis C virus (HCV) NS5A inhibitors with picomolar activity containing 2-[(2S)-pyrrolidin-2-yl]-5-[4-(4-{2-[(2S)-pyrrolidin-2-yl]-1H-imidazol-5-yl}buta-1,3-diyn-1-yl)phenyl]-1H-imidazole core was designed based on molecular modeling study and SAR analysis. The constructed in silico model and docking study provide a deep insight into the binding mode of this type of NS5A inhibitors. Based on the predicted binding interface we have prioritized the most crucial diversity points responsible for improving antiviral activity. The synthesized molecules were tested in a cell-based assay, and compound 1.12 showed an EC₅₀ value in the range of 2.9-34 pM against six genotypes of NS5A HCV, including gT3a, and demonstrated favorable pharmacokinetic profile in rats. This lead compound can be considered as an attractive candidate for further clinical evaluation.

QSAR studies on hepatitis C virus NS5A protein tetracyclic inhibitors in wild type and mutants by CoMFA and CoMSIA

Several 3D-QSAR models were built based on 196 hepatitis C virus (HCV) NS5A protein inhibitors. The bioactivity values EC₉₀ for three types of inhibitors, the wild type (GT1a) and two mutants (GT1a Y93H and GT1a L31V), were collected to build three datasets. The programs OMEGA and ROCS were used for generating conformations and aligning molecules of the dataset, respectively. Each dataset was randomly divided into a training set and a test set three times to reduce the contingency of only one random selection. QSAR models were computed by comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA). For the datasets GT1a, GT1a Y93H, and GT1a L31V, the best models CoMFA-INDX, CoMSIA-SEHA, and CoMSIA-SEHA showed an r² value of 0.682 ± 0.033, 0.779 ± 0.036, and 0.782 ± 0.022 on the test sets, respectively. From the contour maps of the three best models, we summarized the favourable and unfavourable substituents on the tetracyclic core, the Z group, the proline group, and the valine group of inhibitors. We guessed the mutants could change the electrostatic surfaces of the wild type active pocket. In addition, we used ECFP analyses to find important substructures and could intuitively understand the results from QSAR models.

Structure-based screening and validation of potential dengue virus inhibitors through classical and QM/MM affinity estimation

The recurrent outbreaks of dengue virus around the globe represent a huge challenge for governments and public health organizations. With the rapid growth and ease of transportation, dengue disease continues to spread, placing more of the world population under constant threat. Despite decades of research efforts, no effective small molecule antivirals are available against dengue virus. With the efficacy of the recently developed vaccine to be determined, there is an urgent unmet need for small molecule dengue virus treatments. In the current study, we employed state-of-the-art molecular modelling simulations to identify novel inhibitors of the dengue virus envelope protein. The binding modes of all compounds within the conserved β-OctylGlucoside (β-OG) pocket were studied using a combination of docking, molecular dynamics simulations and binding free energy calculations. Here, we describe ten new compounds that significantly reduce production of dengue virus as determined using standard cell-based virological assays. Moreover, we present a comprehensive structural analysis of the identified hits, focusing on their electrostatic and lipophilic binding energy contributions. Finally, we highlight the effect of the desolvation penalty in limiting the activity of some of these compounds. The data presented here paves the way toward rationally designing selective and potent novel inhibitors against dengue virus.

HA1077 displays synergistic activity with daclatasvir against hepatitis C virus and suppresses the emergence of NS5A resistance-associated substitutions in mice

The kinase C-related kinase 2 (PRK2), which phosphorylates hepatitis C virus (HCV) RNA polymerase, is a proviral factor enhancing HCV replication. Here, we report on the in vivo anti-HCV efficacy of HA1077, which inhibits viral genome replication by targeting PRK2 and displays viral entry inhibitory activity by targeting Rho-associated kinase. HA1077 showed synergistic antiviral activity selectively with nonstructural protein 5 A (NS5A) inhibitors including daclatasvir (DCV). HA1077 oral administration substantially reduced serum viral loads in mice bearing HCV genotype 2a-replicating Huh7 xenografts. When administered with DCV, HA1077 potentiated the antiviral efficacy of DCV and suppressed the generation of DCV resistance-associated variants (RAVs). By deep-sequencing analysis, we uncovered an unprecedented DCV-induced polymorphism at the poly-proline motif (PxxPxxP) of NS5A. Coadministration of HA1077 reduced such a polymorphism. Overall, our results demonstrate the potential therapeutic benefit of combination therapy with HA1077 plus DCV for HCV patients carrying emerging or pre-existing RAVs toward NS5A inhibitors.

Applications of computer-aided approaches in the development of hepatitis C antiviral agents

INTRODUCTION

Hepatitis C virus (HCV) is a global health problem that causes several chronic life-threatening liver diseases. The numbers of people affected by HCV are rising annually. Since 2011, the FDA has approved several anti-HCV drugs; while many other promising HCV drugs are currently in late clinical trials. Areas covered: This review discusses the applications of different computational approaches in HCV drug design. Expert opinion: Molecular docking and virtual screening approaches have emerged as a low-cost tool to screen large databases and identify potential small-molecule hits against HCV targets. Ligand-based approaches are useful for filtering-out compounds with rich physicochemical properties to inhibit HCV targets. Molecular dynamics (MD) remains a useful tool in optimizing the ligand-protein complexes and understand the ligand binding modes and drug resistance mechanisms in HCV. Despite their varied roles, the application of in-silico approaches in HCV drug design is still in its infancy. A more mature application should aim at modelling the whole HCV replicon in its active form and help to identify new effective druggable sites within the replicon system. With more technological advancements, the roles of computer-aided methods are only going to increase several folds in the development of next-generation HCV drugs.