Search

Inhibition of HIV-1 protease is a cornerstone of antiretroviral therapy. However, the notorious ability of HIV-1 to develop resistance to protease inhibitors (PIs), particularly darunavir (DRV), poses a major challenge. Using quantum chemistry and computer simulations, this study aims to investigate the interactions between two novel PIs, GRL-004 and GRL-063, as well as a wild-type (WT) HIV strain and a DRV-resistant mutant strain. To do this, we used molecular docking, molecular dynamics simulations, and quantum mechanical calculations to check how well GRL-004 and GRL-063 bound to both WT and DRV-resistant proteases. The results show that GRL-004 and GRL-063 bind very well to ASP29 in the WT strain. ASP29 is an important amino acid in the HIV protease dimer. Remarkably, amino acids such as ILE50 in the WT strains showed substantial binding energies to both drugs. Quantum energy calculations showed a slight reduction in the energy affinity of the interaction between the MUT strain and the ligand GRL-063, compared to the WT strain. GRL-004 showed similar interaction energy for both strains, suggesting that it has greater plasticity than GRL-063 despite its lower interaction affinity. Furthermore, GLY49B demonstrated strong binding energies regardless of mutations. Other relevant residues with strong binding energies include GLY49B, PHE82A, PRO81A, ASP29A, ASP25A and ALA28B. This study improves our understanding of receptor-ligand dynamics and the adaptability of new protease inhibitors (PIs), which has profound implications for the innovation of future antiretroviral drugs.

The human immunodeficiency virus presents a significant global health challenge due to its rapid mutation and the development of resistance mechanisms against antiretroviral drugs. Recent studies demonstrate the impressive performance of machine learning (ML) and deep learning (DL) models in predicting the drug resistance profile of specific FDA-approved inhibitors. However, generalizing ML and DL models to learn not only from isolates but also from inhibitor representations remains challenging for HIV-1 infection. We propose a novel drug-isolate-fold change (DIF) model framework that aims to predict drug resistance score directly from the protein sequence and inhibitor representation. Various ML and DL models, inhibitor representations, and protein representations were analyzed through realistic validation mechanisms. To enhance the molecular learning capacity of DIF models, we employ a transfer learning approach by pretraining a graph neural network (GNN) model for activity prediction on a data set of 4855 HIV-1 protease inhibitors (PIs). By performing various realistic validation strategies on internal and external genotype-phenotype data sets, we statistically show that the learned representations of inhibitors improve the predictive ability of DIF-based ML and DL models. We achieved an accuracy of 0.802, AUROC of 0.874, and r of 0.727 for the unseen external PIs. By comparing the DIF-based models with a null model consisting of isolate-fold change (IF) architecture, it is observed that the DIF models significantly benefit from molecular representations. Combined results from various testing strategies and statistical tests confirm the effectiveness of DIF models in testing novel PIs for drug resistance in the presence of an isolate.

Substituted tetrahydrofuran derivatives were designed and synthesized to serve as the P2 ligand for a series of potent HIV-1 protease inhibitors. Both enantiomers of the tetrahydrofuran derivatives were synthesized stereoselectivity in optically active forms using lipase-PS catalyzed enzymatic resolution as the key step. These tetrahydrofuran derivatives are designed to promote hydrogen bonding and van der Waals interactions with the backbone atoms in the S2 subsite of the HIV-1 protease active site. Several inhibitors displayed very potent HIV-1 protease inhibitory activity. A high-resolution X-ray crystal structure of an inhibitor-bound HIV-1 protease provided important insight into the ligand binding site interactions in the active site.

HIV-1 protease (PR) enzyme is a viable antiretroviral drug target due to its crucial role in HIV maturation. Over many decades, the HIV-1 PR enzyme has exhibited mutations brought on by drug pressure and error-prone nature of HIV-1 reverse transcriptase. Non-active site mutations have played a pivotal role in drug resistance; however, their mechanism of action has not been fully elucidated. We investigated how non-active site mutations affect the conformational stability and drug binding ability of HIV-1 PR. In light of this, we studied a novel HIV-1 subtype C protease variant containing an insertion of valine (↑V) in the hinge region. We analysed the mutations in the presence and absence of ten background mutations. Molecular dynamics simulations revealed that both with and without the background mutations, the PR exhibited increased flexibility of hinge, flaps and fulcrum regions. This allowed the PR to adopt a wider flap conformation when in complex with several inhibitors. Additionally, the simulations revealed that the protease inhibitors (PIs) could not bring the mutated variant proteases into a stable, closed conformation, resulting in increased solvent exposure of the inhibitors. Together, these results suggest that the mutations decrease the favourability of binding by altering the dynamics of the flap regions. Notably, the insertion mutation increased PR hinge flexibility and the introduction of background mutations compensated for this by stabilising the cantilever and hinge regions. Together, these findings provide insight into how non-active site mutations affect PR conformational dynamics in critical areas of the PR thus impacting on drug binding capacity and potentially contributing to drug resistance., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Yasien Sayed reports financial support was provided by University of the WitwatersrandJohannesburg School of Molecular and Cell Biology. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

