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.

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

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

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

Background: In most parts of the world, especially in underdeveloped countries, acquired immunodeficiency syndrome (AIDS) still remains a major cause of death, disability, and unfavorable economic outcomes. This has necessitated intensive research to develop effective therapeutic agents for the treatment of human immunodeficiency virus (HIV) infection, which is responsible for AIDS. Peptide cleavage by HIV-1 protease is an essential step in the replication of HIV-1. Thus, correct and timely prediction of the cleavage site of HIV-1 protease can significantly speed up and optimize the drug discovery process of novel HIV-1 protease inhibitors. In this work, we built and compared the performance of selected machine learning models for the prediction of HIV-1 protease cleavage site utilizing a hybrid of octapeptide sequence information comprising bond composition, amino acid binary profile (AABP), and physicochemical properties as numerical descriptors serving as input variables for some selected machine learning algorithms. Our work differs from antecedent studies exploring the same subject in the combination of octapeptide descriptors and method used. Instead of using various subsets of the dataset for training and testing the models, we combined the dataset, applied a 3-way data split, and then used a "stratified" 10-fold cross-validation technique alongside the testing set to evaluate the models., Results: Among the 8 models evaluated in the "stratified" 10-fold CV experiment, logistic regression, multi-layer perceptron classifier, linear discriminant analysis, gradient boosting classifier, Naive Bayes classifier, and decision tree classifier with AUC, F-score, and B. Acc. scores in the ranges of 0.91-0.96, 0.81-0.88, and 80.1-86.4%, respectively, have the closest predictive performance to the state-of-the-art model (AUC 0.96, F-score 0.80 and B. Acc. ~ 80.0%). Whereas, the perceptron classifier and the K-nearest neighbors had statistically lower performance (AUC 0.77-0.82, F-score 0.53-0.69, and B. Acc. 60.0-68.5%) at p < 0.05. On the other hand, logistic regression, and multi-layer perceptron classifier (AUC of 0.97, F-score > 0.89, and B. Acc. > 90.0%) had the best performance on further evaluation on the testing set, though linear discriminant analysis, gradient boosting classifier, and Naive Bayes classifier equally performed well (AUC > 0.94, F-score > 0.87, and B. Acc. > 86.0%)., Conclusions: Logistic regression and multi-layer perceptron classifiers have comparable predictive performances to the state-of-the-art model when octapeptide sequence descriptors consisting of AABP, bond composition and standard physicochemical properties are used as input variables. In our future work, we hope to develop a standalone software for HIV-1 protease cleavage site prediction utilizing the linear regression algorithm and the aforementioned octapeptide sequence descriptors., (© 2022. The Author(s).)

The global HIV/AIDS epidemic still currently affects approximately 38 million individuals globally. The protease enzyme of the human immunodeficiency virus is a major drug target in antiviral therapy, however, under the influence of reverse transcriptase and in the context of drug pressure, the rapid PR mutation rate contributes significantly to clinical failure. The set of cooperative non-active site mutations, I13V/I62V/V77I, have been associated with reduced inhibitor susceptibility and are the focus of the current study. When compared to the wild-type protease the mutant protease exhibited decreased binding affinities towards ATV and DRV by 64- and 12-fold, respectively, and decreased the overall favourable Gibbs free energy for ATV, DRV, RTV and SQV. Moreover, these mutations decreased the thermal stability of the protease when in complex with ATV and DRV by approximately 6.4 and 4.2 °C, respectively. The crystal structure of the mutant protease revealed that the location of these mutations and their effect on the hydrophobic sliding mechanism may be crucial in their role in resistance., (Copyright © 2022 Elsevier B.V. All rights reserved.)

The emergence of multidrug resistant (MDR) HIV strains severely reduces the effectiveness of antiretroviral therapy. Clinical inhibitor darunavir (1) has picomolar binding affinity for HIV-1 protease (PR), however, drug resistant variants like PR S17 show poor inhibition by 1, despite the presence of only two mutated residues in the inhibitor-binding site. Antiviral inhibitors that target MDR proteases like PR S17 would be valuable as therapeutic agents. Inhibitors 2 and 3 derived from 1 through substitutions at P1, P2 and P2' positions exhibit 3.4- to 500-fold better inhibition than clinical inhibitors for PR S17 with the exception of amprenavir. Crystal structures of PR S17 /2 and PR S17 /3 reveal how these inhibitors target the two active site mutations of PR S17 . The substituted methoxy P2 group of 2 forms new interactions with G48V mutation, while the modified bis-fluoro-benzyl P1 group of 3 forms a halogen interaction with V82S mutation, contributing to improved inhibition of PR S17 ., Competing Interests: Declaration of competing interest None declared., (Copyright © 2021 Elsevier Inc. All rights reserved.)

A novel class of HIV-1 protease inhibitors with flexible piperidine as the P2 ligand was designed with the aim of improving extensive interactions with the active subsites. Many inhibitors exhibited good to excellent inhibitory effect on enzymatic activity and viral infectivity. In particular, inhibitor 3a with (R)-piperidine-3-carboxamide as the P2 ligand and 4-methoxybenzenesulfonamide as the P2' ligand showed an enzyme K i value of 29 pM and antiviral IC 50 value of 0.13 nM, more than six-fold enhancement of activity compared to DRV. Furthermore, there was no significant change in potency against DRV-resistant mutations and HIV-1 NL4-3 variant for 3a. Besides, inhibitor 3a exhibited potent antiviral activity against subtype C variants with low nanomole EC 50 values. In addition, the molecular modeling revealed important hydrogen bonds and other favorable van der Waals interactions with the backbone atoms of the protease and provided insight for designing and optimizing more potent HIV-1 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 © 2021 Elsevier Masson SAS. All rights reserved.)

