Search

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.

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.

Darunavir is a potent HIV protease inhibitor that has been established as an effective tool in the fight against the progression of HIV/AIDS in the global community. The successful application of this drug has spurred the development of derivatives wherein strategic regions (e.g., P1, P1', P2, and P2') of the darunavir framework have been structurally modified. An alternate route for the synthesis of darunavir and three related P1 and P1' derivatives has been developed. This synthetic pathway involves the use of a Crimmins titanium tetrachloride-mediated oxazolidine-2-thione-guided asymmetric glycolate aldol addition reaction. The resultant aldol adduct introduces the P1 fragment of darunavir via an aldehyde. Transamidation with a selected amine (isobutylamine or 2-ethyl-1-butylamine) to cleave the auxiliary yields an amide wherein the P1' component is introduced. From this stage, the amide is reduced to the corresponding β-amino alcohol and the substrate is then bis-nosylated to introduce the requisite p -nitrobenzenesulfonamide component and activate the secondary alcohol for nucleophilic substitution. Treatment with sodium azide yielded the desired azides, and the deprotection of the p -methoxyphenoxy group is achieved with the use of ceric ammonium nitrate. Finally, hydrogenation to reduce both the aniline and azide functionalities with concurrent acylation yields darunavir and its derivatives.

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

In the pursuit of novel antiretroviral therapies for human immunodeficiency virus type-1 (HIV-1) proteases (PRs), recent improvements in drug discovery have embraced machine learning (ML) techniques to guide the design process. This study employs ensemble learning models to identify crucial substructures as significant features for drug development. Using molecular docking techniques, a collection of 160 darunavir (DRV) analogs was designed based on these key substructures and subsequently screened using molecular docking techniques. Chemical structures with high fitness scores were selected, combined, and one-dimensional (1D) screening based on beyond Lipinski's rule of five (bRo5) and ADME (absorption, distribution, metabolism, and excretion) prediction implemented in the Combined Analog generator Tool (CAT) program. A total of 473 screened analogs were subjected to docking analysis through convolutional neural networks scoring function against both the wild-type (WT) and 12 major mutated PRs. DRV analogs with negative changes in binding free energy ( ΔΔ G bind ) compared to DRV could be categorized into four attractive groups based on their interactions with the majority of vital PRs. The analysis of interaction profiles revealed that potent designed analogs, targeting both WT and mutant PRs, exhibited interactions with common key amino acid residues. This observation further confirms that the ML model-guided approach effectively identified the substructures that play a crucial role in potent analogs. It is expected to function as a powerful computational tool, offering valuable guidance in the identification of chemical substructures for synthesis and subsequent experimental testing., (© 2024 Wiley Periodicals LLC.)

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

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-1 protease (PR) plays a crucial role in the treatment of HIV as a key target. The global issue of emerging drug resistance is escalating, and PR mutations pose a substantial challenge to the effectiveness of inhibitors. HIV-1 PR is an ideal model for studying drug resistance to inhibitors. The inhibitor, darunavir (DRV), exhibits a high genetic barrier to viral resistance, but with mutations of residues in the PR, there is also some resistance to DRV. Inhibitors can impede PR in two ways: one involves binding to the active site of the dimerization protease, and the other involves binding to the PR monomer, thereby preventing dimerization. In this study, we aimed to investigate the inhibitory effect of DRV with a modified inhibitor on PR, comparing the differences between wild-type and mutated PR, using molecular dynamics simulations. The inhibitory effect of the inhibitors on PR monomers was subsequently investigated. And molecular mechanics Poisson-Boltzmann surface area evaluated the binding free energy. The energy contribution of individual residues in the complex was accurately calculated by the alanine scanning binding interaction entropy method. The results showed that these inhibitors had strong inhibitory effects against PR mutations, with GRL-142 exhibiting potent inhibition of both the PR monomer and dimer. Improved inhibitors could strengthen hydrogen bonds and interactions with PR, thereby boosting inhibition efficacy. The binding of the inhibitor and mutation of the PR affected the distance between D25 and I50, preventing their dimerization and the development of drug resistance. This study could accelerate research targeting HIV-1 PR inhibitors and help to further facilitate drug design targeting both mechanisms.