The HIV-1 protease is critical for the process of viral maturation and as such, it is one of the most well characterized proteins in the Protein Data Bank. There is some evidence to suggest that the HIV-1 protease is capable of accommodating small molecule fragments at several locations on its surface outside of the active site. However, some pockets on the surface of proteins remain unformed in the apo structure and are termed "cryptic sites." To date, no cryptic sites have been identified in the structure of HIV-1 protease. Here, we characterize a novel cryptic cantilever pocket on the surface of the HIV-1 protease through mixed-solvent molecular dynamics simulations using several probes. Interestingly, we noted that several homologous retroviral proteases exhibit evolutionarily conserved dynamics in the cantilever region and possess a conserved pocket in the cantilever region. Immobilization of the cantilever region of the HIV-1 protease via disulfide cross-linking resulted in curling-in of the flap tips and the propensity for the protease to adopt a semi-open flap conformation. Structure-based analysis and fragment-based screening of the cryptic cantilever pocket suggested that the pocket may be capable of accommodating ligand structures. Furthermore, molecular dynamics simulations of a top scoring fragment bound to the cryptic pocket illustrated altered flap dynamics of the fragment-bound enzyme. Together, these results suggest that the mobility of the cantilever region plays a key role in the global dynamics of retroviral proteases. Therefore, the cryptic cantilever pocket of the HIV-1 protease may represent an interesting target for future in vitro studies., (© 2024 The Author(s). PROTEINS: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Water molecules play various roles in target-ligand binding. For example, they can be replaced by the ligand and leave the surface of the binding pocket or stay conserved in the interface and form bridges with the target. While experimental techniques supply target-ligand complex structures at an increasing rate, they often have limitations in the measurement of a detailed water structure. Moreover, measurements of binding thermodynamics cannot distinguish between the different roles of individual water molecules. However, such a distinction and classification of the role of individual water molecules would be key to their application in drug design at atomic resolution. In this study, we investigate a quantitative approach for the description of the role of water molecules during ligand binding. Starting from complete hydration structures of the free and ligand-bound target molecules, binding enthalpy scores are calculated for each water molecule using quantum mechanical calculations. A statistical evaluation showed that the scores can distinguish between conserved and displaced classes of water molecules. The classification system was calibrated and tested on more than 1000 individual water positions. The practical tests of the enthalpic classification included important cases of antiviral drug research on HIV-1 protease inhibitors and the Influenza A ion channel. The methodology of classification is based on open source program packages, Gromacs, Mopac, and MobyWat, freely available to the scientific community.

This study introduces a novel method named multiple parameter replica exchange Gaussian accelerated molecular dynamics (MP-Rex-GaMD), building on the Gaussian accelerated molecular dynamics (GaMD) algorithm. GaMD enhances sampling and retrieves free energy information for biomolecular systems by adding a harmonic boost potential to smooth the potential energy surface without the need for predefined reaction coordinates. Our innovative approach advances the acceleration power and energetic reweighting accuracy of GaMD by incorporating a replica exchange algorithm that enables the exchange of multiple parameters, including the GaMD boost parameters of force constant and energy threshold, as well as temperature. Applying MP-Rex-GaMD to the three model systems of dialanine, chignolin, and HIV protease, we demonstrate its superior capability over conventional molecular dynamics and GaMD simulations in exploring protein conformations and effectively navigating various biomolecular states across energy barriers. MP-Rex-GaMD allows users to accurately map free energy landscapes through energetic reweighting, capturing the ensemble of biomolecular states from low-energy conformations to rare high-energy transitions within practical computational time scales.

Drug resistance of pathogens, including viruses, is one of the reasons for decreased efficacy of therapy. Considering the impact of HIV type 1 (HIV-1) on the development of progressive immune dysfunction and the rapid development of drug resistance, the analysis of HIV-1 resistance is of high significance. Currently, a substantial amount of data has been accumulated on HIV-1 drug resistance that can be used to build both qualitative and quantitative models of HIV-1 drug resistance. Quantitative models of drug resistance can enrich the information about the efficacy of a particular drug in the scheme of antiretroviral therapy. In our study, we investigated the possibility of developing models for quantitative prediction of HIV-1 resistance to eight protease inhibitors based on the analysis of amino acid sequences of HIV-1 protease for 900 virus variants. We developed random forest regression (RFR), support vector regression (SVR), and self-consistent regression (SCR) models using binary vectors containing values from 0 or 1, depending on the presence of a specific peptide fragment in each amino acid sequence as independent variables, while fold ratio, reflecting the level of resistance, was the predicted variable. The SVR and SCR models showed the highest predictive performances. The models built demonstrate reasonable performances for eight out of nine (R 2 varied from 0.828 to 0.909) protease inhibitors, while R 2 for predicting tipranavir fold ratio was lower (R 2 was 0.642). We believe that the developed approach can be applied to evaluate drug resistance of molecular targets of other viruses where appropriate experimental data are available.

Background & objectives Despite advancements in antiretroviral therapy, drug-resistant strains of HIV (human immunodeficiency virus) remain a global health concern. Natural compounds from medicinal plants offer a promising avenue for developing new HIV-1 PR (protease) inhibitors. This study aimed to explore the potential of compounds derived from Calotropis procera, a medicinal plant, as inhibitors of HIV-1 PR. Methods This in silico study utilized natural compound information and the crystal structure of HIV-1 PR. Molecular docking of 17 steroidal cardenolides from Calotropis procera against HIV-1 PR was performed using AutoDock 4.2 to identify compounds with higher antiviral potential. A dynamic simulation study was performed to provide insights into the stability, binding dynamics, and potential efficacy of the top potential antiviral compound as an HIV-1 therapeutic. Results We found that all tested cardenolides had higher binding affinities than Amprenavir, indicating their potential as potent HIV-1 PR inhibitors. Voruscharin and uscharidin displayed the strongest interactions, forming hydrogen bonds and hydrophobic interactions with HIV-1 PR. Voruscharin showed improved stability with lower RMSD (Root Mean Square Deviation) values and reduced fluctuations in binding site residues but increased flexibility in certain regions. The radius of gyration analysis confirmed a stable binding pose between HIV-1 PR and Voruscharin. Interpretation & conclusions These findings suggest that Calotropis procera could potentially be a source of compounds for developing novel HIV-1 PR inhibitors, contributing to the efforts to combat HIV. Further studies and clinical trials are needed to evaluate the safety and efficacy of these compounds as potential drug candidates for the treatment of HIV-1 infection.

The Gag-Pol polyprotein in human immunodeficiency virus type I (HIV-1) encodes enzymes that are essential for virus replication: protease (PR), reverse transcriptase (RT), and integrase (IN). The mature forms of PR, RT and IN are homodimer, heterodimer and tetramer, respectively. The precise mechanism underlying the formation of dimer or tetramer is not yet understood. Here, to gain insight into the dimerization of PR and RT in the precursor, we prepared a model precursor, PR-RT, incorporating an inactivating mutation at the PR active site, D25A, and including two residues in the p6* region, fused to a SUMO-tag, at the N-terminus of the PR region. We also prepared two mutants of PR-RT containing a dimer dissociation mutation either in the PR region, PR(T26A)-RT, or in the RT region, PR-RT(W401A). Size exclusion chromatography showed both monomer and dimer fractions in PR-RT and PR(T26A)-RT, but only monomer in PR-RT(W401A). SEC experiments of PR-RT in the presence of protease inhibitor, darunavir, significantly enhanced the dimerization. Additionally, SEC results suggest an estimated PR-RT dimer dissociation constant that is higher than that of the mature RT heterodimer, p66/p51, but slightly lower than the premature RT homodimer, p66/p66. Reverse transcriptase assays and RT maturation assays were performed as tools to assess the effects of the PR dimer-interface on these functions. Our results consistently indicate that the RT dimer-interface plays a crucial role in the dimerization in PR-RT, whereas the PR dimer-interface has a lesser role., (© 2024 The Author(s). Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

The battle against viral drug resistance highlights the need for innovative approaches to replace time-consuming and costly traditional methods. Deep generative models offer automation potential, especially in the fight against Human immunodeficiency virus (HIV), as they can synthesize diverse molecules effectively. In this paper, an application of an LSTM-based deep generative model named "LSTM-ProGen" is proposed to be tailored explicitly for the de novo design of drug candidate molecules that interact with a specific target protein (HIV-1 protease). LSTM-ProGen distinguishes itself by employing a long-short-term memory (LSTM) architecture, to generate novel molecules target specificity against the HIV-1 protease. Following a thorough training process involves fine-tuning LSTM-ProGen on a diverse range of compounds sourced from the ChEMBL database. The model was optimized to meet specific requirements, with multiple iterations to enhance its predictive capabilities and ensure it generates molecules that exhibit favorable target interactions. The training process encompasses an array of performance evaluation metrics, such as drug-likeness properties. Our evaluation includes extensive silico analysis using molecular docking and PCA-based visualization to explore the chemical space that the new molecules cover compared to those in the training set. These evaluations reveal that a subset of 12 de novo molecules generated by LSTM-ProGen exhibit a striking ability to interact with the target protein, rivaling or even surpassing the efficacy of native ligands. Extended versions with further refinement of LSTM-ProGen hold promise as versatile tools for designing efficacious and customized drug candidates tailored to specific targets, thus accelerating drug development and facilitating the discovery of new therapies for various diseases., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Albrijawi, Alhajj. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Sabia parviflora ( SP , "xiao hua qing feng teng" in Chinese) was recorded as an important ethnic medicine to be used for treating viral hepatitis. The antiviral activity of four SP extracts and potent antiviral compounds evaluated with cathepsin L protease (Cat L PR) and HIV-1 protease (HIV-1 PR). UPLC-HRMS was used for identifying the bioactive components. In addition, the possible inhibitory mechanism of the identified compounds on viral protease was further discussed by molecular docking. As a result, four extracts of SP exhibited inhibitory activity of HIV-1 PR and Cat L PR with IC 50 range from 0.015 to 0.80 mg/mL. Meanwhile, six compounds inhibited HIV-1 PR with IC 50 range from 0.032 to 0.80 mg/mL. Moreover, procyanidin B2 had good affinity for HIV-1 PR and CatL PR protein, respectively. These findings suggest S. parviflora leaves can be used for treating HIV and procyanidin B2 may play a role in antiviral protease.

HIV-1 protease inhibitors are an essential component of antiretroviral therapy. However, drug resistance is a pervasive issue motivating a persistent search for novel therapies. Recent reports found that when protease activates within the host cell's cytosol, it facilitates the pyroptotic killing of infected cells. This has led to speculation that promoting protease activation, rather than inhibiting it, could help to eradicate infected cells and potentially cure HIV-1 infection. Here, we used a nanoscale flow cytometry-based assay to characterize protease resistance mutations and polymorphisms. We quantified protease activity, viral concentration, and premature protease activation and confirmed previous findings that major resistance mutations generally destabilize the protease structure. Intriguingly, we found evidence that common polymorphisms in the hinge domain of protease can influence its susceptibility to premature activation. This suggests that viral heterogeneity could pose a considerable challenge for therapeutic strategies aimed at inducing premature protease activation in the future.

Protease inhibitors (PIs) remain an important component of antiretroviral therapy for the treatment of HIV-1 infection due to their high genetic barrier to resistance development. Nevertheless, the two most commonly prescribed HIV PIs, atazanavir and darunavir, still require co-administration with a pharmacokinetic boosting agent to maintain sufficient drug plasma levels which can lead to undesirable drug-drug interactions. Herein, we describe GS-9770, a novel investigational non-peptidomimetic HIV PI with unboosted once-daily oral dosing potential due to improvements in its metabolic stability and its pharmacokinetic properties in preclinical animal species. This compound demonstrates potent inhibitory activity and high on-target selectivity for recombinant HIV-1 protease versus other aspartic proteases tested. In cell culture, GS-9770 inhibits Gag polyprotein cleavage and shows nanomolar anti-HIV-1 potency in primary human cells permissive to HIV-1 infection and against a broad range of HIV subtypes. GS-9770 demonstrates an improved resistance profile against a panel of patient-derived HIV-1 isolates with resistance to atazanavir and darunavir. In resistance selection experiments, GS-9770 prevented the emergence of breakthrough HIV-1 variants at all fixed drug concentrations tested and required multiple protease substitutions to enable outgrowth of virus exposed to escalating concentrations of GS-9770. This compound also remained fully active against viruses resistant to drugs from other antiviral classes and showed no in vitro antagonism when combined pairwise with drugs from other antiretroviral classes. Collectively, these preclinical data identify GS-9770 as a potent, non-peptidomimetic once-daily oral HIV PI with potential to overcome the persistent requirement for pharmacological boosting with this class of antiretroviral agents., Competing Interests: All authors were employed and compensated by Gilead Sciences and may be shareholders of the company.

The present work reports the cost-effective, high yielding and environmentally acceptable preparation of unsymmetrical ureas from thiocarbamate salts using sodium percarbonate as an oxidant. Efficacy of the unsymmetrical ureas as potential human immune deficiency virus (HIV-1) protease inhibitors has been evaluated via in silico approach. The results revealed interactions of the urea compounds at the active site of the enzyme with favorable binding affinities causing possible mutations hindering the functioning of the enzyme. Further computational assessment of IC50 using known references satisfactorily authenticated the inhibitory action of the selected compounds against HIV-1 protease. Added to the easy synthesis of the ureas following an environmentally benign protocol, this work may be a valuable addition to the ongoing search for drugs with better efficacy profiles and reduced toxicity against HIV., (© 2023. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

A novel kind of potent HIV-1 protease inhibitors, containing diverse hydroxyphenylacetic acids as the P2-ligands and 4-substituted phenyl sulfonamides as the P2' ligands, were designed, synthesized and evaluated in this work. Majority of the target compounds exhibited good to excellent activity against HIV-1 protease with IC 50 values below 200 nM. In particular, compound 18d with a 2-(3,4-dihydroxyphenyl) acetamide as the P2 ligand and a 4- methoxybenzene sulfonamide P2' ligand exhibited inhibitory activity IC 50 value of 0.54 nM, which was better than that of the positive control darunavir (DRV). More importantly, no significant decline of the potency against HIV-1 DRV R S (DRV-resistant mutation) and HIV-1NL4_3 variant (wild type) for 18d was detected. The molecular docking study of 18d with HIV-1 protease (PDB-ID: 1T3R, www.rcsb.org) revealed possible binding mode with the HIV-1 protease. These results suggested the validity of introducing phenol-derived moieties into the P2 ligand and deserve further optimization which was of great value for future discovery of novel HIV-1 protease., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs., (© 2024. The Author(s).)

HIV infection remains a global health issue plagued by drug resistance and virological failure. Natural polymorphisms (NPs) contained within several African and Brazilian protease (PR) variants have been shown to induce a conformational landscape of more closed conformations compared to the sequence of subtype B prevalent in North America and Western Europe. Here we demonstrate through experimental pulsed EPR distance measurements and molecular dynamic (MD) simulations that the two common NPs D60E and I62V found within subtypes F and H can induce a closed conformation when introduced into HIV-1PR subtype B. Specifically, D60E alters the conformation in subtype B through the formation of a salt bridge with residue K43 contained within the nexus between the flap and hinge region of the HIV-1 PR fold. On the other hand, I62V modulates the packing of the hydrophobic cluster of the cantilever and fulcrum, also resulting in a more closed conformation.

Acquired immunodeficiency syndrome (AIDS) is caused by human immunodeficiency virus (HIV). HIV protease, reverse transcriptase, and integrase are targets of current drugs to treat the disease. However, anti-viral drug-resistant strains have emerged quickly due to the high mutation rate of the virus, leading to the demand for the development of new drugs. One attractive target is Gag-Pol polyprotein, which plays a key role in the life cycle of HIV. Recently, we found that a combination of M50I and V151I mutations in HIV-1 integrase can suppress virus release and inhibit the initiation of Gag-Pol autoprocessing and maturation without interfering with the dimerization of Gag-Pol. Additional mutations in integrase or RNase H domain in reverse transcriptase can compensate for the defect. However, the molecular mechanism is unknown. There is no tertiary structure of the full-length HIV-1 Pol protein available for further study. Therefore, we developed a workflow to predict the tertiary structure of HIV-1 NL4.3 Pol polyprotein. The modeled structure has comparable quality compared with the recently published partial HIV-1 Pol structure (PDB ID: 7SJX). Our HIV-1 NL4.3 Pol dimer model is the first full-length Pol tertiary structure. It can provide a structural platform for studying the autoprocessing mechanism of HIV-1 Pol and for developing new potent drugs. Moreover, the workflow can be used to predict other large protein structures that cannot be resolved via conventional experimental methods.

Proteolytic processing of human immunodeficiency virus type 1 particles mediated by viral protease (PR) is essential for acquiring virus infectivity. Activation of PR embedded in Gag-Pol is triggered by Gag-Pol dimerization during virus assembly. We previously reported that amino acid substitutions at the RT tryptophan repeat motif destabilize virus-associated RT and attenuate the ability of efavirenz (EFV, an RT dimerization enhancer) to increase PR-mediated Gag cleavage efficiency. Furthermore, a single amino acid change at RT significantly reduces virus yields due to enhanced Gag cleavage. These data raise the possibility of the RT domain contributing to PR activation by promoting Gag-Pol dimerization. To test this hypothesis, we investigated the putative involvement of a hydrophobic leucine repeat motif (LRM) spanning RT L282 to L310 in RT/RT interactions. We found that LRM amino acid substitutions led to RT instability and that RT is consequently susceptible to degradation by PR. The LRM mutants exhibited reduced Gag cleavage efficiencies while attenuating the EFV enhancement of Gag cleavage. In addition, an RT dimerization-defective mutant, W401A, reduced enhanced Gag cleavage via a leucine zipper (LZ) motif inserted at the deleted Gag-Pol region. Importantly, the presence of RT and integrase domains failed to counteract the LZ enhancement of Gag cleavage. A combination of the Gag cleavage enhancement factors EFV and W402A markedly impaired Gag cleavage, indicating a disruption of W402A Gag-Pol dimerization following EFV binding to W402A Gag-Pol. Our results support the idea that RT modulates PR activation by affecting Gag-Pol/Gag-Pol interaction. IMPORTANCE A stable reverse transcriptase (RT) p66/51 heterodimer is required for HIV-1 genome replication in host cells following virus entry. The activation of viral protease (PR) to mediate virus particle processing helps viruses acquire infectivity following cell release. RT and PR both appear to be major targets for inhibiting HIV-1 replication. We found a strong correlation between impaired p66/51RT stability and deficient PR-mediated Gag cleavage, suggesting that RT/RT interaction is critical for triggering PR activation via the promotion of adequate Gag-Pol dimerization. Accordingly, RT/RT interaction is a potentially advantageous method for anti-HIV/AIDS therapy if it is found to simultaneously block PR and RT enzymatic activity., Competing Interests: The authors declare no conflict of interest.

Structure-based design, synthesis, X-ray structural studies, and biological evaluation of a new series of potent HIV-1 protease inhibitors are described. These inhibitors contain various pyridyl-pyrimidine, aryl thiazole or alkylthiazole derivatives as the P2 ligands in combination with darunavir-like hydroxyethylamine sulfonamide isosteres. These heterocyclic ligands are inherent to kinase inhibitor drugs, such as nilotinib and imatinib. These ligands are designed to make hydrogen bonding interactions with the backbone atoms in the S2 subsite of HIV-1 protease. Various benzoic acid derivatives have been synthesized and incorporation of these ligands provided potent inhibitors that exhibited subnanomolar level protease inhibitory activity and low nanomolar level antiviral activity. Two high resolution X-ray structures of inhibitor-bound HIV-1 protease were determined. These structures provided important ligand-binding site interactions for further optimization of this class of protease inhibitors., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Masson SAS. All rights reserved.)

The design, synthesis, and biological evaluation of a novel series of HIV-1 protease inhibitors containing pyrrolidines with diverse linkers as the P2 ligands and various aromatic derivatives as the P2' ligands were described. A number of inhibitors demonstrated potent efficacy in both enzyme and cellular assays, as well as relatively low cytotoxicity. In particular, inhibitor 34b with a (R)-pyrrolidine-3-carboxamide P2 ligand and a 4-hydroxyphenyl P2' ligand displayed exceptional enzyme inhibitory activity with an IC 50 value of 0.32 nM. Furthermore, 34b also exhibited robust antiviral activity against both wild-type HIV-1 and drug-resistant variant with low micromolar EC 50 values. In addition, the molecular modelling studies revealed the extensive interactions between inhibitor 34b and the backbone residues of both wild-type and drug-resistant HIV-1 protease. These results suggested the feasibility of utilizing pyrrolidine derivatives as the P2 ligands and provided valuable information for further design and optimization of highly potent HIV-1 protease inhibitors., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Yucheng Wang reports financial support was provided by National Natural Science Foundation of China., (Copyright © 2023 Elsevier Masson SAS. All rights reserved.)

Drug resistance in antiviral treatments is a serious public health problem. Viral proteins mutate very fast, giving them a way to escape drugs by lowering drug binding affinity but with compromised function. Human immunodeficiency virus type I (HIV-1) protease, a critical antiretroviral therapeutic target, represents a model for such viral regulation under inhibition. Drug inhibitors of HIV-1 protease lose effectiveness as the protein evolves through several variants to become more resistant. However, the detailed mechanism of drug resistance in HIV-1 protease is still unclear. Here, we test the hypothesis that mutations throughout the protease alter the protein conformational ensemble to weaken protein-inhibitor binding, resulting in an inefficient protease but still viable virus. Comparing conformational ensembles between variants and the wild type helps detect these function-related dynamical changes. All analyses of over 30 μs simulations converge to the conclusion that conformational dynamics of more drug-resistant variants are more different from that of the wild type. Distinct roles of mutations during viral evolution are discussed, including a mutation predominantly contributing to the increase of drug resistance and a mutation that is responsible (synergistically) for restoring catalytic efficiency. Drug resistance is mainly due to altered flap dynamics that hinder the access to the active site. The mutant variant showing the highest drug resistance has the most ″collapsed″ active-site pocket and hence the largest magnitude of hindrance of drug binding. An enhanced difference contact network community analysis is applied to understand allosteric communications. The method summarizes multiple conformational ensembles in one community network and can be used in future studies to detect function-related dynamics in proteins.

Human immunodeficiency virus 1 (HIV-1) viral protease (PR) is one of the most studied viral enzymes and a crucial antiviral target. Despite its well-characterized role in virion maturation, an increasing body of research is starting to focus on its ability to cleave host cell proteins. Such findings are apparently in contrast with the dogma of HIV-1 PR activity being restricted to the interior of nascent virions and suggest catalytic activity within the host cell environment. Given the limited amount of PR present in the virion at the time of infection, such events mainly occur during late viral gene expression, mediated by newly synthesized Gag-Pol polyprotein precursors, rather than before proviral integration. HIV-1 PR mainly targets proteins involved in three different processes: those involved in translation, those controlling cell survival, and restriction factors responsible for innate/intrinsic antiviral responses. Indeed, by cleaving host cell translation initiation factors, HIV-1 PR can impair cap-dependent translation, thus promoting IRES-mediated translation of late viral transcripts and viral production. By targeting several apoptotic factors, it modulates cell survival, thus promoting immune evasion and viral dissemination. Additionally, HIV-1 PR counteracts restriction factors incorporated in the virion that would otherwise interfere with nascent virus vitality. Thus, HIV-1 PR appears to modulate host cell function at different times and locations during its life cycle, thereby ensuring efficient viral persistency and propagation. However, we are far from having a complete picture of PR-mediated host cell modulation, which is emerging as a field that needs further investigation.

We report here the synthesis and biological evaluation of darunavir derived HIV-1 protease inhibitors and their functional effect on enzyme inhibition and antiviral activity in MT-2 cell lines. The P2' 4-amino functionality was modified to make a number of amide derivatives to interact with residues in the S2' subsite of the HIV-1 protease active site. Several compounds exhibited picomolar enzyme inhibitory and low nanomolar antiviral activity. The X-ray crystal structure of the chloroacetate derivative bound to HIV-1 protease was determined. Interestingly, the active chloroacetate group converted to the acetate functionality during X-ray exposure. The structure revealed that the P2' carboxamide functionality makes enhanced hydrogen bonding interactions with the backbone atoms in the S2'-subsite., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

The influence of the dynamical flexibility of enzymes on reaction mechanisms is a cornerstone in biological sciences. In this study, we aim to 1) study the convergence of the activation free energy by using the first step of the reaction catalysed by HIV-1 protease as a case study, and 2) provide further evidence for a mechanistic divergence in this enzyme, as two different reaction pathways were seen to contribute to this step. We used quantum mechanics/molecular mechanics molecular dynamics simulations, on four different initial conformations that led to different barriers in a previous study. Despite the sampling, the four activation free energies still spanned a range of 5.0 kcal ⋅ mol -1 . Furthermore, the new simulations did confirm the occurrence of an unusual mechanistic divergence, with two different mechanistic pathways displaying equivalent barriers. An active-site water molecule is proposed to influence the mechanistic pathway., (© 2022 Wiley-VCH GmbH.)

