Analysis of the Molecular Structure of TMC114 HIV-1 Protease Inhibitor and the HIV-1 Protease Protein

Molecular Modeling: Jterm 08


John's BLOG

Background

HIV-1 Protease is a protein produced by the HIV virus. It cleaves viral polypeptides into functional proteins that are needed for viral assembly and activity. HIV-1 Protease is a common target in treating AIDS, because the inhibition of this protein would prevent many essential viral proteins from being produced. Therefore the inhibition of protease could significantly slow the progress of the disease. Many drugs have been made to treat AIDS, but many of them have become ineffective. This is because the HIV virus mutates very easily so that the active site or surface receptor the drug targets are no longer recognizable to the drug. If the drug cannot recognize the target protein it can no longer attach to and inhibit the protein.

Most protease inhibitors bind to the active site of the protease protein so that the active site of the protease cannot bind to the substrate (the proteins to be cleaved to produce active viral proteins). When the molecular structure of the active site mutates, the inhibitors no longer can bind to the active site. This is a problem when combatting AIDS because drugs that are produced eventually become ineffective.

The recently released (June 23, 2006) TMC114 protease inhibitor (also known as Darunavir) is unique because unlike many other inhibitors it has conformational flexibility. This means that the molecule is not stiff, the single bonds can rotate so that the molecule has many different conformations. Since the TMC114 inhibitor has so many conformations, it can still bind to the active site of the protease when it mutates.

In this project, I will be investigating the flexibility of this inhibitor. I will try to locate the important areas of the molecule that allow various conformations that can bind to protease mutants. I will also analyze the active site of the protease and determine the amino acid sequences where hydrogen bonding occurs with the inhibitor.



Hydrogen bonding and the HIV-1 Protease active site


This image shows the acidic and basic areas of the HIV-1 Protease. The TMC114inhibitor will form hydrogen bonds with acidic and basic regions when it docks onto the active site. The inhibitor is embedded within the molecule in the image.

This image shows how the inhibitor binds by hydrogen bonding in the active site of the protease.

HIV-1 Protease has many hydrogen bonding sites within the molecule that allow for direct bonding with the active site of the protease. The hydrogen bonding sites can also undergo intramolecular hydrogen bonding, allowing for more conformations of the inhibitor. The intramolecular bonding and flexibility of TMC114 allows for many conformations that can bind to mutated protease active sites. More information on intramolecular bonding will be presented in the next section. This figure highlights tha acidic and basic regions of TMC114 that have the ability to create hydrogen bonds.


The TMC114 inhibitor has a bis-tetrahydrofurane (bis-THF) group that makes it unique from similar inhibitors, such as Amprenavir. Both oxygens in this group hydrogen bond with amino acids ASP 29 and ASP 30 in the protease active site (ASP is the amino acid aspartic acid). This particular pair of hydrogen bonds is apparently essential to the unique stability the TMC114-protease complex has (more about the energetic stablity of TMC114 will be explained in later sections). It is important for an inhibitor to be stable when bound to the protein because it may otherwise fall off, and fail to prevent the protease from acting on its substrate. If this happens the protease will continue to aid in the construction of more viruses.



Intramolecular Hydrogen Bonding and its Impact on Conformations of TMC114

The figure below shows the five types (b-e) of intramolecular hydrogen bonding that occurs in TMC114 to create approximately 43 energetically stable conformations. Theoretically there are 260 conformations, but many of these are energetically unfavored conformations. Figure "a" shows the TMC114 without any intramolecular hydrogen bonding. The following energies (kcal/mol) are relative to the most stable conformation (letter f in the figure) of the 43 conformations tested by Kanda Nivesanond and others: a)7.05 b)0.27 c)3.54 d)2.19 e)6.12

The range of relative energies to most stable conformation is 0 to 8.32 kcal/mol. It is important to note that the conformation with no hydrogen bonding was much higher in energy, and was much less stable then those with intramolecular hydrogen bonding. The intramolecular hydrogen bonding obviously plays an important role in stabilizing the different conformations of the molecule.

