Skip to main content

Advertisement

Pentacycloundecane lactam vs lactone norstatine type protease HIV inhibitors: binding energy calculations and DFT study

Abstract

Background

Novel pentacycloundecane (PCU)-lactone-CO-EAIS peptide inhibitors were designed, synthesized, and evaluated against wild-type C-South African (C-SA) HIV-1 protease. Three compounds are reported herein, two of which displayed IC50 values of less than 1.00 μM. A comparative MM-PB(GB)SA binding free energy of solvation values of PCU-lactam and lactone models and their enantiomers as well as the PCU-lactam-NH-EAIS and lactone-CO-EAIS peptide inhibitors and their corresponding diastereomers complexed with South African HIV protease (C-SA) was performed. This will enable us to rationalize the considerable difference between inhibitory concentration (IC50) of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptides.

Results

The PCU-lactam model exhibited more negative calculated binding free energies of solvation than the PCU-lactone model. The same trend was observed for the PCU-peptide inhibitors, which correspond to the experimental activities for the PCU-lactam-NH-EAIS peptide (IC50 = 0.076 μM) and the PCU-lactone-CO-EAIS peptide inhibitors (IC50 = 0.850 μM). Furthermore, a density functional theory (DFT) study on the natural atomic charges of the nitrogen and oxygen atoms of the three PCU-lactam, PCU-lactim and PCU-lactone models were performed using natural bond orbital (NBO) analysis. Electrostatic potential maps were also used to visualize the electron density around electron-rich regions. The asymmetry parameter (η) and quadrupole coupling constant (χ) values of the nitrogen and oxygen nuclei of the model compounds were calculated at the same level of theory. Electronic molecular properties including polarizability and electric dipole moments were also calculated and compared. The Gibbs theoretical free solvation energies of solvation (∆Gsolv) were also considered.

Conclusions

A general trend is observed that the lactam species appears to have a larger negative charge distribution around the heteroatoms, larger quadrupole constant, dipole moment and better solvation energy, in comparison to the PCU-lactone model. It can be argued that these characteristics will ensure better eletronic interaction between the lactam and the receptor, corresponding to the observed HIV protease activities in terms of experimental IC50 data.

Background

In spite of extensive investigations and clinical efforts made over more than two decades on the extinction of HIV, AIDS is still a substantial threat to global health [1-7]. HIV treatment comprises an amalgam of a wide arsenal of various drugs that target different stages in the viral replication cycle. The majority of synthesized drugs for blocking viral enzymes, act as transcriptase (RT) and protease (PR) inhibitors [8].

Lactams [9,10] and lactones [11] are well known anti-bacterial agents. Lactones exhibit high potency for the treatment of skin and soft tissue infections [11]. We have recently reported pentacycloundecane (PCU) derived cage peptides, functionalised as lactams [12,13] (1) diols [14,15] (3 and 4) and ethers [16] (5) which demonstrated in vitro HIV-protease activities. The basic structures of the various cage compounds are presented in Figure 1. One cage peptoid [15] also showed promising anti HIV PR activity. The synthesis of the PCU-lactam (1) [17-19] and lactone (2)12, [18,20] is well established.

Figure 1
figure1

Structures of different types of cages [12-14,16].

In this paper, we followed our PCU-lactam-NH-EAIS peptide (6) studies [12,13] based on the synthesis and testing of an analogous novel PCU-lactone-CO-EAIS compound (7) (Figure 2). The lactam peptide exhibits an order of magnitude better HIV protease inhibitory activity, in comparison with the lactone analogues. We have recently argued that the activity of PCU-lactam-NH-EAIS peptide stems from its function as a norstatine type transition state analogue [12,13]. Several potential reasons for the discrepancy in HIV PR activities exist. First, the lactam peptide involves C → N amino acid coupling while the lactone peptides consist of N → C coupling; this should induce very different binding energies. Second, the hydrogen bond interaction of the cage lactam group with the protease Asp25/25′ residues may be more advantageous.

Figure 2
figure2

Structures of PCU-lactam-NH-EAIS and PCU-lactone-CO-EAIS peptide inhibitors and their corresponding PCU-models. E: Glutamic acid; A: Alanine; I: Isoleucine; S: Serine (Note that the PCU-models consist of two enantiomers. The cage peptides therefore exist as two diastereomers).

The objective of this study is to investigate MD-based binding affinities and electronic structural features of the proposed compounds to obtain a theoretical explanation for the significant difference in the experimental IC50 data for lactam and lactone inhibitors [12,13].

Extensive molecular modeling techniques on PCU compounds have been applied in our laboratory [12-16,21-30] and we aimed to use these and alternative quantum chemical approaches [31-39] to gain more insight into the discrepancies of the observed bio-activities.

Structures of PCU-lactam-NH-EAIS, and PCU-lactone-CO-EAIS inhibitors and their corresponding PCU-models are demonstrated in Figure 2.

From a computational standpoint, our strategy in this work is summarized as follows:

First, molecular docking and 10 ns MD calculations have been performed for both the PCU-peptide and model compounds as well as their corresponding diasteromers complexed to South African HIV protease (C-SA) in aqueous solution using both MM-PBSA (explicit water as solvent) and MM-GBSA (implicit solvent) methods. The MM-PBSA method includes the calculation of the molecular mechanics gas phase energies, polar continuum electrostatic solvation energies (by solving the linearized Poisson-Boltzmann equation), and non-polar surface area energies. MM-GBSA is a faster method than MM-PBSA and in this approach the GB model approximates the electrostatic contribution to the free energy of solvation. Both of these methods give efficient, reproducible, and reliable binding free energies of solvation [40-43].

A comparative study to determine the MM-PB(GB)SA Gibbs binding free energies of solvation for the PCU-peptide and model compounds was then applied to compare the binding affinities. This will enable us to determine the cause(s) for the significant difference between the IC50 values of PCU-lactam-NH-EAIS in compared with PCU-lactone-CO-EAIS peptides.

Second, we assumed that the considerable difference between experimental IC50 values could be due to a higher charge density on the lactam in comparison to lactone. To verify our assumption, further DFT analysis of these models aimed to explore the electronic structural influences of the lactam versus the lactone moiety. In this regard, NQR studies [44,45] have emerged as a promising theoretical method to study potential correlations between the bioactivity of compounds and their electronic characteristics. Hence, NQR parameters, namely, asymmetry parameters (η) and quadrupole coupling constants (χ) of nitrogen and oxygen nuclei as well as atomic charges derived from natural population analysis (NPA) were calculated for a PCU-model (Figure 2) derived from PCU-lactam, lactim and lactone-CO-EAIS peptides. The tautomeric form of the lactam i.e. lactim [19,20] was studied before.

Finally, to evaluate the natural population analysis (NPA) atomic charges, the electrostatic potentials were calculated and regions rich in electronic charge density were characterized. Polarizability, dipole moment, and Gibbs free energies of solvation (∆Gsolv) values of the title molecules in aqueous solution were also calculated at the same level of theory. Such a theoretical study is useful to address the experimental observations through molecular dynamic and electronic structure features of these inhibitors. The outcome of this work could be helpful for the future design of novel therapeutic agents that inhibit different relevant diseases.

Methods

Synthesis of cage lactone peptides