The retropepsin (PR) of the Bovine leukemia virus (BLV) plays, as in other retroviruses, a crucial role in the transition from the non-infective viral particle to the infective virion by processing the polyprotein Gag. PR is expressed as an immature precursor associated with Gag, after an occasional -1 ribosomal frameshifting event. Self-hydrolysis of PR at specific N- and C-terminal sites releases the monomer that dimerizes giving rise to the active protease. We designed a strategy to express BLV PR in E. coli as a fusion protein with maltose binding protein, with a six-histidine tag at its N-terminal end, and bearing a tobacco etch virus protease hydrolysis site. This allowed us to obtain soluble and mature recombinant PR in relatively good yields, with exactly the same amino acid composition as the native protein. As PR presents relative promiscuity for the hydrolysis sites we designed four fluorogenic peptide substrates based on Förster resonance energy transfer (FRET) in order to characterize the activity of the recombinant enzyme. These substrates opened the way to perform kinetic studies, allowing us to characterize the dimer-monomer equilibrium. Furthermore, we obtained kinetic evidence for the existence of a conformational change that enables the interaction with the substrate. These results constitute a starting point for the elucidation of the kinetic properties of BLV-PR, and may be relevant not only to improve the chemical warfare against this virus but also to better understand other viral PRs., Competing Interests: The authors have declared that no competing interests exist.

