Abstract

In this research, the pyrazoline pyridazine derivative 7-methyl-2-phenyl-4-(3,4,5-trimethoxyphenyl)-2H-pyrazolo[3,4-d]pyridazine (5d) was studied for its interaction with bovine serum albumin (BSA). Various spectroscopic techniques along with molecular docking analysis were utilized to understand the mechanism of interaction. The quenching of BSA fluorescence by using investigational drug 5d was the basic principle for the methodology. Spectrofluorometric methods and UV-absorption studies were conducted for exploration of the 5d and BSA binding mechanism. The fluorescence quenching mechanism involved in BSA and 5d interaction was static quenching, and a complex formation also occurred between them. Both enthalpy and entropy attained positive values suggesting involvement of hydrophobic forces in BSA and 5d interaction. The Förster distance of 2.23 nm was calculated by fluorescence resonance energy transfer (FRET). An alteration in BSA secondary structure was proven from the conformational studies of BSA-5d interaction. This binding interaction study provided a basis to comprehend the binding interaction between 5d and BSA. These results provided information about sites of BSA involved in its interaction with 5d.

1. Introduction

Prostaglandin biosynthesis (cyclooxygenases; COXs) are inhibited by nonsteroidal anti-inflammatory drugs [1, 2]. Arachidonic acid is broken down by the COX enzymes to produce prostaglandins which have a number of important physiological functions. The main use of nonsteroidal anti-inflammatory drugs (NSAIDs) is to alleviate the pain and inflammation due to certain diseased states. COX-1 isoform is mainly induced under all physiological conditions, whereas in carcinogenesis and inflammatory responses, COX-2 is induced [3]. The synthesis of COX-2 inhibitors has aided in the improved management of pain, inflammation, and other associated disorders. In addition to the effective management of certain diseased states, COX-2 inhibitors have established enhanced gastrointestinal tolerability and fewer side effects in comparison with the nonselective NSAIDs [4, 5].

However, a long-term use of either selective COX-inhibitors or nonselective NSAIDs can lead to serious cardiovascular events [6, 7]. In order to overcome the side effects of certain drugs, their structures are either modified or new drugs are discovered with the aim to produce safe and therapeutically competent drugs [8].

Pyrazole scaffolds possess anti-inflammatory [912], analgesic [13], neuroprotective [14], antidepressant, and anticonvulsant activities [15, 16]. Some pyrazoline derivatives that possess analgesic activity and anti-inflammatory activity are phenylbutazone, celecoxib, and deracoxib.

In this context, 5d (Figure 1) was synthesized in-house in the medicinal chemistry laboratory with the aim of finding a new COX-2 inhibitor which might possess better tolerability and fewer side effects [17]. The synthesized compound 5d intensely downregulated COX-2 induction and suppressed its expression up to 42.8% compared to lipopolysaccharide- (LPS-) stimulated controlled cells (100% expression). Also, a strong inhibition of LPS-induced nitric oxide synthase expression was observed with 5d [17].

In the systemic circulation, the synthesized drug has to come in contact with proteins, thus leading to a drug-protein interaction. Serum albumin protein present in blood maintains the osmotic pressure (colloidal) and transports biomolecules such as steroids, fatty acids, and hormones inside the body [1820]. Drug molecules also interact with albumin in serum and are also carried by it in the systemic circulation. The pharmacodynamic and pharmacokinetic (absorption, distribution, metabolism, and excretion) properties of a drug are highly influenced by the protein-drug interaction [1, 2]. A reduced lifetime of a ligand is due to its weak binding to albumin and can also lead to its poor distribution. However, a strong binding interaction with the ligand leads to reduced free ligand concentrations. The study and understanding of binding interaction between ligand and albumin at the molecular level is therefore crucial in evaluating the biological activity of the ligand. Thus, the binding mechanism of drug with albumin, and its distribution needs to be investigated [2023].