The cage lactone (2) [18] was synthesized from Cooksen’s dione [46,47] using a reported procedure utilizing sodium cyanide to generate cyano cage lactam which was converted to cage lactone (2) using hydrochloric acid (Scheme 1).

Scheme 1
scheme1

Synthesis of the cage lactone (i) NaCN, CH3COOH, H2O (ii) 3 M HCl.

The chosen peptide sequence (EAIS, (S)-amino acids were used) is first synthesised on rink resin and the cage lactone is then coupled in the final step. HCTU and DIPEA were used as coupling reagents and piperidine/DMF (2:8) was used to deprotect the Fmoc group from the amino acid. The final cleavage of peptide from the resin was achieved using a mixture of TFA and DCM (95:5% v/v) (Scheme 2).

Scheme 2
scheme2

Synthesis of cage lactone peptides (i) 20:80 piperidine:DMF (ii) HCTU, DIPEA, DMF (ii) 95:5% TFA:DCM.

Over-expression, extraction and purification of the C-SA protease and in vitro HIV-1 protease activity

Over-expression, extraction and purification of the C-SA protease and in vitro HIV-1 protease activity were reported in our earlier publications [12,13]. Genetic diversity among HIV-1 subtypes ranges between 25 and 35% at the nucleotide level and variation within subtypes vary from 15 to 20% [48,49]. Thus, both the high mutation rate within the HIV PR and the genetic variation of HIV is an ongoing challenge toward finding a potent protease inhibitor (PI) with reasonable inhibition of various HIV subtypes [50].

The HIV-1 subtype C PR that predominates in South Africa (C-SA PR) is defined as the consensus subtype C PR [51]. This C-SA PR contains eight polymorphic sites (T12S, I15V, L19I, M36I, R41K, H69K, L89M, and I93L) in relation to the wild-type subtype B HIV-1 PR [52]. Multiple sequence alignment of the C-SA PR, subtype B PR, subtype B-MDR PR, and subtype C-2R8N PR confirmed that subtype C of South African HIV-1 protease is structurally very close to the rest of the HIV-1 proteases [53].

The procedures followed were the same as reported before 12,13.

Molecular dynamic simulation

The geometry optimizations of the considered ligands were performed using the GAUSSIAN03 package [54] at the B3LYP theoretical level and with the 6-311G(d,p) basis set. Afterwards, the partial atomic charges were calculated using the restricted electrostatic potential (RESP) [55,56] method using HF/6-31G(d). The RESP charges were generated based on the calculated electrostatic potential using antechamber implemented in the AMBER 12 package [57,58]. The force field parameters for the ligands were described by General Amber Force Field (GAFF) [59]. The standard AMBER ff99SB force field [60,61] for bio-organic systems was used to describe the protein parameters.

Starting structures for MD simulations were obtained by molecular docking of PCU-models and PCU-EAIS peptide inhibitors as well as their corresponding diastereomers inside the active pocket of HIV PR using Autodock4 program [62]. In order to make sure that all the structures have the same orientation inside the active pocket of HIV PR, all docked structures were aligned using Pymol [63] program. Protons were added to the protein in accordance with pKa calculations using PropKa 3.1 server [64]. The physiological pH conditions and the correct protonation state of ionizable groups in the CSA HIV PR enzyme was considered for molecular dynamic simulations. In this context, one of the aspartates (Asp25) of the catalytic site of the the CSA HIV PR exhibits an increased pKa value of 5.2 in the inhibitor-bound protease [65] while no increased pKa was observed for the free form of the protease (pKa. 4.5) [66]. We have previously reported [14] that different protonation states do not have an adversed effect on the results.

The Sander module implemented in AMBER 12 [57] together with AMBER ff99SB force field [60,61] were used to perform all MD simulations as well as for the minimizations and equilibration protocols. Each system was solvated in a truncated cubic TIP3P water [67] box with 10 Å distance around the complex. Minimizations were carried out in a constant volume (NVT ensemble) by 5000 cycles of steepest descent minimization followed by 5000 cycles of conjugated gradient minimization under harmonic restraints with force constant 500 kcal · mol−1 · Å−2 to all solute atoms. Then, the minimization was followed in a constant volume by 10000 cycles of steepest descent minimization followed by 10000 cycles of conjugated gradient without any harmonic restraints for all atoms of the system. Then, each system was gently annealed from 0 K up to 300 K for 70 ps with an isobaric ensemble. Finally, the 10 ns canonical ensemble (NVT) MD simulations were applied without any restraints. The temperature was regulated at 300 K using the Langevin thermostat [68]. For all MD simulations, 2 fs time step and 12 Å non-bonded cutoff were used. The particle mesh Ewald method [69] was used to treat long-range electrostatics, and bond lengths involving bonds to hydrogen atoms were constrained with the SHAKE algorithm [70]. The root mean square deviation (RMSD) of the average structures over the MD trajectory between a defined starting point of the simulation and all succeeding frames were also calculated.

Since the chances of the inhibitor/complex being trapped as a local minimum structure is quite high, an iterative process was followed where the inhibitor peptide backbone from the lowest energy complex [say PCU-lactam-NH-EAIS(a)] for each of the diastereomeric peptides (see Tables 1 and 2) was imposed on the related analogue [i.e. PCU-lactone-NH-EAIS(a)], for a subsequent MD run. The positions of the pairs of diastereomers (lactam versus lactone) with respect to the active site (particularly with respect to the catalytic Asp25/25′ residues) during the MD were monitored to ensure that the two inhibitors maintained very similar positions during the MD simulations. Several different docked starting structures were used in order to ensure that the lowest possible energy inhibitor/complexes were obtained. This procedure also enabled us to draw meaningful comparisons between the binding energies of the respective lactam/lactone-peptide diastereomers.

Table 1 Binding free energies of solvation and its components for the PCU-based ligands complexed with the HIV protease in in kcal/mol (the 3D structures of these compounds are provided with the supporting information)
Table 2 Binding free energies of solvation and its components using normal mode analysis for the PCU-based ligands complexed with the HIV protease in kcal/mol

In order to obtain a better theoretical understanding about the difference characteristics between the lactam/lactone-peptides, two pairs of theoretical cage peptides were studied, namely PCU-lactone-NH-EAIS(a and b) and PCU-lactam-CO-EAIS(a and b). Due to the current lack of synthetic procedures for the required lactam/lactone starting structures, the synthesis of these peptides are not feasible.

MM-PB (GB) SA binding free energy calculations

MM-PB (GB) SA (Molecular Mechanics-Poisson-Boltzmann or Generalized Born solvent-accessible Surface Area) methods provide an effective computational tool in the analysis of biomolecular interaction [71-73]. This approach is based on the calculation of the average free energies of solvation (ΔGbind) between a target protein and a set of ligands over the trajectory of molecular dynamics (MD) simulation. In this method [74,75], the ΔGbind between a ligand (L) and a receptor (R) to form a complex (RL) is defined as:

