Molecular Dynamics and Free Energy Studies on the Wild-Type and Mutated HIV-1 Protease Complexed with Four Approved Drugs: Mechanism of Binding and Drug Resistance

Functionalized carbon nanotubes as an alternative to traditional anti-HIV-1 protease inhibitors: An understanding towards Nano-medicine development through MD simulations

The Human Immunodeficiency Virus (HIV) has been the source of epidemic infection of AIDS for a longer period. One of the most difficult tasks is identifying novel medications that can help to decrease or control this global health hazard by overcoming drug resistance. In recent decades' nanoparticles are emerging as extremely relevant in drug delivery platforms. In the current study, the pristine (SWCNT) and hydroxyl functionalized (SWCNT-OH) versions of the SWCNT were investigated as inhibitors against the wild-type (WT) and three key mutants of HIV-1 protease (HIV-pr) (I50V, V82A, and I84V). Molecular docking of SWCNT in the catalytic domain and running all-atom MD simulations of all complexes are also part of this project. A thorough inspection of conformational dynamics from 50 ns trajectories reveals that both the pristine and SWCNT-OH can fit right to the pocket region of HIV-pr and govern flap dynamics. The binding affinity of the four HIV-pr-SWCNT/SWCNT-OH complexes was further investigated using MM-PBSA-dependent binding free energy studies. In most mutants and WT systems, SWCNT-OH was reported to bind proportionately many folds (kcal/mol) more than pristine SWCNTs. Hence, SWCNTs are possible HIV-pr inhibitors in terms of their stable existence in the pocket area, stronger binding to the protease, and regulation of flap dynamics in controlling the active site volume, which have vast potential for applications against drug resistance.

Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2

PRS17, a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PRS17 with inhibitor Darunavir (DRV), two complexes of WT PR and PRS17 with substrate analogue CA-p2, two unligand WT PR and PRS17), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PRS17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.

Drug Resistance Mechanism of M46I-Mutation-Induced Saquinavir Resistance in HIV-1 Protease Using Molecular Dynamics Simulation and Binding Energy Calculation

Drug-resistance-associated mutation in essential proteins of the viral life cycle is a major concern in anti-retroviral therapy. M46I, a non-active site mutation in HIV-1 protease has been clinically associated with saquinavir resistance in HIV patients. A 100 ns molecular dynamics (MD) simulation and MM-PBSA calculations were performed to study the molecular mechanism of M46I-mutation-based saquinavir resistance. In order to acquire deeper insight into the drug-resistance mechanism, the flap curling, closed/semi-open/open conformations, and active site compactness were studied. The M46I mutation significantly affects the energetics and conformational stability of HIV-1 protease in terms of RMSD, RMSF, Rg, SASA, and hydrogen formation potential. This mutation significantly decreased van der Waals interaction and binding free energy (∆G) in the M46I–saquinavir complex and induced inward flap curling and a wider opening of the flaps for most of the MD simulation period. The predominant open conformation was reduced, but inward flap curling/active site compactness was increased in the presence of saquinavir in M46I HIV-1 protease. In conclusion, the M46I mutation induced structural dynamics changes that weaken the protease grip on saquinavir without distorting the active site of the protein. The produced information may be utilized for the discovery of inhibitor(s) against drug-resistant HIV-1 protease.

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.

Non-active site mutants of HIV-1 protease influence resistance and sensitisation towards protease inhibitors

Background

HIV-1 can develop resistance to antiretroviral drugs, mainly through mutations within the target regions of the drugs. In HIV-1 protease, a majority of resistance-associated mutations that develop in response to therapy with protease inhibitors are found in the protease’s active site that serves also as a binding pocket for the protease inhibitors, thus directly impacting the protease-inhibitor interactions. Some resistance-associated mutations, however, are found in more distant regions, and the exact mechanisms how these mutations affect protease-inhibitor interactions are unclear. Furthermore, some of these mutations, e.g. N88S and L76V, do not only induce resistance to the currently administered drugs, but contrarily induce sensitivity towards other drugs. In this study, mutations N88S and L76V, along with three other resistance-associated mutations, M46I, I50L, and I84V, are analysed by means of molecular dynamics simulations to investigate their role in complexes of the protease with different inhibitors and in different background sequence contexts.

Results