Protease inhibitors are the most potent antivirals against HIV-1, but they still lose efficacy against resistant variants. Improving the resistance profile is key to developing more robust inhibitors, which may be promising candidates for simplified next-generation antiretroviral therapies. In this study, we explored analogs of darunavir with a P1 phosphonate modification in combination with increasing size of the P1' hydrophobic group and various P2' moieties to improve potency against resistant variants. The phosphonate moiety substantially improved potency against highly mutated and resistant HIV-1 protease variants, but only when combined with more hydrophobic moieties at the P1' and P2' positions. Phosphonate analogs with a larger hydrophobic P1' moiety maintained excellent antiviral potency against a panel of highly resistant HIV-1 variants, with significantly improved resistance profiles. The cocrystal structures indicate that the phosphonate moiety makes extensive hydrophobic interactions with the protease, especially with the flap residues. Many residues involved in these protease-inhibitor interactions are conserved, enabling the inhibitors to maintain potency against highly resistant variants. These results highlight the need to balance inhibitor physicochemical properties by simultaneous modification of chemical groups to further improve resistance profiles., 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.)

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

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.

The availability of drugs capable of blocking the replication of microorganisms has been one of the greatest triumphs in the history of medicine, but the emergence of an ever-increasing number of resistant strains poses a serious problem for the treatment of infectious diseases. The search for new potential ligands for proteins involved in the life cycle of pathogens is, therefore, an extremely important research field today. In this work, we have considered the HIV-1 protease, one of the main targets for AIDS therapy. Several drugs are used today in clinical practice whose mechanism of action is based on the inhibition of this enzyme, but after years of use, even these molecules are beginning to be interested by resistance phenomena. We used a simple artificial intelligence system for the initial screening of a data set of potential ligands. These results were validated by docking and molecular dynamics, leading to the identification of a potential new ligand of the enzyme which does not belong to any known class of HIV-1 protease inhibitors. The computational protocol used in this work is simple and does not require large computational power. Furthermore, the availability of a large number of structural information on viral proteins and the presence of numerous experimental data on their ligands, with which it is possible to compare the results obtained with computational methods, make this research field the ideal terrain for the application of these new computational techniques.

Human Immunodeficiency Virus type 1 protease (HIV-1 PR) is one of the most challenging targets of antiretroviral therapy used in the treatment of AIDS-infected people. The performance of protease inhibitors (PIs) is limited by the development of protease mutations that can promote resistance to the treatment. The current study was carried out using statistics and bioinformatics tools. A series of thirty-three compounds with known enzymatic inhibitory activities against HIV-1 protease was used in this paper to build a mathematical model relating the structure to the biological activity. These compounds were designed by software; their descriptors were computed using various tools, such as Gaussian, Chem3D, ChemSketch and MarvinSketch. Computational methods generated the best model based on its statistical parameters. The model's applicability domain (AD) was elaborated. Furthermore, one compound has been proposed as efficient against HIV-1 protease with comparable biological activity to the existing ones; this drug candidate was evaluated using ADMET properties and Lipinski's rule. Molecular Docking performed on Wild Type, and Mutant Type HIV-1 proteases allowed the investigation of the interaction types displayed between the proteases and the ligands, Darunavir (DRV) and the new drug (ND). Molecular dynamics simulation was also used in order to investigate the complexes' stability allowing a comparative study on the performance of both ligands (DRV & ND). Our study suggested that the new molecule showed comparable results to that of darunavir and maybe used for further experimental studies. Our study may also be used as pipeline to search and design new potential inhibitors of HIV-1 proteases., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Baassi et al. 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.)

