Abstract

In this study, the interaction between the coumarin derivative: N-(diphenylmethyl)-2-[(2-oxo-2H-chromen-4-yl)oxy]acetamide, biologically active drug, and human serum albumin (HSA) was investigated by using various optical spectroscopy techniques along with the computational technique. The results of steady-state fluorescence spectroscopy show that the static quenching occurred while increasing the coumarin drug concentration into HSA. Also, the binding constant (K) and thermodynamical parameters of enthalpy change (ΔH°), entropy change (ΔS°), and Gibbs free energy change (ΔG°) were calculated at different temperatures (293 K, 298 K, and 303 K). The results are in good agreement with those of molecular docking studies, and also, the docking study was carried out to understand the hydrogen bonding and hydrophobic interaction between human serum albumin and coumarin derivative. In addition to the docking, charge distribution analysis was done to understand the internal stability of coumarin derivative active sites of human serum albumin. Further time-resolved emission spectroscopy (TRES) studies were carried out between free HSA and HSA-coumarin complex, and the result confirms the presence of the drug in the protein molecule without cytotoxicity.

1. Introduction

Coumarin and its derivatives have been known for a long time, yet their vital role in plant and animal systems has not been thoroughly analyzed. Many coumarins and their derivatives exert antitumor, antiviral, anticoagulant, anti-inflammatory, antioxidant, vasorelaxant, antimicrobial, and enzyme inhibition properties. These compounds are precisely known to exert an antitumor effect and can cause significant changes in the regulation of the immune response, cell growth, and differentiation. Coumarins have recently attracted much attention due to their broad pharmacological activities as precursor molecules for the synthesis of some synthetic anticoagulant drugs [16]. Coumarin (2-oxo-2H-chromene) is a naturally occurring heterocyclic compound isolated from various plants; it has been found in the microorganisms, animals, essential oils, and fruits [4, 712]. N-(diphenylmethyl)-2-[(2-oxo-2H-chromen-4-yl)oxy]acetamide (Figure 1(b), already reported by Govindhan et al.) [13] which originates from 4-hydroxycoumarin shows the functional biological activity and also the more biologically active site of the coumarin family make more interest on, to understand the pharmacokinetic mechanism into the carrier molecule.

Serum albumins (Figure 1(a)) are the most fruitful proteins in the circulatory route of a wide diversity of various creatures and also blew the lid off the crucial first nature in the binding and transport of various endogenous and exogenous ligands. The molecular interactions of proteins and lots of compounds like medicine and a few organic molecules are investigated successfully. However, the binding of 4-hydroxycoumarin products to serum albumins is now an emerging field to explore the role, in biological systems, and their mode of action in various diseases [14].

In the present study, synthesized N-(diphenylmethyl)-2-[(2-oxo-2H-chromen-4-yl)oxy]acetamide interaction with human serum albumin was carried out using various optical spectroscopy techniques. Also, computational techniques like molecular docking and density functional theory (DFT) were performed further to understand the hydrogen bonding, hydrophobic information with active site residues, internal stability, and charge distribution of the molecule at the atomic level [15].

2. Experimental Section

2.1. Reagents

N-(diphenylmethyl)-2-[(2-oxo-2H-chromen-4-yl)oxy]acetamide is newly synthesized, which is already reported by Govindan et al. [13], and the same drug was dissolved in 0.5% DMSO in the PBS buffer. Human serum albumin was purchased from Sigma-Aldrich, USA, and used without any further purification.

2.2. Fluorescence Measurements

The emission spectroscopy study was performed using a spectrofluorometer (ISA FluoroMax-2; Jobin Yvon-Spex, Edison, NJ); the spectral bandpass was fixed at 5 nm for monochromators of both wavelengths of excitation and emission, and the recorded wavelength region is 300–500 nm at 280 nm excitation. The lifetime fluorescence measurements of the coumarin derivative and human serum albumin complex were detected using Fluorolog (Horiba Jobin Yvon IBH, UK) with a fast response sensitive PMT detector (Hamamatsu Photonics, Japan) for an excitation wavelength of 280 nm to 450 nm. The decays were sliced and normalized each wavelength for HSA and HSA-coumarin complex. The excitation-emission matrix (EEM) was mapped with the following two conditions: the emission wavelength was recorded between 200 and 700 nm and the initial excitation wavelength was set to 200 nm with an increment of 10 nm. In EEM measurements, the spectral bandpasses were kept at 5 nm for both excitation and emission. In addition to that, the Forster resonance energy transfer analysis between the donor and acceptor and further single-point polarization studies between human serum albumin and the coumarin derivative complex were also carried out.

