The Centers for Disease Control and Prevention (CDC) in Healthy People 2020 addressed 10-year national objectives to improve the health of Americans. Of the 600 objectives that cover 42 topic areas, 18 of the objectives focus on HIV and the acquired immunodeficiency syndrome (AIDS). The primary goal of Healthy People 2020 as it relates to HIV/AIDS is to prevent HIV infection and its related illness and death  . As of 2015, there are approximately 39,515 new HIV infections each year in the United States  . Thus, interventions are needed to prevent transmission to those who are uninfected, and optimized treatment is necessary for those who are already infected. The development of new drugs that can prevent transmission as well as disease progression is of paramount public health importance, and targeting a critical viral enzyme such as HIV-1 integrase is a proven approach to meet this need.
Integrases  are enzymes that are packaged into retroviral particles, such as those of the human immunodeficiency virus (HIV-1) along with reverse transcriptase (RT) and the viral ribonucleic acid (RNA) genome. Integrase enzymes catalyze the insertion of the DNA copy of the viral genome synthesized by the RT into a host chromosomal locus to produce a lifelong infection of the host. The clinical development of existing HIV-1 integrase inhibitors or drugs has demonstrated the power of targeting integrase activity as a therapeutic approach     .
A simplified model of a drug binding to the active site of a large molecule such as a protein is shown in Figure 1. When a drug binds to the active site, the estimated binding energy can be computationally calculated from the energy considerations of the different interactions between the drug and the protein.
2. HIV-1 Life Cycle
HIV enters the cell and reverse transcriptase―one of three enzymes involved in the life cycle of HIV―copies single-stranded RNA (HIVRNA) into double stranded viral DNA (DS-cDNA). Integrase, the second of three enzymes, integrates the viral DNA into the host cell’s DNA (HostDNA). The last enzyme, protease, is responsible for cleavage of polyproteins into the structural proteins of HIV such as the matrix, capsid, and nucleocapsid proteins. The mature HIV-1 virion exits the cell. (Figure 2) 
Figure 1. A simple representation of binding.
Figure 2. HIV life cycle. https://cellandbioscience.biomedcentral.com/articles/10.1186/2045-3701-2-32.
3. HIV-1 Integrase Structure and Function
HIV-1 integrase consists of three domains or sections: amino-terminal domain, amino acids 1 - 50 (NTD), core catalytic domain, 50 - 212 (CCD), and the carboxyl-terminal, 212 - 288 (CTD)  .
The catalytic core domain of 1BIS.pdb  was the focus of the research (Figure 3). It contains the active site of the protein that is defined by the amino acids: aspartic acid (D) 64, aspartic acid (D) 116, and glutamic acid (E) 152―the DDE triad  . Using ICM-Pro  , computational software, several pockets were identified but one of the pockets included the active site amino acids, and therefore a binding site for potential drugs. The figure was generated using ICM-Pro, Molsoft, L.L.C.
HIV-1 integrase catalyzes two biologically relevant nicking reactions (Figure 4): The first reaction (1) (referred to as 3' processing) nicks and trims the 3' end of each strand of unintegrated viral DNA at a specific site. The second reaction (2) (referred to as strand transfer) (ST) inserts the trimmed ends of viral DNA into the two strands of cellular DNA  .
4. Drug-Like and Nondrug-Like Molecules/Inhibitors
Many inhibitors have drug-like properties that follow Lipinski’s “Rule of Five” that was introduced for predicting oral bioavailability  .
These rules include the following:
1) Molecular mass ≤500 g/mole,
Figure 3. The catalytic core domain of 1BIS.pdb.
Figure 4. The two nicking reactions of HIV-1 integrase (figure provided by Dr. Michael Katzman, Penn State Hershey Medical Center).
2) Log P (the partition coefficient, defined as the ratio of the concentration of a drug in octanol to its concentration in water) between ±5,
3) Hydrogen bond donors ≤5, and
4) Hydrogen bond acceptors ≤10.
Some integrase inhibitors, however, do not follow these rules. For this research, if an inhibitor failed to follow even one of these rules it was classified as nondrug-like. Using data from ICM-Pro, the inhibitors were divided into two groups―drug-like and nondrug-like.
Figure 5 and Figure 6 (generated using ICM-Pro, Molsoft L.L.C.) represent an example of a drug-like  and nondrug-like molecule  , respectively. The molecules were docked into the active site of 1BIS.pdb  and are shown in the potential binding site pocket identified by ICM-Pro. The parameters for Lipinski’s
Figure 5. Zhao’14, 5o  .
Figure 6. Zhao’14, 5sS  .
“Rule of Five”  from ICM-Pro are shown for each molecule. The nondrug-like molecule, 5sS has an unacceptable number of hydrogen bond donors (HBD), that is, six.
5. Estimated Binding Energy
For the research, the goal was to computationally determine and compare the estimated binding energy (EBE) of a group of drug-like and nondrug-like molecules using 1BIS.pdb and two mutant models of 1BIS.pdb  , Y143R and N155H. An estimated binding energy can be determined using ICM-Pro (Molsoft, L.L.C.) The binding energy function of ICM-Pro takes into account the following terms  :
1) Weighted sum of ligand-target van der Waals interactions and the internal force field energy of the ligand,
2) Free energy changes due to conformational energy loss upon ligand binding,
3) Hydrogen bonding interactions,
4) Hydrogen bond donor-acceptor desolvation energy,
5) Solvation electrostatic energy upon ligand binding,
6) Hydrophobic free energy gain, and
7) Size correction term proportional to the number of ligand atoms.
The more negative the ICM-Pro score, the higher the probability the ligand is a binder and could act as a potential drug to block the action of a protein such as HIV-1 integrase. Therefore, the viral DNA could be prevented from integration into the host’s DNA.
6. The Molecules
From previous research on approximately 200 molecules, initially considered to be potential inhibitors of HIV-1 integrase, a decision was made to address the EBE of a small group of molecules. Molecules with an experimental strand transfer concentration ≤0.01 µM were chosen because strand transfer is the reaction where HIV-1 integrase inserts the trimmed ends of viral DNA into the two strands of cellular DNA  . The smaller the concentration of a drug that is needed to inhibit the integration of viral DNA into host DNA, that is, strand transfer, the fewer the side effects of the drug. Only 30 molecules remained from the 200 considered initially, twenty-four drug-like and six nondrug-like molecules. Figures 7-9 show the drug-like molecules, and Figure 10, the nondrug- like molecules.
Each molecule was docked a minimum of three times within an 8-Å radius of the potential binding site pocket. A 2D ligand interaction of molecule 5sS  is displayed in Figure 11. It represents the interaction of the molecule with the surrounding amino acids of 1BIS.pdb within an 8-Å radius.
Figure 7. Drug-like molecules with an experimental ST ≤ 0.01 µM; Barreca’05  .
Figure 8. Drug-like molecules with an experimental ST ≤ 0.01 µM; Zhao’14  .
An average EBE was recorded. However, after docking and before computational determination of the EBE, the best conformer provided by ICM-Pro was chosen based on using crystal structures from the Cambridge Structural Database  . The exact molecules were not found in the Cambridge Structural Database, but from communication with a computational chemist at Bristol Myers Squibb  , it was suggested that a combination of two to three molecules might be needed to determine the best conformer. For the majority of molecules the first conformer provided by ICM-Pro was the one used to determine the EBE.
Figure 9. Drug-like molecules with an experimental ST ≤ 0.01 µM; Zhuang’03  .
Figure 10. Nondrug-like molecules with an experimental ST ≤ 0.01 µM; Zhao’14 
8. Results of the Estimated Binding Energies with 1BIS.pdb
8.1. Drug-Like Molecules
The estimated binding energy of the drug-like molecules    is shown in Figure 12.
Figure 11. 2D ligand interaction of molecule 5sS  within an 8-Å radius of the active site. The figure was generated by Andy Orry, Molsoft, L.L.C.
Figure 12. The estimated binding energies of the drug-like molecules with 1BIS.pdb. The numbers above the bars correspond to those in Figures 7-9.
8.1.1. Barreca’05  (Figure 7)
Molecules 3 and 21 have a similar framework: aromatic ring, linker, aromatic ring, linker, and two aromatic fused rings  . The only difference is the location of a nitrogen in the ring. However, their EBE values are −8.5 and −6.2, respectively. The placement of the nitrogen may account for the difference in the EBE. Molecule 18 is similar to molecule 3 but instead of two aromatic fused rings there is a sidechain  . They have the same EBE. Molecules 36, 37, 40, 41, 43, 44, and 45 have a similar framework: aromatic ring, linker, indole group, and a sidechain. Molecules 36 and 37 only differ in the presence of a fluorine on the first ring of molecule 37; they have a similar EBE. Molecules 40, 41, and 43 differ in the absence (molecule 40) or presence of fluorine and/or chlorine. They all have a similar EBE. Molecules 44 and 45 have a methoxy group as a sidechain off the benzene ring of the indole group but as in the previous three molecules, the presence of fluorine on molecule 45 does not appear to affect the EBE since they have the same EBE.
8.1.2. Zhao’14  (Figure 8)
All of the Zhao drug-like molecules have a similar framework of two aromatic fused rings, an amide linker, aromatic ring, and a sidechain off the second ring that is different for each of the molecules. There appears to be stronger binding in 5l, 5n, 5o, and 5p; each of these molecules has a larger sidechain off the second ring.
8.1.3. Zhuang’03  (Figure 9)
Molecules 4, 6, and 7 have a similar framework composed of an aromatic ring, a methyl group linker, aromatic ring, a carbonyl group linker, and two aromatic fused rings. Molecule 7 has a large sidechain off the second ring, that is, a methyl group linked to a sulfonamide group that is part of a ring. This group might contribute to its more negative EBE. Molecule 1 had two aromatic linked rings and a large sidechain off the second ring with three carbonyl groups and an hydroxyl group. Its EBE is similar to molecules 4 and 6.
8.2. Nondrug-Like Molecules
The estimated binding energy of the nondrug-like molecules  is shown in Figure 13.
Zhao’14  (Figure 10)
The framework of the Zhao nondrug-like molecules (Figure 10) is two aromatic fused rings with an hydroxyl group off the second ring and a sidechain para to the hydroxyl group, an amide linker, and an aromatic ring with two fluorines
Figure 13. The estimated binding energies of the nondrug-like molecules with 1BIS.pdb. The numbers above the bars correspond to those in Figure 10.
ortho to each other. It is the sidechain, and an R and S configuration of the three pairs of molecules that make them different. The 5rR and 5rS pair had the simplest sidechain and a similar EBE, approximately −10.3 kcal. The other two pairs had a larger sidechain off the second ring than the 5r pair, and a more negative EBE for the R configuration.
9. Results of the Estimated Binding Energies with 1BIS.pdb, Mutation Y143R
As referenced in the HIV Drug Resistance Database  , mutations of HIV-1 integrase in humans have led to increased resistance to HIV-1 integrase drugs such as raltegravir and elvitegravir. One of the mutations is Y143R. Tyrosine (Y) is also one of the amino acids within an 8Å radius of the pocket identified as a potential binding site.
Structurally, tyrosine (Y) has a benzene ring. There are six possible hydrogen bond (HB) interactions for tyrosine and also π − π stacking because of the benzene ring. Arginine (R) can participate in 12 HB interactions but without the phenyl group it loses π − π stacking capability. (Figure 14)
Tyrosine (Y) (Figure 15) and arginine (R) (Figure 16) are shown within the 8Å radius of the active site pocket that includes the DDE triad of HIV-1 integrase. The molecules were generated using ICM-Pro, Molsoft, L.L.C.
Figure 14. Tyrosine (Y) and arginine (R).
Figure 15. The DDE triad and tyrosine in the active site pocket of 1BIS.pdb.
Figure 16. The DDE triad and arginine in the active site pocket of 1BIS.pdb.
9.1. Drug-Like Molecules
Approximately 58% of the drug-like molecules had an EBE that was more negative for the mutant model (shown in black) compared to 1BIS.pdb (shown in green) with a ≥1 kcal/mole difference for 33% of the molecules (Figure 17). There were also 33% of the drug-like molecules with 1BIS.pdb in which the EBE was more negative than for the mutant model.
Figure 17. The difference in the EBE between 1BIS.pdb (shown in green) and the mutated protein Y143R (shown in black) with the drug-like molecules. For those where there is a ≥ 1 kcal/mole difference they are marked with l.
9.2. Nondrug-Like Molecules
For the nondrug-like molecules, there were four molecules with tighter binding with the mutant model―the 5r and 5t molecules―R and S configuration (Figure 18). However, for 5sR and 5sS there was tighter binding wth 1BIS.pdb.
Figure 18. The difference in the EBE between 1BIS.pdb (shown in red) and the mutant model Y143R (shown in black), for the nondrug-like molecules. For those where there is a ≥ 1 kcal/mole difference they are marked with l.
10. Results of the Estimated Binding Energies with 1BIS.pdb, Mutation N155H
This mutation is also referenced in the HIV Drug Resistance Database  and has led to increased resistance of the HIV-1 integrase drugs raltegravir and elvitegravir. Asparagine (N) is also one of the amino acids within an 8Å radius of the pocket identified as a potential binding site.
Structurally, the R group of asparagine (N) has an amide group. There are 11 possible HB interactions for asparagine. The R group of histidine (H) is a polar basic group, an imidazole ring; there are nine possible HB interactions for histidine plus π − π stacking with the imidazole ring (Figure 19).
10.1. Drug-Like Molecules
Approximately 46% of the drug-like molecules had an EBE that was more negative with the mutant model than with 1BIS.pdb (Figure 22). There was a ≥1
Figure 19. Asparagine and histidine.
Figure 20. The DDE triad and asparagine in the active site pocket of 1BIS.pdb.
Figure 21. The DDE triad and histidine in the active site pocket of 1BIS.pdb.
Figure 22. The difference in the EBE between 1BIS.pdb (shown in green) and the mutant model N155H (shown in black) for the drug-like molecules. For those where there is a ≥ 1 kcal/mole difference they are marked with l.
kcal/mole difference for approximately 21% of the molecules. There were also approximately 46% of the drug-like molecules with 1BIS.pdb with a more negative EBE than the mutant model; 29% had a ≥1 kcal/mole difference (Figure 22).
10.2. Nondrug-Like Molecules
In four of the six molecules 1BIS.pdb had a more negative EBE: 5rS, 5sR, 5sS, and 5tR. Two molecules had a more negative EBE with the mutant model: 5rR and 5tS (very slight). There appears to be no structural pattern to tighter binding in the nondrug-like group. In one molecule, 5sR, there was a ≥1 kcal/mole difference with 1BIS.pdb (Figure 23).
Figure 23. The difference in the EBE using 1BIS.pdb (shown in red) and the mutant model N155H (shown in black) for the nondrug-like molecules.
Though drug-like characteristics as defined by Lipinski’s “Rule of Five” were introduced for predicting oral bioavailability, the rules were used to classify 30 potential inhibitors of 1BIS.pdb into two groups. The goal was to examine the estimated binding energy of the two groups with 1BIS.pdb and two mutant models of the protein.
The nondrug-like molecules with 1BIS.pdb generally had a more negative EBE (−12.9 to −10.2), that is, tighter binding than the drug-like molecules. The binding energy function of ICM-Pro addresses van der Waals interactions of the protein/solvent/ligand, free energy changes when the drug/protein complex forms, hydrogen bonding interactions, hydrophobic free energy gain, desolvation energy of hydrogen bond donor-acceptor interactions, and the number of atoms of the drug. The nondrug-like molecules generally have a larger molecular mass, a greater number of hydrogen bond interactions, and a logP value that makes them more hydrophobic, factors that will affect their binding to the protein.
There were seven drug-like molecules with EBE ≥−10 kcal/mole: Barreca 36 and 37, Zhao 5l, 5n, 5o, and 5p, and Zhuang 7. Barreca 36 and 37 were similar in structure except Barreca 37 had a fluorine on the first ring, and a more negative EBE than Barreca 36. Both of these had a sidechain off the indole group with two carbonyl groups and an ethoxide group to contribute to additional interactions. All of the Zhao molecules had the same basic structure except for a sidechain off the second ring. The sidechain of the four Zhao molecules (5l, 5n, 5o, and 5p) was large compared to the other Zhao molecules and could contribute to hydrophobic (5-carbon extension) or hydrogen bond (−NH2 group, 2-carbon sidechain plus an hydroxyl group, 2-carbon extension plus an ester) interactions with the protein and therefore tighter binding. The structure of Zhuang 7 was similar to the other Zhuang molecules except there was a sidechain with a cyclic sulfonamide group. The number of HBA for the Zhuang 7 molecule was 10 and HBD, one, so there were more increased possibilities for hydrogen bond interactions.
11.2. Mutant Model, Y143R
There were structural considerations for mutation of tyrosine (Y) to arginine (R) that could affect the binding of the drug-like molecules in the active site. Tyrosine has six possible HB interactions and a benzene ring that can participate in π − π stacking, a hydrophobic interaction. The change to arginine increased the number of HB possibilities to 12: five HBA and seven HBD, but with loss of π − π stacking without the benzene ring. Approximately 58% of the drug-like molecules had a more negative EBE with the mutant model and there was a ≥1 kcal/mole difference for 33% of the molecules. Narrowing that group of molecules to an EBE ≥ −10.5 kcal/mole there were five drug-like molecules: Barreca 40 and 41, Zhao 5l, Zhuang 1 and 7. There were a total of eight possible HB interactions for Barreca 40, 41, and Zhuang 1, seven for Zhao 5l, and eleven for Zhuang 7. However, there were three other molecules for which there were ten possible HB interactions and two of those, Zhao 5m and 5p, had a ≥1 kcal/mole difference and a more negative EBE with 1BIS.pdb.
All of the nondrug-like molecules had at least 13 possible HB interactions. Four of the molecules had a tighter binding with the mutant model but only two, Zhao 5rS and 5tS, had a ≥1 kcal/mole difference between 1BIS.pdb and the mutant model. There were also two molecules, Zhao 5sR and 5sS, where the binding was tighter with 1BIS.pdb with a ≥1 kcal/mole difference. Again, there did not appear to be a structural connection for tighter binding based on HB interactions.
11.3. Mutant Model, N155H
Asparagine (N) has an amide group and six HBA and five HBD. Histidine’s R group is an imidazole ring that can participate in π − π stacking, a hydrophobic interaction. It also has five HBA and four HBD. 21% of the drug-like molecules had a ≥1 kcal/mole difference between 1BIS.pdb and the mutant model with a more negative EBE for the mutant model. There were only two of those molecules, Barreca 43 and Zhao 5e, where the EBE was ≥−10.8 kcal/mole. Both molecules could participate in hydrophobic interactions such as π − π stacking with the imidazole ring of histidine because of the presence of single and fused rings on the two molecules. Zhao 5e had one more ring than Barreca 43 but this did not contribute to tighter binding with the mutant model. They each had at least one fluorine plus six HBA and three and two HBD for Zhao 5e and Barreca 43, respectively. However 29% of the molecules had a ≥1 kcal/mole difference between 1BIS.pdb and the mutant model with tighter binding for 1BIS.pdb. Two of those molecules had a ≥−10.9 kcal/mole EBE: Barreca 37 and Barreca 41. Both of the molecules had a ring with a fluorine, linker, two fused rings and a sidechain. Barreca 37 also had a chlorine off the fused ring. They differed somewhat in their sidechains but their EBE was similar.
With two exceptions, 5rR and 5tS, the nondrug-like molecules had tighter binding with 1BIS.pdb. Structurally one was an R configuration and the other S and there was only one more HB interaction for 5tS. These differences had a minor impact on the EBE of the drug with the mutant model. There was one molecule, 5sR, with a ≥1 kcal/mole difference between 1BIS.pdb and the mutant model with tighter binding with 1BIS.pdb. The EBE with 1BIS.pdb was −12.9 kcal/mole and −9.8 kcal/mole with the mutant model.
Do the molecules with tighter binding, with either of the two mutant models, share structural characteristics that might be important when designing drugs for HIV-1 integrase? To answer the question, the structure of the drug-like molecules was considered when there was a ≥1 kcal/mole difference between 1BIS. pdb and either of the two mutant models. The two mutations, however, are only two of a number of HIV-1 integrase mutations that have led to increased resistance of HIV-1 integrase drugs such as raltegravir and elvitegravir  .
The drug-like molecules were the following: Barreca 19, 21, 40, 41, and 43; Zhao 5e, 5l, and 5v; and Zhuang 1 and 7. To address the structural characteristics that the 10 molecules share in common, a pharmacophore (Figure 24) was developed using a 3D pharmacophoric potential implemented on a grid  . The IUPAC defines a pharmacophore as “an ensemble of steric and electronic
Figure 24. Pharmacophore developed using atomic property fields. The figure was computationally generated by Andy Orry, Molsoft, L.L.C.
features that is necessary to ensure the optimal supramolecular interactions with a specific biological target and to trigger (or block) its biological response”  . Using the Atomic Property Fields (APF) developed by Molsoft, L.L.C.  seven properties are assigned from empiric physico-chemical components  : hydro- gen bond donors and acceptors, sp2 hybridization, lipophilicity, size, electro- positive/negative, and charge. The structural properties of the pharmacophore that contribute to interactions between a drug and active site amino acids are displayed in Figure 24.
The most obvious feature of the pharmacophore is the aromaticity of the rings (shown in white). All of the molecules had some combination of single aromatic rings and/or two fused rings―a common feature in HIV-1 integrase drugs  . The aromaticity allows for π − π stacking with amino acids in the active site. There were a greater number of hydrogen bond acceptors (shown in red), than hydrogen bond donors (shown in blue) for the 10 molecules. The average was six HBA and 2 HBD. Hydrophobic/lipophilic areas (yellow) are shown by the presence of fluorine atoms on the molecules, also present on the newest HIV-1 integrase drug, Doultegravir  . The aromatic rings are also considered hydrophobic. The hydrophobicity of a molecule is viewed as the driving force for the binding of a protein and a potential drug  because water molecules are displaced (desolvation) around hydrophobic surfaces. Increasing the hydrophobicity of the drug permits better binding in hydrophobic pockets; this results in a free energy gain. Charged areas on a potential drug and surrounding amino acids can be pH dependent with a strong dependence on binding energy because of the changes in proton affinity of ionizable groups  .
When drugs are designed, Lipinski’s Rule of Five should be considered for better bioavailability, that is, the importance of an unchanged drug reaching the systemic circulation. Intermolecular forces, such as hydrogen bonding, hydrophobic interactions, and charge-related bonding between a molecule and receptor can contribute to a more negative estimated binding energy. To minimize the side effects of a molecule, the lowest concentration of the molecule resulting in a physiological response in the active site is paramount. In the case of HIV-1 integrase, the goal of the molecule-receptor interaction is to prevent integration of viral DNA into host DNA.
J.B.E. thanks her administrators for yearly support of the ICM-Pro license. Thank you to Andy Orry, Molsoft, L.L.C., for providing the 2D ligand interaction figure and the pharmacophore.
HIV-1: human immunodeficiency virus
RT: reverse transcriptase
RNA: ribonucleic acid
DNA: deoxyribonucleic acid
D: amino acid, aspartic acid
E: amino acid, glutamic acid
Y143R: amino acid number 143, tyrosine (Y), is mutated to arginine (R)
N155H amino acid number 155 asparagine (N), is mutated to histidine
EBE: estimated binding energy
HBA: hydrogen bond acceptor
HBD: hydrogen bond donor
NTD: amino-terminal domain, amino acids 1-50
CCD: core catalytic domain, amino acids 50-212,
CTD: carboxyl-terminal domain, amino acids 212-288
APF: Atomic Property Fields
Submit or recommend next manuscript to SCIRP and we will provide best service for you:
Accepting pre-submission inquiries through Email, Facebook, LinkedIn, Twitter, etc.
A wide selection of journals (inclusive of 9 subjects, more than 200 journals)
Providing 24-hour high-quality service
User-friendly online submission system
Fair and swift peer-review system
Efficient typesetting and proofreading procedure
Display of the result of downloads and visits, as well as the number of cited articles
Maximum dissemination of your research work
Submit your manuscript at: http://papersubmission.scirp.org/
Or contact firstname.lastname@example.org