A method for predicting HIV drug resistance by using genotypes would greatly assist in selecting appropriate combinations of antiviral drugs. Models reported previously have had two major problems: lack of information on the 3D protein structure and processing of incomplete sequencing data in the modeling procedure. We propose obtaining the 3D structural information of viral proteins by using homology modeling and molecular field mapping, instead of just their primary amino acid sequences. The molecular field potential parameters reflect the physicochemical characteristics associated with the 3D structure of the proteins. We also introduce the Bayesian conditional mutual information theory to estimate the probabilities of occurrence of all possible protein candidates from an incomplete sequencing sample. This approach allows for the effective use of uncertain information for the modeling process. We applied these data analysis techniques to the HIV-1 protease inhibitor dataset and developed drug resistance prediction models with reasonable performance., Competing Interests: The authors have declared that no competing interests exist.

The protonation states of two aspartic acids in the catalytic strands of HIV-1 protease (PR) remarkably affect bindings of inhibitors to PR. It is requisite for the design of potent inhibitors towards PR to investigate the influences of Asp25/Asp25' protonated states on dynamics behaviour of PR and binding mechanism of inhibitors to PR. In this work, molecular dynamics (MD) simulations, MM-GBSA method and principal component (PC) analysis were coupled to explore the effect of Asp25/Asp25' protonation states on conformational changes of PR and bindings of Amprenavir and MKP97 to PR. The results show that the Asp25/Asp25' protonation states exert different impacts on structural fluctuations, flexibility and motion modes of PR. Dynamics analysis verifies that Asp25/Asp25' protonated states highly affect conformational dynamics of two flaps in PR. The binding free energy calculations results suggest that the Asp25/Asp25' protonated states obviously strengthen bindings of inhibitors to PR compared to the non-protonation state. Calculations of residue-based free energy decomposition indicate that the Asp25/Asp25' protonation not only disturbs the interaction network of inhibitors with PR but also stabilizes bindings of inhibitors to PR by cancelling the electrostatic repulsive interaction. Therefore, special attentions should be paid to the Asp25/Asp25' protonation in the design of potent inhibitors towards PR.

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

Drug resistance threatens many critical therapeutics through mutations in the drug target. The molecular mechanisms by which combinations of mutations, especially those remote from the active site, alter drug binding to confer resistance are poorly understood and thus difficult to counteract. A machine learning strategy was developed that coupled parallel molecular dynamics simulations with experimental potency to identify specific conserved mechanisms underlying resistance. Physical features were extracted from the simulations, analyzed, and integrated into one consistent and interpretable elastic network model. To rigorously test this strategy, HIV-1 protease variants with diverse mutations were used, with potencies ranging from picomolar to micromolar to the drug darunavir. Feature reduction resulted in a model with four specific features that predicts for both the training and test sets inhibitor binding free energy within 1 kcal/mol of the experimental value over this entire range of potency. These predictive features are physically interpretable, as they vary specifically with affinity and diagonally transverse across the protease homodimer. This physics-based strategy of parallel molecular dynamics and machine learning captures mechanisms by which complex combinations of mutations confer resistance and identify critical features that serve as bellwethers of affinity, which will be critical in future drug design.

A clinically-relevant, drug-resistant mutant of HIV-1 protease (PR), termed Flap+ (I54V) and containing L10I, G48V, I54V and V82A mutations, is known to produce significant changes in the entropy and enthalpy balance of drug-PR interactions, compared to wild-type PR. A similar mutant, Flap+ (I54A) , which evolves from Flap+ (I54V) and contains the single change at residue 54 relative to Flap+ (I54V) , does not. Yet, how Flap+ (I54A) behaves in solution is not known. To understand the molecular basis of V54A evolution, we compared nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy, isothermal titration calorimetry, and enzymatic assay data from four PR proteins: PR (pWT), Flap+ (I54V) , Flap+ (I54A) , and Flap+ (I54) , a control mutant that contains only L10I, G48V and V82A mutations. Our data consistently show that selection to the smaller side chain at residue 54, not only decreases inhibitor affinity, but also restores the catalytic activity., (© 2020 The Protein Society.)

In the present study, we newly synthesized three types of novel fullerene derivatives: pyridinium-type derivatives trans-3a and 4a-5b, piperidinium-type derivative 9, and proline-type derivatives 10a-12. Among the assessed compounds, 5a, 10e, 10f, 10i, 11a-d, and 12 were found to inhibit both HIV reverse transcriptase and HIV protease (HIV-PR), with IC 50 values in the low micromolar range being observed. Regarding HIV-PR inhibition activity, proline-type derivatives 11a-11d and 12, bearing an alkyl chain between the hydroxylmethylcarbonyl (HMC) moiety and pyrrolidine ring, were more potent than other derivatives. This result might indicate that connecting HMC moieties with proline-type fullerene derivatives through properly sized alkyl chain leads to improved HIV-PR inhibitory activity., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Designing drugs that directly interact with multiple targets is a promising approach for treating complicated diseases. In order to successfully bind to multiple targets of different families and achieve the desired ligand efficiency, multi-target-directed ligands (MTDLs) require a higher level of diversity and complexity. De novo design strategies for creating more diverse chemical entities with desired properties may present an improved approach for developing MTDLs. In this chapter, we describe a computational protocol for developing MTDLs using the first reported multi-target de novo program, LigBuilder 3, which combines a binding site prediction module with de novo drug design and optimization modules. As an illustration of each detailed procedure, we design dual-functional compounds of two well-characterized virus enzymes, HIV protease and reverse transcriptase (PR and RT, respectively), using fragments extracted from known inhibitors. LigBuilder 3 is accessible at http://www.pkumdl.cn/ligbuilder3/ .