Since both the BSA and human serum albumin (HSA) are homologous in structure [2], and due to low procurement cost of BSA, it is a popular approach to substitute HSA with BSA to represent the serum albumins in protein-ligand interaction studies. The binding affinity and mode in the BSA-5d system were determined using fluorescence spectroscopy.

The fluorescence quenching studies including synchronous and three-dimensional fluorescence and UV spectroscopic studies were conducted to comprehend the binding mechanism between BSA and 5d. Binding site of 5d on BSA was established with ligand displacement studies using site probe ligands phenylbutazone for site I and ibuprofen for site II [24, 25]. Molecular docking study for BSA and 5d interaction was also performed to corroborate with experimental results. The results of this study will be helpful in conducting the pharmacokinetic studies of 5d.

2. Materials and Methods

2.1. Materials

Compound 5d was synthesized in-house [17, 21]. The other materials procured were as follows: fatty acid-free BSA (Sigma Aldrich; USA), ibuprofen, and phenylbutazone (National Scientific Company; KSA).

2.2. Sample Preparation

The BSA stock was prepared in phosphate buffer saline (PBS; pH 7.4). The stock of 5d (5 mM) was prepared in acetonitrile and was diluted with PBS pH 7.4.

2.3. UV–Vis Absorption Spectra

The Shimadzu UV-1800 (Shimadzu Japan) spectrophotometer was used to record the UV-Vis absorption spectra. The spectra for BSA both in absence and presence of different concentrations of 5d were obtained. The measurements were performed at room temperature and fixed BSA concentration.

2.4. Fluorescence Spectroscopy

The fluorescence spectrometer (JASCO Model FP-8200) was used to record the fluorescence spectra. The spectra were acquired for BSA (fixed concentration) in absence and presence of different 5d concentrations. The fluorescence spectra were acquired at 298, 303, and 310 K. The spectra were acquired at λex = 280 nm and λem = 300–500 nm.

2.5. Fluorescence Quenching Measurement

The spectra were recorded for BSA in presence of various concentrations of 5d (0–18 μM). The concentration of BSA (1.5 μM) was held constant for all the measurements. The BSA fluorescence is quenched by 5d. Stern–Volmer equation was used for fluorescence data evaluation:where F and F0 are the fluorescence intensity (FI) of BSA in presence and absence of ligand 5d. Ksv represents Stern–Volmer quenching constant; Kq is the bimolecular quenching rate constant. [Q] is the quencher concentration, and “τ0” is fluorophore’s average life time.

Double-log regression plot was used to obtain the binding constant of BSA and 5d. The following equation is used to express the molecular equilibrium:

In (2) “n” is the binding site number and Kb is the binding constant.

The thermodynamic study helps to determine the type of interaction and bonds involved in the interaction. The standard Gibbs free energy change is given as

2.6. Synchronous Fluorescence (SF) and Binding Site Probe Experiment

The SF spectroscopy helps in determination and characterization of microenvironmental changes around the amino acid residues tyrosine (Tyr) and tryptophan (Trp). In the BSA-5d system, these changes were recorded using a constant difference of ∆λ = 60 and 15 nm between excitation and emission wavelengths. These SF spectra were attained at 298 K. Drug displacement studies were conducted with phenylbutazone and ibuprofen to establish BSA binding sites for 5d.

2.7. Three-Dimensional (3D) Fluorescence Spectroscopy

These spectra help in obtaining information about the chromophore. The data aid in finding any conformational alterations in protein (BSA) transpired by interaction with ligand (5d). 3D spectra were acquired for BSA and BSA-5d system.

2.8. Molecular Docking

Molecular docking was done on the molecular operating environment (MOE). The analysis was executed with BSA crystal structure (cocrystallized with naproxen, PDB ID: 4OR0). The 5d structure was drawn within the MOE software. The energy minimization was achieved from default field force parameters of MMFF94X, eps = r, and cutoff (8–10). Default docking parameter triangle matcher was used. The scoring function 1 was put to London dG, and the function 2 was put to GBVI/WSA dG. The most appropriate interaction was designated as per the root mean square values obtained [21].

