Structural insights into the South African HIV-1 subtype C protease: impact of hinge region dynamics and flap flexibility in drug resistance

Mechanism of drug resistance in HIV-1 protease subtype C in the presence of Atazanavir

AIDS is one of the deadliest diseases in the history of humankind caused by HIV. Despite the technological development, curtailing the viral infection inside human host still remains a challenge. Therapies such as HAART uses a combination of drugs to inhibit the viral activity. One of the important targets includes HIV protease and inhibiting its activity will minimize the production of mature structural proteins. However, the genetic diversity and the occurrence of drug resistant mutations adds complexity to effective drug design. In this study, we aimed at understanding the drug binding mechanism of one such subtype, namely subtype C and its insertion variant L38HL. We performed multiple molecular dynamics simulations along with binding free energy analysis of wild-type and L38HL bound to Atazanavir (ATV). From the analysis, we revealed that the insertion alters the hydrogen bond and hydrophobic interaction networks. The alterations in the interaction networks increase flexibility at the hinge-fulcrum interface. Further, the effects of these changes affect flap tip curling. Moreover, the changes in the hinge-fulcrum-cantilever interface alters the concerted motion of the functional regions leading to change in the direction of flap movement thus causing a subtle change in the active site volume. Additionally, formation of intramolecular hydrogen bonds in the ATV docked to L38HL restricted the movement of R1 and R2 groups thereby altering the interactions. Overall, the changes in the flexibility of flap together with the changes in the active site volume and compactness of the ligand provide insights for increased binding affinity of ATV with L38HL.

HIV Protease Hinge Region Insertions at Codon 38 Affect Enzyme Kinetics, Conformational Stability and Dynamics

HIV-1 protease is essential for the production of mature, infectious virions and is a major target in antiretroviral therapy. We successfully purified a HIV-1 subtype C variant, L38↑N↑L- 4, containing an insertion of asparagine and leucine at position 38 without the four background mutations - K20R, E35D, R57K, V82I using a modified purification protocol. Isothermal titration calorimetry indicated that 50% of the variant protease sample was in the active conformation compared to 62% of the wild type protease. The secondary structure composition of the variant protease was unaffected by the double insertion. The specific activity and kcat values of the variant protease were approximately 50% lower than the wild type protease values. The variant protease also exhibited a 1.6-fold increase in kcat/KM when compared to the wild type protease. Differential scanning calorimetry showed a 5 °C increase in Tm of the variant protease, indicating the variant was more stable than the wild type. Molecular dynamics simulations indicated the variant was more stable and compact than the wild type protease. A 3-4% increase in the flexibility of the hinge regions of the variant protease was observed. In addition, increased flexibility of the flaps, cantilever and fulcrum regions of the variant protease B chain was observed. The variant protease sampled only the closed flap conformation indicating a potential mechanism for drug resistance. The present study highlights the direct impact of a double amino acid insertion in hinge region on enzyme kinetics, conformational stability and dynamics of an HIV-1 subtype C variant protease.

Understanding Drug Resistance of Wild-Type and L38HL Insertion Mutant of HIV-1 C Protease to Saquinavir

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.

Non-active site mutations in the HIV protease: Diminished drug binding affinity is achieved through modulating the hydrophobic sliding mechanism

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.

Improved Production of HIV-1 Subtype C Protease from Transgenic E. Coli

Background:

Human Immunodeficiency Virus 1 (HIV-1) subtype C is responsible for the majority of infections of patients in Southern Africa. The HIV protease is a primary target for the development of highly efficient anti-retroviral pharmaceuticals because of its pivotal role in the maturation of the virus in the host cell. For target validation of novel HIV protease inhibitors, there is a need for the availability of an abundance of this protease.

Objective:

This study reports an optimized method to produce HIV-1 protease derived from HIV-1 subtype C.

Methods:

It involves the use of a transgenic E. coli strain that overexpresses the native form of the enzyme via inclusion bodies. A stringent method for the isolation, purification, and renaturation resulted in the production of highly pure active HIV-1 protease. In order to facilitate an increase in protease yields, an optimized growth strategy was developed. In this regard, a chemically defined medium with lower glucose content and devoid of essential amino acids of the TCA cycle was used as an alternative to the widely used nutrient-rich Luria Bertani (LB) medium.