The novel coronavirus (SARS-CoV-2) has infected several million people and caused thousands of deaths worldwide since December 2019. As the disease is spreading rapidly all over the world, it is urgent to find effective drugs to treat the virus. The main protease (Mpro) of SARS-CoV-2 is one of the potential drug targets. Therefore, in this context, we used rigorous computational methods, including molecular docking, fast pulling of ligand (FPL), and free energy perturbation (FEP), to investigate potential inhibitors of SARS-CoV-2 Mpro. We first tested our approach with three reported inhibitors of SARS-CoV-2 Mpro, and our computational results are in good agreement with the respective experimental data. Subsequently, we applied our approach on a database of ∼4600 natural compounds, as well as 8 available HIV-1 protease (PR) inhibitors and an aza-peptide epoxide. Molecular docking resulted in a short list of 35 natural compounds, which was subsequently refined using the FPL scheme. FPL simulations resulted in five potential inhibitors, including three natural compounds and two available HIV-1 PR inhibitors. Finally, FEP, the most accurate and precise method, was used to determine the absolute binding free energy of these five compounds. FEP results indicate that two natural compounds, cannabisin A and isoacteoside, and an HIV-1 PR inhibitor, darunavir, exhibit a large binding free energy to SARS-CoV-2 Mpro, which is larger than that of 13b , the most reliable SARS-CoV-2 Mpro inhibitor recently reported. The binding free energy largely arises from van der Waals interaction. We also found that Glu166 forms H-bonds to all of the inhibitors. Replacing Glu166 by an alanine residue leads to ∼2.0 kcal/mol decreases in the affinity of darunavir to SARS-CoV-2 Mpro. Our results could contribute to the development of potential drugs inhibiting SARS-CoV-2.

HIV infection is a global health epidemic with current FDA-approved HIV-1 Protease inhibitors (PIs) designed against subtype B protease, yet they are used in HIV treatment world-wide regardless of patient HIV classification. In this study, double electron-electron resonance (DEER) electron paramagnetic resonance (EPR) spectroscopy was utilized to gain insights in how natural polymorphisms in several African and Brazilian protease (PR) variants affect the conformational landscape both in the absence and presence of inhibitors. Findings show that Subtypes F and H HIV-1 PR adopt a primarily closed conformation in the unbound state with two secondary mutations, D60E and I62V, postulated to be responsible for the increased probability for closed conformation. In contrast, subtype D, CRF&#95;AG, and CRF&#95;BF HIV-1 PR adopt a primarily semi-open conformation, as observed for PI-naïve-subtype B when unbound by substrate or inhibitor. The impact that inhibitor binding has on shifting the conformational land scape of these variants is also characterized, where analysis provides classification of inhibitor induced shifts away from the semi-open state into weak, moderate and strong effects. The findings are compared to those for prior studies of inhibitor induced conformational shifts in PI-naïve Subtype B, C and CRF&#95;AE., 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 © 2020 Elsevier Inc. All rights reserved.)

Newly designed HIV-1 protease inhibitors that maximize interactions with the protein backbone, especially in the form of hydrogen bonds, may enhance the antiviral potency of these compounds and minimize acquisition of drug-resistant mutations. Herein, we described a series of new HIV-1 PIs containing phenols as the P2 ligands and chiral isopropanol as the P1' ligands, in combination with 4-trifluoromethylphenylsulfonamide or 4-nitrophenylsulfonamide as the P2' ligands. And most of these compounds exhibited nanomolar inhibitory potency. In particular, inhibitors 13c and 13e with 4-trifluoromethylphenylsulfonamide as the P2' ligand and (R) - isopropanol as the P1' ligand, exhibited antiviral IC 50 values of 1.64 nM and 2.33 nM, respectively. Furthermore, they also showed remarkable activity against wild-type and DRV-resistant HIV-1 variants that raised the prospect of designing more effective PIs further., 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 © 2020 Elsevier Ltd. All rights reserved.)

The design, synthesis, and X-ray structural analysis of hybrid HIV-1 protease inhibitors (PIs) containing bis-tetrahydrofuran (bis-THF) in a pseudo-C 2 -symmetric dipeptide isostere are described. A series of PIs were synthesized by incorporating bis-THF of darunavir on either side of the Phe-Phe isostere of lopinavir in combination with hydrophobic amino acids on the opposite P2/P2' position. Structure-activity relationship studies indicated that the bis-THF moiety can be attached at either the P2 or P2' position without significantly affecting potency. However, the group on the opposite P2/P2' position had a dramatic effect on potency depending on the size and shape of the side chain. Cocrystal structures of inhibitors with wild-type HIV-1 protease revealed that the bis-THF moiety retained similar interactions as observed in the darunavir-protease complex regardless of the position on the Phe-Phe isostere. Analyses of cocrystal structures and molecular dynamics simulations provide insights into optimizing HIV-1 PIs containing bis-THF in non-sulfonamide dipeptide isosteres.

