Enthalpic Classification of Water Molecules in Target-Ligand Binding.

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.

Multiple Parameter Replica Exchange Gaussian Accelerated Molecular Dynamics for Enhanced Sampling and Free Energy Calculation of Biomolecular Systems.

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

HIV-1 transmitted drug resistance in newly diagnosed individuals in Italy over the period 2015-21.

BACKGROUND: Transmitted drug resistance (TDR) is still a critical aspect for the management of individuals living with HIV-1. Thus, its evaluation is crucial to optimize HIV care. METHODS: Overall, 2386 HIV-1 protease/reverse transcriptase and 1831 integrase sequences from drug-naïve individuals diagnosed in north and central Italy between 2015 and 2021 were analysed. TDR was evaluated over time. Phylogeny was generated by maximum likelihood. Factors associated with TDR were evaluated by logistic regression. RESULTS: Individuals were mainly male (79.1%) and Italian (56.2%), with a median (IQR) age of 38 (30-48). Non-B infected individuals accounted for 44.6% (N = 1065) of the overall population and increased over time (2015-2021, from 42.1% to 51.0%, P = 0.002). TDR prevalence to any class was 8.0% (B subtype 9.5% versus non-B subtypes 6.1%, P = 0.002) and remained almost constant over time. Overall, 300 transmission clusters (TCs) involving 1155 (48.4%) individuals were identified, with a similar proportion in B and non-infected individuals (49.7% versus 46.8%, P = 0.148). A similar prevalence of TDR among individuals in TCs and those out of TCs was found (8.2% versus 7.8%, P = 0.707).By multivariable analysis, subtypes A, F, and CFR02_AG were negatively associated with TDR. No other factors, including being part of TCs, were significantly associated with TDR. CONCLUSIONS: Between 2015 and 2021, TDR prevalence in Italy was 8% and remained almost stable over time. Resistant strains were found circulating regardless of being in TCs, but less likely in non-B subtypes. These results highlight the importance of a continuous surveillance of newly diagnosed individuals for evidence of TDR to inform clinical practice.

Design of substituted tetrahydrofuran derivatives for HIV-1 protease inhibitors: synthesis, biological evaluation, and X-ray structural studies.

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.

Resistance to protease inhibitors among persons living with HIV in Ghana: a case for viral load and drug resistance monitoring.

OBJECTIVE: Routine viral load and drug resistance testing are well supported in most resource-rich settings and provide valuable benefits in the clinical care of PLWH in these communities. Undoubtedly, there exist financial and political constraints for the scale-up of viral load and drug resistance testing in Sub-Saharan Africa. To achieve the global UNAIDS 95/95/95 targets, there is the need to bridge this inequity in patient care and allow for a universal approach that leaves no community behind. METHODS: Venous blood from 96 PLWH on second-line ART from Korle-Bu Teaching Hospital were collected and processed into plasma for CD4+ T- cell and viral load assessments. Ribonucleic acid (RNA) was extracted from stored plasma and the protease gene amplified, sequenced and analyzed for subtype and drug resistance mutations using the Stanford HIV drug resistance database. RESULTS: Out of the 96 PLWH, 37 experienced virological failure with 8 patients' samples successfully sequenced. The predominant HIV-1 subtype identified was CRF02_AG (6/8, 75.0%) with 12.5% (1/8) each of CFR06_cpx infection and one case unable to subtype. The major PI resistance mutations identified were; M46I, I54V, V82A, I47V, I84V and L90M. CONCLUSIONS: Persons living with HIV who had experienced virologic failure in this study harboured drug resistance mutations to PI, thus compromise the effectiveness of the drugs in the second line. Resistance testing is strongly recommended prior to switching to a new regimen. This will help to inform the choice of drug and to achieve optimum therapeutic outcome among PLWH in Ghana.

Quantitative Prediction of Human Immunodeficiency Virus Drug Resistance.

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

A long-term stable cold-chain-friendly HIV mRNA vaccine encoding multi-epitope viral protease cleavage site immunogens inducing immunogen-specific protective T cell immunity.

