Zahra Beiranvanda,
Farhad Banib,
Ali Kakanejadifarda,
Erik Laurinic,
Maurizio Fermegliac,
Sabrina Pricl*c and
Mohsen Adeli*a
aFaculty of Science, Department of Chemistry, Lorestan University, Khorramabad, Iran. E-mail: mohadeli@yahoo.com
bInstitue of Biochemistry and Biophysic, University of Tehran, Tehran, Iran
cMolecular Simulation Engineering (MOSE) Laboratory, DEA, University of Trieste, Piazzale Europa 1, 34127 Trieste, Italy. E-mail: Sabrina.pricl@di3.units.it
First published on 13th January 2016
Here, we report for the first time a facile method for the preparation of novel drug delivery systems based on supramolecular interactions between β-cyclodextrin-polyglycerol conjugates (β-CD-g-PG) and human serum albumin (HSA). The results obtained by combined experimental/modeling studies showed that the main drivers to the formation of HSA/β-CD-g-PG supramolecular entities are host–guest interactions between β-CD-g-PG and the aromatic side-chains of HSA residues. Due to these interactions, HSA undergoes a confined yet fundamental conformational transition, leading to greater exposure of the hydrophobic protein domains available to hydrophobic drug binding. Next, the binding affinity and loading capacity of the HSA/β-CD-g-PG supramolecular nanovector for doxorubicin (DOX) and paclitaxel (PTX) were investigated. Both drug/HSA binding constants (Kb,HSA/DOX = 3.24 × 103 M−1 and Kb,HSA/PTX = 5.65 × 101 M−1) sensibly increased in the presence of β-CD-g-PG (Kb,HSA/β-CD-g-PG/DOX = 2.78 × 104 M−1 and Kb,HSA/β-CD-g-PG/PTX = 2.82 × 102 M−1). In line, both protein drug loadings increased by about 20% upon HSA/β-CD-g-PG interaction (80% and 71% for DOX and PTX, respectively) with respect to the loading capacity of the bare HSA (60% and 50%). Due to the improved loading capacity with minimal changes in the structure of HSA, this system is a promising vector for future cancer therapy.
Human serum albumin (HSA) plays a key role in maintaining blood osmotic pressure and in acting as one of the main transporter of both endogenous and exogenous ligands, such as bilirubin, fatty acids, ions, and drugs.4 The tertiary structure of this mainly-helical protein resembles a heart shape, with approximate dimensions of 8 × 8 × 3 nm, and consists of three domains: domain I (residues 1–195), II (residues 196–383), and III (384–585). Each domain, in turn, features two subdomains (A and B), and 17 disulphide bridges stabilize the whole protein structure. As a transport protein, HSA possesses several binding sites for its ligands.5 Specifically, aromatic and heterocyclic ligands are found to bind within two hydrophobic pockets in subdomains IIA and IIIA, namely site I and site II.5,6 Contextually, seven binding sites for fatty acids are localized in subdomains IB, IIIA, IIIB and on the subdomain interfaces.5 Taken together, HSA possesses a unique set of properties, which make it well suited for development as drug nanocarrier.
On the other hand, hyperbranched polyglycerols are water soluble, biocompatible, and multivalent materials with many potential applications ranging from drug delivery to catalysis.7 Specifically, the protein-resistance property of polyglycerol is one of the most important feature of this polymer in biomedical applications. Due to this property, nonspecific interactions between polyglycerol and proteins and, consequently, its immunogenicity are very low.8 Cyclodextrins (CDs) are torus-shaped cyclic oligosaccharides made of six, seven, or eight glycosidic units linked by α(1 → 4) bonds and are called α-, β-, and γ-CD, respectively. This molecular structure has the capability to make inclusions (i.e., guest–host) compounds with aliphatic and especially aromatic molecules. Cyclodextrins are used in pharmaceutical applications for numerous purposes, including improving the stability, bioavailability, and delivery of drugs.
Accordingly, in this work we conceived the idea of exploiting all beneficial properties of both polyglycerol and β-cyclodextrin by creating cyclodextrin-polyglycerol conjugates (β-CD-g-PG) having special site for specific interactions.9 In fact, due to the particular supramolecular chemistry for each subdomain of albumin, it interacts with different objects specifically.10 As a result, we expect that β-CD-g-PG conjugates will interact with albumin mainly via cyclodextrin and host–guest interactions. Under this hypothesis, these specific interactions between albumin and β-CD-g-PG should lead to a hybrid system with useful properties originating from both materials.
Thus, aim of the present study is to verify these assertions and demonstrate the potentialities of these novel HSA/β-CD-g-PG host–guest systems for delivery and release of anticancer drugs. Our molecular design is indeed oriented toward the development of low-cost, easy to fabricate, easy to use and robust drug delivery systems. Under this perspective, we designed the present biohybrid-nanomaterials, which can be even produced by in situ administration of the precursors. All supramolecular anticancer drug delivery systems synthesized in this work were also characterized and investigated by various spectroscopic techniques including fluorescence, dynamic light scattering (DLS), ultraviolet-visible (UV-vis), circular dichroism (CD), and Fourier transform infrared (FT-IR) spectroscopy. Molecular simulations were also employed to shed light on the possible mechanism of interactions between all components of these novel supramolecular drug delivery systems. Further, we examined the interaction of two hydrophobic anticancer drugs – doxorubicin and paclitaxel – to gain further information on the effects of drug molecular structure on the HSA/β-CD-g-PG/drug complex formation. Based on the overall experimental and simulation evidences obtained, the present results may provide new insights on the rational design of natural supramolecular drug carriers. To the best of our knowledge, this is the first report of host–guest interactions between albumin and β-CD-g-PG to produce a novel and effective anticancer drug delivery system.
For the preparation of the HSA/β-CD-g-PG supramolecular entities (Fig. 1), β-CD-g-PG was added to a phosphate buffered saline (PBS) solution of HSA, and the resulting mixture was gently shaken for 1 h (Fig. 1). Host–guest interactions between HSA and β-CD-g-PG were investigated using a combination of different spectroscopy methods and molecular simulation techniques.
![]() | ||
Fig. 1 Host–guest interactions between HSA and β-CD-g-PG leading to HSA/β-CD-g-PG supramolecular entities. |
In agreement with these confined spectral perturbations, the calculated α-helical content of HSA decreases only slightly, i.e., from 59.2% for pristine HSA to 52.5% upon formation of the HSA/β-CD-g-PG supramolecular entity (Table 1). Correspondingly, also the other motifs characterizing the secondary structure of HSA manifest only minor changes in their values. Notably, HSA is a dynamic, rapidly moving, “breathing” molecule having a rotational diffusion coefficient of 20 ns (ref. 4) and it is known to have some degree of flexibility, expansion, and contraction, which enables significant conformational transitions; hence, the precise, average shape that HSA takes on in aqueous solution is still under debate.13,14 Therefore, the structural perturbations exerted on HSA by β-CD-g-PG upon binding are indeed negligible.
Structure | HSA | HSA-β-CD-g-PG |
---|---|---|
α-Helix | 59.2% | 52.5% |
β-Antiparallel | 3.8% | 4.7% |
β-Sheet | 4.6% | 5.4% |
β-Turn | 12.9% | 13.8% |
Random coil | 19.5% | 23.5% |
Generally, fluorescence of proteins originates from three residues namely tryptophan, tyrosine, and phenylalanine. Actually, the intrinsic fluorescence of HSA is mainly contributed by tryptophan 214 alone (located in the protein subdomain IIA), whereas the emission of phenylalanine and tyrosine is negligible due to their very low quantum yield.16 The fluorescence of Trp-214 may change when HSA interacts with other molecules, which could be reflected in its fluorescence spectra in the UV region. Here, the fluorescence quenching of HSA in the presence of various concentrations of β-CD-g-PG was measured to shed some light in the quenching mechanism and the relevant, underlying intermolecular interactions.
As shown in Fig. 3, the intensity of the characteristic broad emission band of HSA at 345 nm regularly decreases with increasing β-CD-g-PG, concentration, with a detectable bathochromic shift of the emission maximum of HSA by 5 nm. Accordingly, the quenching and the red shift in the HSA emission maxima not only confirm the onset of intramolecular interactions between β-CD-g-PG and the protein, but also imply that β-CD-g-PG is introducing an increase in the polarity of the micro-environment around the Trp-214 moiety. Eventual dilution effect upon addition of β-CD-g-PG solutions can be ruled out, because the volume of solution containing HSA (4 mL) is far larger than the β-CD-g-PG solution (≤200 μL).
F0/F = 1 + kqτ0[Q] = 1 + KSV[Q] | (1) |
![]() | ||
Fig. 4 Stern–Volmer plots for the quenching of HSA by β-CD-g-PG in PBS solution at 37 °C (A) and at different temperatures (B). [HSA] = 2 × 10−6 M. |
As can be seen in Fig. 4A, the value for kq calculated using eqn (1) is equal to 7.65 × 1012 M−1 s−1. Since this value is much larger than the maximum collisional quenching constant (2.06 × 1010 M−1 s−1), SQ is the dominant quenching mechanism for the HSA/β-CD-g-PG supramolecular complex.21,22 To further verify this assertion, the dependence of kq on temperature was investigated. Accordingly, the Stern–Volmer plots for the quenching of HSA by β-CD-g-PG at three different temperatures (i.e., 20, 30, and 37 °C) are displayed in Fig. 4B. The calculated values of the quenching constant at 20 and 30 °C are 1.35 × 1013 M−1 s−1 and 8.91 × 1012 M−1 s−1, respectively. The inverse correlation of kq with temperature supports the evidence that the fluorescence quenching of HSA by β-CD-g-PG was initiated by complex formation between HSA and β-CD-g-PG rather than by dynamic collision between the two macromolecules.23,24
For stating quenching and under the assumption that there are n equivalent and independent binding sites for the quencher (Q) on the protein, the binding constant Kb for β-CD-g-PG onto HSA can be calculated according to the modified Stern–Volmer equation:16,25
log[(F0 − F)/F] = log![]() ![]() | (2) |
![]() | ||
Fig. 5 Modified Stern–Volmer plot for the quenching of HSA by β-CD-g-PG in PBS solution at different temperatures. [HSA] = 2 × 10−6 M. |
T (°C) | Kb (M−1) | n | R2 | ΔGb (kcal mol−1) |
---|---|---|---|---|
20 | 879 | 1.14 | 0.997 | −3.95 |
30 | 331 | 1.05 | 0.975 | −3.49 |
37 | 244 | 1.02 | 0.997 | −3.39 |
It can be seen from Table 2 that the supramolecular complexes between HSA and β-CD-g-PG can be easily formed at all temperatures; also, the number of binding sites n is always close to 1, indicating that there is only 1 β-CD-g-PG binding site on the protein surface. Interestingly, the calculated Kb values in Table 2 are in excellent agreement with the Kb of 185 M−1 recently reported by Gosh et al.26 for the binding of unmodified β-cyclodextrin (β-CD) to HSA at room temperature. Moreover, the highest affinity of β-CD-g-PG towards the protein with respect to the unmodified β-CD suggests a beneficial effect of the polyglycerol substituent in promoting the formation of the HSA/β-CD-g-PG supramolecular assembly.
ΔGb = ΔHb − TΔSb = −RT![]() ![]() | (3) |
![]() | (4) |
Thus, choosing T1 = 293 K (20 °C) and T2 = 310 K (37 °C), and inserting the corresponding values of Kb,1 and Kb,2 (Table 2) into eqn (4), the relevant value of ΔHb can be estimated. Also, using the corresponding ΔGb values in Table 2, the entropic component TΔSb can be readily obtained from eqn (3) as:
TΔSb = ΔHb − ΔGb | (5) |
The first three columns of Table 3 show the results of this analysis. As can be seen, both the enthalpic and the entropic component of the free energy of binding ΔGb of β-CD-g-PG onto HSA are negative, implying that the binding process between the modified β-CD and the protein is predominantly enthalpic, while the entropic variation opposes binding, as most often seen for ligand/protein complexation.28 Also, the equally negative sign of ΔHb and TΔSb suggests that van der Waals, hydrophobic, and hydrogen bond intermolecular interactions are prevalently leading the HSA/β-CD-g-PG binding process, as it could be expected given the nature of the receptor and its ligand (vide infra).
T (°C) | ΔHb (kcal mol−1) | TΔSb (kcal mol−1) | ΔGb (kcal mol−1) |
---|---|---|---|
20 | −13.61 | −9.66 | −3.95 |
30 | −10.12 | −3.49 | |
37 | −10.22 | −3.39 |
This is in harmony with the small structural perturbation on the overall secondary structure of HSA exerted by β-CD-g-PG upon binding. Contextually, the second, broad peak also increases with a progressive slight shift towards lower wavelengths, suggesting a moderate change in the microenvironment around Trp, in keeping with the results obtained from fluorescence quenching evidences. The entity of these spectral variations is another confirmation of the non-covalent nature of intermolecular interactions underlying the complex formation between β-CD-g-PG and HSA.
One additional method to distinguish static and dynamic quenching is by careful examination of the absorption spectra of the fluorophore. Collisional quenching (or DQ) only affects the excited states of the fluorophore and thus no changes in the absorption spectra are expected. In contrast, ground-state complex formation (i.e., SQ) will frequently result in perturbation of the absorption spectrum of the fluorophore.23 To reconfirm that quenching mechanism of fluorescence of HSA by β-CD-g-PG is initiated by ground-state complex formation, we performed UV-visibile (UV-vis) adsorption experiments of HSA, β-CD-g-PG, and HSA/β-CD-g-PG system (Fig. 7). As seen in Fig. 7, the intensity of both HSA peaks (at 218 nm and 278 nm) slightly increases upon addition of β-CD-g-PG, suggesting the existence of an HSA/β-CD-g-PG ground state complex and, therefore, the static mechanism of HSA fluorescence quenching by the polyglycerol-modified β-CD. Concomitantly, the invariance of the maximum absorption wavelengths implies that the interaction between β-CD-g-PG and HSA is non-covalent in nature.
The value of the binding constant Kb for the formation of the supramolecular complex between HSA and β-CD-g-PG obtained from fluorescence quenching experiments (eqn (2)) can be validated using UV-vis data in Fig. 6 (see ESI and Fig. 6S† for details). Considering the different sensitivity and approximations underlying the two spectroscopic techniques, the corresponding Kb value of 419 M−1 is in agreement with the one obtained by fluorescence spectroscopy (Kb = 244 M−1, Table 2), thereby confirming the moderate strength of the intermolecular, non-covalent interactions between HSA and β-CD-g-PG.
Compound | dh (nm) | ζ (mV) |
---|---|---|
HSA | 6 | −7.27 |
β-CD-g-PG | 3 | −1.23 |
HSA-β-CD-g-PG | 11 | −9.40 |
As seen, the hydrodynamic diameter of native HSA was found to be 6 nm at pH = 7.4 and at 37 °C, in good agreement with previous literature reports.31 With the addition of β-CD-g-PG, the dh value of the resulting complex increases to 11 nm, a value somewhat larger than the sum of the hydrodynamic diameters of the pure components. Given the small structural change of HSA induced by β-CD-g-PG upon binding, the increase in the HSA/β-CD-g-PG supramolecular complex dimension cannot therefore be ascribe to a substantial unfolding of HSA but, rather, to an extended conformation of the polyglycerol chains of β-CD-g-PG into the solvated environment. This result, as discussed later in detail, could be an expected result since the polyglycerol derivative of β-CD can be though of as composed of two molecular moieties: the relatively hydrophobic cyclodextrin cavity and a the polar outer shell, whose hydrophilicity is enhanced by the presence of the PG chains. Accordingly, the former tend to interact with HSA, thereby forming the necessary stabilizing network of hydrophobic, van der Waals and hydrogen bond interactions with the residues lining the binding pocket of HSA. At the same time, the latter can potentiate the binding by shielding these interactions from water.
The zeta potential of a protein is correlated to its electric surface properties, which, in turn, are generated from presence of ionized residues (e.g., aspartic and glutamic acid, arginine and lysine). As such, the effective charge and, hence, the value of ζ on a given protein is affected by pH, ionic strength, and the eventual accumulation of ligand on the protein surface. At pH 7.4, the three homologous, globular domains of native HSA in solution carry a net charge of −9e (domain I, N-terminal), −8e (domain II), and +2e (domain III, C-terminal), giving a total net charge for the protein of −15e (when calculated from the amino acid composition) or −19e if measured experimentally (due to additional bound ions) at pH 7.4.32 In line, the present value of ζ measured for HSA is negative, and equal to −7.27 mV (Table 4). The addition of β-CD-g-PG (ζ = −1.23 mV) results in a ζ value for the supramolecular complex of −9.40 mV. The decrease in the ζ-value upon complex formation supports the main dispersive nature (i.e., hydrophobic, hydrogen bond and van der Waals interactions) of the intermolecular forces leading to the HSA-β-CD-g-PG supramolecular assembly.
Next, by applying the well-validated MM/PBSA (Molecular Mechanics/Poisson Boltzmann Surface Area) approach33 the free energy of binding (ΔGb) between the protein and the cyclodextrin derivative was calculated (see ESI†). The relevant ΔGb value of −3.57 kcal mol−1 at T = 37 °C, in excellent agreement with the corresponding experimental value (Table 3) confirms that β-CD-g-PG binds HSA with a moderate affinity. The deconvolution of ΔGb into its components (see ESI† for details) further reveals that the net electrostatic contribution to the affinity of the complex is significantly unfavorable. Specifically, the mean value of the total electrostatic energy contribution (ΔGele = ΔEele + ΔGPB) is positive and amounts to 22.63 kcal mol−1, while the corresponding value of the dispersive forces (i.e., van der Waals and hydrophobic interaction energies, ΔGdisp = ΔEvdW + ΔGNP) is remarkably favorable, being equal to −36.55 kcal mol−1.
Accordingly, we can conclude that the overall enthalpic contribution governing the spontaneous binding of β-CD-g-PG onto HSA is dominated by dispersive interactions, in agreement with experiments. Performing a per residue decomposition of the enthalpic component of ΔGb shows that the β-CD moiety of β-CD-g-PG affords the most effective contribution to HSA binding, while the polyglycerol branches seem to be almost ineffective. To dissect this aspect in more detail, we performed MD simulation of the complex between the unmodified β-CD and the HSA under the same conditions (Fig. 8B and C). The calculated free energy of binding value for the HSA/β-CD complex, ΔGb = −2.04 kcal mol−1, confirms that the unmodified CD has less affinity for HSA than its PG-modified counterpart (ΔGb = −3.57 kcal mol−1). As can be further inferred from the results displayed in Fig. 14D, the entropic component of the ΔGb, being negative, opposes binding in both HSA/β-CD-g-PG and HSA/β-CD complexes (TΔSb HSA/β-CD-g-PG = −10.35 kcal mol−1 and TΔSb and HSA/β-CD = −6.56 kcal mol−1, respectively). At the same time, however, the unmodified cyclodextrin establishes weaker intermolecular interactions with the HSA since, for this system, non only the overall dispersive component of ΔGb is less favorable (ΔGdisp = −33.77 kcal mol−1) but, at the same time, the β-CD/HSA complex pays higher molecular desolvation penalty upon binding (ΔGPB = 42.16 kcal mol−1) with respect to the polyglycerol-modified β-CD counterpart. In other words, the contribution of the PG branches translates in two meaningful positive effect for the binding with HSA by: (1) rendering the β-CD cavity more hydrophobic thereby improving the interactions with the HSA residues; and (2) providing a hydrophilic shielding of the binding region by exposing the polar moieties to the solvent with a subsequent lower desolvation energy paid by β-CD-g-PG upon binding with HSA.
However, the ellipticity of the protein is significantly reduced compared to that of the pristine HSA in the 200–260 nm wavelength range. Elaboration of CD data indeed reveals that, in the presence of either DOX or PTX, the α-helical content of HSA plummets from 52.5% for HSA/β-CD-g-PG/HSA to 44.2% for HSA/β-CD-g-PG/DOX and 42.1% for HSA/β-CD-g-PG/PTX, respectively (Table 5). The decrease in α-helix content was accompanied by an increase in the β-sheet, β-turn and β-antiparallel motifs; similarly, the unstructured (random coil) portion of the protein increases in both drug complexes (Table 5). This may suggest that both DOX and PTX are able to interact with the residues of the polypeptide chain of HSA, partially destroying the hydrogen bonding networks, thereby evoking some degree of protein destabilization, the adoption of a more open conformation ultimately resulting in better exposure of the hydrophobic cavity available to DOX/PTX binding.
Structure | HSA/β-CD-g-PG/DOX | HSA/β-CD-g-PG/PTX |
---|---|---|
α-Helix | 44.2% | 42.1% |
β-Antiparallel | 6.9% | 6.7% |
β-Sheet | 6.7% | 6.9% |
β-Turn | 15.4% | 15.5% |
Random coil | 26.8% | 28.8% |
Next, the modified Stern–Volmer eqn (2) was used for the determination of the binding constant Kb and number of binding sites n of DOX and PTX onto the HSA/β-CD-g-PG complex. The same calculations were performed also in the case of pristine HSA/drug binding, for comparison. Fig. 11 shows the relevant results, which are numerically summarized in Table 6.
System | Kb (M−1) | n | R2 | ΔGb (kcal mol−1) |
---|---|---|---|---|
HSA/β-CD-g-PG/DOX | 18![]() |
1.03 | 0.993 | −6.07 |
HSA/β-CD-g-PG/PTX | 812 | 0.83 | 0.984 | −4.13 |
HSA/DOX | 3236 | 0.83 | 0.993 | −4.98 |
HSA/PTX | 176 | 0.76 | 0.998 | −3.18 |
The values of Kb listed in Table 6 clearly show the beneficial effect of the presence of β-CD-g-PG on the affinity of both anticancer drugs for HSA. Moreover, DOX is endowed with a more favorable free energy of binding ΔGb with respect to PTX, both towards native HSA and its β-CD-g-PG supramolecular complex, by virtue of a more efficient network of hydrophobic and hydrogen bond intermolecular interactions, as discussed in details in the molecular modeling section below. Of note, the Kb values for HSA/DOX (3256 M−1) and HSA/PTX (176 M−1) derived from the present fluorescence quenching data are in good agreement with previous reports.19,36
Compound | dh (nm) | ζ (mV) |
---|---|---|
HSA/DOX | 22 | −22.47 |
HSA/PTX | 25 | −18.53 |
HSA/β-CD-g-PG/DOX | 16 | −26.96 |
HSA/β-CD-g-PG/PTX | 20 | −24.56 |
These data reveal that the hydrodynamic diameter of HSA increased from 6 to 22 and 25 nm in the presence of DOX and PTX, respectively. On this basis, it can be inferred that, upon drug binding, the resulting HSA/drug complexes aggregate to a certain extent. This inference is confirmed through studies of the size distribution of HSA (Fig. 7S†) and the HSA/drug complexes (Fig. 8S†).
Whilst an almost monodisperse solution of HSA in present in the PSB buffer (pH = 7.4 and 37 °C) (6 nm), the addition of either DOX or PTX causes the formation of aggregates with size increasing to 22 and 25 nm, respectively. These results imply that the binding between HSA and DOX and PTX leads to a conformational change in HSA (Table 5), which contributes to the aggregation of the resulting HSA/drug complexes. The presence of β-CD-g-PG bound on the surface of HSA, on the other hand, seems to disfavor the formation of higher aggregates, likely by hampering optimal protein–protein and/or drug–drug interactions.
The zeta-potential values at physiological pH of the HSA/β-CD-g-PG/DOX and HSA/β-CD-g-PG/PTX complexes are highly negative and approximately equal to −27 and −25 mV, respectively (Table 7). These ζ values confirm the notion that either drug, upon binding onto the HSA/β-CD-g-PG assembly, provokes a perturbation in the secondary structure of the protein (Table 5), which ultimately results in the exposition of negatively charged amino acid side chains on the protein surface. The increased negative protein surface charge – and hence the related enhanced Coulombic repulsions – accounts for the smaller dimensions of the HSA/β-CD-g-PG/DOX and HSA/β-CD-g-PG/PTX bigger aggregates with respect to those of the binary (i.e., drug/protein) complexes (Table 7). This property, coupled with the likely further beneficial effect of intrinsic protein resistance and steric impediment to aggregation exerted by the presence of β-CD-g-PG on the surface of HSA, can be extremely advantageous when considering the polyglycerol-modified β-cyclodextrin–human serum albumin supramolecular assembly as a possible nanovector for drug delivery. In addition, negative surface charge of drug delivery systems will affect their uptake which is a significant factor and should be investigated experimentally.
Finally, the in vitro release of DOX and PTX from the HSA/β-CD-g-PG supramolecular assembly was investigated; the relevant profiles are shown in Fig. 12. It can be seen that the rate of release of both drugs from the modified HSA strongly depends on the pH, the % of drug released being more favorable at acidic pH in both cases. This is an expected result, as it is known that HSA undergoes pH-induced further secondary conformation unfolding which, in turn, might weaken drug binding and, hence, facilitate drug release from its binding site. Contextually, this is a favorable evidence since tumor microenvironment is usually acidic;37 accordingly, all other conditions being equal, the release of both anticancer drugs from the HSA/β-CD-g-PG will be at the tumor site enhanced with respect to the considerable lower one attained at neutral pH.
The drug/protein intermolecular interactions were quantified by carrying out MM/PBSA MD simulations. As shown in Fig. 13B and D, the two anticancer agents established a similar network of interactions with the protein. Specifically, both drug binding is driven by a combination of a weak array of hydrogen bonds (HBs) and hydrophobic forces. In particular, the side chains of the hydrophobic residues V140, M147, Y185, F189, L206, and L209 lean the binding pocket of DOX; more, we confirmed the establishment of HBs among several of the drug heteroatoms both with the polar side chains of the HSA residues R141, H170, R210, and K214 and the backbone NH group of residues R138 and L139 (Fig. 13B and Table 8). At variance, MD simulations carried out with PTX revealed the presence of hydrophobic interactions with the side chains of W238, L262, I314, H312, and A315 while some donor and acceptor atoms of the molecule found their appropriate counterpart in the polar amino acid residues K219, K223, R242, R246, and H266 (Fig. 13D and Table 4). The overall contribution of these favorable interactions calculated by MM/PBSA approach led to ΔGbind of −4.82 kcal mol−1 for HSA/DOX complex and −3.31 kcal mol−1 for HSA/PTX system, respectively. These data confirm the moderate affinity of the two drugs for HSA, in agreement with previously reported results34,35 and with our experimental analysis (see Table 6).
Complex | HSA residue | ADL (Å) |
---|---|---|
HSA/DOX | R138-L139 (bb) | 2.76 ± 0.06 |
R141 | 2.87 ± 0.03 | |
H170 | 3.01 ± 0.03 | |
R210 | 2.67 ± 0.04 | |
K214 | 2.73 ± 0.02 | |
HSA/β-CD-g-PG/DOX | R138-L139 (bb) | 2.66 ± 0.05 |
R141 | 2.49 ± 0.02 | |
H170 | 2.63 ± 0.04 | |
R210 | 2.38 ± 0.05 | |
K214 | 2.25 ± 0.03 | |
HSA/PTX | K219 | 2.93 ± 0.02 |
K223 | 2.99 ± 0.02 | |
R242 | 2.88 ± 0.03 | |
R246 | 3.15 ± 0.04 | |
H266 | 3.19 ± 0.04 | |
HSA/β-CD-g-PG/PTX | K219 | 2.68 ± 0.04 |
K223 | 2.70 ± 0.02 | |
R242 | 2.65 ± 0.05 | |
R246 | 2.92 ± 0.03 | |
H266 | 2.83 ± 0.04 |
The deconvolution of the enthalpic component of ΔGb was next performed for both systems with the aim of quantifying the single contribution of each amino acid involved in the drug binding. The collected results not only substantiate the fact that DOX is a stronger binding of HSA with respect to PTX, as more residues are involved in more efficient protein/drug binding, but also the fact that network of interactions between the protein and the drugs is overall rather weak. Indeed, no residues exhibited an enthalpic gain upon drug binding greater than 0.6 kcal mol−1 (Fig. 14).
![]() | ||
Fig. 14 Comparison of per-residue binding enthalpy decomposition for DOX (A) and PTX (B) in complex with HSA and HSA/β-CD-g-PG. |
In the previous modeling section we described how HSA and β-CD-g-PG interact through the formation of a thermodynamically stable complex. However, to understand the effect of the presence of β-CD-g-PG onto the HSA surface on drug binding, we then carried out MD simulations of the supramolecular HSA/β-CD-g-PG system in complex with DOX and PTX, respectively. Starting from the equilibrated MD structure of the β-CD-g-PG/HSA complex, DOX and PTX were docked in their corresponding HSA binding sites (Fig. 15A and C). Next, we performed extensive MD simulation runs with the aim of checking whether the structural HSA modifications promoted by β-CD-g-PG could eventually affect the binding mode of DOX and PTX onto the protein. Finally, MM/PBSA calculations were applied to estimate the affinity (ΔGb) of DOX and PTX for the supramolecular HSA/β-CD-g-PG complex. In agreement with the experimental results, the calculated affinity against β-CD-g-PG/HSA for both drugs increases, as the corresponding estimated ΔGb values are −6.32 kcal mol−1 for DOX and −4.35 kcal mol−1 for PTX, respectively. The conformational changes induced by the presence of the polyglycerol conjugate indeed perturb both protein drug binding regions, with consequent alteration of the relevant drug binding modes (Fig. 15B and D). This rearrangement reflects into a stronger network of the binding interactions between HSA and DOX and PTX. Indeed, all detected H-bonds became shorter, as revealed by the calculated ADL in Table 8 and, hence, stronger. Furthermore, the per residue enthalpic deconvolution for the HSA/β-CD-g-PG/drug complexes clearly emphasize hydrophobic forces are enhanced for these systems since the involved HSA residues are characterized by a more favorable enthalpic gain upon the binding with the drugs (Fig. 14).
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c5ra25463a |
This journal is © The Royal Society of Chemistry 2016 |