Molecular modeling studies and comparative analysis on structurally similar HTLV and HIV protease using HIV-PR inhibitors

Abstract

Retroviruses are most perilous viral family, which cause much damage to the Homo sapiens. HTLV-1 mechanism found to more similar with HIV-1 and both retroviruses are causative agents of severe and fatal diseases including adult T-cell leukemia (ATL) and the acquired immune deficiency syndrome (AIDS). Both viruses code for a protease (PR) that is essential for replication and therefore represents a key target for drugs interfering with viral infection. In this work, the comparative study of HIV-1 and HTLV-1 PR enzymes through sequence and structural analysis is reported along with approved drugs of HIV-PR. Conformation of each HIV PR drugs have been examined with different parameters of interactions and energy scorings parameters. MD simulations with respect to timescale event of 20?ns favors that, few HIV-PR inhibitors can be more active inside the HTLV-1 PR binding pocket. Overall results suggest that, some of HIV inhibitors like Tipranavir, Indinavir, Darunavir and Amprenavir are having good energy levels with HTLV-1. Due to absence of interactions with MET37, here we report that derivatives of these compounds can be much better inhibitors for targeting HTLV-1 proteolytic activity.     

Litchi chinensis-derived terpenoid as anti-HIV-1 protease agent: structural design from molecular dynamics simulations

The molecular structures of the binding between human immunodeficiency virus-1 protease (HIV-1PR) and various inhibitors including existing extensive natural products extracts have been investigated for anti-HIV drug development. In this study, the binding of HIV-1PR and a terpenoid from Litchi chinensis extracts (3-oxotrirucalla-7,24-dien-21-oic acid) was investigated in order to clarify the inhibition effectiveness of this compound. Molecular dynamics (MD) simulations of HIV-1PR complex with 3-oxotrirucalla-7,24-dien-21-oic acid were performed including water molecules. The MD simulation results indicated the formation of hydrogen bonds between the oxygen atoms of the inhibitor and the catalytic aspartates, which are commonly found in inhibitors–protease complexes. On the other hand, no hydrogen bonding of this particular inhibitor to the flap region was found. In addition, the radial distribution function of water oxygens around the catalytic carboxylate nitrogens of Asp29 and Asp30 suggests that at least one or two water molecules are in the active site region whereas direct interaction of the inhibitor was found for catalytic carboxylate oxygen of Asp25. The results of this simulation, in comparison with the structures of other HIV-PR inhibitor complexes, could lead to a better understanding of the activity of 3-oxotrirucalla-7,24-dien-21-oic acid.     

Comparison of the substrate specificity of the human T-cell leukemia virus and human immunodeficiency virus proteinases.

Human T-cell leukemia virus type-1 (HTLV-1) is associated with a number of human diseases. Based on the therapeutic success of human immunodeficiency virus type 1 (HIV-1) PR inhibitors, the proteinase (PR) of HTLV-1 is a potential target for chemotherapy. To facilitate the design of potent inhibitors, the subsite specificity of HTLV-1 PR was characterized and compared to that of HIV-1 PR. Two sets of substrates were used that contained single amino-acid substitutions in peptides representing naturally occurring cleavage sites in HIV-1 and HTLV-1. The original HIV-1 matrix/capsid cleavage site substrate and most of its substituted peptides were not hydrolyzed by the HTLV-1 enzyme, except for those with hydrophobic residues at the P4 and P2 positions. On the other hand, most of the peptides representing the HTLV-1 capsid/nucleocapsid cleavage site were substrates of both enzymes. A large difference in the specificity of HTLV-1 and HIV-1 proteinases was demonstrated by kinetic measurements, particularly with regard to the S4 and S2 subsites, whereas the S1 subsite appeared to be more conserved. A molecular model of the HTLV-1 PR in complex with this substrate was built, based on the crystal structure of the S9 mutant of Rous sarcoma virus PR, in order to understand the molecular basis of the enzyme specificity. Based on the kinetics of shortened analogs of the HTLV-1 substrate and on analysis of the modeled complex of HTLV-1 PR with substrate, the substrate binding site of the HTLV-1 PR appeared to be more extended than that of HIV-1 PR. Kinetic results also suggested that the cleavage site between the capsid and nucleocapsid protein of HTLV-1 is evolutionarily optimized for rapid hydrolysis.     

Dynamics of Preferential Substrate Recognition in HIV-1 Protease: Redefining the Substrate Envelope

