Crystal structures of inhibitor complexes of human T-cell leukemia virus (HTLV-1) protease

Novel computational and drug design strategies for the inhibition of human T-cell leukemia virus 1-associated lymphoma by Astilbin derivatives

Human T-cell leukemia virus 1 (HTLV-1) associated lymphoma is a devastating malignancy triggered by HTLV-1 infections. We employeda comprehensive drug design and computational strategy in this work to explore the inhibitory activitiesof Astilbin derivatives against HTLV-1-associated lymphoma. We evaluated the stability, binding affinities, and various computational analysis of Astilbin derivatives against target proteins, such as HTLV-1 main protease and HTLV-1 capsid protein. The root mean square deviation (RMSD), root mean square fluctuation, radius of gyration, hydrogen bond analysis, principal component analysis (PCA) and dynamic cross-correlation matrix (DCCM) were applied to characterize these protein-ligand interactions further. Ligand-03 and ligand-04 exhibited notable binding affinity to HTLV-1 capsid protein, while ligand-05 displayed high binding affinity to HTLV-1 protease. MD simulation analysis revealed that ligand-03, bound to HTLV-1 capsid protein, demonstrated enhanced stability with lower RMSD values and fewer conformational changes, suggesting a promising binding orientation. Ligand-04, despite stable binding, exhibited increased structural deviations, making it less suitable. Ligand-05 demonstrated stable binding to HTLV-1 protease throughout the simulation period at 100 nanoseconds. Hydrogen bond analysis indicated that ligand-05 formed persistent hydrogen bonds with significantresidues, contributing to its stability. PCA highlighted ligand-03's more remarkable conformational changes, while DCCM showed ligand-05's distinct dynamics, indicating its different behavior in the complex. Furthermore, binding free energy calculations supported the favorable interactions of ligand-03 and ligand-04 with HTLV-1 capsid protein, while ligand-05 showed weaker interactions with HTLV-1 protease. Molecular electrostatic potential and frontier molecular orbital analyses provided insights into these compounds' charge distribution and stability. In conclusion, this research found Astilbin derivatives as potential inhibitors of HTLV-1-associated lymphoma. Future attempts at drug development will benefit from the steady interaction landscape provided by Ligand-03, Ligand-04 and Ligand-05, which showed the most attractive binding profile with the target protein. These results open up new opportunities for innovative drug development, and more experimental testing should be done between Astilbin derivatives and HTLV-1-associated lymphoma.Communicated by Ramaswamy H. Sarma.

Catalytic Mechanism of Human T-Cell Leukemia Virus Type 1 Protease Investigated by Combined QM/MM Molecular Dynamics Simulations

Combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulations were performed to investigate the catalytic mechanism of human T-cell leukemia virus type 1 (HTLV-1) protease, a retroviral aspartic protease that is a potential therapeutic target for curing HTLV-1-associated diseases. To elucidate the proteolytic cleavage mechanism, we determined the two-dimensional free energy surfaces of the HTLV-1 protease-catalyzed reactions through various possible pathways. The free energy simulations suggest that the catalytic reactions of the HTLV-1 protease occur in the following sequential steps: (1) a proton is transferred from the lytic water to Asp32', followed by the nucleophilic addition of the resulting hydroxyl to the carbonyl carbon of the scissile bond, forming a tetrahedral oxyanion intermediate, and (2) a proton is transferred from Asp32 to the peptide nitrogen of the scissile bond, leading to the spontaneous breakage of the scissile bond. The rate-limiting step of this catalytic process is the proton transfer from Asp32 to the peptide nitrogen of the scissile bond, with a free energy of activation of 21.1 kcal/mol. This free energy barrier is close to the experimentally determined free energy of activation (16.3 kcal/mol) calculated from the measured catalytic rate constant (k_cat). This mechanistic study provides detailed dynamic and structural information that will facilitate the design of mechanism-based inhibitors for the treatment of HTLV-1-associated diseases.

Comparative study of the unbinding process of some HTLV-1 protease inhibitors using unbiased molecular dynamics simulations