$$ \Delta {\mathrm{G}}_{\mathrm{Bind}}=\Delta \mathrm{H}\ \hbox{--}\ \mathrm{T}\mathrm{D}\mathrm{S} \mbox{\fontencoding{T1}\fontfamily{cmr}\selectfont\char"14} \Delta {\mathrm{E}}_{\mathrm{MM}}+\Delta {\mathrm{G}}_{\mathrm{sol}}\hbox{--}\ \mathrm{T}\mathrm{D}\mathrm{S} $$
(1)
$$ \Delta {\mathrm{E}}_{\mathrm{MM}}=\Delta {\mathrm{E}}_{\mathrm{Int}}+\Delta {\mathrm{E}}_{\mathrm{E}\mathrm{le}}+\mathrm{D}{\mathrm{G}}_{\mathrm{VDW}} $$
(2)
$$ \Delta {\mathrm{E}}_{\mathrm{Sol}}=\Delta {\mathrm{G}}_{\mathrm{PB}/\mathrm{GB}}+\Delta {\mathrm{G}}_{\mathrm{SA}} $$
(3)

Where ΔEMM, ΔGSol and -TΔS are the changes of the gas phase MM energy, the free energy of solvation, and the conformational entropy upon binding, respectively. The ΔEMM term includes ΔEInt (bond, angle, and dihedral energies), ΔEEle (electrostatic), and ΔEVDW (van der Waals) energies. The ΔGSol is the solvation binding free energy as the sum of polar (electrostatic solvation energy) and nonpolar free energies of solvation (ΔG PB/GB, and the non-electrostatic solvation component, ΔGSA). The polar contribution is calculated using either the GB or PB model, while the nonpolar energy is estimated by solvent accessible surface area (SASA) [76,77].

DFT study

Gas phase calculations were performed by using GAUSSIAN09 [78] density functional theory (DFT) with the B3LYP functional [79,80]. The 6-311G(d,p) basis set was used for both the geometry optimization and electronic structure calculations. This basis set includes diffuse functions on hydrogen and the heavy atoms and is a suitable option for polar molecules with electron lone-pairs [81]. B3LYP non-local hybrid exchange correlation functional that includes a mixture of Hartree–Fock exchange with DFT exchange–correlation is known as the commonplace substitute approximation. LYP supports the full correlation energy and not only a correction to local spin density approximation (LSDA), its overall performance is sufficiently good for organic molecules [82]. Harmonic vibrational frequencies were also calculated to confirm that all structures were minima on the potential energy surface and to calculate the zero-point vibrational energy (ZPVE) and the Gibbs free energy. Various atomic and molecular electronic structure quantities were calculated from optimized structures. For the calculation of atomic charges, natural bond orbital analysis (NBO) [83-85] was performed as well as the natural population analysis (NPA) as an optimal wave function-based method. To accurately analyze the charge distribution around these atoms, the electrostatic potential energy values were calculated and mapped over an isodensity surface corresponding to 0.002 a.u. This electrostatic potential surface includes the van der Waals volumes of the individual atoms in the molecule and is thus a good representation of the reactive regions around the molecules [86,87].

Moreover, the NQR spectral quantities including asymmetry parameter (η) and quadrupole coupling constant (χ) of nitrogen and oxygen atoms, shine useful insight on biological effectiveness of drug-like compounds [44]. These two quantities can be derived from the electric field gradient tensor using the following equations:

$$ \upchi =\frac{{\mathrm{e}}^2\mathrm{Q}\left\langle \mathrm{q}\mathrm{z}\mathrm{z}\right\rangle }{\mathrm{h}} $$
(4)
$$ \upeta =\frac{\left\langle \mathrm{q}\mathrm{z}\mathrm{z}\right\rangle -\left\langle \mathrm{q}\mathrm{y}\mathrm{y}\right\rangle }{\left\langle \mathrm{q}\mathrm{z}\mathrm{z}\right\rangle } $$
(5)

Where e < qxz>, e < qyy > and e < qzz > are the principal components of the effective electric field gradient tensor defined such that │ < qzz│ < qyy > │ < qxx> > │and Q is the nuclear electric quadrupole moment [88]. The electric quadrupole moments (Q) of 17O and 14 N atoms were taken as 25.58× 10−27 cm2 and 20.44 × 10−27 cm2, respectively [89].

We calculated the asymmetry parameters (η) and quadrupole coupling constants (χ) of nitrogen and oxygen atoms involved in the cage region of PCU-lactam, its tautomer lactim [19,20] and PCU-lactone models in vacuum and in solvent media (water, DMSO and dioxane). The solvation effect (aqueous solution) on the Gibbs free energy of solvation (ΔGsol) were also evaluated utilizing the self-consistent reaction field (SCRF) keyword [90] with the solvation model on density (SMD) [91]. The SMD model computes the electrostatic interaction based on the integral-equation-formalism polarizable continuum model (IEF-PCM) and complements the description of the solute-solvent interaction by adding the missing non-electrostatic term to make accurate predictions of ΔGsol. The latter were calculated in vacuum/gas phase and in solvent (water) for the titled systems by taking the difference of the obtained ΔGsol values in both gas and aqueous media.

Results and discussion

Inhibitory activity of anti-HIV PR compounds

The activity data of HIV PR inhibitors are presented in Table 3. Peptides of the cage lactone furnished significant HIV PR activities. However, the PCU-lactone-CO-EAIS peptide (Figure 2) exhibited an IC50 value of ~0.85 μM, which was eleven times weaker than the corresponding cage lactam peptide (IC50 ~ 0.078 μM). The shorter PCU-lactone-CO-EAI peptide (Figure 3) demonstrated poor activities (IC50 > 2.6 μM) while the PCU-lactone-CO-AIS peptide was a better inhibitor (IC50 ~ 0.83 μM) (Figure 3).

Table 3 Wild type CSA HIV-1 protease inhibition activities of PCU-lactam-NH- and lactone-CO-peptides
Figure 3
figure3

The structures of PCU-lactone-CO-peptides.

MM-PB (GB) SA binding free energy analysis

MM-PB (GB) SA binding free energy method was employed to estimate the average free energies of solvation of the titled complexes obtained over 10 ns MD simulations. The dynamics and stability of the considered complexes were also investigated based on the root mean square deviation (RMSD) graphs. The RMSD graphs of PCU-lactam-NH-EAIS (a,b) and its corresponding lactone, i.e., PCU-lactone-NH-EAIS (a,b) structures over 10 ns MD trajectories are depicted in Figure 4.

Figure 4
figure4

Analysis of the RMSD of energy minimized (1) PCU-lactam-NH-EAIS (a,b) and (2) PCU-lacton-NH-EAIS (a,b) obtained over the course of 10 ns MD trajectories. a and b = diastereomers [MD in explicit solvent model].

According to the plotted RMSD graphs in Figure 4, all the considered PCU-EAIS peptides maintained the same level of stability with an RMSD of 0.80 Å during the entire time of the 10 ns MD simulation.

In order to calculate the theoretical binding energies, we used our MD simulations at 10 ps intervals, as suggested by Stoica and co-workers [92]. The free energies of solvation for all selected complexes averaged over the trajectories of the explicit water 10 ns MD simulations and are listed in Tables 1 and 2. Note that ΔGbind for both MM-PB(GB)SA methods and its components including polar electrostatic (ΔEEle), van der Waals (ΔEVDW) and PB or GB solvation energies [ΔESol(GB) and ΔESol(PB)] are defined in equations 1, 2 and 3.

The first important observation is that for all cases, the PCU-lactam-peptides diastereomers showed better binding energies (more negative ΔGbind) than the corresponding lactone-peptide diastereomers. The theoretical results [both ΔGBind(PBSA) and ΔGBind(GBSA)] for Lactam-NH-EAIS(a) peptide showed better binding affinity against CSA-HIV PR in comparison to PCU-lactone-CO-EAIS(a) peptide inhibitors. On the other hand, the case for (b) has slightly higher binding energy with ΔGBind(PBSA) for the lactam peptide, while with ΔGBind(GBSA) the same pattern as before was observed (better/lower binding energy for the lactam peptide). These results largely support the lower experimental IC50 values for PCU-lactam-NH-EAIS inhibitor (IC50 = 0.076 μm) than the PCU-lactone-CO-EAIS compound (IC50 = 0.850 μm).

The same trend is observed for the calculated binding energies of the isolated lactam/lactone model compounds. The lactam (both enantiomers) give better binding energies than the corresponding lactone. It therefore appears likely that the source of better binding energy for the cage peptides originates from the PCU skeleton.

Since the experimental IC50 values for the peptides were measured for a diastereomeric mixture [resulting from the two cage enantiomers attached to the (S)-peptide] it makes sense to look at the average of the calculated binding free energies of the PCU model compounds as well. In the case of both PB and GB solvation models the same order for binding free energy values was observed for the PCU models (top four models in Tables 1 and 2): Lactam-PCU < Lactone PCU (i.e. PCU-lactam models (either a or b) exhibit stronger ΔGbind in comparison to the PCU-lactone). It is important to note that for both cases (PCU-model and peptide) that the PCU-lactone consistently gave the weakest binding energies.

Amongst the component of binding free energies, the van der Waals (VDW) interactions between the ligands and the HIV protease showed the largest contribution to the ΔGbind. This is the case for both MM-PBSA and MM-GBSA methods. However, the mentioned difference was more evident in MM-GBSA data. Two cases for the electronic energy contributions are worth pointing out. Lower ΔEEle value (−11.47 kcal/mol - smaller difference) was reported for the PCU-lactone-CO-EAIS(a) in comparison to the PCU-lactam-NH-EAIS(a) peptide (−22.36 kcal/mol). Lower ΔEEle value of this compound is attributed to the potential electronic interactions (Coulombic, dipole-dipole interactions, hydrogen bonds etc.). We have previously demonstrated [12,13] that the hydroxyl group (−O3H, see Figure 2) of the PCU-lactam is very important for HIV PR activity since it is presumed to interact with the two catalytic Asp25/25′ residues of HIV PR. Analysis of this distance between the PCU hydroxyl group (−O3H) and Asp revealed a weaker interaction for the PCU-lactam-NH-EAIS (average of about 13 Å) than for the corresponding PCU-lactone peptide (average of about 6.5 Å, see Additional file 1). This can potentially explain the lower ΔEEle value for the PCU-lactone-CO-EAIS (a).

We calculated the entropic contribution using normal mode analysis for the PCU-lactam-NH-EAIS(a) and PCU-lactone-NH-EAIS(a) (Table 2).

The binding free energy has slightly improved, but the same trend was observed.

A more detailed study will be reported later to address the potential segment of the PCU-based compounds causing this change in their binding affinities.

Electronic structural analysis

In this section, we considered the PCU-lactam, lactim and lactone models shown in Figure 5. The two reported sets of 1H NMR signals in various solvents suggested the lactim–lactam tautomeric equilibrium exists [19,20]. In (CD3)2SO the major tautomer is the lactam form, while for dioxane the lactim is dominant. It is therefore most likely that the lactam form is the active form in aqueous solution (also used for the IC50 experiments). The conformation for the amide side chain with the lowest energy was determined. In both cases this conformation involved intramolecular hydrogen bonding (O5 - - HN1 for the lactam and N4H - - O1 for the lactone).

Figure 5
figure5

The electrostatic potential map [B3LYP/6-311G(d,p)] for PCU models in vacuum. (red = electron rich and blue = electron defficient).

First, we have focused on the natural atomic charges obtained from NBO analysis and then we assessed the electrostatic potential map produced with the Molekel [93] program. NBO analysis calculates atomic charges, by summing occupancy of natural atomic orbitals (NAOs). These charge values for nitrogen and oxygen nuclei of the considered models are presented in Table 4.

Table 4 Natural atomic charges (a.u.) on nitrogen and oxygen nuclei of PCU-lactam, PCU-lactim and PCU-lactone models obtained from NBO analysis [B3LYP/6-311G(d,p)]

According to the reported data in Table 4, larger negative atomic charges were observed for most of the hetero-atoms (N1, O2, O3 and N4) of PCU-lactam model in comparison to the corresponding atoms in the lactone case. The less negative atomic charge of O5 for PCU-lactam (−0.63987 a.u) than the corresponding value for the PCU-lactone (−0.65063 a.u) can be explained due to the intramolecular hydrogen bonding interaction in the lactam compound.

Analysis of the charge properties mapped onto the optimized structure will assist with the characterization of how the electronic properties are related to chemical activity [86]. In addition, electrostatic potential map (ESP) exhibits the charge distribution of a molecule based on the properties of the nucleus and nature of electrostatic potential energy. The molecular ESPs of these compounds are illustrated in Figure 5. This assists one to visualize different charged regions of these molecules.

The lactam structure which is supposed to interact with the two aspartic acid residues in the HIV PR (Figure 6, based on general mechanism) [94] is clearly more electron rich in the area (O2 and O5 atoms) (Figure 5). The corresponding regions for the lactone (O2 and N4H) appeared less electron rich for O2 and electron difficient for N4H. With the C = O2 group of the PCU-lactam model, clearly more polarized than the corresponding carbonyl of the lactone. To recognize reactive domains on the compound’s surface in solution, the calculation of ESP charges was also performed in water. According to the obtained ESP values (data are provided in Additional file 1), it is evident that the polarization of these compounds in aqueous medium increase. Again, the same trend is observed. It can be envisaged from the mechanism of HIV-protease that interaction between the aspartic acid groups and the lactam may be more pronounced. This is potentially a reason for better activity (lower IC50) for PCU-lactam inhibitors in comparison to the lactone series.

Figure 6
figure6

Potential reaction of PCU-lactam-NH-AIS peptide inhibitors with HIV protease based on the general HIV PR mechanism [94].

The electronic charge distribution around the nitrogen and oxygen nuclei was further investigated through the calculation of asymmetry parameters (η) and quadrupole coupling constants (χ) generally known as NQR parameters. The calculated parameters for the 14N and 17O atoms of PCU-lactam, lactim and PCU-lactone models in vacuo, water, DMSO, and dioxane are reported in Table 5. These values were obtained from the electric field gradient (EFG) tensor calculations and characterized the local electron density distribution around these atoms. Indeed, the interaction of the nuclear quadrupole moment (Q) with the EFG identifies the degree of the double bond character of chemical bonds in which the quadrupole nuclei is involved [88]. The relationship between asymmetry parameters (η) and quadrupole coupling constants (χ) for oxygen atoms (including carbonyl and alcohol oxygens) was studied [95]. It was reported that the relationship is ordered, but not necessarily linear.

Table 5 Asymmetry parameter (η) and the quadrupole coupling constant (χ) (MHz) of the 14N and 17 O atoms for the PCU-lactam, PCU-lactim and PCU-lactone models [B3LYP/6-311G(d,p)] in different solvent media

The departure of the EFG tensor from the axial symmetry is characterized by the parameter η. The larger the value, the smaller the deviation from axial symmetry. This deviation can be caused by electronic influences such as by hydrogen bonding or packing forces in crystals.

A higher value of the asymmetric parameter (η), around 0.9 was obtained for the oxygen atoms of the lactim hydroxyl group (O2) than that of the lactam carbonyl oxygen atom (0.1). It is notable that larger values for η were observed in more polar solvents, as expected.

From the comparison of the measured IC50 for the synthesized peptides with the obtained ΔGBind(PB/GB-SA) values it was argued that the source of stronger binding affinity for the cage peptides originates from the PCU skeleton. In order to verify our assumption, it serves best to focus on the quadrupole coupling constant data (χ) of the cage (PCU) part of the molecule, which is also expected to interact with the aspartase segments of the PR (Figure 6). Our previous results have shown that O2 and O3 atoms are crucial in this regard for the PCU-lactam-NH-EAIS peptide inhibitors [12,13]. According to the χ values of the hetero-atoms, it is clear that in each of the four cases (vacuum, water, DMSO and dioxane) the χ values for these atoms (O2 and O3) were smaller for the PCU-lactam than for the lactone model. Smaller χ values implied that the quadrupole charge on these atoms were more delocalised. This will enable a stronger hydrogen bond interaction of these atoms (O2 and O3) with the Asp25/Asp25′ groups of the PR. The same observation holded for O3 of the lactim compound. As was the case with the asymmetry parameter, the quadrupole constant of O2 atom for the lactim also reflects the change from carbonyl oxygen to the oxygen atom in hydroxyl group.

The oxygen atom (O5) involved in the lactam side chain carbonyl exhibited less charge density (lower χ values) versus the corresponding lactone oxygen atom. Since the lactam side chain carbonyl (O5) is in fact involved in intramolecular hydrogen bonding with its NH group (N1H), it should be more delocalized resulting in a smaller χ value. This observation is in good agreement with the reported natural atomic charges on the oxygen nuclei of PCU-lactam and lactim (Table 4).

For more information about the electronic character of these compounds, the frontier orbitals (HOMO and LUMO), polarizability and dipole moments of PCU-lactam, lactim and lactone models were calculated next. The HOMO and LUMO frontier molecular orbitals are demonstrated in Figure 7.

Figure 7
figure7

The frontiers orbitals (HOMO and LUMO) of PCU-lactam (a1 and b1), PCU-lactim (a2 and b2) and PCU-lactone (a3 and b3) models [B3LYP/6-311G(d,p)].

It is evident that HOMO orbital of the lactam functional group is larger in comparison to corresponding orbitals located around the lactone and lactim cases. The size of the HOMO orbitals (which is an indication of the electron density) of the PCU-lactam and lactone models can be rationalized in terms of the corresponding experimental binding affinities. The lactam model exhibits larger electron density, which can potentially result into a better electrostatic (i.e. hydrogen bond) interaction between the substrate and the Asp25/Asp25′ groups. The LUMO around lactim functional group appeared to be most prominent of the three compounds.

Polarizability of a molecule is defined as the ease of distortion of the electron cloud of a molecular entity by an electric field (such as when in close proximity of a charged reagent). It is experimentally measured as the ratio of induced dipole moment (μind) to the field E which induces it. The dipole moment is the first derivative of the energy with respect to an applied electric field. It is used as a measure of the asymmetry in the molecular charge distribution. Much effort has been made to look for a correlation between the electric dipole moment of drug-like compounds and their biological activity [96-98] which requires the search for such a correlation between electron density distribution in a molecule of a given compound and its activity [44].

The polarizability and dipole moments for these compounds were calculated and is reported in Table 6.

Table 6 Polarizability (Å 3 ), dipole moment (Debye) and Gibbs free energy of solvation ΔG solv (kcal/mol) values of the PCU-lactam and PCU-lactone models [B3LYP/6-311G(d,p)]

The higher values of both polarizability and dipole moment as well as the more negative Gibbs free energy of solvation, ∆Gsolv for the PCU-lactam in comparison to the lactone model (Table 6), clearly confirms the more polar character of the PCU-lactam model. This characteristic appears to contribute towards a higher binding affinity of the peptide derivatives. The more negative MMPB(GB)SA binding free energy of solvation, ∆Gbind for the PCU-lactam in comparison to the lactone model Tables 1 and 2 is consistent with the more negative Gibbs free energy of solvation, ∆Gsolv for the PCU-lactam in comparison to the lactone model reported in Table 6.

Conclusion

A series of three novel PCU-lactone-CO-EAIS peptides were synthesized and tested for HIV-protease activity. Two of them exhibited significant activities (~1 μM). The most active inhibitor amongst these synthesized PCU-lactone peptides was the PCU-lactone peptide 9 (IC50 ~ 0.80 μM) with the AIS side chain. Comparison of the MM-PB (GB) SA binding free energy data of this compound in lactone series with the one of the most active PCU-lactam-NH-EAIS peptide (IC50 = 0.076 μM) that was synthesized before in our laboratory, reflect a higher binding affinity of the PCU-lactam peptide peptides against South African HIV-protease than the lactone series. This result supports the experimentally observed trend for HIV-PR IC50 values of the PCU-lactam-NH-EAIS inhibitor (IC50 = 0.076 μm) versus the lactone-peptides (IC50 = 0.850 μm).

From an electronic structure standpoint, a relatively more negative atomic charge was observed on the oxygen and nitrogen atoms of cage lactam model in comparison with cage lactone. The NQR results revealed higher delocalisation of charge distribution around oxygen atoms (O2 and O3) of the PCU-lactam, which can potentially react better with the active catalytic aspartic acid residues of HIV PR. It appears that the higher charge density, polarizability and the dipole moment due to the hetero-atoms of cage-lactam plays an essential role in their higher experimental activity and binding affinity than the corresponding cage-lactone peptides.

Supplementary material

NMR spectra and HRMS data of all compounds are provided with the supplementary material. The 3D structures of all calculated compounds are also provided.

Abbreviations

PIs:

Protease inhibitors

HIV-PR:

HIV protease

PCU:

Pentacycloundecane

HCTU:

2-(6-Chloro-1H-benzotriazole-1-yl)-1,1,3,3-tetramethylaminium hexafluorophosphate

MM-PBSA:

Molecular mechanics/poisson-boltzmann surface area

MM-PB(GB)SA methods:

The molecular mechanics/generalized born surface area

References

  1. 1.

    Miller JF, Andrews CW, Brieger M, Furfine ES, Hale MR, Hanlon MH, et al. Ultra-potent P1 modified arylsulfonamide HIV protease inhibitors: the discovery of GW0385. Bioorg Med Chem Lett. 2006;16:1788–94.

  2. 2.

    Beck EJ, Avila C, Gerbase S, Harling G, De Lay P. Counting the cost of not costing HIV health facilities accurately pay now, or pay more later. Pharmacoeconomics. 2012;30:887–902.

  3. 3.

    Deeks SG, Lewin SR, Havlir DV. The end of AIDS: HIV infection as a chronic disease. Lancet. 2013;382:1525–33.

  4. 4.

    El-Sadr W. From HIV to global health: opportunities and challenges. J Int AIDS Soc. 2012;15:1.

  5. 5.

    Martin G, Grant A, D’Agostino M. Global health funding and economic development. Glob Health. 2012;8:8–12.

  6. 6.

    Morens DM, Fauci AS. Emerging infectious diseases: threats to human health and global stability. PLoS Pathog. 2013;9:e1003467- 9.

  7. 7.

    Stephens LL, Swanepoel CC-A, Van Rooyen BA, Abayomi EA. The human immunodeficiency virus, (HIV-1), pandemic: cellular therapies, stem cells and biobanking. Transfus Apher Sci. 2013;49:9–11.

  8. 8.

    Maes M, Loyter A, Friedler A. Peptides that inhibit HIV-1 integrase by blocking its protein-protein interactions. FEBS J. 2012;279:2795–809.

  9. 9.

    Minami NK, Reiner JE, Semple JE. Asymmetric synthesis of novel quaternary alpha-hydroxy-delta-lactam dipeptide surrogates. Bioorg Med Chem Lett. 1999;9:2625–8.

  10. 10.

    Fan G-j, Wang Z, Wee AGH. Regio- and diastereocontrolled C-H insertion of chiral gamma- and delta-lactam diazoacetates. Application to the asymmetric synthesis of (8S,8aS)-8-hydroxyindolizidine. Chem Commun. 2006;35:3732–4.

  11. 11.

    Lorenz-Baath K, Korotkov VS, Lierse von Gostomski C, Sieber SA. In vitro percutaneous absorption of (14) C-labeled beta-lactone promises topical delivery of new bacterial virulence inhibitors. ChemMedChem. 2012;7:1490–5.

  12. 12.

    Makatini MM, Petzold K, Sriharsha SN, Ndlovu N, Soliman MES, Honarparvar B, et al. Synthesis and structural studies of pentacycloundecane-based HIV-1 PR inhibitors: a hybrid 2D NMR and docking/QM/MM/MD approach. Eur J Med Chem. 2011;46:3976–85.

  13. 13.

    Makatini MM, Petzold K, Sriharsha SN, Soliman MES, Honarparvar B, Arvidsson PI, et al. Pentacycloundecane-based inhibitors of wild-type C-South African HIV-protease. Bioorg Med Chem Lett. 2011;21:2274–7.

  14. 14.

    Honarparvar B, Makatini MM, Pawar SA, Petzold K, Soliman MES, Arvidsson PI, et al. Pentacycloundecane-diol-based HIV-1 protease inhibitors: biological screening, 2D NMR, and molecular simulation studies. ChemMedChem. 2012;7:1009–19.

  15. 15.

    Makatini MM, Petzold K, Arvidsson PI, Honarparvar B, Govender T, Maguire GEM, et al. Synthesis, screening and computational investigation of pentacycloundecane-peptoids as potent CSA-HIV PR inhibitors. Eur J Med Chem. 2012;46:3976–85.

  16. 16.

    Karpoormath R, Sayed Y, Govender P, Govender T, Kruger HG, Soliman MES, et al. Pentacycloundecane derived hydroxy acid peptides: a new class of irreversible non-scissile ether bridged type isoster as potential HIV-1 wild type C-SA protease inhibitors. Bioorg Chem. 2012;40:19–29.

  17. 17.

    Martins FJC, Viljoen AM, Kruger HG, Joubert JA. Synthesis Of 8,11-dihydroxy-pentacyclo 5.4.0.02,6.03,10.05,9 undecane-8,11-lactam. Tetrahedron. 1993;49:9573–80.

  18. 18.

    Martins FJC, Viljoen AM, Kruger HG, Joubert JA, Wessels PL. Synthesis of delta-lactams from pentacyclo 5.4.0.0(2,6).0(3,10).0(5,9) undecane-8,11-dione. Tetrahedron. 1994;50:10783–90.

  19. 19.

    Martins FJC, Viljoen AM, Kruger HG, Wessels PL. Structure elucidation of 11-amino-8-hydroxypentacyclo 5.4.0.0(2,6).0(3,10).0(5,9) undecane-8,11-l actam through selective acetylation and complete H-1 and C-13 NMR spectral assignment of the mono-, di- and triacetates. Magn Reson Chem. 2004;42:402–8.

  20. 20.

    Kruger HG, Martins FJC, Viljoen AM. Interconvertions between delta-lactam and delta-lactone derivatives initiated by unique transannular interactions of the rigid cyclohexane boat structure in pentacycloundecane. J Org Chem. 2004;69:4863–6.

  21. 21.

    Kruger HG. Ab initio mechanistic study of the protection of alcohols and amines with anhydrides. J Mol Struct: (THEOCHEM). 2002;577:281–5.

  22. 22.

    Bisetty K, Gomez-Catalan J, Aleman C, Giralt E, Kruger HG, Perez JJ. Computational study of the conformational preferences of the (R)-8-amino-pentacyclo(5.4-0.0(2,6).0(3,10).0(5,9)) undecane-8-carboxylic acid monopeptide. J Pept Sci. 2004;10:274–84.

  23. 23.

    Gokul V, Kruger HG, Govender T, Fourie L, Power TD. An ab initio mechanistic understanding of the regioselective acetylation of 8,11-dihydroxy-pentacyclo[5.4.0.02,6.03,10.05,9] undecane-8,11-lactam. J Mol Struct: THEOCHEM. 2004;672:119–25.

  24. 24.

    Bisetty K, Corcho FJ, Canto J, Kruger HG, Perez JJ. Simulated annealing study of the pentacyclo-undecane cage amino acid tripeptides of the type Ac-X-Y-Z-NHMe. J Mol Struct: (THEOCHEM). 2006;759:145–57.

  25. 25.

    Bisetty K, Corcho FJ, Canto J, Kruger HG, Perez JJ. A theoretical study of pentacyclo-undecane cage peptides of the type (Ac-X-Y-NHMe). J Pept Sci. 2006;12:92–105.

  26. 26.

    Bisetty K, Corcho FJ, Canto J, Kruger HG, Perez JJ. A molecular dynamics study of the pentacyclo-undecane cage amino acid tripeptide. J Mol Struct: (THEOCHEM). 2006;770:221–8.

  27. 27.

    Bisetty K, Kruger HG, Perez JJ. A molecular dynamics study of the pentacyclo-undecane (PCU) cage polypeptides of the type Ac-3Ala-Cage-3Ala-NHMe. Mol Simul. 2007;33:1105–8.

  28. 28.

    Albericio F, Arvidson PI, Bisetty K, Giral E, Govender T, Jali S, et al. Trishomocubane amino acid as a beta-turn scaffold. Chem Biol Drug Des. 2008;71:125–30.

  29. 29.

    Pawar SA, Jabgunde AM, Govender P, Maguire GEM, Kruger HG, Parboosing R, et al. Synthesis and molecular modelling studies of novel carbapeptide analogs for inhibition of HIV-1 protease. Eur J Med Chem. 2012;53:13–21.

  30. 30.

    Singh T, Kruger HG, Bisetty K, Power TD. Theoretical study on the formation of a pentacyclo-undecane cage lactam. Comput Theor Chem. 2012;986:63–70.

  31. 31.

    Monajjemi M, Honarparvar B, Haeri HH, Heshmat M. An ab initio quantum chemical investigation of solvent-induced effect on N-14-NQR parameters of alanine, glycine, valine, and serine using a polarizable continuum model. Russ J Phys Chem. 2006;80:S40–4.

  32. 32.

    Mollaamin F, Baei MT, Monajjemi M, Zhiani R, Honarparvar B. A DFT study of hydrogen chemisorption on V (100) surfaces. Russ J Phys Chem A. 2008;82:2354–61.

  33. 33.

    Monajjemi M, Razavian MH, Mollaamin F, Naderi F, Honarparvar B. A theoretical thermochemical study of solute-solvent dielectric effects in the displacement of codon-anticodon base pairs. Russ J Phys Chem A. 2008;82:2277–85.

  34. 34.

    Monajjemi M, Honarparvar B, Nasseri SM, Khaleghian M. NQR and NMR study of hydrogen bonding interactions in anhydrous and monohydrated guanine cluster model: a computational study. J Struct Chem. 2009;50:67–77.

  35. 35.

    Monajjemi M, Lee VS, Khaleghian M, Honarparvar B, Mollaamin F. Theoretical description of electromagnetic nonbonded Interactions of radical, cationic, and anionic NH2BHNBHNH2 inside of the B18N18 nanoring. J Phys Chem C. 2010;114:15315–30.

  36. 36.

    Monajjemi M, Mahdavian L, Mollaamin F, Honarparvar B. Thermodynamic investigation of EnolKeto tautomerism for alcohol sensors based on carbon nanotubes as chemical sensors. Fullerenes Nanotubes Carbon Nanostruct. 2010;18:45–55.

  37. 37.

    Irani S, Monajjemi M, Honarparvar B, Atyabi SM, Sadeghizadeh M. Investigation of solvent effect and NMR shielding tensors of p53 tumor-suppressor gene in drug design. Int J Nanomedicine. 2011;6:213–8.

  38. 38.

    Monajjemi M, Khosravi M, Honarparvar B, Mollaamin F. Substituent and solvent effects on the structural bioactivity and anticancer characteristic of catechin as a bioactive constituent of green tea. Int J Quantum Chem. 2011;111:2771–7.

  39. 39.

    Sarasia EM, Afsharnezhad S, Honarparvar B, Mollaamin F, Monajjemi M. Theoretical study of solvent effect on NMR shielding tensors of luciferin derivatives. Phys Chem Liq. 2011;49:561–71.

  40. 40.

    Zuo Z, Gandhi NS, Mancera RL. Calculations of the free energy of interaction of the c-Fos-c-Jun coiled coil: effects of the solvation model and the inclusion of polarization effects. J Chem Inf Model. 2010;50:2201–12.

  41. 41.

    Jiao D, Zhang J, Duke RE, Li G, Schnieders MJ, Ren P. Trypsin-ligand binding free energies from explicit and implicit solvent simulations with polarizable potential. J Comput Chem. 2009;30:1701–11.

  42. 42.

    Vorontsov II, Miyashita O. Crystal molecular dynamics simulations to speed up MM/PB(GB)SA evaluation of binding free energies of di-mannose deoxy analogs with P51G-m4-Cyanovirin-N. J Comput Chem. 2011;32:1043–53.

  43. 43.

    Slynko I, Scharfe M, Rumpf T, Eib J, Metzger E, Schuele R, et al. Virtual screening of PRK1 inhibitors: ensemble docking, rescoring using binding free energy calculation and QSAR model development. J Chem Inf Model. 2014;54:138–50.

  44. 44.

    Latosinska JN. Structure-activity study of thiazides by magnetic resonance methods (NQR, NMR, EPR) and DFT calculations. J Mol Graph Model. 2005;23:329–37.

  45. 45.

    Latosinska JN. Applications of nuclear quadrupole resonance spectroscopy in drug development. Expert Opin Drug Discovery. 2007;2:225–48.

  46. 46.

    Cookson RC, Crundwell E, Hudec J. Synthesis of cage-like molecules by irradiation of Diels-Alder adducts. Chem Indust. 1958;32:1003–4.

  47. 47.

    Cookson RC, Hill RR, Hudec J. Stereochemistry of adducts of P-benzoquinone with 2 molecules of cyclopentadiene. Charge-transfer from olefinic double bonds to P-benzoquinone + ene-1.4-dione groups. J Chem Soc. 1964;0:3043–8.

  48. 48.

    Hemelaar J, Gouws E, Ghys PD, Osmanov S. Global and regional distribution of HIV-1 genetic subtypes and recombinants in 2004. AIDS. 2006;20:W13–23.

  49. 49.

    McCutchan FE. Global epidemiology of HIV. J Med Virol. 2006;78:S7–12.

  50. 50.

    Kantor R, Katzenstein D. Polymorphimsm in HIV-1 non-subtype B protease and reverse transcriptase and its potential impact on drug susceptibility and drug resistance evolution. AIDS Rev. 2003;5:25–35.

  51. 51.

    Los Alamos HIV sequence database. http://www.HIV.lanl.gov.

  52. 52.

    Mosebi S, Morris L, Dirr HW, Sayed Y. Active-site mutations in the South African human immunodeficiency virus type 1 subtype C protease have a significant impact on clinical inhibitor binding: kinetic and thermodynamic study. J Virol. 2008;82:11476–9.

  53. 53.

    Naicker P, Achilonu I, Fanucchi S, Fernandes M, Ibrahim MAA, Dirr HW, et al. Structural insights into the South African HIV-1 subtype C protease: impact of hinge region dynamics and flap flexibility in drug resistance. J Biomol Struct Dyn. 2013;31:1370–80.

  54. 54.

    Gaussian 03 RC, Frisch MJT, Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, et al. Gaussian, Inc., Wallingford CT. 2004.

  55. 55.

    Cornell WD, Cieplak P, Bayly CI, Kollman PA. Application of Resp charges to calculate conformational energies, hydrogen-bond energies, and free-energies of solvation. J Am Chem Soc. 1993;115:9620–31.

  56. 56.

    Bayly CI, Cieplak P, Cornell WD, Kollman PA. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges - the Resp model. J Phys Chem. 1993;97:10269–80.

  57. 57.

    Case DA, Cheatham TE, Darden T, Gohlke H, Luo R, Merz KM, et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26:1668–88.

  58. 58.

    Salomon-Ferrer R, Case DA, Walker RC. An overview of the Amber biomolecular simulation package. Wiley Interdiscip Rev: Comput Mol Sci. 2013;3:198–210.

  59. 59.

    Wang JM, Wolf RM, Caldwell JW, Kollman PA, Case DA. Development and testing of a general amber force field. J Comput Chem. 2004;25:1157–74.

  60. 60.

    Cornell WD, Cieplak P, Bayly CI, Gould IR, Merz KM, Ferguson DM, et al. A second generation force field for the simulation of proteins, nucleic acids, and organic molecules (vol 117, pg 5179, 1995). J Am Chem Soc. 1996;118:2309.

  61. 61.

    Hornak V, Abel R, Okur A, Strockbine B, Roitberg A, Simmerling C. Comparison of multiple amber force fields and development of improved protein backbone parameters. Proteins: Struct Funct Bioinf. 2006;65:712–25.

  62. 62.

    Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, et al. Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function. J Comput Chem. 1998;19:1639–62.

  63. 63.

    The PyMOL Molecular Graphics System. Version 1.5.0.4 Schrödinger, LLC.

  64. 64.

    Li H, Robertson AD, Jensen JH. Very fast empirical prediction and rationalization of protein pK(a) values. Proteins: Struct Funct Bioinf. 2005;61:704–21.

  65. 65.

    Shen C-H, Wang Y-F, Kovalevsky AY, Harrison RW, Weber IT. Amprenavir complexes with HIV-1 protease and its drug-resistant mutants altering hydrophobic clusters. FEBS J. 2010;277:3699–714.

  66. 66.

    Smith R, Brereton IM, Chai RY, Kent SBH. Ionization states of the catalytic residues in HIV-1 protease. Nat Struct Biol. 1996;3:946–50.

  67. 67.

    Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79:926–35.

  68. 68.

    Grest GS, Kremer K. Molecular-dynamics simulation for polymers in the presence of a heat bath. Phys Rev A. 1986;33:3628–31.

  69. 69.

    Toukmaji A, Sagui C, Board J, Darden T. Efficient particle-mesh Ewald based approach to fixed and induced dipolar interactions. J Chem Phys. 2000;113:10913–27.

  70. 70.

    Ryckaert JP, Ciccotti G, Berendsen HJC. Numerical-integration of cartesian equations of motion of a system with constraints - molecular-dynamics of N-alkanes. J Comput Phys. 1977;23:327–41.

  71. 71.

    Hou TJ, Guo SL, Xu XJ. Predictions of binding of a diverse set of ligands to gelatinase-A by a combination of molecular dynamics and continuum solvent models. J Chem Phys B. 2002;106:5527–35.

  72. 72.

    Bonnet P, Bryce RA. Scoring binding affinity of multiple ligands using implicit solvent and a single molecular dynamics trajectory: application to influenza neuraminidase. J Mol Graph Model. 2005;24:147–56.

  73. 73.

    Hou T, Wang J, Li Y, Wang W. Assessing the performance of the MM/PBSA and MM/GBSA Methods. 1. The accuracy of binding free energy calculations based on molecular dynamics simulations. J Chem Inf Model. 2011;51:69–82.

  74. 74.

    Kollman PA, Massova I, Reyes C, Kuhn B, Huo SH, Chong L, et al. Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc Chem Res. 2000;33:889–97.

  75. 75.

    Miller III BR, McGee Jr TD, Swails JM, Homeyer N, Gohlke H, Roitberg AE. MMPBSA.py: an efficient program for End-state free energy calculations. J Chem Theory Comput. 2012;8:3314–21.

  76. 76.

    Srinivasan J, Cheatham TE, Cieplak P, Kollman PA, Case DA. Continuum solvent studies of the stability of DNA, RNA, and phosphoramidate - DNA helices. J Am Chem Soc. 1998;120:9401–9.

  77. 77.

    Rastelli G, Del Rio A, Degliesposti G, Sgobba M. Fast and accurate predictions of binding free energies using MM-PBSA and MM-GBSA. J Comput Chem. 2010;31:797–810.

  78. 78.

    Gaussian 09 RA, Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, et al. Gaussian, Inc., Wallingford CT. 2009.

  79. 79.

    Becke AD. Density-functional thermochemistry.3. The role of exact exchange. J Chem Phys. 1993;98:5648–52.

  80. 80.

    Lee CT, Yang WT, Parr RG. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron-density. Phys Rev B. 1988;37:785–9.

  81. 81.

    Witanowski M, Sicinska W, Biedrzycka Z, Webb GA. Hydrogen-bond and solvent polarity effects on the nitrogen NMR shielding of urea systems. J Mol Struct. 1999;476:133–8.

  82. 82.

    Chai J-D, Head-Gordon M. Systematic optimization of long-range corrected hybrid density functionals. J Chem Phys. 2008;128:084106.

  83. 83.

    Reed AE, Weinhold F. Natural bond orbital analysis of near-Hartree-Fock water dimer. J Chem Phys. 1983;78:4066–73.

  84. 84.

    Reed AE, Weinstock RB, Weinhold F. Natural-population analysis. J Chem Phys. 1985;83:735–46.

  85. 85.

    Reed AE, Weinhold F. Natural bond orbital analysis of internal-rotation barriers and related phenomena. Israel J Chem. 1991;31:277–85.

  86. 86.

    Flurchick K, Bartolotti L. Visualizing properties of atomic and molecular-systems. J Mol Graph. 1995;13:10–3.

  87. 87.

    Bartolotti LJ, Ayers PW. An example where orbital relaxation is an important contribution to the Fukui function. J Phys Chem A. 2005;109:1146–51.

  88. 88.

    Lucken EAC. Nuclear Quadrupole Coupling Constants. London: Academic Press; 1992.

  89. 89.

    Pyykko P. Spectroscopic nuclear quadrupole moments. Mol Phys. 2001;99:1617–29.

  90. 90.

    Curutchet C, Orozco M, Luque FJ, Mennucci B, Tomasi J. Dispersion and repulsion contributions to the solvation free energy: comparison of quantum mechanical and classical approaches in the polarizable continuum model. J Comput Chem. 2006;27:1769–80.

  91. 91.

    Marenich AV, Cramer CJ, Truhlar DG. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J Phys Chem B. 2009;113:6378–96.

  92. 92.

    Stoica I, Sadiq SK, Coveney PV. Rapid and accurate prediction of binding free energies for saquinavir-bound HIV-1 proteases. J Am Chem Soc. 2008;130:2639–48.

  93. 93.

    Portmann S, Luthi HP. MOLEKEL: an interactive molecular graphics tool. Chimia. 2000;54:766–70.

  94. 94.

    Brik A, Wong CH. HIV-1 protease: mechanism and drug discovery. Org Biomol Chem. 2003;1:5–14.

  95. 95.

    Gready JE. The relationship between nuclear-quadrupole coupling-constants and the asymmetry parameter - The interplay of theory and experiment. J Am Chem Soc. 1981;103:3682–91.

  96. 96.

    Schnarre R. Drug Design. New York: Academic Press; 1971.

  97. 97.

    Wohl A. Drug Design. New York: Academic Press; 1971.

  98. 98.

    Kier L. Orbital Studies in Chemical Pharmacology. New York: Springer; 1997.

Download references

Acknowledgements

We thank the National Research Foundation, South Africa for financial support, UKZN and the CHPC (www.chpc.ac.za) for computational resources.

Author information

Correspondence to Hendrik G Kruger.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SAP performed the experimental experiments under supervision of TG, GEMM and HGK. BH executed the computational work with assistance by JRAS and with support by CAN, JL, GEMM and HGK. All authors read and approved the final manuscript.

Additional file

Additional file 1:

Supporting Information.

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Pentacyclo-undecane (PCU)
  • Lactone
  • Lactam and lactim peptides
  • Human immunodeficiency virus protease (HIV PR)
  • Inhibitory concentration (IC50)
  • Molecular dynamics (MD)
  • MM-PB (GB) SA method