Darunavir (DRV) is exceptional among potent HIV-1 protease inhibitors (PIs) in high drug concentrations that are achieved in vivo. Little is known about the de novo resistance pathway for DRV. We selected for resistance to high drug concentrations against 10 PIs and their structural precursor DRV. Mutations accumulated through two pathways (anchored by protease mutations I50V or I84V). Small changes in the inhibitor P1'-equivalent position led to preferential use of one pathway over the other. Changes in the inhibitor P2'-equivalent position determined differences in potency that were retained in the resistant viruses and that impacted the selected mutations. Viral variants from the two pathways showed differential selection of compensatory mutations in Gag cleavage sites. These results reveal the high level of selective pressure that is attainable with fifth-generation PIs and how features of the inhibitor affect both the resistance pathway and the residual potency in the face of resistance., Competing Interests: ES, SL, SZ, GL, MH, AS, KK, EN, AA, NY, CS, RS No competing interests declared, (© 2023, Spielvogel et al.)

Acquired immunodeficiency syndrome (AIDS) is one of the most challenging infectious diseases to treat on a global scale. Understanding the mechanisms underlying the development of drug resistance is necessary for novel therapeutics. HIV subtype C is known to harbor mutations at critical positions of HIV aspartic protease compared to HIV subtype B, which affects the binding affinity. Recently, a novel double-insertion mutation at codon 38 (L38HL) was characterized in HIV subtype C protease, whose effects on the interaction with protease inhibitors are hitherto unknown. In this study, the potential of L38HL double-insertion in HIV subtype C protease to induce a drug resistance phenotype towards the protease inhibitor, Saquinavir (SQV), was probed using various computational techniques, such as molecular dynamics simulations, binding free energy calculations, local conformational changes and principal component analysis. The results indicate that the L38HL mutation exhibits an increase in flexibility at the hinge and flap regions with a decrease in the binding affinity of SQV in comparison with wild-type HIV protease C. Further, we observed a wide opening at the binding site in the L38HL variant due to an alteration in flap dynamics, leading to a decrease in interactions with the binding site of the mutant protease. It is supported by an altered direction of motion of flap residues in the L38HL variant compared with the wild-type. These results provide deep insights into understanding the potential drug resistance phenotype in infected individuals., Competing Interests: The authors declare no conflicts of interest.

Human immunodeficiency virus-1 (HIV-1) protease is one of the important targets in AIDS therapy. The majority of HIV infections are caused due to non-B subtypes in developing countries. The co-occurrence of mutations along with naturally occurring polymorphisms in HIV-1 protease cause resistance to the FDA approved drugs, thereby posing a major challenge in the treatment of antiretroviral therapy. In this work, the resistance mechanism against SQV due to active site mutations G48V and V82F in CRF01_AE (AE) protease was explored. The binding free energy calculations showed that the direct substitution of valine at position 48 introduces a bulkier side chain, directly impairing the interaction with SQV in the binding pocket. Also, the intramolecular hydrogen bonding network of the neighboring residues is altered, indirectly affecting the binding of SQV. Interestingly, the substitution of phenylalanine at position 82 induces conformational changes in the 80's loop and the flap region, thereby favoring the binding of SQV. The V82F mutant structure also maintains similar intramolecular hydrogen bond interactions as observed in AE-WT.Communicated by Ramaswamy H. Sarma.

The development of dual inhibitors of HIV-1 protease and reverse transcriptase is an attractive strategy for multi-target therapeutic of AIDS, which may be privileged in delaying the occurrence of drug resistance. We herein designed a novel kind of dual inhibitors with benzofuran or indole cores. Biological results showed that a number of inhibitors displayed significant activity against both HIV-1 protease and reverse transcriptase. Among which, inhibitor 10f exhibited a good correlation with an approximate ratio of 1: 2 between the two enzymes. Furthermore, the dual inhibitors illustrated similar potency against both the wild-type virus and drug-resistant mutant. In addition, the molecular dynamic simulation studies verified the dual actions of such 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 © 2022 Elsevier Masson SAS. All rights reserved.)