Results:

Results indicated an increase in protease yield up to twice the amount, thereby making this medium an attractive alternative for increasing biomass and HIV protease production for future research.

Conclusion:

An optimized method for HIV-1 protease derived from HIV-1 subtype C production using chemically defined media was established. This was achieved using a known method to isolate and purify the enzyme with the use of a specialized feeding strategy.

Investigating potency of TMC-126 against wild-type and mutant variants of HIV-1 protease: a molecular dynamics and free energy study

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.

Elasticity-Associated Functionality and Inhibition of the HIV Protease

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.

Cantilever-centric mechanism of cooperative non-active site mutations in HIV protease: Implications for flap dynamics

The HIV-1 protease is an important drug target in antiretroviral therapy due to the crucial role it plays in viral maturation. A greater understanding of the dynamics of the protease as a result of drug-induced mutations has been successfully elucidated using computational models in the past. We performed induced-fit docking studies and molecular dynamics simulations on the wild-type South African HIV-1 subtype C protease and two non-active site mutation-containing protease variants; HP3 PR and HP4 PR. The HP3 PR contained the I13V, I62V, and V77I mutations while HP4 PR contained the same mutations with the addition of the L33F mutation. The simulations were initiated in a cubic cell universe containing explicit solvent, with the protease variants beginning in the fully closed conformation. The trajectory for each simulation totalled 50 ns. The results indicate that the mutations increase the dynamics of the flap, hinge, fulcrum and cantilever regions when compared to the wild-type protease while in complex with protease inhibitors. Specifically, these mutations result in the protease favouring the semi-open conformation when in complex with inhibitors. Moreover, the HP4 PR adopted curled flap tip conformers which coordinated several water molecules into the active site in a manner that may reduce inhibitor binding affinity. The mutations affected the thermodynamic landscape of inhibitor binding as there were fewer observable chemical contacts between the mutated variants and saquinavir, atazanavir and darunavir. These data help to elucidate the biophysical basis for the selection of cooperative non-active site mutations by the HI virus.

Rutin and flavone analogs as prospective SARS-CoV-2 main protease inhibitors: In silico drug discovery study

Coronavirus disease 2019 (COVID-19) is a new pandemic characterized by quick spreading and illness of the respiratory system. To date, there is no specific therapy for Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2). Flavonoids, especially rutin, have attracted considerable interest as a prospective SARS-CoV-2 main protease (M^pro) inhibitor. In this study, a database containing 2017 flavone analogs was prepared and screened against SARS-CoV-2 M^pro using the molecular docking technique. According to the results, 371 flavone analogs exhibited good potency towards M^pro with docking scores less than −9.0 kcal/mol. Molecular dynamics (MD) simulations, followed by molecular mechanics-generalized Born surface area (MM/GBSA) binding energy calculations, were performed for the top potent analogs in complex with M^pro. Compared to rutin, PubChem-129-716-607 and PubChem-885-071-27 showed better binding affinities against SARS-CoV-2 M^pro over 150 ns MD course with ΔG_binding values of −69.0 and −68.1 kcal/mol, respectively. Structural and energetic analyses demonstrated high stability of the identified analogs inside the SARS-CoV-2 M^pro active site over 150 ns MD simulations. The oral bioavailabilities of probable SARS-CoV-2 M^pro inhibitors were underpinned using drug-likeness parameters. A comparison of the binding affinities demonstrated that the MM/GBSA binding energies of the identified flavone analogs were approximately three and two times less than those of lopinavir and baicalein, respectively. In conclusion, PubChem-129-716-607 and PubChem-885-071-27 are promising anti-COVID-19 drug candidates that warrant further clinical investigations.

Exploring the concerted mechanistic pathway for HIV-1 PR-substrate revealed by umbrella sampling simulation

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.

Drug susceptibility and replication capacity of a rare HIV-1 subtype C protease hinge region variant

BACKGROUND