Human immunodeficiency virus type 1 (HIV-1) protease (PR) permits viral maturation by processing the gag and gag-pro-pol polyproteins. HIV-1 PR inhibitors (PIs) are used in combination antiviral therapy but the emergence of drug resistance has limited their efficacy. The rapid evolution of HIV-1 necessitates consideration of drug resistance in novel drug design. Drug-resistant HIV-1 PR variants no longer inhibited efficiently, continue to hydrolyze the natural viral substrates. Though highly diverse in sequence, the HIV-1 PR substrates bind in a conserved three-dimensional shape we termed the substrate envelope. Earlier, we showed that resistance mutations arise where PIs protrude beyond the substrate envelope, because these regions are crucial for drug binding but not for substrate recognition. We extend this model by considering the role of protein dynamics in the interaction of HIV-1 PR with its substrates. We simulated the molecular dynamics of seven PR-substrate complexes to estimate the conformational flexibility of the bound substrates. Interdependence of substrate-protease interactions might compensate for variations in cleavage-site sequences and explain how a diverse set of sequences are recognized as substrates by the same enzyme. This diversity might be essential for regulating sequential processing of substrates. We define a dynamic substrate envelope as a more accurate representation of PR-substrate interactions. This dynamic substrate envelope, described by a probability distribution function, is a powerful tool for drug design efforts targeting ensembles of resistant HIV-1 PR variants with the aim of developing drugs that are less susceptible to resistance.     

Combinatorial design of nonsymmetrical cyclic urea inhibitors of aspartic protease of HIV-1

The aspartic protease (PR) of the human immunodeficiency virus type 1 (HIV-1) is an important target for the design of specific antiviral agents dedicated to treatment of HIV-1 infection. We have employed computer-assisted combinatorial chemistry methods to design a small focused virtual library of nonsymmetrically substituted cyclic urea inhibitors of the PR. Nonsymmetrical compounds with decreased peptidic character were namely found to inhibit the PR with comparable inhibition potencies as their C2-pseudosymmetric counterparts and to possess superior pharmacokinetic properties. To generate the virtual library of fully nonsymmetrical cyclic urea analogs, diverse reagents were selected from databases of available chemicals with characteristics similar to those of the building blocks of known potent PR inhibitors. The X-ray structure of the protease-inhibitor complex PR-XV-638 was used as the receptor model in the structure-based focusing and in silico screening of the virtual library. A target-specific LUDI-type scoring function, parameterized for a QSAR training set of known cyclic urea inhibitors and validated on a set of compounds not included into the training set, was used to predict the inhibition constants (Ki) of the generated analogs toward the HIV-1 PR. The fragments most frequently occurring in the analogs with the highest predicted inhibition potencies (Ki*<10 pM) were then selected to constitute a highly focused library subset containing novel nonsymmetrical cyclic ureas with predicted Ki*s 1 order of magnitude lower than the most potent known cyclic urea inhibitors. ADME properties calculated for the most promising analogs suggested that the cyclic ureas are endowed with a wide range of favorable pharmacokinetic properties, which may favor the discovery of a potent orally administrable antiviral drug.     

Molecular basis for the relative substrate specificity of human immunodeficiency virus type 1 and feline immunodeficiency virus proteases

We have used a random hexamer phage library to delineate similarities and differences between the substrate specificities of human immunodeficiency virus type 1 (HIV-1) and feline immunodeficiency virus (FIV) proteases (PRs). Peptide sequences were identified that were specifically cleaved by each protease, as well as sequences cleaved equally well by both enzymes. Based on amino acid distinctions within the P3-P3' region of substrates that appeared to correlate with these cleavage specificities, we prepared a series of synthetic peptides within the framework of a peptide sequence cleaved with essentially the same efficiency by both HIV-1 and FIV PRs, Ac-KSGVF/VVNGLVK-NH(2) (arrow denotes cleavage site). We used the resultant peptide set to assess the influence of specific amino acid substitutions on the cleavage characteristics of the two proteases. The findings show that when Asn is substituted for Val at the P2 position, HIV-1 PR cleaves the substrate at a much greater rate than does FIV PR. Likewise, Glu or Gln substituted for Val at the P2' position also yields peptides specifically susceptible to HIV-1 PR. In contrast, when Ser is substituted for Val at P1', FIV PR cleaves the substrate at a much higher rate than does HIV-1 PR. In addition, Asn or Gln at the P1 position, in combination with an appropriate P3 amino acid, Arg, also strongly favors cleavage by FIV PR over HIV PR. Structural analysis identified several protease residues likely to dictate the observed specificity differences. Interestingly, HIV PR Asp30 (Ile-35 in FIV PR), which influences specificity at the S2 and S2' subsites, and HIV-1 PR Pro-81 and Val-82 (Ile-98 and Gln-99 in FIV PR), which influence specificity at the S1 and S1' subsites, are residues which are often involved in development of drug resistance in HIV-1 protease. The peptide substrate KSGVF/VVNGK, cleaved by both PRs, was used as a template for the design of a reduced amide inhibitor, Ac-GSGVF Psi(CH(2)NH)VVNGL-NH(2.) This compound inhibited both FIV and HIV-1 PRs with approximately equal efficiency. These findings establish a molecular basis for distinctions in substrate specificity between human and feline lentivirus PRs and offer a framework for development of efficient broad-based inhibitors.     