Antiretroviral drug resistance is a therapeutic obstacle for people with HIV. HIV protease inhibitors darunavir and lopinavir are recommended for resistant infections. We characterized a protease mutant (PR10x) derived from a highly resistant clinical isolate including 10 mutations associated with resistance to lopinavir and darunavir. Compared to the wild-type protease, PR10x exhibits ∼3-fold decrease in catalytic efficiency and K i values of 2-3 orders of magnitude worse for darunavir, lopinavir, and potent investigational inhibitor GRL-519. Crystal structures of the mutant were solved in a ligand-free form and in complex with GRL-519. The structures show altered interactions in the active site, flap-core interface, hydrophobic core, hinge region, and 80s loop compared to the corresponding wild-type protease structures. The ligand-free crystal structure exhibits a highly curled flap conformation which may amplify drug resistance. Molecular dynamics simulations performed for 1 μs on ligand-free dimers showed extremely large fluctuations in the flaps for PR10x compared to equivalent simulations on PR with a single L76V mutation or wild-type protease. This analysis offers insight about the synergistic effects of mutations in highly resistant variants., 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 Inc. All rights reserved.)

The Human Immunodeficiency Virus (HIV) has been the source of epidemic infection of AIDS for a longer period. One of the most difficult tasks is identifying novel medications that can help to decrease or control this global health hazard by overcoming drug resistance. In recent decades' nanoparticles are emerging as extremely relevant in drug delivery platforms. In the current study, the pristine (SWCNT) and hydroxyl functionalized (SWCNT-OH) versions of the SWCNT were investigated as inhibitors against the wild-type (WT) and three key mutants of HIV-1 protease (HIV-pr) (I50V, V82A, and I84V). Molecular docking of SWCNT in the catalytic domain and running all-atom MD simulations of all complexes are also part of this project. A thorough inspection of conformational dynamics from 50 ns trajectories reveals that both the pristine and SWCNT-OH can fit right to the pocket region of HIV-pr and govern flap dynamics. The binding affinity of the four HIV-pr-SWCNT/SWCNT-OH complexes was further investigated using MM-PBSA-dependent binding free energy studies. In most mutants and WT systems, SWCNT-OH was reported to bind proportionately many folds (kcal/mol) more than pristine SWCNTs. Hence, SWCNTs are possible HIV-pr inhibitors in terms of their stable existence in the pocket area, stronger binding to the protease, and regulation of flap dynamics in controlling the active site volume, which have vast potential for applications against drug resistance., 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 Inc. All rights reserved.)

With the increasing prevalence of drug-resistant variants, novel potent HIV-1 protease inhibitors with broad-spectrum antiviral activity against multidrug-resistant causative viruses are urgently needed. Herein, we designed and synthesized a new series of HIV-1 protease inhibitors with phenols or polyphenols as the P2 ligands and a variety of sulfonamide analogs as the P2' ligands. A number of these new inhibitors showed superb enzymatic inhibitory activity and antiviral activity. In particular, inhibitors 15d and 15f exhibited potent enzymatic inhibitory activity in the low picomolar range, and the latter showed excellent activity against the Darunavir-resistant HIV-1 variant. Furthermore, the molecular modeling studies provided insight into the ligand-binding site interactions between inhibitors and the enzyme cavity, and they sparked inspiration for the further optimization of potent inhibitors.

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 human immunodeficiency virus type 1 (HIV-1) has continued to be a global concern. With the new HIV incidence, the emergence of multi-drug resistance and the untoward side effects of currently used anti-HIV drugs, there is an urgent need to discover more efficient anti-HIV drugs. Modern computational tools have played vital roles in facilitating the drug discovery process. This research focuses on a pharmacophore-based similarity search to screen 111,566,735 unique compounds in the PubChem database to discover novel HIV-1 protease inhibitors (PIs). We used an in silico approach involving a 3D-similarity search, physicochemical and ADMET evaluations, HIV protease-inhibitor prediction (IC 50 /percent inhibition), rigid receptor-molecular docking studies, binding free energy calculations and molecular dynamics (MD) simulations. The 10 FDA-approved HIV PIs (saquinavir, lopinavir, ritonavir, amprenavir, fosamprenavir, atazanavir, nelfinavir, darunavir, tipranavir and indinavir) were used as reference. The in silico analysis revealed that fourteen out of the twenty-eight selected optimized hit molecules were within the acceptable range of all the parameters investigated. The hit molecules demonstrated significant binding affinity to the HIV protease (PR) when compared to the reference drugs. The important amino acid residues involved in hydrogen bonding and п-п stacked interactions include ASP25, GLY27, ASP29, ASP30 and ILE50. These interactions help to stabilize the optimized hit molecules in the active binding site of the HIV-1 PR (PDB ID: 2Q5K). HPS/002 and HPS/004 have been found to be most promising in terms of IC 50 /percent inhibition (90.15%) of HIV-1 PR, in addition to their drug metabolism and safety profile. These hit candidates should be investigated further as possible HIV-1 PIs with improved efficacy and low toxicity through in vitro experiments and clinical trial investigations.

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

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