2.3. FTIR and CD Spectroscopy Measurements

Fourier-transform infrared spectra (FTIR) of the human serum albumin solutions were administered at temperature by using an FTIR spectrometer (Interspec 2020; Estonia) in the attenuated total reflection (ATR) mode with a resolution of 4 cm−1 and 32 scans. The secondary structural analysis was carried out using circular dichroism (CD) for free human albumin and complex using a spectropolarimeter (Jasco J715; MD, USA) employing a quartz cell of 0.1 cm beam bath.

2.4. Docking Studies

The human serum albumin structure was taken from the Protein Data Bank (PDB ID: 4L8U), and the coumarin derivative is drawn using the ChemDraw software. Proteins and small molecules are energetically minimized using the OPLS 2005 force field (Schrödinger LLC) used for molecular docking [15].

2.5. DFT Methods

Energy minimization has been performed using the density functional theory (DFT) methodology with B3LYP functional and 6-311G∗∗ basis sets for both coumarin molecule and cofactor with active site amino acids of HSA using the Gaussian 03 software. The electrostatic potential countermaps were pictured using the GView software [16].

2.6. Cell Viability Studies

The cytotoxic assay for the drug was obtained using the Cell Counting Kit-8 (CCK-8) (Dojindo Laboratories, Osaka, Japan). The cell density in the 96-well plate was 5000 cells/well, 24 hours later; the drug was added to the culture plate and incubated for another 24 h. After drug exposure, the cells were incubated with 10 µL CCK-8 for 2 h, and the absorbance at 450 nm was measured using a microplate reader (MTP-880Lab; Corona, Hitachinaka, Japan). The cytotoxicity was represented as a percentage of untreated control cells.

3. Results and Discussion

3.1. Fluorescence Spectroscopy Studies

The steady-state fluorescence spectroscopy is the most effective technique to understand the further binding mechanism of the HSA-coumarin derivative complex. In general, the human serum albumin emission will be at 350 nm due to the presence of tryptophan residue and other phenylalanines; tyrosine residues emission is considered to be less compared to tryptophan residues emission. The emission spectrum of HSA (1 µM) and various concentrations of the coumarin derivative (0–3.6 µM) is shown in Figure 2(a), and it is indicated that there is a decrease without a shift in emission maximum of the HSA-coumarin derivative complex. Hence, we may conclude that fluorescence quenching occurs due to the presence of the coumarin derivative [17].

To know more about the type of quenching, we performed Stern–Volmer (Equation (1)) calculation at a different temperature as shown in Figure 2(b), and the calculated results of the quenching constant are tabulated in Table 1 [18, 19]:where and is the free HSA and HSA complex (coumarin drug), is the bimolecular quenching constant (), is the lifetime of the fluorophore without the drug molecule (10−8 s), is the Stern–Volmer quenching constant (0.79 × 10−5, 1.00 × 10−5, and 1.33 × 10−5 mol−1·s−1·L−1), and [Q] denotes the concentration of the drug. From Table 1, the calculation result of indicates that static quenching occurs due to the complex formation of the coumarin derivative. The binding constant and number of binding sites were calculated using the following equation [20]:where are the total concentration of HSA and coumarin drug, respectively, and the binding constant and number of binding sites were obtained and tabulated in Table 1 using the least-square fitting algorithm (Figure 2(c)).

3.2. Thermodynamic Parameter Analysis

The thermodynamic study is important for the drug-binding mechanism to calculate the binding stability of the drug with the macromolecule complex. Thermodynamics parameters ΔH°, ΔS°, and ΔG° are calculated from the following equations:where is the gas constant and is the temperature, and here, three different temperatures are used: 293 K, 298 K, and 303 K.