Protease inhibitors form the main component of second-line antiretroviral treatment in South Africa. Despite their efficacy, mutations arising within the HIV-1 gag and protease coding regions contribute to the development of resistance against this class of drug. In this paper we investigate a South African HIV-1 subtype C Gag-protease that contains a hinge region mutation and insertion (N37T↑V).

METHODS

In vitro single-cycle drug susceptibility and viral replication capacity assays were performed on W1201i, a wild-type reference isolate (MJ4) and a chimeric construct (MJ4GagN37T↑VPR). Additionally, enzyme assays were performed on the N37T↑V protease and a wild-type reference protease.

RESULTS

W1201i showed a small (threefold), but significant (P<0.0001) reduction in drug susceptibility to darunavir compared with MJ4. Substitution of W1201i-Gag with MJ4-Gag resulted in an additional small (twofold), but significant (P<0.01) reduction in susceptibility to lopinavir and atazanavir. The W1201i pseudovirus had a significantly (P<0.01) reduced replication capacity (16.4%) compared with the MJ4. However, this was dramatically increased to 164% (P<0.05) when W1201i-Gag was substituted with MJ4-Gag. Furthermore, the N37T↑V protease displayed reduced catalytic processing compared with the SK154 protease.

CONCLUSIONS

Collectively, these data suggest that the N37T↑V mutation and insertion increases viral infectivity and decreases drug susceptibility. These variations are classified as secondary mutations, and indirectly impact inhibitor binding, enzyme fitness and enzyme stability. Additionally, polymorphisms arising in Gag can modify the impact of protease with regards to viral replication and susceptibility to protease inhibitors.

Two-Step Preparation of Highly Pure, Soluble HIV Protease from Inclusion Bodies Recombinantly Expressed in Escherichia coli

Heterologous expression of exogenous proteases in Escherichia coli often results in the formation of insoluble inclusion bodies. When sequestered into inclusion bodies, the functionality of the proteases is minimized. To be characterized structurally and functionally, however, proteases must be obtained in their native conformation. HIV protease is readily expressed as inclusion bodies, but must be recovered from the inclusion bodies. This protocol describes an efficient method for recovering HIV protease from inclusion bodies, as well as refolding and purifying the protein. HIV protease-containing inclusion bodies are treated with 8 M urea and purified via cation-exchange chromatography. Subsequent refolding by buffer exchange via dialysis and further purification by anion-exchange chromatography produces highly pure HIV protease that is functionally active. © 2020 by John Wiley & Sons, Inc. Basic Protocol: Recovery, refolding, and purification of HIV protease from inclusion bodies Support Protocol 1: Expression and extraction of inclusion bodies containing HIV protease expressed in Escherichia coli Support Protocol 2: Determination of the active site concentration of HIV protease via isothermal titration calorimetry.

Concerted hydrolysis mechanism of HIV-1 natural substrate against subtypes B and C-SA PR: insight through molecular dynamics and hybrid QM/MM studies

It is well known that understanding the catalytic mechanism of HIV-1 PR is the rationale on which its inhibitors were developed; therefore, a better understanding of the mechanism of natural substrate hydrolysis is important. Herein, the reaction mechanism of HIV-1 natural substrates with subtypes B and common mutant in South Africa (subtype C-SA) protease were studied through transition state modelling, using a general acid-general base (GA-GB) one-step concerted process. The activation free energies of enzyme-substrate complexes were compared based on their rate of hydrolysis using a two-layered ONIOM (B3LYP/6-31++G(d,p):AMBER) method. We expanded our computational model to obtain a better understanding of the mechanism of hydrolysis as well as how the enzyme recognises or chooses the cleavage site of the scissile bonds. Using this model, a potential substrate-based inhibitor could be developed with better potency. The calculated activation energies of natural substrates in our previous study correlated well with experimental data. A similar trend was observed for the Gag and Gag-Pol natural substrates in the present work for both enzyme complexes except for the PR-RT substrate. Natural bond orbital (NBO) analysis was also applied to determine the extent of charge transfer within the QM part of both enzymes considered and the PR-RT natural substrate. The result of this study shows that the method can be utilized as a dependable computational technique to rationalize lead compounds against specific targets.

Kinetic and thermodynamic characterisation of HIV-protease inhibitors against E35D↑G↑S mutant in the South African HIV-1 subtype C protease

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