Developing novel drug molecules against HIV is a scientific quest necessitated by development of drug resistance against used drugs. We report comparative results of molecular dynamics simulation studies on 11 structural analogues of Saquinavir (SQV) - against HIV-protease that were earlier examined for pharmacodynamic and pharmacokinetic properties. We reported analogues S1, S5 and S8 to qualify the ADMET criterion and may be considered as potential lead molecules. In this study the designed molecules were successively docked with native HIV-protease at AutoDock. Docking scores established relative goodness of the 11 analogues against the benchmark for Saquinavir. The docked complexes were subjected to molecular dynamics simulation studies using GROMACS 4.6.2. Four parameters viz. H-bonding, RMSD, Binding energy and Protein-Ligand Distance were used for comparative analyses of the analogues relative to Saquinavir. The comparison and analysis of the results are indicative that analogues S8, S9 and S1 are promising candidates among all the analogues studied. From our earlier work and present study it is evident that among the three S8 and S1 qualify the ADMET criterion and between S1 and S8, the analogue S8 shows more target efficacy and specificity over S1 and have better molecular dynamics simulation results. Thus, of the 11 de novo Saquinavir analogues, the S8 appears to be the most promising candidate as lead molecule for HIV-protease inhibitor and is best suited for testing under biological system. Further validation of the proposed lead molecules through wet lab studies involving antiviral assays however is required.Communicated by Ramaswamy H. Sarma.

Based upon the preliminary design of enhancing genetic barrier to drug-resistant viral mutants by maximizing hydrogen-bonding or other van der Waals contacts, we have designed, synthesized and biologically evaluated a new class of HIV-1 protease inhibitors with phenol derived P2 ligands and nitro or halogens in P2' ligands. Results indicate that a majority of inhibitors exhibit robust enzyme inhibitory with IC 50 values in picomolar or single digit nanomolar ranges. Among which, compound 17d displays potency with IC 50 value of 21 pM and high protease selectivity. Of note, 17d exhibits greater antiviral activity against the DRV-resistant variant than the efficacy against the wild type virus. Furthermore, the molecular modeling studies demonstrate important interactions between 17d and the active sites of both the wild-type and DRV-resistant HIV-1 protease, as well as furnish insights for further optimization of new inhibitors., (Copyright © 2022 Elsevier Ltd. All rights reserved.)

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

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.

PR S17 , a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PR S17 with inhibitor Darunavir (DRV), two complexes of WT PR and PR S17 with substrate analogue CA-p2, two unligand WT PR and PR S17 ), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PR S17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.

The binding behavior of ritonavir (RTV), a HIV/AIDS protease inhibitor, with ct-DNA was characterized through multiple testing technologies and theoretical calculation. The findings revealed that the RTV-DNA complex was formed through the noncovalent interaction mainly including conventional hydrogen bonds and carbon hydrogen bonds as well as hydrophobic interactions (pi-alkyl interactions). The stoichiometry and binding constant of the RTV-DNA complex were 1:1 and 1.87 × 10 3  M -1 at 298 K, respectively, indicating that RTV has moderate affinity with ct-DNA. The findings confirmed that RTV binds to the minor groove of DNA. The outcomes of CD experiments showed that the binding with RTV changed the conformation of DNA slightly. However, the conformation of RTV had obvious changes after binding to DNA, meaning that the flexibility of RTV molecule played an important role in stabilizing the RTV-DNA complex. Meanwhile, the results of DFT calculation revealed that the RTV and DNA interaction caused the changes in the frontier molecular orbitals, dipole moment and atomic charge distribution of RTV, altering the chemical properties of RTV when it bound to DNA. Communicated by Ramaswamy H. Sarma.

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.