Structural evidence for effectiveness of darunavir and two related antiviral inhibitors against HIV-2 protease

No drug has been targeted specifically for HIV-2 (human immunodeficiency virus type 2) infection despite its increasing prevalence worldwide. The antiviral HIV-1 (human immunodeficiency virus type 1) protease (PR) inhibitor darunavir and the chemically related GRL98065 and GRL06579A were designed with the same chemical scaffold and different substituents at P2 and P2′ to optimize polar interactions for HIV-1 PR (PR1). These inhibitors are also effective antiviral agents for HIV-2-infected cells. Therefore, crystal structures of HIV-2 PR (PR2) complexes with the three inhibitors have been solved at 1.2-Å resolution to analyze the molecular basis for their antiviral potency. Unusually, the crystals were grown in imidazole and zinc acetate buffer, which formed interactions with the PR2 and the inhibitors. Overall, the structures were very similar to the corresponding inhibitor complexes of PR1 with an RMSD of 1.1 Å on main-chain atoms. Most hydrogen-bond and weaker C-H…O interactions with inhibitors were conserved in the PR2 and PR1 complexes, except for small changes in interactions with water or disordered side chains. Small differences were observed in the hydrophobic contacts for the darunavir complexes, in agreement with relative inhibition of the two PRs. These near-atomic-resolution crystal structures verify the inhibitor potency for PR1 and PR2 and will provide the basis for the development of antiviral inhibitors targeting PR2.     

Atomistic simulations of the HIV-1 protease folding inhibition

Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini β-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.     

Molecular dynamics simulations of HIV-1 protease monomer: Assembly of N-terminus and C-terminus into beta-sheet in water solution

HIV-1 protease (HIV-PR) consists of two identical subunits that are united together through a four-stranded antiparallel beta-sheet formed of the peptide termini of each monomer. Since the active site exists only in the dimer, a strategy that is attracting more and more attention in inhibitor design and which may overcome the serious drug resistance caused by competitive inhibitors is to block the peptide termini of the monomer, thereby interfering with formation of the active dimer. In the present work, we performed several extensive molecular dynamics (MD) simulations of the HIV-PR monomer in water to illustrate its solvated conformation and dynamics behavior. We found that the peptide termini usually assembled into beta-sheet after several nanoseconds' simulation, and became much less flexible. This beta-sheet is stabilized by intramolecular interactions and is not easily disaggregated under the present MD simulation conditions. This transformation may be an important transition during the relaxing and equilibrating of the HIV-PR monomer in aqueous solution, and the terminal beta-sheet may be one of the major conformations of the solvated HIV-PR monomer termini in water. This work may provide new insights into the dynamics behavior and dimerization mechanism of HIV-PR, and more significantly, offer a more rational receptor model for the design and discovery of novel dimerization inhibitors than crystalline structures.     

Molecular basis for substrate recognition and drug resistance from 1.1 to 1.6 angstroms resolution crystal structures of HIV-1 protease mutants with substrate analogs

HIV-1 protease (PR) and two drug-resistant variants--PR with the V82A mutation (PR(V82A)) and PR with the I84V mutation (PR(I84V))--were studied using reduced peptide analogs of five natural cleavage sites (CA-p2, p2-NC, p6pol-PR, p1-p6 and NC-p1) to understand the structural and kinetic changes. The common drug-resistant mutations V82A and I84V alter residues forming the substrate-binding site. Eight crystal structures were refined at resolutions of 1.10-1.60 A. Differences in the PR-analog interactions depended on the peptide sequence and were consistent with the relative inhibition. Analog p6(pol)-PR formed more hydrogen bonds of P2 Asn with PR and fewer van der Waals contacts at P1' Pro compared with those formed by CA-p2 or p2-NC in PR complexes. The P3 Gly in p1-p6 provided fewer van der Waals contacts and hydrogen bonds at P2-P3 and more water-mediated interactions. PR(I84V) showed reduced van der Waals interactions with inhibitor compared with PR, which was consistent with kinetic data. The structures suggest that the binding affinity for mutants is modulated by the conformational flexibility of the substrate analogs. The complexes of PR(V82A) showed smaller shifts of the main chain atoms of Ala82 relative to PR, but more movement of the peptide analog, compared to complexes with clinical inhibitors. PR(V82A) was able to compensate for the loss of interaction with inhibitor caused by mutation, in agreement with kinetic data, but substrate analogs have more flexibility than the drugs to accommodate the structural changes caused by mutation. Hence, these structures help to explain how HIV can develop drug resistance while retaining the ability of PR to hydrolyze natural substrates.     