The infectivity of HIV-1 requires its protease (PR) cleave multiple cut-sites with low sequence similarity. The diversity of cleavage sites has made it challenging to investigate the underlying sequence properties that determine binding and turnover of substrates by PR. We engineered a mutational scanning approach utilizing yeast display, flow cytometry, and deep sequencing to systematically measure the impacts of all individual amino acid changes at 12 positions in three different cut-sites (MA/CA, NC/p1, and p1/p6). The resulting fitness landscapes revealed common physical features that underlie cutting of all three cut-sites at the amino acid positions closest to the scissile bond. In contrast, positions more than two amino acids away from the scissile bond exhibited a strong dependence on the sequence background of the rest of the cut-site. We observed multiple amino acid changes in cut-sites that led to faster cleavage rates, including a preference for negative charge five and six amino acids away from the scissile bond at locations where the surface of protease is positively charged. Analysis of individual cut sites using full-length matrix-capsid proteins indicate that long-distance sequence context can contribute to cutting efficiency such that analyses of peptides or shorter engineered constructs including those in this work should be considered carefully. This work provides a framework for understanding how diverse substrates interact with HIV-1 PR and can be extended to investigate other viral PRs with similar properties., (© 2022 The Protein Society.)

Starting three decades ago and spreading rapidly around the world, acquired immunodeficiency syndrome (AIDS) is an infectious disease distinct from other contagious diseases by its unique ways of transmission. Over the past few decades, research into new drug compounds has been accompanied by extensive advances, and the design and manufacture of drugs that inhibit virus enzymes is one way to combat the AIDS virus. Since blocking enzyme activity can kill a pathogen or correct a metabolic imbalance, the design and use of enzyme inhibitors is a new approach against viruses. We carried out an in-depth analysis of the efficacy of atazanavir and its newly designed analogs as human immunodeficiency virus (HIV) protease inhibitors using molecular docking. The best-designed analogs were then compared with atazanavir by the molecular dynamics simulation. The most promising results were ultimately found based on the docking analysis for HIV protease. Several exhibited an estimated free binding energy lower than -9.45 kcal/mol, indicating better prediction results than the atazanavir. ATV7 inhibitor with antiviral action may be more beneficial for infected patients with HIV. Molecular dynamics analysis and binding energy also showed that the ATV7 drug had more inhibitory ability than the atazanavir drug., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