HIV protease (HIVPR) is a key target in AIDS therapeutics. All ten FDA-approved drugs that compete with substrates in binding to this dimeric enzyme's active site have become ineffective due to the emergence of drug resistant mutants. Blocking the dimerization interface of HIVPR is thus being explored as an alternate strategy. The latest drug, darunavir (DRV), which exhibited a high genetic barrier to viral resistance, is said to have a dual mode of action - (i) binding to the dimeric active site, and (ii) preventing the dimerization by binding to the HIVPR monomer. Despite several reports on DRV complexation with dimeric HIVPR, the mode and mechanism of the binding of DRV to the HIVPR monomer are poorly understood. In this study, we utilized all-atomic MD simulations and umbrella sampling techniques to identify the best possible binding mode of DRV to the monomeric HIVPR and its mechanism of association. The results suggest that DRV binds between the active site and the flap of the monomer, and the flap plays a crucial role in directing the drug to bind and driving the other protein domains to undergo induced fit changes for stronger complexation. The obtained binding mode of DRV was validated by comparing with various mutational data from clinical isolates to reported in vitro mutations. The identified binding pose was also able to successfully reproduce the experimental K i value in the picomolar range. The residue-level information extracted from this study could accelerate the structure-based drug designing approaches targeting HIVPR dimerization.

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.

HIV-1 protease (HIV-1 PR) is an essential enzyme for the replication process of its virus, and therefore considered an important target for the development of drugs against the acquired immunodeficiency syndrome (AIDS). Our previous study shows that the catalytic mechanism of subtype B/C-SA HIV-1 PR follows a one-step concerted acyclic hydrolysis reaction process using a two-layered ONIOM B3LYP/6-31++G(d,p) method. This present work is aimed at exploring the proposed mechanism of the proteolysis catalyzed by HIV-1 PR and to ensure our proposed mechanism is not an artefact of a single theoretical technique. Hence, we present umbrella sampling method that is suitable for calculating potential mean force (PMF) for non-covalent ligand/substrate-enzyme association/dissociation interactions which provide thermodynamic details for molecular recognition. The free activation energy results were computed in terms of PMF analysis within the hybrid QM(DFTB)/MM approach. The theoretical findings suggest that the proposed mechanism corresponds in principle with experimental data. Given our observations, we suggest that the QM/MM MD method can be used as a reliable computational technique to rationalize lead compounds against specific targets such as the HIV-1 protease.

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.

Purpose: To understand the role of different surfactants, incorporated into amorphous solid dispersions (ASDs) of ritonavir and copovidone, in terms of their impact on release, phase behavior and stabilization of amorphous precipitates formed following drug release., Methods: Ternary ASDs with ritonavir, copovidone and surfactants (30:70:5 w/w/w) were prepared by rotary evaporation. ASD release performance was tested using Wood's intrinsic dissolution rate apparatus and compared to the binary drug-polymer ASD with 30% drug loading. Size measurement of amorphous droplets was performed using dynamic light scattering. Solid state characterization was performed using attenuated total reflectance-infrared spectroscopy, differential scanning calorimetry and scanning electron microscopy., Results: All surfactant-containing ASDs showed improvement over the binary ASD. Span 85 and D-α-tocopheryl polyethylene glycol succinate (TPGS) showed complete release with no evidence of AAPS or crystallization whereas Span 20 and Tween 80 showed < 50% release with amorphous amorphous phase separation (AAPS). Span 20 also induced solution crystallization. Sodium dodecyl sulfate (SDS) showed very rapid, albeit incomplete (~ 80%) release. AAPS was not observed with SDS. However, crystallization on the dissolving solid surface was noted. Span 20 and TPGS formed the smallest and most size-stable droplets with ~ 1 µm size whereas coalescence was noted with other surfactants., Conclusions: Surfactants improved the release performance relative to the binary ASD. Different surfactant types impacted overall performance to varying extents and affected different attributes. Overall, Span 85 showed best performance (complete release, no crystallization/AAPS and small droplet size). Correlation between physicochemical properties and surfactant performance was not observed., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

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

