Project description: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.
Project description:BackgroundDrug resistance in HIV is the major problem limiting effective antiviral therapy. Computational techniques for predicting drug resistance profiles from genomic data can accelerate the appropriate choice of therapy. These techniques can also be used to select protease mutants for experimental studies of resistance and thereby assist in the development of next-generation therapies.ResultsThe machine learning produced highly accurate and robust classification of HIV protease resistance. Genotype data were mapped to the enzyme structure and encoded using Delaunay triangulation. Generative machine learning models trained on one inhibitor could classify resistance from other inhibitors with varying levels of accuracy. Generally, the accuracy was best when the inhibitors were chemically similar.ConclusionsRestricted Boltzmann Machines are an effective machine learning tool for classification of genomic and structural data. They can also be used to compare resistance profiles of different protease inhibitors.
Project description: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.
Project description:The reaction of HIV protease to inhibitor therapy is characterized by the emergence of complex mutational patterns which confer drug resistance. The response of HIV protease to drugs often involves both primary mutations that directly inhibit the action of the drug, and a host of accessory resistance mutations that may occur far from the active site but may contribute to restoring the fitness or stability of the enzyme. Here we develop a probabilistic approach based on connected information that allows us to study residue, pair level and higher-order correlations within the same framework.We apply our methodology to a database of approximately 13,000 sequences which have been annotated by the treatment history of the patients from which the samples were obtained. We show that including pair interactions is essential for agreement with the mutational data, since neglect of these interactions results in order-of-magnitude errors in the probabilities of the simultaneous occurence of many mutations. The magnitude of these pair correlations changes dramatically between sequences obtained from patients that were or were not exposed to drugs. Higher-order effects make a contribution of as much as 10% for residues taken three at a time, but increase to more than twice that for 10 to 15-residue groups. The sequence data is insufficient to determine the higher-order effects for larger groups. We find that higher-order interactions have a significant effect on the predicted frequencies of sequences with large numbers of mutations. While relatively rare, such sequences are more prevalent after multi-drug therapy. The relative importance of these higher-order interactions increases with the number of drugs the patient had been exposed to.Correlations are critical for the understanding of mutation patterns in HIV protease. Pair interactions have substantial qualitative effects, while higher-order interactions are individually smaller but may have a collective effect. Together they lead to correlations which could have an important impact on the dynamics of the evolution of cross-resistance, by allowing the virus to pass through otherwise unlikely mutational states. These findings also indicate that pairwise and possibly higher-order effects should be included in the models of protein evolution, instead of assuming that all residues mutate independently of one another.
Project description:The protease from type 1 human immunodeficiency virus (HIV-1) is a critical drug target against which many therapeutically useful inhibitors have been developed; however, the set of viral strains in the population has been shifting to become more drug-resistant. Because indirect effects are contributing to drug resistance, an examination of the dynamic structures of a wild-type and a mutant could be insightful. Consequently, this study examined structural properties sampled during 22 nsec, all atom molecular dynamics (MD) simulations (in explicit water) of both a wild-type and the drug-resistant V82F/I84V mutant of HIV-1 protease. The V82F/I84V mutation significantly decreases the binding affinity of all HIV-1 protease inhibitors currently used clinically. Simulations have shown that the curling of the tips of the active site flaps immediately results in flap opening. In the 22-nsec MD simulations presented here, more frequent and more rapid curling of the mutant's active site flap tips was observed. The mutant protease's flaps also opened farther than the wild-type's flaps did and displayed more flexibility. This suggests that the effect of the mutations on the equilibrium between the semiopen and closed conformations could be one aspect of the mechanism of drug resistance for this mutant. In addition, correlated fluctuations in the active site and periphery were noted that point to a possible binding site for allosteric inhibitors.
Project description:In this preliminary study we show that in 2008, 3 years after antiretroviral therapy was introduced into the Karonga District, Malawi, a greater than expected number of drug-naive individuals have been infected with HIV-1 subtype C virus harboring major and minor drug resistance mutations (DRMs). From a sample size of 40 reverse transcriptase (RT) consensus sequences from drug-naive individuals we found five showing NRTI and four showing NNRTI mutations with one individual showing both. From 29 protease consensus sequences, again from drug-naive individuals, we found evidence of minor DRMs in three. Additional major and minor DRMs were found in clonal sequences from a number of individuals that were not present in the original consensus sequences. This clearly illustrates the importance of sequencing multiple HIV-1 variants from individuals to fully assess drug resistance.
Project description:BackgroundDrug resistance is a critical problem limiting effective antiviral therapy for HIV/AIDS. Computational techniques for predicting drug resistance profiles from genomic data can accelerate the appropriate choice of therapy. These techniques can also be used to identify protease mutants for experimental studies of resistance and thereby assist in the development of next-generation therapies. Few studies, however, have assessed the evolution of resistance from genotype-phenotype data.ResultsThe machine learning produced highly accurate and robust classification of resistance to HIV protease inhibitors. Genotype data were mapped to the enzyme structure and encoded using Delaunay triangulation. Estimates of evolutionary relationships, based on this encoding, and using Minimum Spanning Trees, showed clusters of mutations that closely resemble the wild type. These clusters appear to evolve uniquely to more resistant phenotypes.ConclusionsUsing the triangulation metric and spanning trees results in paths that are consistent with evolutionary theory. The majority of the paths show bifurcation, namely they switch once from non-resistant to resistant or from resistant to non-resistant. Paths that lose resistance almost uniformly have far lower levels of resistance than those which either gain resistance or are stable. This strongly suggests that selection for stability in the face of a rapid rate of mutation is as important as selection for resistance in retroviral systems.
Project description:The diversity of human immunodeficiency virus type 1 (HIV-1) has given rise to multiple subtypes and recombinant strains. The majority of research into antiretroviral agents and drug resistance has been performed on subtype B viruses, yet non-subtype B strains are responsible for 90% of global infections. Although it seems that combination antiretroviral regimens are effective against all HIV-1 subtypes, there is emerging evidence of subtype differences in drug resistance, relevant to antiretroviral strategies in different parts of the world. For this purpose, extensive sampling of HIV genetic diversity, curation and analyses are required to inform antiretroviral strategies in different parts of the world.
Project description:Antiviral inhibitors of HIV-1 protease are a notable success of structure-based drug design and have dramatically improved AIDS therapy. Analysis of the structures and activities of drug resistant protease variants has revealed novel molecular mechanisms of drug resistance and guided the design of tight-binding inhibitors for resistant variants. The plethora of structures reveals distinct molecular mechanisms associated with resistance: mutations that alter the protease interactions with inhibitors or substrates; mutations that alter dimer stability; and distal mutations that transmit changes to the active site. These insights will inform the continuing design of novel antiviral inhibitors targeting resistant strains of HIV.
Project description:NMR and MD simulations have demonstrated that the flaps of HIV-1 protease (HIV-1p) adopt a range of conformations that are coupled with its enzymatic activity. Previously, a model was created for an allosteric site located between the flap and the core of HIV-1p, called the Eye site (Biopolymers 2008, 89, 643-652). Here, results from our first study were combined with a ligand-based, lead-hopping method to identify a novel compound (NIT). NIT inhibits HIV-1p, independent of the presence of an active-site inhibitor such as pepstatin A. Assays showed that NIT acts on an allosteric site other than the dimerization interface. MD simulations of the ligand-protein complex show that NIT stably binds in the Eye site and restricts the flaps. That bound state of NIT is consistent with a crystal structure of similar fragments bound in the Eye site (Chem. Biol. Drug Des. 2010, 75, 257-268). Most importantly, NIT is equally potent against wild-type and a multidrug-resistant mutant of HIV-1p, which highlights the promise of allosteric inhibitors circumventing existing clinical resistance.