A major challenge in computer-aided drug design is the accurate estimation of ligand binding affinity. Here, a new approach that combines the adaptive steered molecular dynamics (ASMD) and partial atomic charges calculated by semi-empirical quantum mechanics (SQMPC), namely ASMD-SQMPC, is suggested to predict the ligand binding affinities, with 24 HIV-1 protease inhibitors as testing examples. In the ASMD-SQMPC, the relative binding free energy (ΔG) is reflected by the average maximum potential of mean force ( max ) between bound and unbound states. The correlation coefficient (R 2 ) between the max and experimentally determined ΔG is 0.86, showing a significant improvement compared with the conventional ASMD (R 2 = 0.52). Therefore, this study provides an efficient approach to predict the relative ΔG and reveals the significance of precise partial atomic charges in the theoretical simulations., (© 2020 Wiley Periodicals, Inc.)

HIV-1 nucleocapsid (NC) becomes an attractive target for the development of novel anti-HIV-1 agents. Discovering of non-antibody scaffolds that disrupt the function of NC will be a potential aspect for disturbing viral maturation process. Correspondingly, we explored the specific binding site of the thermoresistant-scaffold protein, αRep9A8 which formerly demonstrated the inhibitory effect on HIV-1 replication. The portion of Gag, CA21-SP1-NC has been used as a template for designing nine overlapping peptides (P4-P12). The P9 peptide showed the strongest binding activity followed by P8 and P12 respectively. The amino acid sequences on those peptides resemble the N-terminal domain of the NC proximity to the SP1-NC initial cleavage site and across the conserved CCHC zinc finger 1 (ZF1) of NC. The interaction KD between αRep9A8 with its target was 224.9 ± 57.4 nM. Consequently, αRep9A8 demonstrated the interference of the HIV-1 protease function by hindering a protease cleavage site. The released NC product from CA21-SP1-NC was diminished. The present study provided an additional information of αRep9A8 function in interfering of viral maturation processes resulting in the decremental efficiency of viral infectivity., (© 2020 The Author(s).)

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

Aim: HIV-1 protease inhibitors regimens suffered from a number of drawbacks, among which, the most egregious issue was the growing emergence of drug-resistant strains. Materials & methods: The design strategy of maximizing the protease active site interactions with the inhibitor, especially promoting extensive hydrogen bonding with the protein backbone atoms, might be in favor of combating drug resistance. A series of HIV-1 protease inhibitors that incorporated enantiomeric isopropanols as the P1' ligands in combination with phenols as the P2 ligands were reported herein. Results: A number of inhibitors displayed potent protease enzyme inhibition activity. In particular, inhibitor 14c showed comparable potency as darunavir with IC 50 value of 1.91 nM and activity against darunavir-resistant HIV-1 variants. Conclusion: The new kind of HIV-1 protease inhibitors deserves further study.

Obtaining reactivity information from the molecular electronic structure of a chemical system is a computationally intensive process. As a way of probing reactivity information around that, there exist electron density response variables, such as the Fukui functions (FFs), which are well-established descriptors that summarize the local susceptibility to react. These properties only require few single-point quantum chemical calculations, but even then, the intrinsic high cost and unfavorable computational complexity with respect to the number of atoms in the system makes this approach available only to small fragments and systems. In this study, we explore the computation of FFs, showing that semiempirical quantum chemical methods can be used to obtain the reactivity information equivalent to that of a Density Functional Theory (DFT) functional, for the eight entire polypeptide chains. The combination of semiempirical methods with the frozen orbital approximation allows for the obtention of these reactivity descriptors for biological systems with reasonable accuracy and speed, unlocking the utilization of these methods for such systems. These results for the frozen orbital approximation can be additionally improved when other molecular orbitals from the frontier band are employed in the computation. We also show the potential of this computational protocol in the ligand-protein complexes of HIV-1 protease, predicting which of those ligands are active inhibitors., (© 2020 Wiley Periodicals, Inc.)

A series of novel HIV-1 protease inhibitors has been designed and synthesized, which contained morpholine derivatives as the P2 ligands and hydrophobic cyclopropyl as the P1' ligand at the meantime in this study, with the aim of improving the interactions between the active sites of HIV-1 protease and the inhibitors. Twenty-eight compounds were synthesized and assessed, among which inhibitors m18 and m1 exhibited excellent inhibitory effect on the activity of HIV-1 protease with IC 50 value of 47 nM and 53 nM, respectively. The molecular modeling of m1 revealed possible hydrogen bondings or van der Waals between the inhibitor and the protease, worthy of in-depth study., 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 © 2020 Elsevier Ltd. All rights reserved.)

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations., (© 2019 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals, Inc.)

Over the past several decades, atomistic simulations of biomolecules, whether carried out using molecular dynamics or Monte Carlo techniques, have provided detailed insights into their function. Comparing the results of such simulations for a few closely related systems has guided our understanding of the mechanisms by which changes such as ligand binding or mutation can alter the function. The general problem of detecting and interpreting such mechanisms from simulations of many related systems, however, remains a challenge. This problem is addressed here by applying supervised and unsupervised machine learning techniques to a variety of thermodynamic observables extracted from molecular dynamics simulations of different systems. As an important test case, these methods are applied to understand the evasion by human immunodeficiency virus type-1 (HIV-1) protease of darunavir, a potent inhibitor to which resistance can develop via the simultaneous mutation of multiple amino acids. Complex mutational patterns have been observed among resistant strains, presenting a challenge to developing a mechanistic picture of resistance in the protease. In order to dissect these patterns and gain mechanistic insight into the role of specific mutations, molecular dynamics simulations were carried out on a collection of HIV-1 protease variants, chosen to include highly resistant strains and susceptible controls, in complex with darunavir. Using a machine learning approach that takes advantage of the hierarchical nature in the relationships among the sequence, structure, and function, an integrative analysis of these trajectories reveals key details of the resistance mechanism, including changes in the protein structure, hydrogen bonding, and protein-ligand contacts.