3. Results and Discussion

3.1. Fluorescence Quenching Studies

Protein and ligand interaction can be understood with fluorescence quenching studies of proteins by ligands. The fluorescent behavior of BSA is due to the amino acid residues, Trp, Tyr and Phe residues, that are embedded in it. The ratio of fluorescence intensity of the three amino acids present in BSA is 100 : 9 : 0.5 for Trp, Tyr, and Phe, respectively [26]. Thus, the two Trp residues present in BSA are mainly responsible for its fluorescence. The residue Trp-134 is on the subdomain IB surface, whereas Trp-213 is in the subdomain IIA hydrophobic pocket. The BSA fluorescence because of Trp and Tyr residues has been used to study the conformation changes and ligand interactions. The quantum yield of Tyr is almost the same as Trp; however, the indole group of Trp in proteins is mainly liable for UV absorbance (∼280 nm) and emission (∼350 nm). The excitation spectra measured at λem = 340 nm (Figure 2(a)) for BSA and BSA-5d were compared, and they provided information of the effect of 5d concentration increase on the Trp behavior [27].

The BSA concentration was held constant. 5d concentration was varied in different samples to obtain the fluorescence quenching spectra (Figure 2(b)). The FI of BSA reduced as the concentrations of 5d increased. A small red shift of 1 nm in the λ(emissions) was detected which implied that the amino acids were present in a polar environment and had a higher exposure to the solvents [28]. The quenching mechanisms are classified as static and dynamic quenching. The dynamic quenching and static quenching are differentiated on the basis of effect of temperature and viscosity on them. In dynamic quenching, the quenching constant rises with temperature due to high diffusion and collision quenching. In contrast, higher temperature causes instability of the complex and lower quenching constants in the static quenching. Thus, lower quenching constant values at higher temperatures indicate formation of a nonfluorescent complex and static quenching mechanism [29].

Stern–Volmer equation (1) is used to interpret the quenching mechanism in the protein-ligand interaction. Low Ksv was observed at high temperatures (Figure 3(a)), thus inferring formation of BSA-5d ground-state complex (Table 1). Similarly, Kq (collision quenching constant) also helps in elucidation of quenching mechanism and is given as

The maximum value of Kq for a dynamic quenching process is 2 × 1010 L·mol−1·s−1 [25, 30], and the values for the BSA-5d system were quite greater than the maximum value (Table 1). These results suggest that static quenching was involved in BSA-5d interaction.

3.2. Binding Constant and Number of Binding Sites

From equation (2), the plot log[(F0/F) − 1] vs log[Q] yielded a straight line, where the binding stoichiometry (n) equals the slope and Kb (binding constant) from the intercept (Figure 3(b)) [31]. The binding constant at different temperatures and “n” values are given in Table 2 [32].

Site specific markers were used for identification of binding site studies using competitive binding experiments. The specific site markers phenylbutazone (PBT) for site I and ibuprofen (IBF) for site II of BSA were used to study the binding site involved in the competition with 5d (Figure 4) [24, 25]. The spectra were recorded for BSA in presence of site marker holding the concentration of both constants and increasing 5d concentrations at room temperature (Figure 3(d)). A huge reduction of the binding constant of the BSA-5d system with addition of either of the site markers PBT and IBF was observed (Table 3). As the binding constants decreased in either of the PBT and IBF and the decrease was more in presence of PBT than IBF, it was inferred that 5d might be binding to both the BSA at sites I and II and more preferably to site I.

3.3. Thermodynamic Parameters