The figure above also shows the regions on the molecule that undergo intramolecular bonding. Over half of the conformations have hydrogen bonding at the -OH and O=S site (H1 in the figure). This intramolecular hydrogen bond is probably the most common because the bonding regions are very close to each other and are in a flexible region (area with numerous single bonds) of the molecule.

Intramolecular hydrogen bonding and flexibility give TMC114 the ability to have several stable conformations. As a result TMC114 has the ability to adapt to the mutating active site of the protease, therefore making it an effective drug to combat the rapidly evolving HIV virus.


The Energetic Stability of Various TMC114 Conformations when Bound to the HIV-1 Protease

The stability of TMC114 conformations is important, but this means nothing if the complex it forms with protease is not stable. The table below shows the binding energies Nivesanond and others obtained of various TMC114 conformations bound to wild and mutant protease. It is important to note that sometimes stable conformations do not have low binding energies, such as conf17 with the wild protease. Particular conformations have lower energies depending on if the molecule is bound to the wild or mutant protease.

Relative Conformational Energies (kcal/mol) and Binding Energies (Wild and Mutant)
Relative energy Binding energy with wild type HIV-1 Protease Binding energy with mutant type HIV-1 Protease
confwt 7.05 -13.99 -13.29
conf21 3.08 -5.27 -7.64
conf20 2.78 -2.98 -7.21
conf04 1.92 -0.62 -6.09
conf17 0.00 +6.49 -8.06

These energy values show that particular conformations form very stable complexes with wild and mutant proteases. These values demonstrate how TMC114 can be an effective inhibitor even to a mutated protease; conformations of TMC114 have significantly negative energy values for both wild and mutant types.

Comparison of TMC114 and Amprenavir

TMC114 and Amprenavir are very similar HIV-1 protease inhibitors. They only differ in that TMC114 has two THF groups on the the molecule, and Amprenavir has only one THF group. This small difference, however, has an important effect whenthe inhibitor binds to the active site of the protease. Since TMC114 has two THF groups, it has more places for hydrogen bonding to occur. Specifically, the two THF groups interact with the amino groups of the amino acids ASP 29 and ASP 30. Kanda Nivesanond and others emphasize that this particular hydrogen bonding interaction is essential to TMC114's stability when it is bound to the protease. The extra THF group on TMC114 ultimately makes the protease-TMC114 complex more stable than a protease-Amprenavir complex.


Energy Calculations for TMC114 and Amprenavir

The following is a table of energy calculations I conducted through WebMO for TMC114 and Amprenavir.

AM1 Energy Calculations(kcal/mol) using WebMO
Energy (kcal/mol)
TMC114 (rcsb) -255.0295
TMC114 (binding database) -251.4993
TMC114 (manually constructed) -255.7708
Amprenavir (rcsb) -214.0717
Amprenavir (binding database) -212.8114

Both TMC114 and Amprenavir are stable molecules by themselves since their energies are very negative. TMC114 is more stable than Amprenavir, an inhibitor very similar to TMC114. The TMC114 molecule I constructed manually on WebMO seems to have an accurate energy when compared the other molecules downloaded from databases.

Conclusions

The conformational flexibility and the strong hydrogen bonds of the THF groups allows TMC114 to form a very stable complex with HIV-1 protease. These properties make TMC114 more stable than other inhibitors, and able to bind to mutant forms of protease.

Things I would have liked to have done...

  • Use AMBER to get my own calculations on conformational energies.
  • Use the RMSD tool to compare the different conformations of TMC114, and to compare TMC114 to Amprenavir.
  • Make some kind of movie, because that would have been really cool.
  • Have a more detailed comparison of the binding energies of TMC114 and Amprenavir
  • References

    Nivesanond K.; Peeters A.;Lamoen D.; Van Alsenoy C. Conformational Analysis of TMC114, a Novel HIV-1 Protease Inhibitor. J. Chem. Inf. Model. 2008, 48, 99-108.

    King N. and others; Structural and Thermodynamic Basis for the Binding of TMC114, a Next-Generation Human Immunodeficiency Virus Type 1 Protease Inhibitor. J Virol. 2004 November; 78(21): 12012.12021.