According to Ross and Subramanian [17, 21, 22], the enthalpy change (ΔH°) and entropy change (ΔS°) tabulated in Table 1 indicate that (ΔH° < 0 and ΔS° > 0) type of interaction between HSA and the coumarin derivative complex is maybe due to the electrostatic interaction.

3.3. Time-Resolved Emission Spectroscopy Studies (TRES) and Fluorescence Lifetime Studies

Time-resolved emission spectroscopy is one of the highly noninvasive techniques to understand the real-time dynamics of biomolecules. Figure 3(a) shows the emission decay profile of free HSA and HSA-coumarin derivative complex; each point consists of average photons of respective wavelengths. All the recorded emission profiles of HSA and HSA complex followed the multiexponential decay kinetics, and the average lifetime decay for free HSA was found to be 4.83 ns which is almost equal to previously reported values [21, 23]. However, in the case of HSA-coumarin derivative complex, the decay is slightly faster compared to that of free HSA, and the average estimation of the HSA complex is 4.67 ns (Table 2). From Figure 3(b), the observed result indicates that quenching occurs due to the binding of the coumarin derivative with HSA. It is important to mention that the emission decay profile of the HSA complex is faster than that of free HSA in all wavelengths due to carrying of the coumarin derivative. This present result elucidates that protein unfolding is possible in the HSA complex because of the presence of the drug, and also, the tryptophan region gets more exposed to water. These results suggest that free HSA is more labile and has faster hydration dynamics, and the similar process occurs in the HSA-coumarin derivative complex without any structural disturbance. Figure 3(b) represents the decay parameter-fitted values τ0 of free HSA and HSA-coumarin derivative complex for easy understanding.

3.4. Excitation-Emission Matrix (EEM) Analysis

The study is used to further understand the conformational change of HSA more carefully during the interaction with the coumarin derivative. Figure 4 shows the contour map of free HSA (a) and coumarin derivative-HSA complex (b) with the corresponding characteristic parameters presented. There are two major fluorescence peaks observed in the figures: peak 1 and peak 2. Peak 1 (λex = 220 nm and λem = 325 nm) is mainly due to polypeptide backbone structures, which correlated with the secondary structure of the protein, and peak 2 arises mainly (λex = 278 nm and λem = 344 nm) due to the presence of tryptophan and tyrosine residues. From Figure 4, it is observed that there is a decrease in intensity maxima while adding the coumarin derivative into free HSA, and the quenching in both peak 1 and peak 2 may be due to the presence of the coumarin derivative, suggesting that the compound may be nearby the HSA chromophore. These results conclude that there may be slight confirmation changes in HSA due to the presence of the coumarin derivative [23, 24].

3.5. FTIR Spectroscopy Studies

In general, the FTIR spectral technique is performed to understand the biochemical changes; here, the FTIR spectra of HSA were recorded to investigate the conformational changes to understand the binding mechanism of the coumarin derivative. Figure 5(a) shows the spectrum of free HSA and HSA-coumarin derivative complex, and the amide I peak of HSA is lying between 1600 and 1700 cm−1 and amide II region is almost in 1558 cm−1, which is almost similar to the secondary structure of proteins. Comparing these two bands, amide I is strong (1644 cm−1) as shown in Figure 5(a), and also while adding the coumarin derivative into free HSA, there is a decrease observed in the absorption of amide I maxima, which may be due to the coumarin derivative disturbance of the HSA amide region [21, 23, 24].

3.6. Circular Dichroism (CD) Spectral Analysis

To understand the conformational changes of the protein secondary structure, CD spectral analysis was carried out. From Figure 5(b), CD spectra of HSA show that there are two negative bands in the ultraviolet region, at 208 nm and 222 nm. When the coumarin derivative is added into the free HSA molecule, there are notable increases observed in both 208 and 222 nm intensities while increasing the concentration and also particularly at 208 nm region of HSA (1 µM) and HSA-coumarin derivative (1.6 µM) show a slight shift towards the lower wavelength region, indicating that the α-helical structure may lose stability or rearrange the orientation of total structural confirmation due to presence of the coumarin derivative.

