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

Fluorinated derivatives of tetrahydroaltersolanol molecule on COVID-19, HIV, and HTLV protease by DFT and molecular docking approaches

Structural, optoelectronic, and biological properties of tetrahydroaltersolanol (C₁₆H₂₀O₇) and fluorinated derivatives are calculated using density functional theory (DFT) and molecular docking approaches. It is shown that the pure C₁₆H₂₀O₇ molecule has a direct HOMO-LUMO energy gap about 3.1 eV. The substitution of F atom at A category decreases the electronic energy gap, while it is constant at B category. In A category, the behavior of the pure molecule changed from insulator to semiconductor with various substitution of F atom. The electronic properties were depended on the F sites in the pure molecule. The molecular electrical transport properties and charge-transfer possibilities increase with decreasing energy gap. The pure C₁₆H₂₀O₇ molecule with high energy gap has low chemical reactivity and substitution of F atom at considered molecule increases chemical reactivity. Obtained results show that F-O bonds in trifurcation bonds of C₁₆H₁₉O₇(F14), C₁₆H₁₉O₇(F16), and C₁₆H₁₉O₇(F17) molecules play a key role in confronting with COVID-19, HIV, and HTLV proteases, respectively. Optical spectra, such as the dielectric functions, electron energy-loss spectroscopy, refractive index, extinction coefficient, and reflection spectra show that fluorinated derivatives of C₁₆H₂₀O₇ at B category can be used in the new drugs.

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.

Repurposing antiviral drugs against HTLV-1 protease by molecular docking and molecular dynamics simulation