Using these simulations for alchemical calculations to estimate the effects of mutations M46I, I50L, I84V, N88S, and L76V on binding free energies shows they are in general in line with the mutations’ effect on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$IC_{50}$$\end{document}IC50 values. For the primary mutation L76V, however, the presence of a background mutation M46I in our analysis influences whether the unfavourable effect of L76V on inhibitor binding is sufficient to outweigh the accompanying reduction in catalytic activity of the protease. Finally, we show that L76V and N88S changes the hydrogen bond stability of these residues with residues D30/K45 and D30/T31/T74, respectively.

Conclusions

We demonstrate that estimating the effect of both binding pocket and distant mutations on inhibitor binding free energy using alchemical calculations can reproduce their effect on the experimentally measured \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$IC_{50}$$\end{document}IC50 values. We show that distant site mutations L76V and N88S affect the hydrogen bond network in the protease’s active site, which offers an explanation for the indirect effect of these mutations on inhibitor binding. This work thus provides valuable insights on interplay between primary and background mutations and mechanisms how they affect inhibitor binding.

Drug Resistance Mutation L76V Alters Nonpolar Interactions at the Flap-Core Interface of HIV-1 Protease

Four HIV-1 protease (PR) inhibitors, clinical inhibitors lopinavir and tipranavir, and two investigational compounds 4 and 5, were studied for their effect on the structure and activity of PR with drug-resistant mutation L76V (PR_L76V). Compound 5 exhibited the best K _i value of 1.9 nM for PR_L76V, whereas the other three inhibitors had K _i values of 4.5-7.6 nM, 2-3 orders of magnitude worse than for wild-type enzymes. Crystal structures showed only minor differences in interactions of inhibitors with PR_L76V compared to wild-type complexes. The shorter side chain of Val76 in the mutant lost hydrophobic interactions with Lys45 and Ile47 in the flap, and with Asp30 and Thr74 in the protein core, consistent with decreased stability. Inhibitors forming additional polar interactions with the flaps or dimer interface of PR_L76V were unable to compensate for the decrease in internal hydrophobic contacts. These structures provide insights for inhibitor design.

Systematic molecular dynamics, MM-PBSA, and ab initio approaches to the saquinavir resistance mechanism in HIV-1 PR due to 11 double and multiple mutations

Mutations in the human immunodeficiency virus (HIV) enable virus replication even when appropriate antiretroviral therapy is followed, thus leading to the emergence of drug resistance. In a previous work, we systematically examined seven single mutations that are associated with saquinavir (SQV) resistance in HIV-1 protease (Tzoupis, H.; Leonis, G.; Mavromoustakos, T.; Papadopoulos, M. G. J. Chem. Theory Comput. 2013, 9, 1754-1764). Herein, we extend our analysis, which includes seven double (G48V-V82A, L10I-G48V, G48V-L90M, I84V-L90M, L10I-V82A, L10I-L63P, A71V-G73S) and four multiple (L10I-L63P-A71V, L10I-G48V-V82A, G73S-I84V-L90M, L10I-L63P-A71V-G73S-I84V-L90M) SQV-HIV-1 PR mutant complexes, in an attempt to generalize our findings and formulate the main elements of the SQV resistance mechanism in the protease. On the basis of molecular dynamics (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA), and ab initio computational approaches, we identified specific features that constitute the HIV-1 PR mechanism of resistance at the molecular level: the low flexibility of SQV in the binding cavity and the preservation of hydrogen bonding (HB) and van der Waals interactions between SQV and several active-site (Gly27/27', Asp29/29'/30/30', especially Asp25/25') and flap (Ile50/50', Gly48/48') residues of the protease contribute significantly to efficient binding. The total enthalpy loss in all mutants is mostly due to the loss in enthalpy of the active-site region. Furthermore, it was observed that mutation accumulation may induce stabilization to SQV and to the flaps through enhanced HB interactions that lead to improved inhibition (e.g., accumulation of mutations in complexes containing L10I, G48V, L63P, I84V, or L90M single mutations). It was also concluded that permanent flap closure is obtained independently of mutations and SQV binding is mostly driven by van der Waals, nonpolar, and exchange-energy contributions. Importantly, it was indicated that the optimal positioning of SQV and the structure of the binding cavity are tightly coupled, since small changes in geometry may affect the binding energy greatly. The results of our theoretical approaches are in agreement with experimental evidence and provide a reliable description of SQV resistance in HIV-1 PR.

Molecular modeling and residue interaction network studies on the mechanism of binding and resistance of the HCV NS5B polymerase mutants to VX-222 and ANA598