The lack of success in clinical trials for HIV vaccines highlights the need to explore novel strategies for vaccine development. Research on highly exposed seronegative (HESN) HIV-resistant Kenyan female sex workers revealed naturally protective immunity is correlated with a focused immune response mediated by virus-specific CD8 T cells. Further studies indicated that the immune response is unconventionally focused on highly conserved sequences around HIV viral protease cleavage sites (VPCS). Thus, taking an unconventional approach to HIV vaccine development, we designed lipid nanoparticles loaded with mRNA that encodes multi-epitopes of VPCS (MEVPCS-mRNA LNP), a strategic design to boost antigen presentation by dendritic cells, promoting effective cellular immunity. Furthermore, we developed a novel cold-chain compatible mRNA LNP formulation, ensuring long-term stability and compatibility with cold-chain storage/transport, widening accessibility of mRNA LNP vaccine in low-income countries. The in-vivo mouse study demonstrated that the vaccinated group generated VPCS-specific CD8 memory T cells, both systemically and at mucosal sites of viral entry. The MEVPCS-mRNA LNP vaccine-induced CD8 T cell immunity closely resembled that of the HESN group and displayed a polyfunctional profile. Notably, it induced minimal to no activation of CD4 T cells. This proof-of-concept study underscores the potential of the MEVPCS-mRNA LNP vaccine in eliciting CD8 T cell memory specific to the highly conserved multiple VPCS, consequently having a broad coverage in human populations and limiting viral escape mutation. The MEVPCS-mRNA LNP vaccine holds promise as a candidate for an effective prophylactic HIV vaccine.

Premature Activation of the HIV-1 Protease Is Influenced by Polymorphisms in the Hinge Region.

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

LSTM-driven drug design using SELFIES for target-focused de novo generation of HIV-1 protease inhibitor candidates for AIDS treatment.

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.

Interplay between protease and reverse transcriptase dimerization in a model HIV-1 polyprotein.

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

Performance of the Applied Biosystems HIV-1 Genotyping Kit with Integrase.

HIV genotyping is used to assess HIV susceptibility to antiretroviral drugs. The Applied Biosystems HIV-1 Genotyping Kit with Integrase (AB kit, Thermo Fisher Scientific) detects resistance-associated mutations (RAMs) in HIV protease (PR), reverse transcriptase (RT), and integrase (IN). We compared results from the AB kit with results obtained previously with the ViroSeq HIV-1 Genotyping System. DNA amplicons from the AB kit were also analyzed using next-generation sequencing (NGS). HIV RNA was extracted using the MagNA Pure 24 instrument (Roche Diagnostics; 96 plasma samples, HIV subtype B, viral load range: 530-737,741 copies/mL). FASTA files were generated from AB kit data using Exatype (Hyrax Biosciences). DNA amplicons from the AB kit were also analyzed by NGS using the Nextera XT kit (Illumina). Drug resistance was predicted using the Stanford HIV Drug Resistance Database. The mean genetic distance for sequences from ViroSeq and the AB kit was 0.02% for PR/RT and 0.04% for IN; 103 major RAMs were detected by both methods. Four additional major RAMs were detected by the AB kit only. These four major RAMs were also detected by NGS (detected in 18.1%-38.2% of NGS reads). NGS detected 27 major RAMs that were not detected with either of the Sanger sequencing-based kits. All major RAMs detected with ViroSeq were detected with the AB kit; additional RAMs were detected with the AB kit only. DNA amplicons from the AB kit can be used for NGS for more sensitive detection of RAMs.

Integrating Explicit and Implicit Fullerene Models into UNRES Force Field for Protein Interaction Studies.