The HTLV-1 protease is one of the major antiviral targets to overwhelm this virus. Several research groups have developed protease inhibitors, but none has been successful. In this regard, developing new HTLV-1 protease inhibitors to fix the defects in previous inhibitors may overcome the lack of curative treatment for this oncovirus. Thus, we decided to study the unbinding pathways of the most potent (compound 10, PDB ID 4YDF, Ki = 15 nM) and one of the weakest (compound 9, PDB ID 4YDG, Ki = 7900 nM) protease inhibitors, which are very structurally similar. We conducted 12 successful short and long simulations (totaling 14.8 μs) to unbind the compounds from two monoprotonated (mp) forms of protease using the Supervised Molecular Dynamics (SuMD) without applying any biasing force. The results revealed that Asp32 or Asp32′ in the two forms of mp state similarly exert powerful effects on maintaining both potent and weak inhibitors in the binding pocket of HTLV-1 protease. In the potent inhibitor’s unbinding process, His66′ was a great supporter that was absent in the weak inhibitor’s unbinding pathway. In contrast, in the weak inhibitor’s unbinding process, Trp98/Trp98′ by pi-pi stacking interactions were unfavorable for the stability of the inhibitor in the binding site. In our opinion, these results will assist in designing more potent and effective inhibitors for the HTLV-1 protease.

The Effects of Side-Chain Configurations of a Retro-Inverso-Type Inhibitor on the Human T-Cell Leukemia Virus (HTLV)-1 Protease

In this study, the effects of side-chain configurations of D-Ile residues of a retro-inverso (RI)-type inhibitor on the human T-cell leukemia virus type 1 (HTLV-1) protease containing a hydroxyethylamine dipeptide isostere were clarified. Prior to evaluation using the RI-type inhibitor, the effects of side-chain configurations of Ile residues of the substrate peptide on the HTLV-1 protease were examined to estimate the influence of side-chain configurations on enzyme activity. Based on the estimation of the influence of side-chain configurations on protease affinity, the RI-type inhibitors containing a D-allo-Ile residue in the corresponding substrate sequence, instead of a D-Ile residue, were synthesized via 9-fluorenylmethoxycarbonyl-based solid-phase peptide synthesis. Refolded recombinant HTLV-1 protease (1-116, L40I) was used for the simple and short evaluation of the inhibitory activities of the synthesized RI-type inhibitors. The results clearly indicated that mimicking the whole topology, comprising both the main- and side-chain structures of the parent inhibitor, is effective for the design of potent RI-modified protease inhibitors.

Viral proteases: Structure, mechanism and inhibition

Viral proteases are diverse in structure, oligomeric state, catalytic mechanism, and substrate specificity. This chapter focuses on proteases from viruses that are relevant to human health: human immunodeficiency virus subtype 1 (HIV-1), hepatitis C (HCV), human T-cell leukemia virus type 1 (HTLV-1), flaviviruses, enteroviruses, and coronaviruses. The proteases of HIV-1 and HCV have been successfully targeted for therapeutics, with picomolar FDA-approved drugs currently used in the clinic. The proteases of HTLV-1 and the other virus families remain emerging therapeutic targets at different stages of the drug development process. This chapter provides an overview of the current knowledge on viral protease structure, mechanism, substrate recognition, and inhibition. Particular focus is placed on recent advances in understanding the molecular basis of diverse substrate recognition and resistance, which is essential toward designing novel protease inhibitors as antivirals.

Combining Molecular Dynamic Information and an Aspherical-Atom Data Bank in the Evaluation of the Electrostatic Interaction Energy in Multimeric Protein-Ligand Complex: A Case Study for HIV-1 Protease

Computational analysis of protein–ligand interactions is of crucial importance for drug discovery. Assessment of ligand binding energy allows us to have a glimpse of the potential of a small organic molecule to be a ligand to the binding site of a protein target. Available scoring functions, such as in docking programs, all rely on equations that sum each type of protein–ligand interactions in order to predict the binding affinity. Most of the scoring functions consider electrostatic interactions involving the protein and the ligand. Electrostatic interactions constitute one of the most important part of total interactions between macromolecules. Unlike dispersion forces, they are highly directional and therefore dominate the nature of molecular packing in crystals and in biological complexes and contribute significantly to differences in inhibition strength among related enzyme inhibitors. In this study, complexes of HIV-1 protease with inhibitor molecules (JE-2147 and darunavir) were analyzed by using charge densities from the transferable aspherical-atom University at Buffalo Databank (UBDB). Moreover, we analyzed the electrostatic interaction energy for an ensemble of structures, using molecular dynamic simulations to highlight the main features of electrostatic interactions important for binding affinity.