Four main forces, namely, hydrophobic interactions, hydrogen bonding, and van der Waals and electrostatic interaction forces are involved in the ligand-protein interaction [3335]. The thermodynamic parameters are calculated from the Van’t Hoff equation (Figure 3(c)):where Kb and R represent the binding and the gas constant, respectively. The ∆H° and ∆S° are the enthalpy and entropy change, respectively; ∆G° represents Gibbs free energy. The calculated values of these parameters are given in Table 2. The value of ∆H° > 0 and ∆S° > 0 suggests a hydrophobic interaction. The van der Waals force or hydrogen bond formation is implied if ∆H° < 0 and ∆S° < 0, whereas ∆H° ≈ 0 and ∆S° > 0 indicate an electrostatic force. The results (negative ∆G°) clearly indicate a spontaneous binding interaction between 5d and BSA. The positive ∆H° and ∆S° infer major influence of hydrophobic forces in the interaction of BSA with 5d.

3.4. UV-Vis Spectroscopy

UV-Vis spectroscopy is widely applied for investigating protein-ligand interaction and structural changes that might have occurred because of it [3638]. The protein molecules present in the complex absorb light at λmax = 280 nm, primarily due to amino acid residues Phe, Trp, and Tyr present in BSA, and the other absorption peak in BSA occurs at 205 nm and is attributed to the polypeptide backbone of BSA.

The microenvironmental alterations within the vicinity of amino acids due to protein-ligand interaction provide information regarding the binding mechanism [39, 40]. An increase in the absorbance of BSA was observed at 280 nm and 205 nm [41]. This increase was attributed to the interaction of 5d and BSA, and the maximum absorbance shifted (1 nm) towards a higher wavelength (Figure 5(a)) at both the studied wavelengths. These shifts at 280 nm and at 205 nm suggest change in the microenvironment of tyrosine and peptide bonds in the protein. Thus, collective shifts at 205 nm and 280 nm suggest conformational changes in the protein structure [42]. A plot between the mixture of BSA and 5d versus concentration of 5d at 205 and 280 nm revealed that the plots obtained were nonlinear and deviated from the straight line (Figure 5(c)). This behavior of UV data suggests complex formation conformational changes in the BSA on its interaction with 5d. The normalized fluorescence excitation spectra (normalized for quenching effect) were also obtained for BSA and 5d mixtures and are given in Figure 5(d). After normalization, no change in excitation spectra were observed at ∼280 nm; however, the peak at ∼226 nm which represents the lowest excited state of aromatic amino acids (Trp and Tyr) of BSA showed a decrease in the intensity upon addition of 5d and suggests microenvironmental changes in the vicinity of Trp and Tyr residues. The involvement of dynamic quenching in the interaction of BSA and 5d was not considered as the absorption spectrum of BSA in presence of 5d might have been unaffected in this case. However, an increase in the absorption intensity of BSA on its interaction with 5d suggests complex formation between them (Figure 5(b)) and involvement of static quenching [43]. Moreover, the no-overlap was observed in the absorption spectra of BSA, 5d, BSA-5d, and (BSA-5d)-5d. These results in Figure 5(b) also confirm that complex formation took place between BSA and 5d.

3.5. Energy Transfer

The binding distance in BSA and 5d was calculated using fluorescence resonance energy transfer (FRET) [34, 35] (Figure 6). FRET occurs due to an overlap in the emission spectrum (donor) and absorption spectrum (acceptor). The FRET efficiency is dependent on the degree of overlap and distance between the emission and absorption spectrum and is calculated as [44]where energy transfer efficiency is given as E. “r” is the donor and acceptor distance, and R0 is the critical binding distance when energy transfers (50%) and is represented bywhere K2 denotes the dipole orientation; medium refractive index is given as “n”; “ϕ” is the donor’s quantum yield in absence of acceptor; and J denotes the spectral overlap extent in emission and absorption spectra and is given aswhere F(λ) is the emission spectrum fluorescent donor. ϵ(λ) denotes molar absorptivity coefficient of acceptor. The value for K2 = 2/3 (random orientation of fluid solution). The refractive index of the medium was taken as 1.336 and ϕD = 0.118 for BSA [43, 44]. The calculated values for E, R0, and J(λ) are presented in Table 4.