Fullerenes, particularly C60, exhibit unique properties that make them promising candidates for various applications, including drug delivery and nanomedicine. However, their interactions with biomolecules, especially proteins, remain not fully understood. This study implements both explicit and implicit C60 models into the UNRES coarse-grained force field, enabling the investigation of fullerene-protein interactions without the need for restraints to stabilize protein structures. The UNRES force field offers computational efficiency, allowing for longer timescale simulations while maintaining accuracy. Five model proteins were studied: FK506 binding protein, HIV-1 protease, intestinal fatty acid binding protein, PCB-binding protein, and hen egg-white lysozyme. Molecular dynamics simulations were performed with and without C60 to assess protein stability and investigate the impact of fullerene interactions. Analysis of contact probabilities reveals distinct interaction patterns for each protein. FK506 binding protein (1FKF) shows specific binding sites, while intestinal fatty acid binding protein (1ICN) and uteroglobin (1UTR) exhibit more generalized interactions. The explicit C60 model shows good agreement with all-atom simulations in predicting protein flexibility, the position of C60 in the binding pocket, and the estimation of effective binding energies. The integration of explicit and implicit C60 models into the UNRES force field, coupled with recent advances in coarse-grained modeling and multiscale approaches, provides a powerful framework for investigating protein-nanoparticle interactions at biologically relevant scales without the need to use restraints stabilizing the protein, thus allowing for large conformational changes to occur. These computational tools, in synergy with experimental techniques, can aid in understanding the mechanisms and consequences of nanoparticle-biomolecule interactions, guiding the design of nanomaterials for biomedical applications.

A validated in-house assay for HIV drug resistance mutation surveillance from dried blood spot specimens.

Despite increasing scale-up of antiretroviral therapy (ART) coverage, challenges related to adherence and HIV drug resistance (HIVDR) remain. The high cost of HIVDR surveillance is a persistent challenge with implementation in resource-constrained settings. Dried blood spot (DBS) specimens have been demonstrated to be a feasible alternative to plasma or serum for HIVDR genotyping and are more suitable for lower resource settings. There is a need for affordable HIVDR genotyping assays which can amplify HIV-1 sequences from DBS specimens, particularly those with low viral loads, at a low cost. Here, we present an in-house assay capable of reliably amplifying HIV-1 protease and partial reverse transcriptase genes from DBS specimens, which covers the complete World Health Organization 2009 list of drug resistance mutations under surveillance. DBS specimens were prepared using whole blood spiked with HIV-1 at concentrations of 10,000, 5000, 1000, and 500 copies/mL (n=30 for each concentration). Specimens were tested in triplicate. A two-step approach was used consisting of cDNA synthesis followed by nested PCR. The limit of detection of the assay was calculated to be approximately 5000 (95% CI: 3200-10,700) copies/mL for the protease gene and 3600 (95% CI: 2200-10,000) copies/mL for reverse transcriptase. The assay was observed to be most sensitive with higher viral load specimens (97.8% [95% CI: 92.2-99.7]) for both protease and reverse transcriptase at 10,000 copies/mL with performance decreasing with the use of specimens with lower viral loads (46.7% [36.1-57.5] and 60.0% [49.1-70.2] at 500 copies/mL for protease and reverse transcriptase, respectively). Ultimately, this assay presents a promising opportunity for use in resource-constrained settings. Future work should involve validation under field conditions including sub-optimal storage conditions and preparation of DBS with fingerprick blood in order to accurately reflect real-world collection scenarios.

HIV protease resistance mutations in patients receiving second-line antiretroviral therapy in Libreville, Gabon.