Inhibiting HTLV-1 Protease: A Viable Antiviral Target

Human T-cell lymphotropic virus type 1 (HTLV-1) is a retrovirus that can cause severe paralytic neurologic disease and immune disorders as well as cancer. An estimated 20 million people worldwide are infected with HTLV-1, with prevalence reaching 30% in some parts of the world. In stark contrast to HIV-1, no direct acting antivirals (DAAs) exist against HTLV-1. The aspartyl protease of HTLV-1 is a dimer similar to that of HIV-1 and processes the viral polyprotein to permit viral maturation. We report that the FDA-approved HIV-1 protease inhibitor darunavir (DRV) inhibits the enzyme with 0.8 μM potency and provides a scaffold for drug design against HTLV-1. Analogs of DRV that we designed and synthesized achieved submicromolar inhibition against HTLV-1 protease and inhibited Gag processing in viral maturation assays and in a chronically HTLV-1 infected cell line. Cocrystal structures of these inhibitors with HTLV-1 protease highlight opportunities for future inhibitor design. Our results show promise toward developing highly potent HTLV-1 protease inhibitors as therapeutic agents against HTLV-1 infections.

Biochemical Characterization, Specificity and Inhibition Studies of HTLV-1, HTLV-2, and HTLV-3 Proteases

The human T-lymphotropic viruses (HTLVs) are causative agents of severe diseases including adult T-cell leukemia. Similar to human immunodeficiency viruses (HIVs), the viral protease (PR) plays a crucial role in the viral life-cycle via the processing of the viral polyproteins. Thus, it is a potential target of anti-retroviral therapies. In this study, we performed in vitro comparative analysis of human T-cell leukemia virus type 1, 2, and 3 (HTLV-1, -2, and -3) proteases. Amino acid preferences of S4 to S1′ subsites were studied by using a series of synthetic oligopeptide substrates representing the natural and modified cleavage site sequences of the proteases. Biochemical characteristics of the different PRs were also determined, including catalytic efficiencies and dependence of activity on pH, temperature, and ionic strength. We investigated the effects of different HIV-1 PR inhibitors (atazanavir, darunavir, DMP-323, indinavir, ritonavir, and saquinavir) on enzyme activities, and inhibitory potentials of IB-268 and IB-269 inhibitors that were previously designed against HTLV-1 PR. Comparative biochemical analysis of HTLV-1, -2, and -3 PRs may help understand the characteristic similarities and differences between these enzymes in order to estimate the potential of the appearance of drug-resistance against specific HTLV-1 PR inhibitors.

Dimer Interface Organization is a Main Determinant of Intermonomeric Interactions and Correlates with Evolutionary Relationships of Retroviral and Retroviral-Like Ddi1 and Ddi2 Proteases

The life cycles of retroviruses rely on the limited proteolysis catalyzed by the viral protease. Numerous eukaryotic organisms also express endogenously such proteases, which originate from retrotransposons or retroviruses, including DNA damage-inducible 1 and 2 (Ddi1 and Ddi2, respectively) proteins. In this study, we performed a comparative analysis based on the structural data currently available in Protein Data Bank (PDB) and Structural summaries of PDB entries (PDBsum) databases, with a special emphasis on the regions involved in dimerization of retroviral and retroviral-like Ddi proteases. In addition to Ddi1 and Ddi2, at least one member of all seven genera of the Retroviridae family was included in this comparison. We found that the studied retroviral and non-viral proteases show differences in the mode of dimerization and density of intermonomeric contacts, and distribution of the structural characteristics is in agreement with their evolutionary relationships. Multiple sequence and structure alignments revealed that the interactions between the subunits depend mainly on the overall organization of the dimer interface. We think that better understanding of the general and specific features of proteases may support the characterization of retroviral-like proteases.

Comparison of a retroviral protease in monomeric and dimeric states