Aim: This letter investigates the role of radial distribution function-based descriptors for in silico design of new drugs. Methodology: The multiple linear regression models for HIV-1 protease and its complexes with a series of inhibitors were constructed. A detailed analysis of major atomic contributions to the radial distribution function descriptor weighted by the number of valence shell electrons identified residues Arg8, Asp29 and residues of the catalytic triad as crucial for the correlation with the inhibition constant, together with residues Asp30 and Ile50, whose mutations are known to cause an emergence of drug resistant variants. Conclusion: This study demonstrates an easy and fast assessment of the activity of potential drugs and the derivation of structural information of their complexes with the receptor or enzyme.

Since dual inhibitors may yield lower toxicity and reduce the likelihood of drug resistance, as well as inhibitors of HIV-1 PR and RT constitute the core of chemotherapy for AIDS treatment, we herein designed and synthesized new coumarin derivatives characterized by various linkers that exhibited excellent potency against PR and a weak inhibition of RT. Among which, compounds 6f and 7c inhibited PR with IC 50 values of 15.5 and 62.1 nM, respectively, and weakly affected also RT with IC 50 values of 241.8 and 188.7 μM, respectively, showing the possibility in the future of developing dual HIV-1 PR/RT inhibitors. Creative stimulation for further research of more potent dual HIV-1 inhibitors was got according to the molecular docking studies., (Copyright © 2019 Elsevier Masson SAS. All rights reserved.)

In this study, we isolated nine compounds from the acid hydrolysate of the flower buds of Lonicera fulvotomentosa Hsu et S. C. Cheng and characterized their chemical structures using 1 H-NMR, 13 C-NMR, and electron ionization mass spectroscopy (EI-MS). These compounds were identified as β-sitosterol ( 1 ), 5,5'-dibutoxy-2,2'-bifuran ( 2 ), nonacosane-10-ol ( 3 ), ethyl (3β)-3,23-dihydroxyolean-12-en-28-oate ( 4 ), oleanolic acid ( 5 ), ethyl caffeate ( 6 ), caffeic acid ( 7 ), isovanillin ( 8 ), and hederagenin ( 9 ), with 4 as a new triterpene compound. Inhibitory activity against human immunodeficiency virus (HIV) protease was also evaluated for the compounds, and only ethyl caffeate, caffeic acid, and isovanillin ( 6 , 7 , and 8 ) exhibited inhibitory effects, with IC 50 values of 1.0 μM, 1.5 μM, and 3.5 μM, respectively. Molecular docking with energy minimization and subsequent molecular dynamic (MD) simulation showed that ethyl caffeate and caffeic acid bound to the active site of HIV protease, while isovanillin drifted out from the active site and dissociated into bulk water during MD simulations, and most of the binding residues of HIV protease have been previously identified for HIV protease inhibitors. These results suggest that caffeic acid derivatives may possess inhibitory activities towards HIV protease other than previously reported inhibitory activities against HIV integrase, and thus ethyl caffeate and caffeic acid could be used as lead compounds in developing potential HIV protease inhibitors, and possibly even dual-function inhibitors against HIV.

Background: Acquired immunodeficiency syndrome (AIDS) is caused by the Human immunodeficiency virus type-1 (HIV-1). HIV-1 protease (HIV-1 PR) is an essential enzyme for the HIV replication, and therefore, it is an important target for antiretroviral drugs development, particularly from natural products. Auricularia polytricha (AP) is an edible mushroom with several important therapeutic properties. These properties will be investigated as HIV-1 PR inhibitors., Methods: The sequential hexane (APH), ethanol (APE) and water (APW) extracts from AP were screened for inhibitory activity against HIV-1 PR. The extract that consistently showed the strong HIV-1 PR inhibition was further investigated for its phytochemical constituents. The compounds were purified by column chromatography. The isolated compounds were structurally elucidated using 1D and 2D NMR, HRMS, FTIR, and GC/MS techniques. Each compound was screened against HIV-1 PR to determine its inhibitory activity and to provide an explanation for the activity found in the extract., Results: Hexane crude extract of AP (APH) exhibited significant inhibition on HIV-1 PR activity. Four major compounds isolated from APH fraction were identified to be two triacylglycerols, linoleic acid and ergosterol. Moreover, all four compounds showed significant inhibition of HIV-1 PR activity., Conclusion: The findings from this study suggest that AP is a good source of fatty esters, fatty acids and ergosterol. These natural products exhibit anti-HIV-1 properties by blocking HIV-1 PR. These important biological results warrant further development of AP as an alternative antiretroviral drug.