Hepatitis C virus (HCV) NS5B protein is an RNA-dependent RNA polymerase (RdRp) with essential functions in viral genome replication and represents a promising therapeutic target to develop direct-acting antivirals (DAAs). Multiple nonnucleoside inhibitors (NNIs) binding sites have been identified within the polymerase. VX-222 and ANA598 are two NNIs targeting thumb II site and palm I site of HCV NS5B polymerase, respectively. These two molecules have been shown to be very effective in phase II clinical trials. However, the emergence of resistant HCV replicon variants (L419M, M423T, I482L mutants to VX-222 and M414T, M414L, G554D mutants to ANA598) has significantly decreased their efficacy. To elucidate the molecular mechanism about how these mutations influenced the drug binding mode and decreased drug efficacy, we studied the binding modes of VX-222 and ANA598 to wild-type and mutant polymerase by molecular modeling approach. Molecular dynamics (MD) simulations results combined with binding free energy calculations indicated that the mutations significantly altered the binding free energy and the interaction for the drugs to polymerase. The further per-residue binding free energy decomposition analysis revealed that the mutations decreased the interactions with several key residues, such as L419, M423, L474, S476, I482, L497, for VX-222 and L384, N411, M414, Y415, Q446, S556, G557 for ANA598. These were the major origins for the resistance to these two drugs. In addition, by analyzing the residue interaction network (RIN) of the complexes between the drugs with wild-type and the mutant polymerase, we found that the mutation residues in the networks involved in the drug resistance possessed a relatively lower size of topology centralities. The shift of betweenness and closeness values of binding site residues in the mutant polymerase is relevant to the mechanism of drug resistance of VX-222 and ANA598. These results can provide an atomic-level understanding about the mechanisms of drug resistance conferred by the studied mutations and will be helpful to design more potent inhibitors which could effectively overcome drug resistance of antivirus agents.

Computational study on the drug resistance mechanism against HCV NS3/4A protease inhibitors vaniprevir and MK-5172 by the combination use of molecular dynamics simulation, residue interaction network, and substrate envelope analysis

Hepatitis C virus (HCV) NS3/4A protease is an important and attractive target for anti-HCV drug development and discovery. Vaniprevir (phase III clinical trials) and MK-5172 (phase II clinical trials) are two potent antiviral compounds that target NS3/4A protease. However, the emergence of resistance to these two inhibitors reduced the effectiveness of vaniprevir and MK-5172 against viral replication. Among the drug resistance mutations, three single-site mutations at residues Arg155, Ala156, and Asp168 in NS3/4A protease are especially important due to their resistance to nearly all inhibitors in clinical development. A detailed understanding of drug resistance mechanism to vaniprevir and MK-5172 is therefore very crucial for the design of novel potent agents targeting viral variants. In this work, molecular dynamics (MD) simulation, binding free energy calculation, free energy decomposition, residue interaction network (RIN), and substrate envelope analysis were used to study the detailed drug resistance mechanism of the three mutants R155K, A156T, and D168A to vaniprevir and MK-5172. MD simulation was used to investigate the binding mode for these two inhibitors to wild-type and resistant mutants of HCV NS3/4A protease. Binding free energy calculation and free energy decomposition analysis reveal that drug resistance mutations reduced the interactions between the active site residues and substituent in the P2 to P4 linker of vaniprevir and MK-5172. Furthermore, RIN and substrate envelope analysis indicate that the studied mutations of the residues are located outside the substrate (4B5A) binding site and selectively decrease the affinity of inhibitors but not the activity of the enzyme and consequently help NS3/4A protease escape from the effect of the inhibitors without influencing the affinity of substrate binding. These findings can provide useful information for understanding the drug resistance mechanism against vaniprevir and MK-5172. The results can also provide some potential clues for further design of novel inhibitors that are less susceptible to drug resistance.

State-of-the-art and dissemination of computational tools for drug-design purposes: a survey among Italian academics and industrial institutions