Retroviral proteases (RPs) are of high interest owing to their crucial role in the maturation process of retroviral particles. RPs are obligatory homodimers, with a pepsin-like active site built around two aspartates (in DTG triads) that activate a water molecule, as the nucleophile, under two flap loops. Mason-Pfizer monkey virus (M-PMV) is unique among retroviruses as its protease is also stable in the monomeric form, as confirmed by an existing crystal structure of a 13 kDa variant of the protein (M-PMV PR) and its previous biochemical characterization. In the present work, two mutants of M-PMV PR, D26N and C7A/D26N/C106A, were crystallized in complex with a peptidomimetic inhibitor and one mutant (D26N) was crystallized without the inhibitor. The crystal structures were solved at resolutions of 1.6, 1.9 and 2.0 Å, respectively. At variance with the previous study, all of the new structures have the canonical dimeric form of retroviral proteases. The protomers within a dimer differ mainly in the flap-loop region, with the most extreme case observed in the apo structure, in which one flap loop is well defined while the other flap loop is not defined by electron density. The presence of the inhibitor molecules in the complex structures was assessed using polder maps, but some details of their conformations remain ambiguous. In all of the presented structures the active site contains a water molecule buried deeply between the Asn26-Thr27-Gly28 triads of the protomers. Such a water molecule is completely unique not only in retropepsins but also in aspartic proteases in general. The C7A and C106A mutations do not influence the conformation of the protein. The Cys106 residue is properly placed at the homodimer interface area for a disulfide cross-link, but the reducing conditions of the crystallization experiment prevented S-S bond formation. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:Acta_Cryst_D:S2059798319011355.

New directions for protease inhibitors directed drug discovery

Proteases play crucial roles in various biological processes, and their activities are essential for all living organisms-from viruses to humans. Since their functions are closely associated with many pathogenic mechanisms, their inhibitors or activators are important molecular targets for developing treatments for various diseases. Here, we describe drugs/drug candidates that target proteases, such as malarial plasmepsins, β-secretase, virus proteases, and dipeptidyl peptidase-4. Previously, we reported inhibitors of aspartic proteases, such as renin, human immunodeficiency virus type 1 protease, human T-lymphotropic virus type I protease, plasmepsins, and β-secretase, as drug candidates for hypertension, adult T-cell leukaemia, human T-lymphotropic virus type I-associated myelopathy, malaria, and Alzheimer's disease. Our inhibitors are also described in this review article as examples of drugs that target proteases. © 2015 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 106: 563-579, 2016.

Design of new potent HTLV-1 protease inhibitors: in silico study

HTLV-1 and HIV-1 are two major causes for severe T-cell leukemia disease and acquired immune deficiency syndrome (AIDS). HTLV-1 protease, a member of aspartic acid protease family, plays important roles in maturation during virus replication cycle. The impairment of these proteases results in uninfectious HTLV-1virions.Similar to HIV-1protease deliberate mutations that confer drug resistance on HTLV-1 are frequently seen in this protease. Therefore, inhibition of HTLV-1 protease activity is expected to disrupt HTLV-1's ability to replicate and infect additional cells. In this study, we initially designed fifteen inhibitory compounds based on the conformations of a class of HIV-1 aspartyl protease inhibitors, sulfonamid-peptoid. Five compounds were chosen based on the goodness of their Drug-Likeness scoreusing "Lipinsk's rule of five". Here, using protein-ligand docking approach we compared the inhibitory constants of these compounds to those available in literatures and observed significantly higher inhibition for two compounds, SP-4 and SP-5. Our data suggest that the addition of two cyclic hydrocarbons to both ends of sulfonamide peptoids leads to the formation of new hydrophobic interactions due to the semi-circular form of these compounds, connecting the first chain of protease to the two ends of tested ligands via Hydrophobic interactions. We conclude that hydrophobic force plays an important role in suppressing protease activity especially for HTLV-1 protease, which in turn prevents the virus maturity. Therefore, designing and development of new ligands based on aromatic hydrocarbons in both ends of inhibitors is very promising for efficient treatment.

Characterizing the protonation states of the catalytic residues in apo and substrate-bound human T-cell leukemia virus type 1 protease