HIV-1 protease (PR) is a viral enzyme that cleaves the Gag and Gag-Pol polyprotein precursors to convert them into their functional forms, a process which is essential to generate infectious viral particles. Due to its broad substrate specificity, HIV-1 PR can also cleave certain host cell proteins. Several studies have identified host cell substrates of HIV-1 PR and described the potential impact of their cleavage on HIV-1-infected cells. Of particular interest is the interaction between PR and the caspase recruitment domain-containing protein 8 (CARD8) inflammasome. A recent study demonstrated that CARD8 can sense HIV-1 PR activity and induce cell death. While PR typically has low levels of intracellular activity prior to viral budding, premature PR activation can be achieved using certain non-nucleoside reverse transcriptase inhibitors (NNRTIs), resulting in CARD8 cleavage and downstream pyroptosis. Used together with latency reversal agents, the induction of premature PR activation to trigger CARD8-mediated cell killing may help eliminate latent reservoirs in people living with HIV. This represents a novel strategy of utilizing PR as an antiviral target through premature activation rather than inhibition. In this review, we discuss the viral and host substrates of HIV-1 protease and highlight potential applications and advantages of targeting CARD8 sensing of HIV-1 PR.

HIV-1 encodes a viral protease that is essential for the maturation of infectious viral particles. While protease inhibitors are effective antiretroviral agents, recent studies have shown that prematurely activating, rather than inhibiting, protease function leads to the pyroptotic death of infected cells, with exciting implications for efforts to eradicate viral reservoirs. Despite 40 years of research into the kinetics of protease activation, it remains unclear exactly when protease becomes activated. Recent reports have estimated that protease activation occurs minutes to hours after viral release, suggesting that premature protease activation is challenging to induce efficiently. Here, monitoring viral protease activity with sensitive techniques, including nanoscale flow cytometry and instant structured illumination microscopy, we demonstrate that the viral protease is activated within cells prior to the release of free virions. Using genetic mutants that lock protease into a precursor conformation, we further show that both the precursor and mature protease have rapid activation kinetics and that the activity of the precursor protease is sufficient for viral fusion with target cells. Our finding that HIV-1 protease is activated within producer cells prior to release of free virions helps resolve a long-standing question of when protease is activated and suggests that only a modest acceleration of protease activation kinetics is required to induce potent and specific elimination of HIV-infected cells. IMPORTANCE HIV-1 protease inhibitors have been a mainstay of antiretroviral therapy for more than 2 decades. Although antiretroviral therapy is effective at controlling HIV-1 replication, persistent reservoirs of latently infected cells quickly reestablish replication if therapy is halted. A promising new strategy to eradicate the latent reservoir involves prematurely activating the viral protease, which leads to the pyroptotic killing of infected cells. Here, we use highly sensitive techniques to examine the kinetics of protease activation during and shortly after particle formation. We found that protease is fully activated before virus is released from the cell membrane, which is hours earlier than recent estimates. Our findings help resolve a long-standing debate as to when the viral protease is initially activated during viral assembly and confirm that prematurely activating HIV-1 protease is a viable strategy to eradicate infected cells following latency reversal.