Nelfinavir is one of the FDA-approved HIV-1 protease inhibitors and a part of highly active anti-retroviral therapy (HAART) for the treatment of HIV-AIDS. Nelfinavir was the first HIV-1 protease inhibitor to be approved as a paediatric formulation. The application of HAART had resulted in significant improvement in the lives of AIDS patients. However, the emergence of drug resistance in HIV-1 protease has limited the use of many of these drugs including nelfinavir. A unique mutation observed frequently in patients treated with nelfinavir is D30N as it is selected exclusively by nelfinavir. The D30N mutation imparts very high resistance to nelfinavir but unlike other primary mutations does not give cross-resistance to the majority of other drugs. D30N mutation also significantly reduces cleavage activity of HIV-1 protease and affects viral fitness. Here, we have determined crystal structures of D30N HIV-1 protease in unliganded form and in complex with nelfinavir. These structures provide the rationale for reduced cleavage activity and the molecular basis of drug resistance induced by D30N mutation. The loss of coulombic interaction part of a crucial hydrogen bond between the drug and the protease is likely to play a major role in reduced affinity and resistance towards nelfinavir. The decreased catalytic activity of D30N HIV-1 protease due to altered interaction with the substrates and reduced stability of folding core may be the reason for the reduced replicative capacity of the virus harboring mutant HIV-1 protease.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 is highly homologous to SARS-CoV. To date, the main protease (Mpro) of SARS-CoV-2 is regarded as an important drug target for the treatment of Coronavirus Disease 2019 (COVID-19). Some experiments confirmed that several HIV protease inhibitors present the inhibitory effects on the replication of SARS-CoV-2 by inhibiting Mpro. However, the mechanism of action has still not been studied very clearly. In this work, the interaction mechanism of four HIV protease inhibitors Darunavir (DRV), Lopinavir (LPV), Nelfinavir (NFV), and Ritonavire (RTV) targeting SARS-CoV-2 Mpro was explored by applying docking, molecular dynamics (MD) simulations, and MM-GBSA methods using the broad-spectrum antiviral drug Ribavirin (RBV) as the negative and nonspecific control. Our results revealed that LPV, RTV, and NFV have higher binding affinities with Mpro, and they all interact with catalytic residues His41 and the other two key amino acids Met49 and Met165. Pharmacophore model analysis further revealed that the aromatic ring, hydrogen bond donor, and hydrophobic group are the essential infrastructure of Mpro inhibitors. Overall, this study applied computational simulation methods to study the interaction mechanism of HIV-1 protease inhibitors with SARS-CoV-2 Mpro, and the findings provide useful insights for the development of novel anti-SARS-CoV-2 agents for the treatment of COVID-19.

Acquired immunodeficiency syndrome (AIDS) caused by the human immunodeficiency virus (HIV) continues to be a public health problem. In 2020, 680,000 people died from HIV-related causes, and 1.5 million people were infected. Antiretrovirals are a way to control HIV infection but not to cure AIDS. As such, effective treatment must be developed to control AIDS. Developing a drug is not an easy task, and there is an enormous amount of work and economic resources invested. For this reason, it is highly convenient to employ computer-aided drug design methods, which can help generate and identify novel molecules. Using the de novo design, novel molecules can be developed using fragments as building blocks. In this work, we develop a virtual focused compound library of HIV-1 viral protease inhibitors from natural product fragments. Natural products are characterized by a large diversity of functional groups, many sp 3 atoms, and chiral centers. Pseudo-natural products are a combination of natural products fragments that keep the desired structural characteristics from different natural products. An interactive version of chemical space visualization of virtual compounds focused on HIV-1 viral protease inhibitors from natural product fragments is freely available in the supplementary material.