Low Frequency of Protease Inhibitor Resistance Mutations and Insertions in HIV-1 Subtype C Protease Inhibitor-Naïve Sequences

Human immunodeficiency virus-1 (HIV-1) protease sequences from 2,225 protease inhibitor (PI)-naïve HIV-1 subtype C-infected individuals collected over a 14-year period were analyzed for polymorphisms. Over 50% of sequences differed from an HIV-1 subtype B consensus sequence at 8 of the 99 amino acids at residues 12, 15, 19, 36, 41, 69, 89, and 93, but not in the functionally important regions. The frequency of primary resistance and accessory mutations occurred in <1% of the sequences. Of note, 11 sequences (0.5%) harbored amino acid insertions between residues 36 and 39, located in the elbow of the flap region. The insertions were found throughout the 13-year period. Occurrence of insertions in subtype C viruses is rare and viruses remain sensitive to currently used PIs (lopinavir/r, atazanavir/r, and darunavir/r). However, ongoing characterization of isolates is required to identify changes that may impact PI treatment since PIs are part of standard SA regimens.

Exploring the flap dynamics of the South African HIV subtype C protease in presence of FDA-approved inhibitors: MD study

HIV-1 protease (HIV PR) is considered as one of the most attractive targets for the treatment of HIV and the impact of flap dynamics of HIV PR on the binding affinities of protease inhibitors (PIs) is a crucial ongoing research field. Recently, our research group evaluated the binding affinities of different FDA approved PIs against the South African HIV-1 subtype C (C-SA) protease (PR). The CSA-HIV PR displayed weaker binding affinity for most of the clinical PIs compared to HIV-1 B subtype for West and Central Europe, the Americas. In the current work, the flap dynamics of four different systems of HIV-1 C-SA PR complexed to FDA approved second generation PIs and its impact on binding was explored over the molecular dynamic trajectories. It was observed that the interactions of the selected drugs with the binding site residues of the protease may not be the major contributor for affinity towards PIs. Various post-MD analyses were performed, also entropic contributions, solvation free energies and hydrophobic core formation interactions were studied to assess how the flap dynamics of C-SA PR which is affected by such factors. From these contributions, large van der Waals interactions and low solvation free energies were found to be major factors for the higher activity of ATV against C-SA HIV PR. Furthermore, a comparatively stable hydrophobic core may be responsible for higher stability of the PR flaps of the ATV complex. The outcome of this study provides significant guidance to how the flap dynamics of C-SA PR is affected by various factors as a result of the binding affinity of various protease inhibitors. It will also assist with the design of potent inhibitors against C-SA HIV PR that apart from binding in the active site of PR can interacts with the flaps to prevent opening of the flaps resulting in inactivation of the protease.

Overexpression, Purification and Functional Characterisation of Wild-Type HIV-1 Subtype C Protease and Two Variants Using a Thioredoxin and His-Tag Protein Fusion System

In recent years, various strategies have been used to overexpress and purify HIV-1 protease because it is an essential drug target in anti-retroviral therapy. Obtaining sufficient quantities of the enzyme, however, remains challenging. Overexpression of large quantities is prevented due to the enzyme's autolytic nature and its inherent cytotoxicity in Escherichia coli cells. Here, we describe a novel HIV-1 protease purification method using a thioredoxin-hexahistidine fusion system for the wild-type and two variant proteases. The fusion proteases were overexpressed in E. coli and recovered by immobilised metal ion affinity chromatography. The proteases were cleaved from the fusion constructs using thrombin. When compared to the standard overexpression and purification protocol in use in our laboratory, the expression of the fusion-derived wild-type protease was increased from 0.83 to 2.5 mg/l of culture medium. The expression levels of the two variant proteases ranged from 1.5 to 2 mg/l of culture medium. The fusion wild-type and variant proteases were inactive before the cleavage of the thioredoxin-hexahistidine fusion tag as no enzymatic activity was observed. The proteases were, however, active after cleavage of the tag. The novel thioredoxin-hexahistidine fusion system, therefore, enables the successful overexpression and purification of catalytically active HIV-1 proteases.