The L33F darunavir resistance mutation acts as a molecular anchor reducing the flexibility of the HIV-1 protease 30s and 80s loops

HIV-1 protease (PR) is a 99 amino acid protein responsible for proteolytic processing of the viral polyprotein – an essential step in the HIV-1 life cycle. Drug resistance mutations in PR that are selected during antiretroviral therapy lead to reduced efficacy of protease inhibitors (PI) including darunavir (DRV). To identify the structural mechanisms associated with the DRV resistance mutation L33F, we performed X-ray crystallographic studies with a multi-drug resistant HIV-1 protease isolate that contains the L33F mutation (MDR769 L33F). In contrast to other PR L33F DRV complexes, the structure of MDR769 L33F complexed with DRV reported here displays the protease flaps in an open conformation. The L33F mutation increases noncovalent interactions in the hydrophobic pocket of the PR compared to the wild-type (WT) structure. As a result, L33F appears to act as a molecular anchor, reducing the flexibility of the 30s loop (residues 29–35) and the 80s loop (residues 79–84). Molecular anchoring of the 30s and 80s loops leaves an open S1/S1′ subsite and distorts the conserved hydrogen-bonding network of DRV. These findings are consistent with previous reports despite structural differences with regards to flap conformation.     

Inhibition of XMRV and HIV-1 proteases by pepstatin A and acetyl-pepstatin

The kinetic properties of two classical inhibitors of aspartic proteases (PRs), pepstatin A and acetyl-pepstatin, were compared in their interactions with HIV-1 and xenotropic murine leukemia virus related virus (XMRV) PRs. Both compounds are substantially weaker inhibitors of XMRV PR than of HIV-1 PR. Previous kinetic and structural studies characterized HIV-1 PR-acetyl-pepstatin and XMRV PR-pepstatin A complexes and suggested dramatically different binding modes. Interaction energies were calculated for the possible binding modes and suggested a strong preference for the one-inhibitor binding mode for HIV-1 PR-acetyl-pepstatin and the two-inhibitor binding mode for XMRV PR-pepstatin A interactions. Comparison of the molecular models suggested that in the case of XMRV PR the relatively unfavorable interactions at S3' and the favorable interactions at S4 and S4' sites with the statine residues may shift the ground state binding towards the two-inhibitor binding mode, whereas the single molecule ground state binding of statines to the HIV-1 PR appear to be more favorable. The preferred single molecular binding to HIV-1 PR allows the formation of the transition state complex, represented by substantially better binding constants. Intriguingly, the crystal structure of the complex of acetyl-pepstatin with XMRV PR has shown a mixed type of binding: the unusual binding mode of two molecules of the inhibitor to the enzyme, in a mode very similar to the previously determined complex with pepstatin A, together with the classical binding mode found for HIV-1 PR. The structure is thus in good agreement with the very similar interaction energies calculated for the two types of binding. Database The final coordinates of the crystal structure of XMRV protease complexed with acetyl-pepstatin are available in the Protein Data Bank under the accession number 4EXH Structured digital abstract ? HIV-1 PR?and?HIV-1 PR?bind?by?biochemical?(View interaction) ? XMRV PR?cleaves?MLV Gag?by?enzymatic study?(View interaction) ? XMRV PR?and?XMRV PR?bind?by?biochemical?(View interaction) ? XMRV PR?and?XMRV PR?bind?by?x-ray crystallography?(View interaction).     

Binding mode prediction of biologically active compounds from plant Salvia Miltiorrhiza as integrase inhibitor