The design, synthesis, X-ray structural, and biological evaluation of a series of highly potent HIV-1 protease inhibitors are reported herein. These inhibitors incorporate novel cyclohexane-fused tricyclic bis-tetrahydrofuran as P2 ligands in combination with a variety of P1 and P2' ligands. The inhibitor with a difluoromethylphenyl P1 ligand and a cyclopropylaminobenzothiazole P2' ligand exhibited the most potent antiviral activity. Also, it maintained potent antiviral activity against a panel of highly multidrug-resistant HIV-1 variants. The corresponding inhibitor with an enantiomeric ligand was significantly less potent in these antiviral assays. The new P2 ligands were synthesized in optically active form using enzymatic desymmetrization of meso-diols as the key step. To obtain molecular insight, two high-resolution X-ray structures of inhibitor-bound HIV-1 protease were determined and structural analyses have been highlighted., (© 2022 Wiley-VCH GmbH.)

A novel strategy for lead identification that we have dubbed the "Pocket-to-Lead" strategy is demonstrated using HIV-1 protease as a model target. Sometimes, it is difficult to obtain hit compounds because of the difficulties in satisfying the complex pharmacophoric features. In this study, a virtual fragment hit which does not match all of the pharmacophore features but has key interactions and vectors that could grow into remaining pharmacophore features was optimized in silico . The designed compound 9 demonstrated weak but evident inhibitory activity (IC 50 = 54 μM), and the design concept was proven by the co-crystal structure. Then, structure-based drug design promptly gave compound 14 (IC 50 = 0.0071 μM, EC 50 = 0.86 μM), an almost 10,000-fold improvement in activity from 9 . The structure of the designed molecules proved to be novel with high synthetic feasibility, indicating the usefulness of this strategy to tackle tough targets with complex pharmacophore.

Motivation: Human immunodeficiency virus (HIV) drug resistance is a global healthcare issue. The emergence of drug resistance influenced the efficacy of treatment regimens, thus stressing the importance of treatment adaptation. Computational methods predicting the drug resistance profile from genomic data of HIV isolates are advantageous for monitoring drug resistance in patients. However, existing computational methods for drug resistance prediction are either not suitable for emerging HIV strains with complex mutational patterns or lack interpretability, which is of paramount importance in clinical practice. The approach reported here overcomes these limitations and combines high accuracy of predictions and interpretability of the models., Results: In this work, a new methodology based on generative topographic mapping (GTM) for biological sequence space representation and quantitative genotype-phenotype relationships prediction purposes was introduced. The GTM-based resistance landscapes allowed us to predict the resistance of HIV strains based on sequencing and drug resistance data for three viral proteins [integrase (IN), protease (PR) and reverse transcriptase (RT)] from Stanford HIV drug resistance database. The average balanced accuracy for PR inhibitors was 0.89 ± 0.01, for IN inhibitors 0.85 ± 0.01, for non-nucleoside RT inhibitors 0.73 ± 0.01 and for nucleoside RT inhibitors 0.84 ± 0.01. We have demonstrated in several case studies that GTM-based resistance landscapes are useful for visualization and analysis of sequence space as well as for treatment optimization purposes. Here, GTMs were applied for the in-depth analysis of the relationships between mutation pattern and drug resistance using mutation landscapes. This allowed us to predict retrospectively the importance of the presence of particular mutations (e.g. V32I, L10F and L33F in HIV PR) for the resistance development. This study highlights some perspectives of GTM applications in clinical informatics and particularly in the field of sequence space exploration., Availability and Implementation: https://github.com/karinapikalyova/ISIDASeq., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

By following up on the design vector of optimizing amine-based HIV-1 protease inhibitors, we have designed and biologically evaluated a novel class of inhibitors with the free nitrogen or sulphone in morpholine cores as P2 ligands in combination with diverse substituted phenylsulfonamide P2' ligands. As it turns out, a majority of these inhibitors exhibit prominent enzymatic inhibitory activity in low nanomolar ranges with relatively low cytotoxicity. Particularly, inhibitor 1e containing a morpholine carboxamide P2 ligand and a 4-hydroxyphenylsulfonamide P2' ligand illustrates a robust enzyme inhibitory IC 50 value of 90 pM. Furthermore, 1e demonstrates impressive in vivo antiviral activity with EC 50 value of 89 nM and a degree of inhibitory potency against the DRV-resistant variant. More importantly, 1e exhibits remarkable activity with EC 50 values of 13.59 nM and 8.23 nM against subtype C HIV-1 strains ZM246 and Indie, respectively. Furthermore, the in silico studies provide molecular insights into binding features of inhibitors with HIV-1 protease, and furnish a valuable forecast on further process., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Masson SAS. All rights reserved.)

Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I-saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.

A human immunodeficiency virus-1 (HIV-1) protease is a homodimeric aspartic protease essential for the replication of HIV. The HIV-1 protease is a target protein in drug discovery for antiretroviral therapy, and various inhibitor molecules of transition state analogues have been developed. However, serious drug-resistant mutants have emerged. For understanding the molecular mechanism of the drug resistance, an accurate examination of the impacts of the mutations on ligand binding and enzymatic activity is necessary. Here, we present a molecular simulation study on the ligand binding of indinavir, a potent transition state analogue inhibitor, to the wild-type protein and a V82T/I84V drug-resistant mutant of the HIV-1 protease. We employed a hybrid ab initio quantum mechanical/molecular mechanical (QM/MM) free-energy optimization technique which combines a highly accurate QM description of the ligand molecule and its interaction with statistically ample conformational sampling of the MM protein environment by long-time molecular dynamics simulations. Through the free-energy calculations of protonation states of catalytic groups at the binding pocket and of the ligand-binding affinity changes upon the mutations, we successfully reproduced the experimentally observed significant reduction of the binding affinity upon the drug-resistant mutations and elucidated the underlying molecular mechanism. The present study opens the way for understanding the molecular mechanism of drug resistance through the direct quantitative comparison of ligand binding and enzymatic reaction with the same accuracy.