Investigation of the binding free energies of FDA approved drugs against subtype B and C-SA HIV PR: ONIOM approach

Human immune virus subtype C is the most widely spread HIV subtype in Sub-Sahara Africa and South Africa. A profound structural insight on finding potential lead compounds is therefore necessary for drug discovery. The focus of this study is to rationalize the nine Food and Drugs Administration (FDA) HIV antiviral drugs complexed to subtype B and C-SA PR using ONIOM approach. To achieve this, an integrated two-layered ONIOM model was used to optimize the geometrics of the FDA approved HIV-1 PR inhibitors for subtype B. In our hybrid ONIOM model, the HIV-1 PR inhibitors as well as the ASP 25/25' catalytic active residues were treated at high level quantum mechanics (QM) theory using B3LYP/6-31G(d), and the remaining HIV PR residues were considered using the AMBER force field. The experimental binding energies of the PR inhibitors were compared to the ONIOM calculated results. The theoretical binding free energies (?G_bind) for subtype B follow a similar trend to the experimental results, with one exemption. The computational model was less suitable for C-SA PR. Analysis of the results provided valuable information about the shortcomings of this approach. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide much improved binding energies for complex enzyme drug interactions.

Flap Dynamics in Aspartic Proteases: A Computational Perspective

Recent advances in biochemistry and drug design have placed proteases as one of the critical target groups for developing novel small-molecule inhibitors. Among all proteases, aspartic proteases have gained significant attention due to their role in HIV/AIDS, malaria, Alzheimer's disease, etc. The binding cleft is covered by one or two β-hairpins (flaps) which need to be opened before a ligand can bind. After binding, the flaps close to retain the ligand in the active site. Development of computational tools has improved our understanding of flap dynamics and its role in ligand recognition. In the past decade, several computational approaches, for example molecular dynamics (MD) simulations, coarse-grained simulations, replica-exchange molecular dynamics (REMD) and metadynamics, have been used to understand flap dynamics and conformational motions associated with flap movements. This review is intended to summarize the computational progress towards understanding the flap dynamics of proteases and to be a reference for future studies in this field.

Virological failure in patients with HIV-1 subtype C receiving antiretroviral therapy: an analysis of a prospective national cohort in Sweden

BACKGROUND

People with HIV-1 in low-income and middle-income countries increasingly need second-line regimens with boosted protease inhibitors. However, data are scarce for treatment response in patients with HIV-1 subtype C (HIV-1C), which is predominant in these regions. We aimed to examine factors associated with virological failure in patients in a standardised national health-care setting.

METHODS

We analysed data for participants in InfCare HIV, a prospective national cohort that includes more than 99% of people with HIV in Sweden. We extracted data for the cohort from the InfCare HIV database on Jan 14, 2015. Baseline was initiation of antiretroviral therapy. We used logistic regression to assess factors associated with primary virological failure (failure to suppress HIV-1 within 9 months) in patients with HIV-1B and HIV-1C and calculated odds ratios (OR) for failure. We also used Cox regression models to calculate hazard ratios (HR) for time-to-secondary virological failure (detectable viral load after initial virological suppression). We did homology-based molecular modelling to assess docking.

FINDINGS

We included 1077 patients with HIV-1B and 596 with HIV-1C. In multivariate regression analysis, pre-therapy higher viral load (OR 1·82, 95% CI 1·49-2·21; p<0·0001), subtype C infection (1·75, 1·06-2·88; p=0·028), and boosted protease inhibitor-based regimens (1·55, 1·45-2·11; p=0·004) were associated with increased risk of primary virological failure. Individuals with HIV-1C who were given therapy with boosted protease inhibitors had earlier time-to-secondary virological failure than did those with HIV-1B given similar regimens (adjusted HR 1·92, 95% CI 1·30-2·83; p=0·002). Molecular modelling suggested lower affinity for protease inhibitors to HIV-1C protease than to HIV-1B.

INTERPRETATION

Our findings suggest an increased risk of virological failure in patients with HIV-1C, especially in those on boosted protease inhibitor-based regimens. Future studies should further dissect the biochemical and viral mechanisms of resistance to protease inhibitors in patients with non-B subtypes of HIV-1, including clinical studies to assess the efficacy of boosted protease inhibitor-based regimens in low-income and middle income countries.