Amorphous solid dispersions (ASDs) are regarded as one of the most promising techniques for poorly-soluble active pharmaceutical ingredients (API). However, the thermodynamic instability of ASDs at supersaturated state makes them easy to recrystallize in aqueous media. In this study, ritonavir (RTV) was selected as a model drug for evaluating the solubility enhancement and recrystallization inhibition effect of various cellulose derivatives and the combinations of them with typical surfactants. Combination of HPMCAS-HF/SLS was filtrated for preparing ternary RTV solid dispersions (RTV SD) via solvent evaporation method. RTV SD exhibited enhanced dissolution manner, while the oral bioavailability of RTV SD was equivalent with the Reference Standard Norvir® but increased significantly compared to the ternary physical mixture. Thus, the ternary SD system might be promisingly employed as efficient drug delivery system for RTV, while the HPMCAS-HF/SLS combination could be recommended as effective excipient for fabricating steady solid dispersions loading poorly soluble API., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

A photocatalytic decarboxylative radical conjugate addition-elimination-oxa-Michael reaction of hydroxyalkylated carboxylic acids with cyclopentenones is developed to construct diverse cyclopentanonyl-fused functionalized 5-7 membered cyclic ethers. The stereoselective synthetic strategy is amenable to substructural variation, establishing a direct total synthetic route to two diastereomers of C3-amino cyclopentyltetrahydrofuranyl-derived potent HIV-1 protease inhibitors with low nanomolar IC 50 values.

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

Upon the basis of both possible ligand-binding site interactions and the uniformity of key residues in active sites, a novel class of HIV-1 PR/RT dual inhibitors was designed and evaluated. Cinnamic acids or phenylpropionic acids with more flexible chain and smaller steric hindrance were introduced into the inhibitors, giving rise to significant improvement in HIV-1 RT inhibitory activity by one or two orders of magnitude, with comparable or even improved potency against PR at the same time, compared with coumarin anologues in our previous studies. Among these inhibitors, 38d displayed a 19-fold improvement in anti-PR activity with IC 50 value of 0.081 nM compared to the control DRV. In addition, inhibitor 38c exhibited an excellent anti-RT IC 50 value of 0.43 μM, only a 4.7-fold less potent activity than the control EFV. More significantly, the disparate ratio between HIV-1 PR and RT inhibition became more reasonable with ratio of 1: 10.4, just as 37b. Furthermore, the assays on HIV-1 late stage and early stage supported the rationality of designing dual inhibitors. The SAR data as well as molecular modeling studies provided new insight for further optimization of more potent HIV-1 PR/RT dual 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 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.

Purpose: Application of multi-scale modelling workflows to characterise polymorphism in ritonavir with regard to its stability, bioavailability and processing., Methods: Molecular conformation, polarizability and stability are examined using quantum mechanics (QM). Intermolecular synthons, hydrogen bonding, crystal morphology and surface chemistry are modelled using empirical force fields., Results: The form I conformation is more stable and polarized with more efficient intermolecular packing, lower void space and higher density, however its shielded hydroxyl is only a hydrogen bond donor. In contrast, the hydroxyl in the more open but less stable and polarized form II conformation is both a donor and acceptor resulting in stronger hydrogen bonding and a more stable crystal structure but one that is less dense. Both forms have strong 1D networks of hydrogen bonds and the differences in packing energies are partially offset in form II by its conformational deformation energy difference with respect to form I. The lattice energies converge at shorter distances for form I, consistent with its preferential crystallization at high supersaturation. Both forms exhibit a needle/lath-like crystal habit with slower growing hydrophobic side and faster growing hydrophilic capping habit faces with aspect ratios increasing from polar-protic, polar-aprotic and non-polar solvents, respectively. Surface energies are higher for form II than form I and increase with solvent polarity. The higher deformation, lattice and surface energies of form II are consistent with its lower solubility and hence bioavailability., Conclusion: Inter-relationship between molecular, solid-state and surface structures of the polymorphic forms of ritonavir are quantified in relation to their physical-chemical properties.