A genetic selection system for activity of HIV protease is described that is based on a synthetic substrate constructed as a modified AraC regulatory protein that when cleaved stimulate l-arabinose metabolism in an Escherichia coli araC strain. Growth stimulation on selective plates was shown to depend on active HIV protease and the scissile bond in the substrate. In addition, the growth of cells correlated well with the established cleavage efficiency of the sites in the viral polyprotein, Gag, when these sites were individually introduced into the synthetic substrate of the selection system. Plasmids encoding protease variants selected based on stimulation of cell growth in the presence of saquinavir or cleavage of a site not cleaved by wild-type protease, were indistinguishable with respect to both phenotypes. Also, both groups of selected plasmids encoded side chain substitutions known from clinical isolates or displayed different side chain substitutions but at identical positions. One highly frequent side chain substitution, E34V, not regarded as a major drug resistance substitution was found in variants obtained under both selective conditions and is suggested to improve protease processing of the synthetic substrate. This substitution is away from the substrate-binding cavity and together with other substitutions in the selected reading frames supports the previous suggestion of a substrate-binding site extended from the active site binding pocket itself., (© 2022 The Author(s). Published by Portland Press Limited on behalf of the Biochemical Society.)

To date, there are no specific treatment regimens for HIV-1-related central nervous system (CNS) complications, such as HIV-1-associated neurocognitive disorders (HAND). Here, we report that two newly generated CNS-targeting HIV-1 protease (PR) inhibitors (PIs), GRL-08513 and GRL-08613, which have a P1-3,5- bis -fluorophenyl or P1- para -monofluorophenyl ring and P2-tetrahydropyrano-tetrahydrofuran ( Tp -THF) with a sulfonamide isostere, are potent against wild-type HIV-1 strains and multiple clinically isolated HIV-1 strains (50% effective concentration [EC 50 ]: 0.0001 to ∼0.0032 μM). As assessed with HIV-1 variants that had been selected in vitro to propagate at a 5 μM concentration of each HIV-1 PI (atazanavir, lopinavir, or amprenavir), GRL-08513 and GRL-08613 efficiently inhibited the replication of these highly PI-resistant variants (EC 50 : 0.003 to ∼0.006 μM). GRL-08513 and GRL-08613 also maintained their antiviral activities against HIV-2 ROD as well as severely multidrug-resistant clinical HIV-1 variants. Additionally, when we assessed with the in vitro blood-brain barrier (BBB) reconstruction system, GRL-08513 and GRL-08613 showed the most promising properties of CNS penetration among the evaluated compounds, including the majority of FDA-approved combination antiretroviral therapy (cART) drugs. In the crystallographic analysis of compound-PR complexes, it was demonstrated that the Tp -THF rings at the P2 moiety of GRL-08513 and GRL-08613 form robust hydrogen bond interactions with the active site of HIV-1 PR. Furthermore, both the P1-3,5- bis -fluorophenyl- and P1- para -monofluorophenyl rings sustain greater contact surfaces and form stronger van der Waals interactions with PR than is the case with darunavir-PR complex. Taken together, these results strongly suggest that GRL-08513 and GRL-08613 have favorable features for patients infected with wild-type/multidrug-resistant HIV-1 strains and might serve as candidates for a preventive and/or therapeutic agent for HAND and other CNS complications.

HIV protease plays a critical role in the life cycle of the virus through the generation of mature and infectious virions. Detailed knowledge of the structure of the enzyme and its substrate has led to the development of protease inhibitors. However, the development of resistance to all currently available protease inhibitors has contributed greatly to the decreased success of antiretroviral therapy. When therapy failure occurs, multiple mutations are found within the protease sequence starting with primary mutations, which directly impact inhibitor binding, which can also negatively impact viral fitness and replicative capacity by decreasing the binding affinity of the natural substrates to the protease. As such, secondary mutations which are located outside of the active site region accumulate to compensate for the recurrently deleterious effects of primary mutations. However, the resistance mechanism of these secondary mutations is not well understood, but what is known is that these secondary mutations contribute to resistance in one of two ways, either through increasing the energetic penalty associated with bringing the protease into the closed conformation, or, through decreasing the stability of the protein/drug complex in a manner that increases the dissociation rate of the drug, leading to diminished inhibition. As a result, the elasticity of the enzyme-substrate complex has been implicated in the successful recognition and catalysis of the substrates which may be inferred to suggest that the elasticity of the enzyme/drug complex plays a role in resistance. A realistic representation of the dynamic nature of the protease may provide a more powerful tool in structure-based drug design algorithms., (© 2021. Springer Nature Switzerland AG.)

