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.

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.

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

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.

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

A series of potent HIV-1 protease inhibitors, containing diverse piperidine analogues as the P2-ligands, 4-substituted phenylsulfonamides as the P2\'-ligands and a hydrophobic cyclopropyl group as the P1\'-ligand, were designed, synthesized and evaluated in this work. Among these twenty-four target compounds, many of them exhibited excellent activity against HIV-1 protease with half maximal inhibitory concentration (IC50) values below 20 nM. Particularly, compound 22a containing a (R)-piperidine-3-carboxamide as the P2-ligand and a 4-methoxylphenylsulfonamide as the P2\'-ligand exhibited the most effective inhibitory activity with an IC50 value of 3.61 nM. More importantly, 22a exhibited activity with inhibition of 42% and 26% against wild-type and Darunavir (DRV)-resistant HIV-1 variants, respectively. Additionally, the molecular docking of 22a with HIV-1 protease provided insight into the ligand-binding properties, which was of great value for further study.

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

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

GSK3532795 (formerly BMS955176) is a second-generation maturation inhibitor (MI) that progressed through a Phase 2b study for treatment of HIV-1 infection. Resistance development to GSK3532795 was evaluated through in vitro methods and was correlated with information obtained in a Phase 2a proof-of-concept study in HIV-1 infected participants. Both low and high concentrations of GSK3532795 were used for selections in vitro, and reduced susceptibility to GSK3532795 mapped specifically to amino acids near the capsid/ spacer peptide 1 (SP1) junction, the cleavage of which is blocked by MIs. Two key substitutions, A364V or V362I, were selected, the latter requiring secondary substitutions to reduce susceptibility to GSK3532795. Three main types of secondary substitutions were observed, none of which reduced GSK3532795 susceptibility in isolation. The first type was in the capsid C-terminal domain and downstream SP1 region (including (Gag numbering) R286K, A326T, T332S/N, I333V and V370A/M). The second, was an R41G substitution in viral protease that occurred with V362I. The third was seen in the capsid N-terminal domain, within the cyclophilin A binding domain (V218A/M, H219Q and G221E). H219Q increased viral replication capacity and reduced susceptibility of poorly growing viruses. In the Phase 2a study, a subset of these substitutions was also observed at baseline and some were selected following GSK35323795 treatment in HIV-1-infected participants.

Aspartate proteases are potential targets for various diseases, and many of their inhibitors are FDA-approved drugs. However, these peptidomimetic and reversibly bound drugs become ineffective upon prolonged use. Attempts have been made to design and synthesize various nonpeptidic epoxide-based irreversible inhibitors to combat the drug-resistance enigma. Here, we study the mechanism of epoxide ring opening in two widely studied aspartate proteases, HIV-1 protease and pepsin. Our results from QM/MM molecular dynamics show that the epoxide ring opening in aspartate proteases follow a two-step mechanism with the formation of an oxyanion intermediate, stabilized by a set of water molecules in the protein active site. These water molecules by virtue of \"low-barrier hydrogen bonds\" with the epoxide ring reduce the intrinsic reaction barrier while remaining structurally unperturbed, thus playing a cocatalytic role. We validated our results by reproducing the experimentally observed protease/pepsin-epoxide covalent complexes as end products. The observed stability of our oxyanion intermediate in a four-water-coordinated state is also consistent with the reported stable state of the hydroxide ion in water as OH-(H2O)4. Our study could pave the way for the design of new class \"HIV protease irreversible inhibitors\" from the acquired knowledge of the structures of intermediate and transition states traced during the explored reaction mechanism.