Integrase (IN), an essential enzyme for HIV-1 replication, has been targeted in antiretroviral drug therapy. The emergence of HIV-1 variants clinically resistant to antiretroviral agents has lead to the development of alternative IN inhibitors. In the present work, binding modes of a high potent IN inhibitor, M₅22 and M₅32, within the catalytic binding site of wild type (WT) IN were determined using molecular docking calculation. Both M₅22 and M₅32 displayed similar modes of binding within the IN putative binding pocket and exhibited favorable interactions with the catalytic Mg²⁺ ions, the nearby amino acids and viral DNA through metal-ligand chelation, hydrogen bonding and π-π stacking interactions. Furthermore, the modes of action of these two compounds against the mutated Y212R, N224H and S217H PFV IN were also predicted. Although the replacement of amino acid could somehow disturb inhibitor binding mode, almost key interactions which detected in the WT complexes were fairly conserved. Detailed information could highlight the application of M₅22 and M₅32 as candidate IN inhibitors for drug development against drug resistant strains.     

Solvation effects are responsible for the reduced inhibitor affinity of some HIV-1 PR mutants.

The formulation of HIV-1 PR inhibitors as anti-viral drugs has been hindered by the appearance of protease strains that present drug resistance to these compounds. The mechanism by which the HIV-1 PR mutants lower their affinity for the inhibitor is not yet fully understood. We have applied a modified Poisson-Boltzmann method to the evaluation of the molecular interactions that contribute to the lowering of the inhibitor affinity to some polar mutants at position 82. These strains present drug resistance behavior and hence are ideally suited for these studies. Our results indicate that the reduction in binding affinity is due to the solvation effects that penalize the binding to the more polar mutants. The inhibitor binding ranking of the different mutants can be explained from the analysis of the different components of our free energy scoring function.     

Molecular docking and 3D-QSAR studies of HIV-1 protease inhibitors

HIV-1 protease is an obligatory enzyme in the replication process of the HIV virus. The abundance of structural information on HIV-1PR has made the enzyme an attractive target for computer-aided drug design strategies. The daunting ability of the virus to rapidly generate resistant mutants suggests that there is an ongoing need for new HIV-1PR inhibitors with better efficacy profiles and reduced toxicity. In the present investigation, molecular modeling studies were performed on a series of 54 cyclic urea analogs with symmetric P2/P2′ substituents. The binding modes of these inhibitors were determined by docking. The docking results also provided a reliable conformational superimposition scheme for the 3D-QSAR studies. To gain insight into the steric, electrostatic, hydrophobic and hydrogen-bonding properties of these molecules and their influence on the inhibitory activity, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed. Two different alignment schemes viz. receptor-based and atom-fit alignment, were used in this study to build the QSAR models. The derived 3D-QSAR models were found to be robust with statistically significant r ² and r ² _pred values and have led to the identification of regions important for steric, hydrophobic and electronic interactions. The predictive ability of the models was assessed on a set of molecules that were not included in the training set. Superimposition of the 3D-contour maps generated from these models onto the active site of enzyme provided additional insight into the structural requirements of these inhibitors. The CoMFA and CoMSIA models were used to design some new inhibitors with improved binding affinity. Pharmacokinetic and toxicity predictions were also carried out for these molecules to gauge their ADME and safety profile. The computational results may open up new avenues for synthesis of potent HIV-1 protease inhibitors.     

Structural characterization of B and non-B subtypes of HIV-protease: insights into the natural susceptibility to drug resistance development

Although a majority of HIV-1 infections in Brazil are caused by the subtype B virus (also prevalent in the United States and Western Europe), viral subtypes F and C are also found very frequently. Genomic differences between the subtypes give rise to sequence variations in the encoded proteins, including the HIV-1 protease. The current anti-HIV drugs have been developed primarily against subtype B and the effects arising from the combination of drug-resistance mutations with the naturally existing polymorphisms in non-B HIV-1 subtypes are only beginning to be elucidated. To gain more insights into the structure and function of different variants of HIV proteases, we have determined a 2.1 A structure of the native subtype F HIV-1 protease (PR) in complex with the protease inhibitor TL-3. We have also solved crystal structures of two multi-drug resistant mutant HIV PRs in complex with TL-3, from subtype B (Bmut) carrying the primary mutations V82A and L90M, and from subtype F (Fmut) carrying the primary mutation V82A plus the secondary mutation M36I, at 1.75 A and 2.8 A resolution, respectively. The proteases Bmut, Fwt and Fmut exhibit sevenfold, threefold, and 54-fold resistance to TL-3, respectively. In addition, the structure of subtype B wild type HIV-PR in complex with TL-3 has been redetermined in space group P6(1), consistent with the other three structures. Our results show that the primary mutation V82A causes the known effect of collapsing the S1/S1' pockets that ultimately lead to the reduced inhibitory effect of TL-3. Our results further indicate that two naturally occurring polymorphic substitutions in subtype F and other non-B HIV proteases, M36I and L89M, may lead to early development of drug resistance in patients infected with non-B HIV subtypes.