The denatured state of several proteins has been shown to display transient structures that are relevant for folding, stability, and aggregation. To detect them by nuclear magnetic resonance (NMR) spectroscopy, the denatured state must be stabilized by chemical agents or changes in temperature. This makes the environment different from that experienced in biologically relevant processes. Using high-resolution heteronuclear NMR spectroscopy, we have characterized several denatured states of a monomeric variant of HIV-1 protease, which is natively structured in water, induced by different concentrations of urea, guanidinium chloride, and acetic acid. We have extrapolated the chemical shifts and the relaxation parameters to the denaturant-free denatured state at native conditions, showing that they converge to the same values. Subsequently, we characterized the conformational properties of this biologically relevant denatured state under native conditions by advanced molecular dynamics simulations and validated the results by comparison to experimental data. We show that the denatured state of HIV-1 protease under native conditions displays rich patterns of transient native and non-native structures, which could be of relevance to its guidance through a complex folding process., (© 2021 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

Diverse phylogenetic methods require a substitution model of evolution that should mimic, as accurately as possible, the real substitution process. At the protein level, empirical substitution models have traditionally been based on a large number of different proteins from particular taxonomic levels. However, these models assume that all of the proteins of a taxonomic level evolve under the same substitution patterns. We believe that this assumption is highly unrealistic and should be relaxed by considering protein-specific substitution models that account for protein-specific selection processes. In order to test this hypothesis, we inferred and evaluated four new empirical substitution models for the protease and integrase of HIV and other viruses. We found that these models more accurately fit, compared with any of the currently available empirical substitution models, the evolutionary process of these proteins. We conclude that evolutionary inferences from protein sequences are more accurate if they are based on protein-specific substitution models rather than taxonomic-specific (generalist) substitution models. We also present four new empirical substitution models of protein evolution that could be useful for phylogenetic inferences of viral protease and integrase.

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

HIV-1 protease (PR) is thought to be efficient targets of anti-AIDS drug design. Molecular dynamics (MD) simulations and multiple post-processing analysis technologies were applied to decipher molecular mechanism underlying binding of three drugs Lopinavir (LPV), Nelfinavir (NFV) and Atazanavir (ATV) to the PR. Binding free energies calculated by molecular mechanics generalized Born surface area (MM-GBSA) suggest that compensation between binding enthalpy and entropy plays a vital role in binding of drugs to PR. Dynamics analyses show that binding of LPV, NFV and ATV highly affects structural flexibility, motion modes and dynamics behaviour of the PR, especially for two flaps. Computational alanine scanning and interaction network analysis verify that although three drugs have structural difference, they share similar binding modes to the PR and common interaction clusters with the PR. The current findings also confirm that residues located interaction clusters, such as Asp25/Asp25', Gly27/Gly27', Ala28/Ala28', Asp29, Ile47/Ile47', Gly49/Gly49', Ile50/Ile50', Val82/Val82' and Ile84/Ile84, can be used as efficient targets of clinically available inhibitors towards the PR.

Drug resistance is a serious problem for controlling the HIV/AIDS pandemic. Current antiviral drugs show several orders of magnitude worse inhibition of highly resistant clinical variant PRS17 of HIV-1 protease compared with wild-type protease. We have analyzed the effects of a common resistance mutation G48V in the flexible flaps of the protease by assessing the revertant PRS17 V48G for changes in enzyme kinetics, inhibition, structure, and dynamics. Both PRS17 and the revertant showed about 10-fold poorer catalytic efficiency than wild-type enzyme (0.55 and 0.39 μM -1 min -1 compared to 6.3 μM -1 min -1 ). Clinical inhibitors, amprenavir and darunavir, showed 2-fold and 8-fold better inhibition, respectively, of the revertant than of PRS17, although the inhibition constants for PRS17 V48G were still 25 to 1,200-fold worse than for wild-type protease. Crystal structures of inhibitor-free revertant and amprenavir complexes with revertant and PRS17 were solved at 1.3-1.5 Å resolution. The amprenavir complexes of PRS17 V48G and PRS17 showed no significant differences in the interactions with inhibitor, although changes were observed in the conformation of Phe53 and the interactions of the flaps. The inhibitor-free structure of the revertant showed flaps in an open conformation, however, the flap tips do not have the unusual curled conformation seen in inhibitor-free PRS17. Molecular dynamics simulations were run for 1 μs on the two inhibitor-free mutants and wild-type protease. PRS17 exhibited higher conformational fluctuations than the revertant, while the wild-type protease adopted the closed conformation and showed the least variation. The second half of the simulations captured the transition of the flaps of PRS17 from a closed to a semi-open state, whereas the flaps of PRS17 V48G tucked into the active site and the wild-type protease retained the closed conformation. These results suggest that mutation G48V contributes to drug resistance by altering the conformational dynamics of the flaps., (Copyright © 2021 Elsevier Inc. All rights reserved.)

A detailed computational study was performed to investigate the conformational changes of flap region and the mechanism underlying the binding of the inhibitor TMC-126 to HIV-1 protease (PR1) and its mutant variants through molecular dynamics simulations in conjunction with the molecular mechanics Poisson-Boltzmann (MM-PBSA) free energy calculation. Further, we have studied the effectiveness of the inhibitor against HIV-2 protease (PR2). The MM-PBSA calculation suggests that TMC-126 loses its potency against mutant variants and PR2 compared to wild-type PR1 mainly due to the loss in intermolecular electrostatic interactions. The potency of the inhibitor decreases in the order: wild type PR1 > M46L > MDR20 > I50V > PR2 > V32I > A28S. Our study reveals that the flap of PR1 adopts a semi-open conformation due to the mutation I50V or MDR20. The dissimilar nature of the movement of the flap tip of both monomers is evident from the dynamic cross-correlation map. The protein structural network analysis displays that mutation causes structural rearrangements and changes the communication path between residues. Overall, we believe our study may help explore and accelerate the development of novel HIV-1/HIV-2 protease inhibitors with better potency.

Non-B strains cause nearly 90% of the worldwide human immunodeficiency virus (HIV) infections. At the same time, the protease inhibitors (PIs) were designed for subtype B. Therefore, the use of PIs in the non-B subtype context requires further investigation. Herein, we have investigated the effectiveness of currently used four PIs, namely atazanavir, darunavir, lopinavir, and tipranavir against subtype CRF01_AE (PR CRF ) by employing the MD/MMPBSA (molecular dynamics/molecular mechanics Poisson-Boltzmann surface area) scheme. Our investigation reveals that tipranavir is the most potent inhibitor against PR CRF while the other three PIs display a similar binding affinity. The energetic penalty arises due to the desolvation of polar groups always disfavor the association between PR CRF and PI, and this contribution is the least in the case of tipranavir/PR CRF compared to the other three PI-PR CRF complexes resulting in a better binding affinity for tipranavir. Further, it is revealed that the primary interaction controlling the binding of inhibitors with PR CRF is the van der Waals forces. The dynamic cross-correlation map and principal component analysis show that the anti-correlated motion at the flap region of PR CRF is diminished after the ligand binding. Further, our studies indicate that D25' forms a stable H-bond with darunavir, lopinavir, and tipranavir, while D25 forms a stable H-bond with atazanavir. The per-residue based decomposition of free energy reveals the actual residual origin of the binding free energy and identify the hotspot residues. Overall, the data presented in this study can guide the computer-aided rational design of more potent drugs targetting HIV-1 PR CRF .Communicated by Ramaswamy H. Sarma.

In the present study, we newly synthesized four types of novel fullerene derivatives: pyridinium/ethyl ester-type derivatives 3b-3l, pyridinium/carboxylic acid-type derivatives 4a, 4e, 4f, pyridinium/amide-type derivative 5a, and pyridinium/2-morpholinone-type derivative 6a. Among the assessed compounds, cis-3c, cis-3d, trans-3e, trans-3h, cis-3l, cis-4e, cis-4f, trans-4f, and cis-5a were found to inhibit HIV-1 reverse transcriptase (HIV-RT), HIV-1 protease (HIV-PR), and HCV NS5B polymerase (HCV NS5B), with IC 50 values observed in the micromolar range. Cellular uptake of pyridinium/ethyl ester-type derivatives was higher than that of corresponding pyridinium/carboxylic acid-type derivatives and pyridinium/amide-type derivatives. This result might indicate that pyridinium/ethyl ester-type derivatives are expected to be lead compounds for multitargeting drugs to treat HIV/HCV coinfection., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

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., Competing Interests: The authors have declared that no competing interests exist.

Initially, the SARS-CoV-2 virus was emerged from Wuhan, China and rapidly spreading across the world and urges the scientific community to develop antiviral therapeutic agents. Among several strategies, drug repurposing will help to react immediately to overcome the COVID-19 pandemic. In the present study, we have chosen two clinical trial drugs against HIV-1 protease namely, TMB607 and TMC310911 to use as the inhibitors of SARS-CoV-2 main protease (M pro ) enzyme. To make use of these two inhibitors as the repurposed drugs for COVID-19, it is essential to know the molecular basis of the binding mechanism of these two molecules with the SARS-CoV-2 M pro . To understand the binding mechanism, we have performed molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations against the SARS-CoV-2 M pro . The docking results indicate that both molecules form intermolecular interactions with the active site amino acids of M pro enzyme. However, during the MD simulations, TMB607 forms strong interaction with the key amino acids of M pro , and remains intact. The RMSD and RMSF values of both complexes were stable throughout the MD simulations. The MM-GBSA binding free energy values of both complexes are -43.7 and -34.9   kcal/mol, respectively. This in silico study proves that the TMB607 molecule binds strongly with the SARS-CoV-2 M pro enzyme and it may be suitable for the drug repurposing of COVID-19 and further drug designing.Communicated by Ramaswamy H. Sarma.