Herein, we report the effect of nine FDA approved protease inhibitor drugs against a new HIV-1 subtype C mutant protease, E35D↑G↑S. The mutant has five mutations, E35D, two insertions, position 36 (G and S), and D60E. Kinetics, inhibition constants, vitality, Gibbs free binding energies are reported. The variant showed a decreased affinity for substrate and low catalytic efficiency compared to the wild type. There was a significant decrease in the binding of seven FDA approved protease inhibitors against the mutant ( p  < .0001). Amprenavir and ritonavir showed the least decrease, but still significant reduced activity in comparison to the wildtype (4 and 5 folds, respectively, p  = .0021 and .003, respectively). Nelfinavir and atazanavir were the worst inhibitors against the variant as seen from the IC 50 , with values of 1401 ± 3.0 and 685 ± 3.0 nM, respectively. Thermodynamics data showed less favourable Gibbs free binding energies for the protease inhibitors to the mutant.

We report the synthesis and biological evaluation of phenylcarboxylic acid and phenylboronic acid containing HIV-1 protease inhibitors and their functional effect on enzyme inhibition and antiviral activity in MT-2 cell lines. Inhibitors bearing bis-THF ligand as P2 ligand and phenylcarboxylic acids and carboxamide as the P2' ligands, showed very potent HIV-1 protease inhibitory activity. However, carboxylic acid containing inhibitors showed very poor antiviral activity relative to carboxamide-derived inhibitors which showed good antiviral IC 50 value. Boronic acid derived inhibitor with bis-THF as the P2 ligand showed very potent enzyme inhibitory activity, but it showed lower antiviral activity than darunavir in the same assay. Boronic acid containing inhibitor with a P2-Crn-THF ligand also showed potent enzyme K i but significantly decreased antiviral activity. We have evaluated antiviral activity against a panel of highly drug-resistant HIV-1 variants. One of the inhibitors maintained good antiviral activity against HIV DRV R P20 and HIV DRV R P30 viruses. We have determined high resolution X-ray structures of two synthetic inhibitors bound to HIV-1 protease and obtained molecular insight into the ligand-binding site interactions., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

HIV-1 protease (HIV-1-pr) plays an important role in viral replication and maturation, making it one of the most attractive targets for anti-retroviral therapy. To design new effective inhibitors able to combat drug resistance in mutant HIV-1-pr variants, it is essential to gain further understanding about the mechanisms by which the recently proposed inhibitors deactivate the mutant HIV-1-pr variants. In the present work, we explored the interactions between two P2-ligands (DRV, and one new derivative, 4UY) with wild type (WT) and two multiple mutant HIV-1-pr variants (p20 and p51) with all atom molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA). The trajectories of MD simulations show that both 4UY and DRV primarily bind with the active sites, flap and 80s loop regions of HIV-1-pr variants through either hydrogen bonds or hydrophobic interactions. More hydrogen bonds and hydrophobic interactions were located for 4UY/HIV-1-pr complexes than for DRV/HIV-1-pr counterparts. More importantly, 4UY was found to have an extra hydrogen bond with the backbone of Gly48' in the flap region of the HIV-1-prs. The flap tip-tip distance (I50-I50') and flap tip-active site distance (I50-D25 and I50'-D25') indicate that the flaps turn more closed in 4UY bound HIV-1-prs than DRV bound ones, and the former also have more compact hydrophobic cavities than the latter. Further, the vector projections from PCA indicate that 4UY/DRV inhibitor binding projects the closing of flap in HIV-1-pr variants. In line with the above trajectory analysis, the thermodynamics calculation with MM-PBSA method suggests much stronger binding affinity for 4UY/HIV-1-pr than DRV/HIV-1-pr by 4.3-6.4 kcal/mol. Although p20 and p51 also induce weaker binding due to multiple mutants for 4UY inhibitor by 1.9-1.8 kcal/mol, their bindings to the new P2 ligand (4UY) are indeed significantly enhanced as compared to DRV. The thermodynamic components responsible for the binding differences and the contribution from key residues to the binding were also discussed in detail., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Drug-resistance threatens effective treatment of HIV/AIDS. Clinical inhibitors, including darunavir (1), are ineffective for highly resistant protease mutant PR20, however, antiviral compound 2 derived from 1 with fused tricyclic group at P2, extended amino-benzothiazole P2' ligand and two fluorine atoms on P1 shows 16-fold better inhibition of PR20 enzyme activity. Crystal structures of PR20 and wild-type PR complexes reveal how the extra groups of 2 counteract the expanded ligand-binding pocket, dynamic flaps, and faster dimer dissociation of PR20., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Discovery and optimization of small molecule inhibitors as therapeutic drugs have immensely benefited from rational structure-based drug design. With recent advances in high-resolution structure determination, computational power, and machine learning methodology, it is becoming more tractable to elucidate the structural basis of drug potency. However, the applicability of machine learning models to drug design is limited by the interpretability of the resulting models in terms of feature importance. Here, we take advantage of the large number of available inhibitor-bound HIV-1 protease structures and associated potencies to evaluate inhibitor diversity and machine learning models to predict ligand affinity. First, using a hierarchical clustering approach, we grouped HIV-1 protease inhibitors and identified distinct core structures. Explicit features including protein-ligand interactions were extracted from high-resolution cocrystal structures as 3D-based fingerprints. We found that a gradient boosting machine learning model with this explicit feature attribution can predict binding affinity with high accuracy. Finally, Shapley values were derived to explain local feature importance. We found specific van der Waals (vdW) interactions of key protein residues are pivotal for the predicted potency. Protein-specific and interpretable prediction models can guide the optimization of many small molecule drugs for improved potency.