The binding distance r = 3.89 infers close proximity of the donor and acceptor to one another. Further, the distance amongst donor and acceptor (r < 8 nm) suggests a high possibility of nonradiative energy transfer from BSA to 5d. In addition, static fluorescence quenching mechanism can also be inferred between BSA and 5d.

4. Conformational Studies

4.1. Synchronous Fluorescence Studies

Synchronous fluorescence studies provide information regarding the microenvironmental changes in the amino acids responsible for fluorescence of protein [4447]. These studies measure the shift in the λem, which might have occurred due to variations near the chromophores of a molecule. The fluorescent residues Trp and Tyr of BSA can be analyzed for any changes in their microenvironment at ∆λ = 15 and 60 nm, respectively. The polarity near the chromophore molecule relates to the shift in the emission of maximum intensity. The shift in the position of fluorescence emission maximum corresponds to changes of the polarity around the chromophore molecule. Therefore, a blue shift in the maximum emission suggests a more hydrophobic microenvironment for amino acid residues, and a red shift infers a polar environment for them and more solvent exposure.

No shift in the spectrum at Δλ = 15 nm suggests insignificant effect of 5d on the Tyr microenvironment (Figure 7(a)). Moreover, insignificant red shift (1 nm) was also observed at ∆λ = 60 nm was observed (Figure 7(b)). These nonsignificant shifts in the emission wavelengths suggest that 5d interaction with BSA has no affect over the conformation of the regions where Tyr and Trp are located [42].

4.2. Three-Dimensional Fluorescence Spectroscopy

The 3D fluorescence spectroscopy has gained quite popularity to access the protein structural changes. It gives a detailed account of the conformation modifications in protein due to protein-ligand interaction [21, 43]. The changes in the conformation and the microenvironment of BSA are explored by relating the changes in the spectrum of BSA in absence and after addition of 5d (Figure 7). The identified peaks (Figure 8(a) and 8(b)) were as follows: peak a (Rayleigh scattering peak; λex = λem) and peak b (2nd order scattering peak; λem = 2ex) [46]. The intensity of peak a increased with the increasing concentration of 5d. The increased peak intensity is attributed to complex formation, and thus, a bigger macromolecule size leads to higher scattering effect. Peak 1 (λex/λem = 278/340 nm) represents the Trp and Tyr residues, and both are excited at ∼280 nm and are responsible for the fluorescence of BSA. The fluorescence of Phe is negligible and is not taken into account. Peak 2 (λex/λem = 226/340 nm) represents the fluorescence from the excited electronic states of Trp and Tyr residues present in BSA. The 3D-contour plots for fluorescence studies are presented in Figures 8(c) and 8(d). The ratio of intensities of peak 1 for the BSA: BSA-5d system was 1 : 0.77 and of peak 2 was 1 : 0.69 (Figure 8). Both the peaks are from different residues in BSA, and alterations in the FI upon interaction with 5d suggest some sort of microenvironmental changes in the BSA. Peak 2 which was believed to be from the peptide backbone has been proved wrong [48]. Peak 2 is due to the aromatic residues that get absorbed in the far UV-region, swiftly relax to the lowest excited state, and then emit in the range 300–350 nm. Therefore, Peak 2 also provides information about the microenvironmental changes that might have occurred during the Trp and Tyr residues as in case of excitation at ∼280 nm. The synchronous fluorescence results indicated no change in the protein conformation near Trp and Tyr residues. However, the excitation spectra in Figure 5(d) showed changes in the excitation of BSA upon addition of 5d. Figure 9 represents the Tyr and Trp residues corresponding to ∆λ = 15 and 60 nm, respectively, and 3D contour plot of fluorescence. The figure confirms that the value for Δλ = 60 nm leads to direct probe of tryptophan emission.

4.3. Molecular Docking