INTRODUCTION: In 2022, the WHO reported that 29.8 million people around the world were living with HIV (PLHIV) and receiving antiretroviral treatment (ART), including 25| 375 people in Gabon (54% of all those living with HIV in the country). The literature reports a frequency of therapeutic failure with first-line antiretrovirals (ARVs) of between 20% and 82%. Unfortunately, data relating to the failure of second-line ARVs are scarce in Gabon. This study aims to determine the profiles of HIV drug resistance mutations related to protease inhibitors in Gabon. METHODOLOGY: Plasma from 84 PLHIV receiving ARVs was collected from 2019 to 2021, followed by RNA extraction, amplification, and sequencing of the protease gene. ARV resistance profiles were generated using the Stanford interpretation algorithm version 8.9-1 ( https://hivdb.stanford.edu ) and statistical analyses were performed using EpiInfo software version 7.2.1.0 (CDC, USA). RESULTS: Of 84 HIV plasma samples collected from 45 men and 39 women, 342 mutations were detected. Of these, 43.3% (148/342) were associated with nucleoside reverse transcriptase inhibitors (NRTIs), 30.4% (104/342) with non-nucleoside reverse transcriptase inhibitors (NNRTIs), and 26.3% (90/342) with protease inhibitors (PIs). Most NRTI mutations were associated with thymidine analogues (TAMs) (50.7%; 75/148), including T215F/V (14.9%; 22/148), D67DN/E/G/N/T (10.1%; 15/148), M41L (9.5%; 14/148), and K70E/KN/S/R (9.5%; 14/148). Resistance mutations related to non-TAM NRTIs (33.1%; 49/148) were M184V (29.1%; 43/148), and L74I/V (8.1%; 12/148). NNRTI mutations were predominantly K103N/S (32.7%; 34/104), V108I (10.6%; 11/104), A98G (10.6%; 11/104), and P225H (9.6%; 10/104). Minor mutations associated with PIs (60.0%; 54/90) were predominantly K20I (15.6%; 14/90) and L10F/I/V (14.5%; 13/90). The major mutations associated with PIs (40.0%; 36/90) were M41L (12.2%; 11/90), I84V (6.7%; 06/90), and V82A (6.7%; 06/90). The four most prescribed therapeutic regimens were TDF + 3TC + LPV/r (20.3%; 17/84), ABC + DDI + LPV/r (17.9%; 15/84), TDF + FTC + LPV/r (11.9%; 10/84), and ABC + 3TC + LPV/r (11.9%; 10/84). CONCLUSION: This study revealed that HIV drug resistance mutations are common in Gabon. The major mutations associated with PIs were M41L, I84V, and V82A. There is a need for access to new NRTIs, NNRTIs, and PIs for a better therapeutic management of PLHIV in Gabon.

Natural Polymorphisms D60E and I62V Stabilize a Closed Conformation in HIV-1 Protease in the Absence of an Inhibitor or Substrate.

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

Rational design, synthesis and biological evaluation of novel HIV-1 protease inhibitors containing 2-phenylacetamide derivatives as P2 ligands with potent activity against DRV-Resistant HIV-1 variants.

A novel kind of potent HIV-1 protease inhibitors, containing diverse hydroxyphenylacetic acids as the P2-ligands and 4-substituted phenyl sulfonamides as the P2' ligands, were designed, synthesized and evaluated in this work. Majority of the target compounds exhibited good to excellent activity against HIV-1 protease with IC50 values below 200 nM. In particular, compound 18d with a 2-(3,4-dihydroxyphenyl) acetamide as the P2 ligand and a 4- methoxybenzene sulfonamide P2' ligand exhibited inhibitory activity IC50 value of 0.54 nM, which was better than that of the positive control darunavir (DRV). More importantly, no significant decline of the potency against HIV-1DRVRS (DRV-resistant mutation) and HIV-1NL4_3 variant (wild type) for 18d was detected. The molecular docking study of 18d with HIV-1 protease (PDB-ID: 1T3R, www.rcsb.org) revealed possible binding mode with the HIV-1 protease. These results suggested the validity of introducing phenol-derived moieties into the P2 ligand and deserve further optimization which was of great value for future discovery of novel HIV-1 protease.

Modeling and Analysis of HIV-1 Pol Polyprotein as a Case Study for Predicting Large Polyprotein Structures.

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

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors.

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.

Preclinical characterization of a non-peptidomimetic HIV protease inhibitor with improved metabolic stability.

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

Advanced molecular mechanisms of modified DRV compounds in targeting HIV-1 protease mutations and interrupting monomer dimerization.

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.