During the first edition of the Computationally Driven Drug Discovery meeting, held in November 2011 at Dompé Pharma (L'Aquila, Italy), a questionnaire regarding the diffusion and the use of computational tools for drug-design purposes in both academia and industry was distributed among all participants. This is a follow-up of a previously reported investigation carried out among a few companies in 2007. The new questionnaire implemented five sections dedicated to: research group identification and classification; 18 different computational techniques; software information; hardware data; and economical business considerations. In this article, together with a detailed history of the different computational methods, a statistical analysis of the survey results that enabled the identification of the prevalent computational techniques adopted in drug-design projects is reported and a profile of the computational medicinal chemist currently working in academia and pharmaceutical companies in Italy is highlighted.

Molecular dynamics investigation on a series of HIV protease inhibitors: assessing the performance of MM-PBSA and MM-GBSA approaches

The binding free energies (ΔG(Bind)) obtained from molecular mechanics with Poisson-Boltzmann surface area (MM-PBSA) or molecular mechanics with Generalized Born surface area (MM-GBSA) calculations using molecular dynamics (MD) trajectories are the most popular procedures to measure the strength of interactions between a ligand and its receptor. Several attempts have been made to correlate the ΔG(Bind) and experimental IC(50) values in order to observe the relationship between binding strength of a ligand (with its receptor) and its inhibitory activity. The duration of MD simulations seems very important for getting acceptable correlation. Here, we are presenting a systematic study to estimate the reasonable MD simulation time for acceptable correlation between ΔG(Bind) and experimental IC(50) values. A comparison between MM-PBSA and MM-GBSA approaches is also presented at various time scales. MD simulations (10 ns) for 14 HIV protease inhibitors have been carried out by using the Amber program. MM-PBSA/GBSA based ΔG(Bind) have been calculated and correlated with experimental IC(50) values at different time scales (0-1 to 0-10 ns). This study clearly demonstrates that the MM-PBSA based ΔG(Bind) (ΔG(Bind)-PB) values provide very good correlation with experimental IC(50) values (quantitative and qualitative) when MD simulation is carried out for a longer time; however, MM-GBSA based ΔG(Bind) (ΔG(Bind)-GB) values show acceptable correlation for shorter time of simulation also. The accuracy of ΔG(Bind)-PB increases and ΔG(Bind)-GB remains almost constant with the increasing time of simulation.

Binding of single walled carbon nanotube to WT and mutant HIV-1 proteases: analysis of flap dynamics and binding mechanism

Most of the currently treated HIV-1 protease (HIV-PR) inhibitors have been prone to suffer from the mutations associated drug resistance. Therefore, it is necessary to search for potent alternatives against the drug resistance. In the current study we have tested the single-walled carbon nanotube (SWCNT) as an inhibitor in wild type (WT) as well as in three primary mutants (I50V(PR), V82A(PR) and I84V(PR)) of the HIV-1-PR through docking the SWCNT in the active site region, and then performed all-atom MD simulations for the complexes. The conformational dynamics of HIV-PR with a 20 ns trajectory reveals that the SWCNT can effectively bind to the HIV-1-PR active site and regulate the flap dynamics such as maintaining the flap-flap closed. To gain an insight into the binding affinity, we also performed the MM-PBSA based binding free energy calculations for the four HIV-PR/SWCNT complexes. It was observed that, although the binding between the SWCNT and the HIV-PR decreases due to the mutations, the SWCNTs bind to the HIV-PRs 3-5 folds stronger than the most potent HIV-1-PR inhibitor, TMC114. Remarkably, the significant interactions with binding energy higher than 1kcal/mol focus on the flap and active regions, which favors closing flap-flap and deactivating the active residues of the HIV-PR. The flap dynamics and binding strength information for HIV-PR and SWCNTs can help design SWCNT-based HIV-1-PR inhibitors.

Tackling the challenges posed by target flexibility in drug design

IMPORTANCE OF THE FIELD

Current computational docking methods are often effective in predicting accurate drug-binding geometries in cases of relatively rigid target structures. However, binding of drug-like ligands to protein receptor molecules frequently involves or even requires conformational adaptation. Realistic prediction of ligand-receptor binding geometries and complex stability needs in many cases an appropriate inclusion of conformational changes, not only for the ligand, but also for the receptor molecule.

AREAS COVERED IN THIS REVIEW

Recent approaches to efficiently account for target receptor flexibility during docking simulations are reviewed.

WHAT THE READER WILL GAIN

The reader will gain insights into methods to efficiently treat protein side-chain flexibility and approaches for continuous adaptation of backbone conformations in pre-calculated essential or soft collective degrees of freedom. In addition, molecular dynamics or Monte Carlo based methods providing simultaneous inclusion of receptor and ligand flexibility are discussed as well as promising new developments to generate conformationally diverse ensembles of a protein structure. The large variety of possible conformational changes in proteins on ligand binding is illustrated for the enzyme reverse transcriptase of HIV-1, which is an important drug target.

TAKE HOME MESSAGE

If the backbone conformation of the binding site does not change, current docking programs can perform well by taking side-chain reorientations into account only. Future progress to account for full target flexibility in docking requires both accurate prediction of the essential modes of backbone motion and improvements in scoring to enhance selectivity. Thus, the scoring function should realistically cover energetic and particularly entropic contributions to binding, which would allow more realistic estimates of binding free energies.

From base pair to bedside: molecular simulation and the translation of genomics to personalized medicine

Despite the promises made that genomic sequencing would transform therapy by introducing a new era of personalized medicine, relatively few tangible breakthroughs have been made. This has led to the recognition that complex interactions at multiple spatial, temporal, and organizational levels may often combine to produce disease. Understanding this complexity requires that existing and future models are used and interpreted within a framework that incorporates knowledge derived from investigations at multiple levels of biological function. It also requires a computational infrastructure capable of dealing with the vast quantities of data generated by genomic approaches. In this review, we discuss the use of molecular modeling to generate quantitative and qualitative insights at the smallest scales of the systems biology hierarchy, how it can play an important role in the development of a systems understanding of disease and in the application of such knowledge to help discover new therapies and target existing ones on a personal level.

The L76V mutation in HIV-1 protease is potentially associated with hypersusceptibility to protease inhibitors Atazanavir and Saquinavir: is there a clinical advantage?

Background

Although being considered as a rarely observed HIV-1 protease mutation in clinical isolates, the L76V-prevalence increased 1998-2008 in some European countries most likely due to the approval of Lopinavir, Amprenavir and Darunavir which can select L76V. Beside an enhancement of resistance, L76V is also discussed to confer hypersusceptibility to the drugs Atazanavir and Saquinavir which might enable new treatment strategies by trying to take advantage of particular mutations.

Results

Based on a cohort of 47 L76V-positive patients, we examined if there might exist a clinical advantage for L76V-positive patients concerning long-term success of PI-containing regimens in patients with limited therapy options.

Genotypic- and phenotypic HIV-resistance tests from 47 mostly multi-resistant, L76V-positive patients throughout Germany were accomplished retrospectively 1999-2009. Five genotype-based drug-susceptibility predictions received from online interpretation-tools for Atazanavir, Saquinavir, Amprenavir and Lopinavir, were compared to phenotype-based predictions that were determined by using a recombinant virus assay along with a Virtual Phenotype™(Virco). The clinical outcome of the L76V-adapted follow-up therapy was determined by monitoring viral load for 96 weeks.

Conclusions

In this analysis, the mostly used interpretation systems overestimated the L76V-mutation concerning Atazanavir- and SQV resistance. In fact, a clear benefit in drug susceptibility for these drugs was observed in phenotype analysis after establishment of L76V. More importantly, long-term therapy success was significantly higher in patients receiving Atazanavir and/or Saquinavir plus one L76V-selecting drug compared to patients without L76V-selecting agents (p = 0.002).

In case of L76V-occurrence ATV and/or SQV may represent encouraging options for patients in deep salvage situations.

A comparative study of HIV-1 and HTLV-I protease structure and dynamics reveals a conserved residue interaction network

The two retroviruses human T-lymphotropic virus type I (HTLV-I) and human immunodeficiency virus type 1 (HIV-1) are the causative agents of severe and fatal diseases including adult T-cell leukemia and the acquired immune deficiency syndrome (AIDS). Both viruses code for a protease that is essential for replication and therefore represents a key target for drugs interfering with viral infection. The retroviral proteases from HIV-1 and HTLV-I share 31% sequence identity and high structural similarities. Yet, their substrate specificities and inhibition profiles differ substantially. In this study, we performed all-atom molecular dynamics (MD) simulations for both enzymes in their ligand-free states and in complex with model substrates in order to compare their dynamic behaviors and enhance our understanding of the correlation between sequence, structure, and dynamics in this protein family. We found extensive similarities in both local and overall protein dynamics, as well as in the energetics of their interactions with model substrates. Interestingly, those residues that are important for strong ligand binding are frequently not conserved in sequence, thereby offering an explanation for the differences in binding specificity. Moreover, we identified an interaction network of contacts between conserved residues that interconnects secondary structure elements and serves as a scaffold for the protein fold. This interaction network is conformationally stable over time and may provide an explanation for the highly similar dynamic behavior of the two retroviral proteases, even in the light of their rather low overall sequence identity.