Molecular docking helps to understand and interpret the interaction mechanism between BSA and 5d [47, 49, 50]. Molecular docking acts as an in-silico tool to analyze the protein-ligand interaction and authenticate experimental interpretations. As is already known, the hydrophobic cavities of BSA in subdomains IIA and IIIA represent sites I and II, respectively. The site probe displacement studies suggested subdomain IIA/IIIA or both of BSA were involved in binding to 5d. The conformations which had least binding energy for site I and site II are given in Figure 10. The ligand 5d embedded in the hydrophobic cavity of subdomains IIA and IIIA with a free binding energy of −31.04 kJ·mol−1 and −30.12 kJ·mol−1. The ligand 5d is surrounded by the amino acid residues Ser 343, Ser 201, Ser 453, Arg 217, Arg 198, Arg 194, Leu 454, Leu 480, Leu 346, Leu 210, Leu 197, Ala 290, Ala 209, Val 342, Val 481, Tyr 451, and Asp 450 in the hydrophobic cavity of IIA (Figures 10(a) and 10(c)). Tryptophan residue Trp 213 was also found in close vicinity of 5d at site I. Further, one hydrogen bond with Arg 194 (N-Arg194, 2.88 Å) was observed. The hydrophobic residues involved in the interaction were Ala 290, Val 342, Leu 197, Leu 454, Leu 346, Leu 480, and Ala 209.

Similarly, the hydrophobic cavity of subdomain IIIA has Asn 390, Ser 488, Lys 413, Arg 409, Gln 389, Leu 386, Leu 490, Leu 429, Leu 456, Leu 452, Thr 448, Val 432, Phe 487, and Phe 402, surrounding residues IIA (Figures 10(b) and 10(d)). Tyrosine residue (Tyr 410) was found in the vicinity of 5d site II, and a hydrogen bond was found between 5d and Tyr 410 (N-Tyr410, 2.35 Å). Another hydrogen bond was found with Lys 413(N-Lys413, 3.15 Å). The hydrophobic residues in the vicinity of the site II were Leu 386, Leu 452, and Gln 389 [51, 52].

These results are in confirmation of the thermodynamic studies and site probe studies where it was inferred that hydrophobic forces contributed in the BSA-5d interaction as well as both sites I and II were also involved. The rigid and flexible docking were performed for the site I of BSA. In both the flexible and rigid docking, only one hydrogen bond was found with slight difference in the bond length between Arg 194 and 5d in the flexible (3.12 Å) and rigid (2.88 Å) docking. The amino acid residue His 241 contains an imidazole ring and can quench the Trp 213 fluorescence. The distance between His 241 and Trp 213 is 3.9 Å in the rigid docking, and when compared to the flexible docking, the distance was reduced to 3.7 Å. Since His 241 is placed closer to the Trp 213 residue, influence of His 241 on Trp 213 cannot be ruled out. Further, the distances between Trp-213 and 5d were 3.3 Å in the rigid and 3.5 Å in the flexible docking studies, and these results were consistent with the results of the fluorescence studies, where no significant shift was observed on interaction between BSA and 5d [53].

5. Conclusion

The study presents the binding interaction between BSA and 5d using different spectrometric and molecular docking techniques. The BSA-5d interaction followed a static quenching mechanism. The binding sites, site I and/or site II of BSA, participated in the interaction with 5d. The site probe and molecular docking study results corroborated one another. Three-dimensional fluorescence and UV-spectrophotometric studies indicated conformational changes in BSA due to its interaction with 5d. Positive enthalpy and entropy of the BSA-5d system suggested hydrophobic interaction. The present study can be useful to understand pharmacological and biochemical behavior of 5d since strong or weak binding to these proteins plays a major rule in the therapeutic capability of drugs. Binding of the drug to serum albumin also influences their distribution and elimination, and 5d has a moderate binding affinity to BSA.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to extend their sincere appreciation to the Deanship of Scientific Research. This work was supported by the Deanship of Scientific Research, King Saud University (research group no. RG-1435-073).