Secondary structure calculations were carried out using the DichroWeb (http://dichroweb.cryst.bbk.ac.uk/html/home.shtml) online tool for the free HSA that consists of ∼54.2% α-helix, ∼22.1% β sheets, and 16.3% of the random coil, which are almost comparable with previously published data [17, 24]. For the coumarin derivative at 1.6 µM concentration, the secondary structure parameters increased with respect to ∼58.1% α-helix, ∼24.7% β sheets, and 19.2% of the random coil, respectively. These results show that the secondary structure of HSA may be partially unfolded or reorientation occurs in the overall structure due to the presence of the coumarin derivative.

3.7. Molecular Docking Studies

The molecular docking method is useful to understand the binding insight mechanism between the macromolecule and the small molecule at the atomic level [25, 26]. In this work, the molecular docking between HSA and coumarin derivative was performed using AutoDock 4.2 to clarify further the binding mode of coumarin derivative and HSA. Here, both small molecules and macromolecules are allowed for flexible docking, and the free energy values and hydrogen bonding and hydrophobic information results are shown in Table 3. Docking results show that the binding interaction between HSA and the coumarin derivative is in good agreement regarding energy and interaction. Figures 6(a) and 6(b) show the hydrogen bonding and hydrophobic information of the HSA-coumarin derivative complex further to understand more about internal stability of the coumarin derivative, and the charge distribution analysis was carried out at active residues of the HSA-coumarin derivative complex.

3.8. Charge Distribution Analysis

The Mulliken atomic charges were analyzed for both coumarin derivative and cofactor with interaction amino acids of HSA. The MPA charges are presented in Table 4 (Figure 1(c)). The nitrogen atoms N8, N10, and N11 in the amino acid arginine possess high negative charges; simultaneously, the oxygen atoms O48 and O56 possess low negative charges. Because of these nitrogen atoms, interaction occurs with O48 and O56 atoms (N8: −0.44e, N10: −0.45e, N11: −0.46e, O48: −0.16e, and O56: −0.38e), respectively; these charges are low compared with those in the cofactor. The oxygen atom O69 in tyrosine has a very low negative charge, which interacts with O91. At the same time, the ligand contains nitrogen atom N45, which also has a low negative charge (N45: −0.39e, O69: −0.16e, and O91: −0.30e). From this, the MPA charges decreased while increasing the distance between the protein and ligand. The carbon C9 atom possesses a high positive charge due to the carbon atom surrounded by nitrogen atoms [27].

The global reactivity descriptors [2834] of the coumarin derivative and cofactor molecule such as electronegativity, electrophilicity, and chemical hardness have been calculated from HOMO and LUMO energies. The calculated values of ionization potential, electronegativity, and electron affinity of those molecules are 6.5/7.2, 8.1/9.6, and 1.6/2.4 eV (Table 5), and these values show that the coumarin derivative molecule has high tendency to accept electrons than the cofactor. The calculated electrophilicity of both molecules is 6.8/9.6 eV; this low value indicates that the molecule has low toxicity [35]. The chemical hardness of the molecule is 4.85/4.82 eV, which reveals that the molecule exhibits high kinetic stability and reactivity [16, 3638]. The energy gap between the molecular orbitals is also calculated, and the HLG value is 4.86/4.82 eV, which indicates that both molecules have high internal stability. Figure 7 shows the electrostatic potential countermap of the coumarin derivative and cofactor with active site amino acids of HSA.

3.9. Steady-State Fluorescence Anisotropy Studies

Fluorescence anisotropy is widely used to measure the binding constants and kinetics of reactions that cause a change in the rotational time of the fluorescent molecules [24]. This technique will give information about changing the orientation of our coumarin derivative in human serum albumin. Here, the fluorescence anisotropy measurements were carried out between human serum albumin and a coumarin derivative of concentration ranging from 0 to 100 µM. Most reports say that anisotropy of small molecules may be zero and there is no rotation when it increases to the maximum value of 0.4. From Figure 8, coumarin derivative versus various increasing concentrations of human serum albumin shows that the anisotropy value of the coumarin derivative increases up to 20 µM and then reaches the steady-state maximum value (r = 0.13) after 40 µM of the HSA environment. This saturated value may be due to the coumarin derivative restriction in rotation motion of the HSA molecule, which also may be due to the formation of a hydrogen bond and hydrophobic interaction in active sites of HSA with the coumarin derivative.

3.10. Fluorescence Resonance Energy Transfer

The Forster energy transfer method is used to understand the energy distance between the donor and acceptor of biomolecules. The overlap absorbance spectrum of the drug and emission spectrum of human serum albumin are shown in Figure 9. The energy transfer efficiency (E) is used to calculate the distance between the coumarin derivative and human serum albumin. The rate of energy transfer depends on the integral overlap (J) of the donor (HSA) and acceptor (a coumarin derivative). Transition dipoles and distance between the drug and protein molecules also play a role in deciding the rate of energy transfer [17, 24]. The distance from the tryptophan residue which is located in the donor (HSA) and bound drug acceptor (a coumarin derivative) in HSA can be calculated using Forster’s theory.

The energy transfer efficiency (E) was calculated using the following equation [39]:where is the distance between the donor and acceptor and is the critical distance between the donor and acceptor when their transfer efficiency is 50%:where is the orientation factor of the dipole, φD is the fluorescence quantum yield of the donor in the absence of the acceptor, and is the refractive index of the complex medium:

J is the integral overlap of the emission spectrum of the HSA (donor) with the absorption spectrum of the coumarin derivative (acceptor). FD(λ) is the fluorescence intensity of the donor at wavelength λ to λ + Δλ, with the total intensity normalized to unity, and εA(λ) is the molar extinction coefficient of the acceptor at the wavelength (λ). It is possible for us to calculate the J value from the spectrum in which it is estimated that J = 1.3182 × 10−20 M−1·cm3. Under the experimental conditions, we have found R0= 2.5 nm from Equation (5). Using the following values, E = 0.018, k2 = 2/3, and ϕD = 0.118, r = 4.88 nm has been calculated. Both R0 and r values were found to be lower than the maximum values of R0 (<10 nm), and the distance (r) between the donor and acceptor lies <7 nm.

3.11. Cytotoxicity of the Coumarin Derivative

The cytotoxicity of a coumarin derivative, which was used to evaluate its biocompatibility, was studied using the HeLa cell as a model cell by the CCK-8. Figure 10 shows the cell viability test results conducted for the coumarin derivative up to 30 µg/mL. As shown in Figure 10, the amount of the coumarin derivative increases and the cell viability decreases in a tiny percentage even at the maximum level of concentration at 30 µg/mL. However, coumarin derivative is a safe molecule to bind with the serum albumin protein and deliver the targeted tissue region.

4. Conclusion

In this context, the synthesized coumarin derivative molecules have been subjected to various spectroscopy and computational techniques for understanding the interaction mechanism in the HSA molecule. From the UV-Vis absorbance result, it can be seen that the HSA-coumarin derivative complex may be due to the ground state complex; further fluorescence technique reveals that the coumarin derivative in the HSA complex is in the static quenching mode, and also, the thermal parameters calculated show good binding results. Time-resolved emission spectroscopy results show that the HSA molecule could be able to carry the coumarin derivative, and the calculated average lifetime value found for free HSA is 4.83 ns and for the HSA-coumarin derivative complex is 4.67 ns. Energy transfer analysis also calculated between donor and acceptor and r values is found that 4.88 nm, which is in academic value. The anisotropy results conclude that the rotation motion of HSA is saturated at r = 0.13 after the concentration of HSA 40 µM which is maybe due to the presence of the coumarin derivative. The excitation-emission matrix results suggest that the confirmation of HSA changed due to the presence of the coumarin derivative, and following this, the FTIR and CD spectral analysis also confirmed those secondary structure changes. Computational analysis elucidate the interaction mechanism of the HSA-coumarin derivative complex at the atomic level and information on bonding and internal stability; from this overall result, the binding nature of coumarin derivatives with HSA is good, and it could carry the molecule to a particular target without cytotoxicity.

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.