Human T-cell leukemia virus type I (HTLV-1) belongs to the delta retrovirus family and the etiological agent of adult T-cell leukemia (ATL(. While the current HTLV-1 therapy, relies on using Zidovudine plus IFN-γ, there is no FDA approved drugs against it. In silico drug repurposing is a fast and accurate way for screening US-FDA approved drugs to find a therapeutic option for the HTLV-1 infection. So that, this research aims to analyze a dataset of approved antiviral drugs as a potential prospect for an anti-viral drug against HTLV-1 infection. Molecular docking simulation was performed to identify interactions of the antiviral drugs with the key residues in the HTLV-1 protease binding site. Then, molecular dynamics simulation was also performed for the potential protein-ligand complexes to confirm the stable behavior of the ligands inside the binding pocket. The best docking scores with the target was found to be Simeprevir, Atazanavir, and Saquinavir compounds which indicate that these drugs can firmly bind to the HTLV-1 protease. The MD simulation confirmed the stability of Simeprevir-protease, Atazanavir-Protease, and Saquinavir-Protease interactions. Clearly, these compounds should be further evaluated in experimental assays and clinical trials to confirm their actual activity against HTLV-1 infection.Communicated by Ramaswamy H. Sarma.

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.

Reconstruction of the binding pathway of an anti-HIV drug, Indinavir, in complex with the HTLV-1 protease using unaggregated unbiased molecular dynamics simulation

Retroviruses are a growing concern for the health of human beings, and one of the dangerous members of this family is the Human T-cell Leukemia Virus 1 (HTLV-1) virus. It has affected more than 20 million people so far, and since there are no registered treatments against it yet, urgent treatment solutions are needed. One of the most promising drug targets to fight this virus is the protease enzyme of the virus's protein machinery. In this study, by utilizing a computational method called Unaggregated Unbiased Molecular Dynamics (UUMD), we reconstructed the binding pathway of a HTLV-1 protease inhibitor, Indinavir, to find the details of the binding pathway, the influential residues, and also the stable states of the binding pathway. We achieved the native conformation of the inhibitor in 6 rounds, 360 replicas by performing over 4 micro-seconds of UMD simulations. We found 3 Intermediate states between the solvated state and the native conformation state in the binding pathway. We also discovered that aromatic residues such as Trp98 and Trp98', catalytic residues Asp32 and Asp32', and the flap region's residues have the most influential roles in the binding pathway and also have the most contribution to the total interaction energies. We believe that the details found in this study would be a great guide for developing new treatment solutions against the HTLV-1 virus by inhibiting the HTLV-1 protease.

Comparative analysis of the unbinding pathways of antiviral drug Indinavir from HIV and HTLV1 proteases by supervised molecular dynamics simulation

Determining the unbinding pathways of potential small molecule compounds from their target proteins is of great significance for designing efficacious treatment solutions. One of these potential compounds is the approved HIV-1 protease inhibitor, Indinavir, which has a weak effect on the HTLV-1 protease. In this work, by employing the SuMD method, we reconstructed the unbinding pathways of Indinavir from HIV and HTLV-1 proteases to compare and understand the mechanism of the unbinding and to discover the reasons for the lack of inhibitory activity of Indinavir against the HTLV-1 protease. We achieved multiple unbinding events from both HIV and HTLV-1 proteases in which the RMSD values of Indinavir reached over 40 Å. Also, we found that the mobility and fluctuations of the flap region are higher in the HTLV-1 protease, making the drug less stable. We realized that critically positioned aromatic residues such as Trp98/Trp98' and Phe67/Phe67' in the HTLV-1 protease could make strong π-Stacking interactions with Indinavir in the unbinding pathway, which are unfavorable for the stability of Indinavir in the active site. The details found in this study can make a reasonable explanation for the lack of inhibitory activity of this drug against HTLV-1 protease. We believe the details discovered in this work can help design more effective and selective inhibitors for the HTLV-1 protease.

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.

Antiviral Drug Targets of Single-Stranded RNA Viruses Causing Chronic Human Diseases

Ribonucleic acid (RNA) viruses associated with chronic diseases in humans are major threats to public health causing high mortality globally. The high mutation rate of RNA viruses helps them to escape the immune response and also is responsible for the development of drug resistance. Chronic infections caused by human immunodeficiency virus (HIV) and hepatitis viruses (HBV and HCV) lead to acquired immunodeficiency syndrome (AIDS) and hepatocellular carcinoma respectively, which are one of the major causes of human deaths. Effective preventative measures to limit chronic and re-emerging viral infections are absolutely necessary. Each class of antiviral agents targets a specific stage in the viral life cycle and inhibits them from its development and proliferation. Most often, antiviral drugs target a specific viral protein, therefore only a few broad-spectrum drugs are available. This review will be focused on the selected viral target proteins of pathogenic viruses containing single-stranded (ss) RNA genome that causes chronic infections in humans (e.g. HIV, HCV, Flaviviruses). In the recent past, an exponential increase in the number of available three-dimensional protein structures (>150000 in Protein Data Bank), allowed us to better understand the molecular mechanism of action of protein targets and antivirals. Advancements in the in silico approaches paved the way to design and develop several novels, highly specific small-molecule inhibitors targeting the viral proteins.

Portrait of the Intrinsically Disordered Side of the HTLV-1 Proteome

Intrinsically disordered proteins (IDPs) lack an ordered 3D structure. These proteins contain one or more intrinsically disordered protein regions (IDPRs). IDPRs interact promiscuously with other proteins, which leads to their structural transition from a disordered to an ordered state. Such interaction-prone regions of IDPs are known as molecular recognition features. Recent studies suggest that IDPs provide structural plasticity and functional diversity to viral proteins that are involved in rapid replication and immune evasion within the host cells. In the present study, we evaluated the prevalence of IDPs and IDPRs in human T lymphotropic virus type 1 (HTLV-1) proteome. We also investigated the presence of MoRF regions in the structural and nonstructural proteins of HTLV-1. We found abundant IDPRs in HTLV-1 bZIP factor, p30, Rex, and structural nucleocapsid p15 proteins, which are involved in diverse functions such as virus proliferation, mRNA export, and genomic RNA binding. Our study analyzed the HTLV-1 proteome with the perspective of intrinsic disorder identification. We propose that the intrinsic disorder analysis of HTLV-1 proteins may form the basis for the development of protein disorder-based drugs.

HIV co-infection in HTLV-1 carriers in Spain

BACKGROUND

Human retroviruses HIV and HTLV share transmission routes. HIV widely spread in Spain during the 80 s through injection drug use and sex, and nowadays HIV rates in Spain account for one of the largest in Europe. In contrast, HTLV-1 is not endemic in Spain, despite hosting huge numbers of migrants from highly endemic regions. Herein, we report the rate and main features of the HIV-HTLV co-infected population in Spain.

METHODS

A national registry exists in Spain for HTLV since year 1989. Data from standardized case report forms and one centralized lab repository were reviewed, especially for the subset with HTLV-HIV co-infection.

RESULTS

Up to December 2018, a total of 369 individuals with HTLV-1 had been diagnosed in Spain. 64% of the population were females, and Latin American individuals accounted for 64.5%. Classical HTLV-associated illnesses were found in 12.7% (myelopathy) and 7.6% (leukemia). HIV coinfection was found in 12 (3.2%). Of those, 3 patients (25%) were female and 39 (75%) were of non-Spanish origin. All but two harbored HIV-1 subtype B, being non-B variants found in the two West Africans. Exposure had been sexual in most cases, being 4 homosexual men. Seven HTLV-HIV co-infected patients had developed AIDS and two had developed myelopathy. There was no evidence for increased HTLV-1 clinical pathogenicity due to HIV coinfection.

CONCLUSION

HIV coinfection is infrequent (<5%) among HTLV-1 carriers in Spain. More than half of co-infected patients come from Latin America. Sexual contact is the most frequent risk behavior, being MSM one third of cases. Late diagnosis explains the high rate (9/12) of clinical manifestations in our HIV-HTLV co-infected population.

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

Characterizing Drug-Target Residence Time with Metadynamics: How To Achieve Dissociation Rate Efficiently without Losing Accuracy against Time-Consuming Approaches

Drug-target residence time plays a vital role in drug efficacy. However, there is still no effective strategy to predict drug residence time. Here, we propose to use the optimized (or minimized) structures derived from holo-state proteins to calculate drug residence time, which could give a comparable or even better prediction accuracy compared with those calculated utilizing a large number of molecular dynamics (MD) structures based on the Poisson process. Besides, in addition to the Poisson process, one may use fewer samples for predicting residence time due to the reason that, in a large extent, the calculated drug residence time is stable and independent of the number of samples used for the prediction. With remarkably reduced computational load, the proposed strategy may be promising for large-scale drug residence time prediction, such as post-processing in virtual screening (VS) and lead compound optimization.

Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach

Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?G_bind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.

Privileged Structures Meet Human T-Cell Leukemia Virus-1 (HTLV-1): C2-Symmetric 3,4-Disubstituted Pyrrolidines as Nonpeptidic HTLV-1 Protease Inhibitors

3,4-disubstituted pyrrolidines originally designed to inhibit the closely related HIV-1 protease were evaluated as privileged structures against HTLV-1 protease (HTLV-1 PR). The most potent inhibitor of this series exhibits two-digit nanomolar affinity and represents, to the best of our knowledge, the most potent nonpeptidic inhibitor of HTLV-1 PR described so far. The X-ray structures of two representatives bound to HTLV-1 PR were determined, and the structural basis of their affinity is discussed.