FUNDING

Karolinska Institutet Research Foundation, Swedish Research Council, Stockholm County Council, Swedish Physicians against AIDS, US National Institutes of Health, University of Missouri.

Purification and characterization of naturally occurring HIV-1 (South African subtype C) protease mutants from inclusion bodies

Human immunodeficiency virus (HIV) infections in sub-Saharan Africa represent about 56% of global infections. Many studies have targeted HIV-1 protease for the development of drugs against AIDS. Recombinant HIV-1 protease is used to screen new drugs from synthetic compounds or natural substances. Along with the wild type (C-SA) we also over-expressed and characterized two mutant forms from patients that had shown resistance to protease inhibitors. Using recombinant DNA technology, we constructed three recombinant plasmids in pGEX-6P-1 and expressed them containing a sequence encoding wild type HIV protease and two mutants (I36T↑T contains 100 amino acids and L38L↑N↑L contains 101 amino acids). These recombinant proteins were isolated from inclusion bodies by using QFF anion exchange and GST trap columns. In SDS-PAGE, we obtained these HIV proteases as single bands of approximately 11.5, 11.6 and 11.7 kDa for the wild type, I36T↑Tand L38L↑N↑L mutants, respectively. The enzyme was recovered efficiently (0.25 mg protein/L of Escherichia coli culture) and had high specific activity of 2.02, 2.20 and 1.33 μmol min(-1) mg(-1) at an optimal pH of 5 and temperature of 37 °C for the wild type, I36T↑T and L38L↑N↑L, respectively. The method employed here provides an easy and rapid purification of the HIV-1(C-SA) protease from the inclusion bodies, with high yield and high specific activities.

Binding Free Energy Calculations of Nine FDA-approved Protease Inhibitors Against HIV-1 Subtype C I36T↑T Containing 100 Amino Acids Per Monomer

In this work, have investigated the binding affinities of nine FDA-approved protease inhibitor drugs against a new HIV-1 subtype C mutated protease, I36T↑T. Without an X-ray crystal structure, homology modelling was used to generate a three-dimensional model of the protease. This and the inhibitor models were employed to generate the inhibitor/I36T↑T complexes, with the relative positions of the inhibitors being superimposed and aligned using the X-ray crystal structures of the inhibitors/HIV-1 subtype B complexes as a reference. Molecular dynamics simulations were carried out on the complexes to calculate the average binding free energies for each inhibitor using the molecular mechanics generalized Born surface area (MM-GBSA) method. When compared to the binding free energies of the HIV-1 subtype B and subtype C proteases (calculated previously by our group using the same method), it was clear that the I36T↑T proteases mutations and insertion had a significant negative effect on the binding energies of the non-pepditic inhibitors nelfinavir, darunavir and tipranavir. On the other hand, ritonavir, amprenavir and indinavir show improved calculated binding energies in comparison with the corresponding data for wild-type C-SA protease. The computational model used in this study can be used to investigate new mutations of the HIV protease and help in establishing effective HIV drug regimes and may also aid in future protease drug design.

Multi-drug resistance profile of PR20 HIV-1 protease is attributed to distorted conformational and drug binding landscape: molecular dynamics insights