We describe the design, synthesis, and biological evaluation of novel HIV-1 protease inhibitors containing a squaramide-derived scaffold as the P2 ligand in combination with a (R)-hydroxyethylamine sulfonamide isostere. Inhibitor 3h with an N-methyl-3-(R)-aminotetrahydrofuranyl squaramide P2-ligand displayed an HIV-1 protease inhibitory K i value of 0.51 nM. An energy minimized model of 3h revealed the major molecular interactions between HIV-1 protease active site and the tetrahydrofuranyl squaramide scaffold that may be responsible for its potent activity., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

A structure-guided design strategy was used to improve the resistance profile of HIV-1 protease inhibitors by optimizing hydrogen bonding and van der Waals interactions with the protease while staying within the substrate envelope. Stereoisomers of 4-(1-hydroxyethyl)benzene and 4-(1,2-dihydroxyethyl)benzene moieties were explored as P2' ligands providing pairs of diastereoisomers epimeric at P2', which exhibited distinct potency profiles depending on the configuration of the hydroxyl group and size of the P1' group. While compounds with the 4-(1-hydroxyethyl)benzene P2' moiety maintained excellent antiviral potency against a panel of multidrug-resistant HIV-1 strains, analogues with the polar 4-(1,2-dihydroxyethyl)benzene moiety were less potent, and only the ( R )-epimer incorporating a larger 2-ethylbutyl P1' group showed improved potency. Crystal structures of protease-inhibitor complexes revealed strong hydrogen bonding interactions of both ( R )- and ( S )-stereoisomers of the hydroxyethyl group with Asp30'. Notably, the ( R )-dihydroxyethyl group was involved in a unique pattern of direct hydrogen bonding interactions with the backbone amides of Asp29' and Asp30'. The SAR data and analysis of crystal structures provide insights for optimizing these promising HIV-1 protease inhibitors.

We describe here the syntheses of optically pure (3a S ,4 S ,7a R )-hexahydro-4 H -furo[2,3- b ]pyran-4-ol and (3a R ,4 R ,7a S )-hexahydro-4 H -furo[2,3- b ]pyran-4-ol. These stereochemically defined heterocycles are important high-affinity P2 ligands for a variety of highly potent HIV-1 protease inhibitors. The key steps involve an efficient Paternò-Büchi [2 + 2] photocycloaddition, catalytic hydrogenation, acid-catalyzed cyclization to form the racemic ligand alcohol, and an enzymatic resolution with immobilized Amano Lipase PS-30. Optically active ligands (-)- 6 and (+)- 6 were obtained with high enantiomeric purity. Enantiomer (-)- 6 has been converted to potent HIV-1 protease inhibitor 3 .

Over the last two decades, both the sensitivity of NMR and the time scale of molecular dynamics (MD) simulation have increased tremendously and have advanced the field of protein dynamics. HIV-1 protease has been extensively studied using these two methods, and has presented a framework for cross-evaluation of structural ensembles and internal dynamics by integrating the two methods. Here, we review studies from our laboratories over the last several years, to understand the mechanistic basis of protease drug-resistance mutations and inhibitor responses, using NMR and crystal structure-based parallel MD simulations. Our studies demonstrate that NMR relaxation experiments, together with crystal structures and MD simulations, significantly contributed to the current understanding of structural/dynamic changes due to HIV-1 protease drug resistance mutations.

Human immunodeficiency virus type 1 protease is essential for virus replication and maturation and has been considered as one of the important drug target for the antiretroviral treatment of HIV infection. The majority of HIV infections are caused due to non-B subtypes in developing countries. Subtype AE is spreading rapidly and infecting huge population worldwide. Understanding the interdependence of active and non-active site mutations in conferring drug resistance is crucial for the development effective inhibitors in subtype AE protease. In this work, we have investigated the mechanism of resistance against indinavir (IDV) due to therapy selected active site mutation V82F, non-active site mutations PF82V and their cooperative effects PV82F in subtype AE-protease using molecular dynamics simulations and binding free energy calculations. The simulations suggested all the three complexes lead to decrease in binding affinity of IDV, whereas the PF82V complex resulted in an enhanced binding affinity compared to V82F and PV82F complexes. Large positional deviation of IDV was observed in V82F complex. The preservation of hydrogen bonds of IDV with active site Asp25/Asp25' and flap residue Ile50/50' via a water molecule is crucial for effective binding. Owing to the close contact of 80s loop with Ile50' and Asp25, the alteration between residues Thr80 and Val82, further induces conformational change thereby resulting in loss of interactions between IDV and the residues in the active site cavity, leading to drug resistance. Our present study shed light on the effect of active, non-active site mutations and their cooperative effects in AE protease. Communicated by Ramaswamy H. Sarma.

HIV-1 protease is one of the prime targets of agents used in antiretroviral therapy against HIV. However, under selective pressure of protease inhibitors, primary mutations at the active site weaken inhibitor binding to confer resistance. Darunavir (DRV) is the most potent HIV-1 protease inhibitor in clinic; resistance is limited, as DRV fits well within the substrate envelope. Nevertheless, resistance is observed due to hydrophobic changes at residues including I50, V82, and I84 that line the S1/S1' pocket within the active site. Through enzyme inhibition assays and a series of 12 crystal structures, we interrogated susceptibility of DRV and two potent analogues to primary S1' mutations. The analogues had modifications at the hydrophobic P1' moiety compared to DRV to better occupy the unexploited space in the S1' pocket where the primary mutations were located. Considerable losses of potency were observed against protease variants with I84V and I50V mutations for all three inhibitors. The crystal structures revealed an unexpected conformational change in the flap region of I50V protease bound to the analogue with the largest P1' moiety, indicating interdependency between the S1' subsite and the flap region. Collective analysis of protease-inhibitor interactions in the crystal structures using principle component analysis was able to distinguish inhibitor identity and relative potency solely based on van der Waals contacts. Our results reveal the complexity of the interplay between inhibitor P1' moiety and S1' mutations and validate principle component analyses as a useful tool for distinguishing resistance and inhibitor potency.