Human T-cell leukemia virus type 1 (HTLV-1) protease is an attractive target when developing inhibitors to treat HTLV-1 associated diseases. To study the catalytic mechanism and design novel HTLV-1 protease inhibitors, the protonation states of the two catalytic aspartic acid residues must be determined. Free energy simulations have been conducted to study the proton transfer reaction between the catalytic residues of HTLV-1 protease using a combined quantum mechanical and molecular mechanical (QM/MM) molecular dynamics simulation. The free energy profiles for the reaction in the apo-enzyme and in an enzyme - substrate complex have been obtained. In the apo-enzyme, the two catalytic residues are chemically equivalent and are expected to be both unprotonated. Upon substrate binding, the catalytic residues of HTLV-1 protease evolve to a singly protonated state, in which the OD1 of Asp32 is protonated and forms a hydrogen bond with the OD1 of Asp32', which is unprotonated. The HTLV-1 protease-substrate complex structure obtained from this simulation can serve as the Michaelis complex structure for further mechanistic studies of HTLV-1 protease while providing a receptor structure with the correct protonation states for the active site residues toward the design of novel HTLV-1 protease inhibitors through virtual screening.

Structural basis for HTLV-1 protease inhibition by the HIV-1 protease inhibitor indinavir

HTLV-1 protease (HTLV-1 PR) is an aspartic protease which represents a promising drug target for the discovery of novel anti-HTLV-1 drugs. The X-ray structure of HTLV-1 PR in complex with the well-known and approved HIV-1 PR inhibitor Indinavir was determined at 2.40 Å resolution. In this contribution, we describe the first crystal structure in complex with a nonpeptidic inhibitor that accounts for rationalizing the rather moderate affinity of Indinavir against HTLV-1 PR and provides the basis for further structure-guided optimization strategies.

The application of bioisosteres in drug design for novel drug discovery: focusing on acid protease inhibitors

INTRODUCTION

A bioisostere is a powerful concept for medicinal chemistry. It allows the improvement of the stability; oral absorption; membrane permeability; and absorption, distribution, metabolism and excretion (ADME) of drug candidate, while retaining their biological properties. The term 'bioisostere' is derived from 'isostere', whose physical and chemical properties, such as steric size, hydrophobicity, and electronegativity, are similar to those of a functional or atomic group, and is considered to possess biological properties. Here, the authors highlight the recent applications of bioisosteres in drug design, mainly based on our drug discovery studies.

AREAS COVERED

This review discusses the application of bioisosteres for novel drug discovery with focus on the authors' drug discovery studies such as renin, HIV-protease, and β-secretase inhibitors. The authors highlight that some bioisosteres can form the scaffolding for drug candidates, namely substrate transition state, amide/ester, and carboxylic acid bioisosteres. Moreover, the authors propose the new terms 'electron-donor bioisostere' and 'conformational bioisostere' for drug discovery.

EXPERT OPINION

The authors discuss the importance of bioisostere's design concept based on specific interaction with the corresponding biomolecule. In addition, some strategies for drug discovery based on the bioisostere concept are introduced. Many bioisosteres, which are recognized by corresponding target biomolecules as exhibiting similar biological properties, have been reported to date; most of the recently developed bioisosteres were designed by cheminformatics approaches. Some molecular design softwares and databases are introduced.

Design and synthesis of several small-size HTLV-I protease inhibitors with different hydrophilicity profiles

The human T cell leukemia/lymphotropic virus type 1 (HTLV-I) is clinically associated with adult T cell leukemia/lymphoma, HTLV-I associated myelopathy/tropical spastic paraparesis, and a number of other chronic inflammatory diseases. To stop the replication of the virus, we developed highly potent tetrapeptidic HTLV-I protease inhibitors. In a recent X-ray crystallography study, several of our inhibitors could not form co-crystal complexes with the protease due to their high hydrophobicity. In the current study, we designed, synthesized and evaluated the HTLV-I protease inhibition potency of compounds with hydrophilic end-capping moieties with the aim of improving pharmaceutic and pharmacokinetic properties.

Maintaining potent HTLV-I protease inhibition without the P3-cap moiety in small tetrapeptidic inhibitors

The human T cell lymphotropic/leukemia virus type 1 (HTLV-I) causes adult T cell lymphoma/leukemia. The virus is also responsible for chronic progressive myelopathy and several inflammatory diseases. To stop the manufacturing of new viral components, in our previous reports, we derived small tetrapeptidic HTLV-I protease inhibitors with an important amide-capping moiety at the P(3) residue. In the current study, we removed the P(3)-cap moiety and, with great difficulty, optimized the P(3) residue for HTLV-I protease inhibition potency. We discovered a very potent and small tetrapeptidic HTLV-I protease inhibitor (KNI-10774a, IC(50)=13 nM).