The PR20 HIV-1 protease, a variant with 20 mutations, exhibits high levels of multi-drug resistance; however, to date, there has been no report detailing the impact of these 20 mutations on the conformational and drug binding landscape at a molecular level. In this report, we demonstrate the first account of a comprehensive study designed to elaborate on the impact of these mutations on the dynamic features as well as drug binding and resistance profile, using extensive molecular dynamics analyses. Comparative MD simulations for the wild-type and PR20 HIV proteases, starting from bound and unbound conformations in each case, were performed. Results showed that the apo conformation of the PR20 variant of the HIV protease displayed a tendency to remain in the open conformation for a longer period of time when compared to the wild type. This led to a phenomena in which the inhibitor seated at the active site of PR20 tends to diffuse away from the binding site leading to a significant change in inhibitor-protein association. Calculating the per-residue fluctuation (RMSF) and radius of gyration, further validated these findings. MM/GBSA showed that the occurrence of 20 mutations led to a drop in the calculated binding free energies (ΔGbind) by ~25.17 kcal/mol and ~5 kcal/mol for p2-NC, a natural peptide substrate, and darunavir, respectively, when compared to wild type. Furthermore, the residue interaction network showed a diminished inter-residue hydrogen bond network and changes in inter-residue connections as a result of these mutations. The increased conformational flexibility in PR20 as a result of loss of intra- and inter-molecular hydrogen bond interactions and other prominent binding forces led to a loss of protease grip on ligand. It is interesting to note that the difference in conformational flexibility between PR20 and WT conformations was much higher in the case of substrate-bound conformation as compared to DRV. Thus, developing analogues of DRV by retaining its key pharmacophore features will be the way forward in the search for novel protease inhibitors against multi-drug resistant strains.

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC₅₀ values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC₅₀) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC₅₀ = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC₅₀ = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆G_solv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC₅₀ data.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-015-0115-5) contains supplementary material, which is available to authorized users.

The impact of active site mutations of South African HIV PR on drug resistance: Insight from molecular dynamics simulations, binding free energy and per-residue footprints

Molecular dynamics simulations and binding free energy calculations were used to provide an understanding of the impact of active site drug-resistant mutations of the South African HIV protease subtype C (C-SA HIV PR), V82A and V82F/I84V on drug resistance. Unique per-residue interaction energy 'footprints' were developed to map the overall drug-binding profiles for the wild type and mutants. Results confirmed that these mutations altered the overall binding landscape of the amino acid residues not only in the active site region but also in the flaps as well. Four FDA-approved drugs were investigated in this study; these include ritonavir (RTV), saquinavir (SQV), indinavir (IDV), and nelfinavir (NFV). Computational results compared against experimental findings were found to be complementary. Against the V82F/I84V variant, saquinavir, indinavir, and nelfinavir lose remarkable entropic contributions relative to both wild-type and V82A C-SA HIV PRs. The per-residue energy 'footprints' and the analysis of ligand-receptor interactions for the drug complexes with the wild type and mutants have also highlighted the nature of drug interactions. The data presented in this study will prove useful in the design of more potent inhibitors effective against drug-resistant HIV strains.

The role of select subtype polymorphisms on HIV-1 protease conformational sampling and dynamics

HIV-1 protease is an essential enzyme for viral particle maturation and is a target in the fight against HIV-1 infection worldwide. Several natural polymorphisms are also associated with drug resistance. Here, we utilized both pulsed electron double resonance, also called double electron-electron resonance, and NMR (15)N relaxation measurements to characterize equilibrium conformational sampling and backbone dynamics of an HIV-1 protease construct containing four specific natural polymorphisms commonly found in subtypes A, F, and CRF_01 A/E. Results show enhanced backbone dynamics, particularly in the flap region, and the persistence of a novel conformational ensemble that we hypothesize is an alternative flap orientation of a curled open state or an asymmetric configuration when interacting with inhibitors.

New insights into the in silico prediction of HIV protease resistance to nelfinavir

The Human Immunodeficiency Virus type 1 protease enzyme (HIV-1 PR) is one of the most important targets of antiretroviral therapy used in the treatment of AIDS patients. The success of protease-inhibitors (PIs), however, is often limited by the emergence of protease mutations that can confer resistance to a specific drug, or even to multiple PIs. In the present study, we used bioinformatics tools to evaluate the impact of the unusual mutations D30V and V32E over the dynamics of the PR-Nelfinavir complex, considering that codons involved in these mutations were previously related to major drug resistance to Nelfinavir. Both studied mutations presented structural features that indicate resistance to Nelfinavir, each one with a different impact over the interaction with the drug. The D30V mutation triggered a subtle change in the PR structure, which was also observed for the well-known Nelfinavir resistance mutation D30N, while the V32E exchange presented a much more dramatic impact over the PR flap dynamics. Moreover, our in silico approach was also able to describe different binding modes of the drug when bound to different proteases, identifying specific features of HIV-1 subtype B and subtype C proteases.