The design, synthesis and SAR study of a new series of HIV-1 protease inhibitors with pentacyclic triterpenoids as P2 ligands and phenylsulfonamide as P2' ligands were discussed. These compounds exhibited micromolar inhibitory potency, among which compound T1c displayed HIV-1 protease inhibition with IC 50 values of 0.12 μM, which was 67 times the inhibitory activity of its raw material Ursolic acid (8.0 μM)., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

HIV protease is essential for processing the Gag polyprotein to produce infectious virions and is a major target in antiretroviral therapy. We have identified an unusual HIV-1 subtype C variant that contains insertions of leucine and asparagine (L38↑N↑L) in the hinge region of protease at position 38. This was isolated from a protease inhibitor naïve infant. Isothermal titration calorimetry showed that 10% less of L38↑N↑L protease was in the active conformation as compared with a reference strain. L38↑N↑L protease displayed a ±50% reduction in K M and k cat The catalytic efficiency ( k cat / K M ) of L38↑N↑L protease was not significantly different from that of wild type although there was a 42% reduction in specific activity for the variant. An in vitro phenotypic assay showed the L38↑N↑L protease to be susceptible to lopinavir (LPV), atazanavir (ATV) and darunavir in the context of an unrelated Gag. However, in the presence of the related Gag, L38↑N↑L showed reduced susceptibility to darunavir while remaining susceptible to LPV and ATV. Furthermore, a reduction in viral replication capacity (RC) was observed in combination with the related Gag. The reduced susceptibility to darunavir and decrease in RC may be due to PTAPP duplication in the related Gag. The present study shows the importance of considering the Gag region when looking at drug susceptibility of HIV-1 protease variants., (© 2019 The Author(s). Published by Portland Press Limited on behalf of the Biochemical Society.)

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

The use of antiretrovirals (ARVs) has drastically improved the life quality and expectancy of HIV patients since their introduction in health care. Several millions are still afflicted worldwide by HIV and ARV resistance is a constant concern for both healthcare practitioners and patients, as while treatment options are finite, the virus constantly adapts via complex mutation patterns to select for resistant strains under the pressure of drug treatment. The HIV protease is a crucial enzyme for viral maturation and has been a game changing drug target since the first application. Due to similarities in protease inhibitor designs, drug cross-resistance is not uncommon across ARVs of the same class. It is known that resistance against protease inhibitors is associated with a wider active site, but results from our large scale molecular dynamics simulations combined with statistical tests and network analysis further show, for the first time, that there are regions of local expansions and compactions associated with high levels of resistance conserved across eight different protease inhibitors visible in their complexed form within closed receptor conformations. The observed conserved expansion sites may provide an alternative drug-targeting site. Further, the method developed here is novel, supplementary to methods of variation analysis at sequence level, and should be applicable in analysing the structural consequences of mutations in other contexts using molecular ensembles.

We describe the design, synthesis, and biological evaluation of a series of novel HIV-1 protease inhibitors with carboxamide derivatives as the P2 ligands. We have specifically designed aminothiochromane and aminotetrahydronaphthalene-based carboxamide ligands to promote hydrogen bonding and van der Waals interactions in the active site of HIV-1 protease. Inhibitors 4e and 4j have shown potent enzyme inhibitory and antiviral activity. High resolution X-ray crystal structures of 4d- and 4k-bound HIV-1 protease revealed molecular insights into the ligand-binding site interactions., (Copyright © 2018 Elsevier Masson SAS. All rights reserved.)

Recently, allosteric regulations of HIV-1 protease (PR) are suggested as a promising approach to relieve drug resistance of mutations toward inhibitors targeting the active site of PR. Replica-exchange molecular dynamics (REMD) simulations and normal mode analysis (NMA) are integrated to enhance conformational sampling of PR. Molecular mechanics generalized Born surface area (MM-GBSA) method was applied to calculate binding free energies of three inhibitors APV, DRV, and NIT to the wild-type (WT) and multidrug resistance (MDR) PRs. The results suggest that binding free energies of APV and DRV are decreased in the MDR PR relative to the WT PR, suggesting drug resistance of mutations on these two inhibitors. However, the binding ability of the allosteric inhibitor NIT is not impaired in the MDR PR. In addition, internal dynamics analysis based on REMD simulations proves that mutations hardly produce obvious effect on the conformation of the MDR PR in comparison to the WT PR. Scanning of hydrophobic contacts and hydrogen bond contacts of inhibitors with residues of PRs on the concatenated trajectories of REMD demonstrates that mutations change the symmetric interaction networks of APV and DRV with PR, but do not generate obvious influence on the asymmetric interaction network of NIT with PR. In summary, allosteric inhibitor NIT can adapt the MDR PR better than those inhibitors toward the active site of PR, thus allosteric inhibitors of PR may be a possible channel to overcome drug resistance of PR., (© 2018 Wiley Periodicals, Inc.)