Abstract

Diorganotin(IV) antitumor compound bis-[2,6-difluoro-N-(hydroxyl-<κ>O)benzamidato-<κ>O] (DBDF2,6T) was one of the novel patent organotin compounds with high antitumor activity and relatively low toxicity. In this study, several methods were used to study the interaction between DBDF2,6T and hPPARγ protein, including fluorescence quenching, three-dimensional (3D) fluorescence, drug affinity responsive target stability (DARTS), ultrafiltration-LC, and molecular docking. According to the experimental results, the quenching process of the hPPARγ protein was induced by static quenching mode to form a nonradiative ground-state complex with DBDF2,6T spontaneously, mainly through the hydrophobic force. DBDF2,6T could bind to the hPPARγ protein directly and give the protein the ability of antienzymatic hydrolysis. And the binding mode of DBDF2,6T into hPPARγ protein appeared to have an orientation towards residues of SER342 and GLY284. In conclusion, these methods could comprehensively reveal the interaction details of DBDF2,6T and the hPPARγ protein and established a feasible way to preliminarily identify the agonist compounds for the hPPARγ protein.

1. Introduction

Organotin compounds have many uses in our life, which could act as stabilizers in plastics, fungicides, industrial catalysts, and so on. [1]. Our research group had synthesized a series of organotin patent compounds which possessed high anticancer activity with low toxicity and devoted to clarify its mechanism of action [2]. From the results of proteomics data, these compounds might play the physiological role through the PPAR (peroxisome proliferator-activated receptor) signaling pathway, which was consistent with the reports that organotin compounds may function as endocrine-disrupting chemicals by affecting the function of the protein PPARγ [3]. Consequently, a reasonable hypothesis was made to assume that these biologically active compounds might function through the PPAR signaling pathway as an agonist to the protein PPARγ and further influence the expressions of the target genes.

PPARs proteins belong to the most important members of the nuclear receptor superfamily and can act as the ligand-activated transcription factors [4]. When the PPARs proteins bind to a specific ligand, the ligand-binding domains of PPARs will encounter the conformational change followed by promoting the recruitment of nuclear receptor coregulators such as steroid rector coactivator-1 (SRC-1) and eventually influence the transcription of downstream target genes [5]. The PPARs proteins have three isotypes which had been identified as PPARα, PPARβ/δ, and PPARγ. These three subtypes exhibit distinct tissue distributions and have unique biological functions [6, 7]. In particular, PPARγ has received much focus these years for the important physiological functions played by its ligands. For example, thiazolidinediones (TZDs), a class of PPARγ agonist compounds, had been used as a therapeutic compound for metabolic disorders such as type 2 diabetes and obesity [8], and it was also reported that the agonists to PPARγ protein had a potential to be used as a new therapeutic approach to cancers, immune disorders, and so on. [9, 10]. Therefore, the experiments established to find ligands which could interact with the PPARγ protein are promising works nowadays.

In this study, a patent organotin compound DBDF2,6T (bis-[2,6-difluoro-N-(hydroxyl-<κ>O)benzamidato-<κ>O]dibutylitin) (patent number: CN200910074795.X and ZL01135148.9 (P)) which showed a high antitumor activity was assumed as a potential agonist. Several different methods were adopted to test and verify the interaction between DBDF2,6T and the hPPARγ protein. Spectroscopic study, one of the most widely used methods for analyzing the interaction between small molecule and protein, was applied to provide parameters such as binding constants and types of interaction forces. DARTS and ultrafiltration-LC were used to verify such interaction while molecular docking was used to evaluate affinity between receptors and ligands in a theoretical way. These methods meet the requirements of low cost and high feasibility and perfectly supplement and verify each other, which could be used to find the new agonists of the hPPARγ protein preliminarily and offered references for the interaction analysis between synthesized compounds and proteins. The structures of DBDF2,6T and the hPPARγ protein are shown in Figure 1.

2. Materials and Methods

2.1. Reagent

DBDF2,6T was synthesized by Shanxi Medical University with purity over 99%. The protein human PPARγ (hPPARγ) was purchased from Flarebio Company (Flarebio Biotech LLC, Wu Han, China) and stored at −20°C. Pronase was purchased from Roche Diagnostic GmbH (Mannheim, Germany) and stored at 4°C. The Coomassie Blue R-250 was purchased from Sigma Company (Shanghai, China), and the PageRuler Prestained Protein Ladder was made by Thermo Fisher Scientific (Massachusetts, America). The solvent of organotin(IV) was configured by propanediol (Tianjin Fengchuan Chemical Reagent Science and Technology Co., Tianjin, China), ethylenediamine, and normal saline (Shijiazhuang Pharmaceutical, Shijiazhuang, China) (90 : 9 : 1, v/v/v). The disodium hydrogen phosphate dodecahydrate was purchased from Tianjin Fengchuan Chemical Reagent Company (Tianjin, China). The potassium phosphate monobasic was purchased from Tianjin Beichen Fangzheng Company (Tianjin, China).

All the other reagents used in this study were of analytical grade and were obtained commercially.

2.2. Fluorescence Quenching Spectrum

Amino acid residues such as tryptophan, tyrosine, and phenylalanine could empower the proteins with the ability to generate endogenous fluorescence. The fluorescence peaks of those three amino acids were located at 348 nm, 303 nm, and 282 nm, respectively. Actually, 95% of protein fluorescence was contributed to the tryptophan residue [11, 12]. Compared with other methods, fluorescence spectroscopy had many superior advantages including high sensitivity, selectivity, and easy operation [13]. Therefore, in this paper, the fluorescence quenching method was used to analyze the interaction between DBDF2,6T and the hPPARγ protein.

The experiments were performed at 293 K and 310 K on a U-3900 spectrofluorophotometer (BaHens Instrument Co. Ltd., China). Protein hPPARγ (10 μg) was dissolved in a 2 mL PBS buffer. Several concentrations of DBDF2,6T (0.5 × 10−6, 1.0 × 10−6, 1.5 × 10−6, 2.5 × 10−6, 3.0 × 10−6, and 3.5 × 10−6 mol/L) were, respectively, incubated with the certain concentration of the hPPARγ protein. Samples of protein hPPARγ and hPPARγ-DBDF2,6T complexes were measured in a 1 cm2 quartz cuvette. And the excitation and emission spectral slit widths were set as 10 nm. The emission spectra were recorded for light-scattering effects from 300 nm to 450 nm while the exciting wavelength was set as 280 nm.

2.3. Three-Dimensional Fluorescence Spectrum

The coordinate axes of the three-dimensional (3D) fluorescence spectrum were excitation wavelength, emission wavelength, and fluorescence intensity. It had been proved that the 3D fluorescence spectrum was an effective analytical technique to analyze the conformation changes of a protein in its solution state [14]. And this method could not only test the molecular structure change with much selectivity and sensitivity but also display fluorescent information of the sample solution comprehensively [15].

Experiments were performed at the temperature of 293 K. Protein hPPARγ (10 μg) was dissolved in a 2 mL PBS buffer and incubated with 2.0 × 10−6 mol/L DBDF2,6T for two minutes. Samples were tested on a U-3900 spectrofluorophotometer (BaHens Instrument Co. Ltd., China) with the parameters set as follows: excitation wavelength was from 200 nm to 300 nm; emission wavelength was from 320 nm to 450 nm; spectral slit width was 10 nm; and the gain value was 2.

2.4. DARTS with Pure hPPARγ Protein

DARTS had been proved to be an efficient approach to efficiently verify drug-protein interactions when the protein was available in relatively pure form [16]. The basic principle of DARTS was that compounds were proposed to stabilize the combined protein globally or locally by reducing protease sensitivity of the target protein. This phenomenon was attributed to a specific conformational change caused by such a binding process, which would further induce protease recognition sites of the protein to be masked [17]. In this study, protein hPPARγ regarded as target protein and pure hPPARγ protein generated from recombinant plasmid were used in this experiment. Whether the presence of DBDF2,6T could reduce the proteolysis of the protein to validate the interaction between protein hPPARγ and DBDF2,6T should be observed after incubating the hPPARγ protein with DBDF2,6T.

Eight sample groups were set and divided into blank group, negative control group, and test group; each sample contained 0.5 μg hPPARγ protein. Except for the blank group, other samples were incubated with 2 μL DMSO or 2 μL DBDF2,6T with the concentration ranging from 1.0 × 10−2 mol/L to 1.0 × 10−4 mol/L for 60 min at 4°C and then digested with pronase (1 : 100) at room temperature for 30 min. The digestion was stopped by adding 5× SDS-PAGE sample loading buffer and boiling at 100°C for 10 min immediately. Samples were then subjected to electrophoresis on 8% SDS-PAGE. After electrophoresis, the gel was stained with Coomassie Brilliant Blue for 1 h and was observed after eluted overnight.

2.5. Ultrafiltration-Liquid Chromatography Experiment

Ultrafiltration-liquid chromatography (ultrafiltration-LC) was developed to verify the agonists of protein hPPARγ. The main principle of ultrafiltration-LC was that the agonists of hPPARγ had the ability to bind to the protein and would not be filtered out through the membrane of the ultrafiltration centrifuge tube, while after the protein denatured by dissolving in organic solvents, the compounds would be unbound to the protein and could be washed out through the membrane of the ultrafiltration centrifuge tube. Based on the ultraviolet absorption of the compound, the DBDF2,6T which were bound to hPPARγ could be detected through analyzing the washed solution by liquid chromatography. This method was first used in such confirmatory experiment and marked by its simplicity, generality, and applicability [18].

The recombination protein hPPARγ (20 μg) was incubated with compound DBDF2,6T (2 μL·10−3 mol/L) for 24 h at 4°C. After being filtered through a 10000 Da molecular weight cutoff ultrafiltration membrane (Millipore, UFC500396) by centrifugation at 13000 r/min for 8 min at 4°C, the sample was washed three times with 150 μL PBS buffer (pH 7.4) and centrifuged at 13000 r/min for 12 min at 4°C to remove the unbond compounds. The washed solution was transferred to a new 10000 Da molecular weight cutoff ultrafiltration centrifuge tube and dissolved it in 400 μL methanol. Centrifugation at 13000 r/min for 12 min was performed to wash out the compounds which were combined with the protein, and the washed solution was collected. Following reconstitution in 100 μL of 50% aqueous methanol, the compound was analyzed using HPLC (Agilent Technologies). Denatured protein was used as a negative control, and in this experiment, the protein was heated at 98°C for 15 min to make it denatured.

The HPLC analysis was carried out using mobile phase methanol/0.5% phosphoric acid (28 : 72, v/v, pH 3.0) on a C18 column (Agilent TC-C18, 4.6 × 250 mm i.d., 5 μm) at a flow rate of 0.8 mL/min and at 25°C. The detection wavelength was set at 264 nm.

2.6. Molecular Docking

Surflex-Dock, docking module in SYBYL software (UCSF), was performed to determine the binding model of protein hPPARγ (4a4w.pdb) and DBDF2,6T. It used prototype molecule (protomal) to represent protein binding pocket, utilizing probe to test the qualities of protein pocket such as surface hydrophobicity and could generate the invert transform of an active protein pocket. This method had high docking accuracy and could be used to research on the interaction between biomacromolecules and small molecular ligands [19].

Sybly × 2.0 was used to draw two-dimensional structure of DBDF2,6T with standard bonds and angles. In the process of optimizing the compound structure, minimize details and parameters of modify were set as follows. Minimize details: the iterations were set as 10000, and the color option was set as force. Parameters of modify: the force field was set as Tripos and the charges were set as Gasteiger–Marsili. In the docking process, A/YFB99 was chosen as extracted ligand structure and hydrogen molecules were added to the protein, and the modify details were set as follows: the force field was AMBER7 FF99 and the charge was AMBER. All other parameters were used the default value of SYBYL during the protein pocket generation and the molecular docking.

3. Results and Discussions

3.1. Fluorescence Spectrum of Protein hPPARγ and hPPARγ-DBDF2,6T Complexes

Fluorescence curves of the protein hPPARγ and hPPARγ-DBDF2,6T complexes are shown in Figure 2. The experiments were carried out in two temperatures with other experimental conditions unchanged. It can be observed from Figure 2 that, in both the experimental system of 273 K and 310 K, the intensity of fluorescence of the hPPARγ protein was decreased while increasing the concentration of DBDF2,6T which was incubated with the hPPARγ protein. And the highest fluorescent intensity values of hPPARγ were recorded for next step of theoretical calculation.

3.2. Mechanism of Fluorescence Quenching

While the pH, temperature, and ionic strength were kept as constants, the types of fluorescence quenching could be classified into two categories: dynamic quenching and static quenching [20]. Dynamic quenching was caused by the fluorescent chromophore interacted with a quencher in excitation state, while the causes of static quenching were of three types: the first one, the fluorescent chromophore interacted with a quencher in ground state and came into being a nonfluorescent compound; the second one, the medium near the fluorescent chromophore had a polarity change, which caused by the conformational change of the protein attributing to the combination with the quencher; and the third one, a radiationless energy transfer between the fluorescent chromophore and the quencher [21].

The dynamic quenching obeyed the Stern–Volmer equation, and the formulas are shown as follows:where is the Stern–Volmer quenching constant, is the bimolecular quenching constant, and is the average lifetime of the molecule which always be considered as 1.0 × 10−8 s [22, 23]. is the concentration of the quencher. and correspondingly represent the intensity of the fluorescence of the protein added with the quencher or not.

The Stern–Volmer quenching curves were drawn according to the data obtained from fluorescence quenching spectra, and the corresponding linear regression equations and correlation coefficients are shown in Table 1.

When a quencher interacted with a biomacromolecule, the maximum value of diffusion collision rate constant was considered as 2.0 × 1010 L/mol/s. According to the computing results shown in Table 1, the dynamic quenching constant () between DBDF2,6T and the hPPARγ protein was of the order of magnitude of 1012, which was much bigger than the maximum value of diffusion collision rate constant. Consequently, the type of fluorescence quenching of the hPPARγ protein induced by DBDF2,6T was preliminary defined as a kind of static quenching [24]. In dynamic quenching, which was associated with diffusion, quenching constant of fluorescent material was increasing as the temperature increased. But from the Figure 3, the slope of the Stern–Volmer lines was decreased while increasing the temperature of the experimental system, which further confirmed the quenching mechanism of DBDF2,6T with the hPPARγ protein was static quenching.

3.3. Binding Constants and Binding Site Numbers

Equation (2) is the Lineweaver–Burk double-reciprocal equation, and (3) was deduced by (2) [25]: where KA is the binding constant and n is the number of independent binding sites. When were plotted against , a straight line could be drawn and is shown in Figure 4. The corresponding computing results are shown in Table 2.

In both the experimental temperatures (273 K and 310 K), the computing binding site numbers were near to 1, which meant that the hPPARγ-DBDF2,6T complexes were formed by protein hPPARγ and DBDF2,6T at the ratio approximately to 1 : 1. And the binding constants were of the order of magnitude of 103, which meant that the binding ability between them was pretty strong.

3.4. Thermodynamic Parameters and Interaction Forces

The interaction forces between small molecules and biomacromolecules were belonged to noncovalent force including hydrogen bond, van der Waals force, electrostatic attraction, and so on. The main acting force between hPPARγ protein and DBDF2,6T could be judged according to the thermodynamic parameters which were calculated based on the Van’t Hoff equation [26]. From the previous researches on interaction abilities, it was assumed that different proteins and compounds had different main acting force [27]. The thermodynamic parameters were calculated according to following equations:where is the gas constant, is the Gibbs free energy change, is the entropy change, is the enthalpy change, and is the Stern–Volmer quenching constant. could be considered as a constant when the temperature changed in small range. And the Ross law indicated that, if  > 0 and  > 0, the main acting force between small molecules and biomacromolecules would be hydrophobic force; if  < 0 and  < 0, it would be hydrogen bond and van der Waals force; and if  ≈ 0 and  > 0, it would be electrostatic force [28].

According to the computing results shown in Table 3, was 19.03 kJ/mol, was 3.92 J/mol·K and 3.90 J/mol·K correspondingly at 293 K and 310 K, and was −20.18 kJ/mol and −20.24 kJ/mol correspondingly at 293 K and 310 K. Based on the Ross law, the main acting force between hPPARγ protein and DBDF2,6T was hydrophobic force. In addition, the regulator effect of several kinds of interaction forces and relevant microenvironments were both responsible for the macroscopic consequence [29].

3.5. Conformational Change of hPPARγ Protein

The 3D fluorescence spectrum is shown in Figure 5 in the form of intensive contour map. The related data are shown in Table 4.

It can be observed from Figure 5 that two typical fluorescence peaks of proteins were located approximately at λem = 340 nm. In order to observe the peaks of fluorescent groups more clearly, the excitation wavelength range was set smaller than emission wavelength range, so the spectra of Rayleigh scattering were not available on the picture [19]. After the hPPARγ protein was incubated with DBDF2,6T, the location of both fluorescence peaks did not show a significant change, but the intensity of each peak was reduced at different degrees. From the 3D fluorescence spectra of the hPPARγ protein (Figure 5(a)), the intensity ratio of the big peak to the small one was 7.97 : 1, while after the protein was incubated with DBDF2,6T (Figure 5(b)), the value was changed to 8.00 : 1, and the DBDF2,6T showed a strong quenching effect on the big peak which was located at about 290/340 (λex/λem). The 3D fluorescence spectrum indicated a conformational change of the specific structures of the hPPARγ protein, which could further validate the interaction between the hPPARγ protein and DBDF2,6T [30].

3.6. Confirmation of the Interaction Ability Using DARTS Technique

To identify the binding targets for small molecules, the key advantage of DARTS method was no sample pretreatments such as labeling the ligand [17]. And the method was particularly useful when a compound had a lower affinity with the target, even the binding constant was in micromolar range [31]. To confirm the feasibility of the DARTS applying to the hPPARγ protein, a preliminary experiment had been performed to research the digestion effects of the protease on the hPPARγ protein. And the pronase was used for digestion because it had been proved to be more useful for DARTS than any other protease [16]. In preliminary experiment, the digestion effects of time of enzymolysis and the concentration of pronase had been investigated.

After electrophoresis and staining, the protein bands of each sample are shown in Figure 6. Compared with the DMSO control, under certain conditions, the antienzymatic hydrolysis ability of hPPARγ protein did exist and was closely related to the concentration of DBDF2,6T incubated with the protein. In 0.5 × 10−4 mol/L DBDF2,6T, the strongest antienzymatic hydrolysis ability of the hPPARγ protein would appear, and such ability could be weakened with the change in the concentration. Although DBDF2,6T had shown hydrolysis ability to the hPPARγ protein at relatively high concentration, the protective functions of DBDF2,6T to the hPPARγ protein still could be observed and existed concentration-effect relationships in some extent. Consequently, the interaction between hPPARγ protein and DBDF2,6T could be indirectly verified [32].

3.7. Confirmation of the Binding Ability Using Ultrafiltration-LC Technique

The results of ultrafiltration-LC experiment are shown in Figure 7. It can be observed from the chromatograms that both the sample and the negative control had an obvious peak at location about 8.75 min which belongs to compound DBDF2,6T. Significant signal enhancement of the peak of compound DBDF2,6T between the sample and the negative control indicated a specific binding between DBDF2,6T and the recombinant protein hPPARγ, while the signal of DBDF2,6T in the negative control was attributed to the nonspecific binding [33]. The big impurity peak was located at about 7 min attributed to the solution of recombination protein hPPARγ. Therefore, the experiment of ultrafiltration-LC did verify that compound DBDF2,6T could bind to pure protein hPPARγ directly in physiological environment [34].

3.8. Exploration of the Theoretical Binding Details Using Molecular Docking

In order to further understand the interaction between DBDF2,6T and hPPARγ protein, molecular docking was used to explore the theoretical binding details of them [35]. Among the docking of 12 conformers of DBDF2,6T to generate pocket of the hPPARγ protein, the highest total score was 7.14 and the corresponding crash score and polar score were −1.86 and 0.00, respectively, which meant that DBDF2,6T had a pretty strong affinity to hPPARγ protein, and such docking process was under a relatively comfortable level of molecules [36]. The generated pocket of hPPARγ is shown in Figure 8(a), and the hydrogen bond graph is shown in Figure 8(b). In conclusion, DBDF2,6T could theoretically bind to hPPARγ protein with pretty strong binding strength, and it could directly interact with SER342 and GLY284 of hPPARγ protein by hydrogen bond. The hydrogen bond lengths between DBDF2,6T and SER342 were 2.50 Å and 2.42 Å, and that between DBDF2,6T and GLY284 was 2.74 Å.

4. Conclusions

This study analyzed the interaction between the novel patent organotin compound DBDF2,6T and the hPPARγ protein under physiological condition with the methods of fluorescence quenching, 3D fluorescence, DARTS, ultrafiltration-LC, and computer molecular docking. According to the spectroscopic experimental data, DBDF2,6T could interact with the hPPARγ protein and formed a nonradiative ground-state complex of hPPARγ-DBDF2,6T, mainly through hydrophobic force. Such a reaction was spontaneous and could cause a conformational change of the hPPARγ protein. And the experiments of DARTS and ultrafiltration-LC preliminarily proved the possibility of DBDF2,6T to be an agonist compound to hPPARγ protein. Considering the anticancer activity of DBDF2,6T and various physiological functions performed by agonists of PPARγ protein, the conclusion could be drawn that DBDF2,6T had a possibility to interact with the hPPARγ protein as an agonist and finally inducing physiological effects such as anticancer activity. This work successfully revealed the interaction of DBDF2,6T with hPPARγ protein and established a feasible way to validate the agonist compounds for hPPARγ protein.

Abbreviations

DBDF2,6T:Bis-[2,6-difluoro-N-(hydroxyl-<κ>O)benzamidato-<κ>O]dibutylitin
hPPARγ:Human peroxisome proliferator-activated receptor gamma
3D fluorescence:Three-dimensional fluorescence
DARTS:Drug affinity responsive target stability
LC:Liquid chromatography.

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

Financial supports funded by the Shanxi Key Laboratory of Innovative Drug for the Treatment of Serious Diseases Basing on the Chronic Inflammation Shanxi University of Chinese Medicine, Program for the Natural Science Foundation of Shanxi Province (no. 2014011027-1), the Fund for Shanxi Key Subjects Construction (2016), the Program for Science and Technology of Shanxi Province (no. 20130313021-9), the Top Science and Technology Innovation Teams of Higher Learning Institutions of Shanxi Province, the Program for the Top Young and Middle-aged Innovative Talents of Higher Learning Institutions of Shanxi Province (2015), the Foundation of Young Academic Leader in Shanxi Province, the Project of Shanxi Key Laboratory for Innovative Drugs on Inflammation-Based Major Disease (2018), and the Project of Center of Comprehensive Development, Utilization and Innovation of Shanxi Medicine (2017-JYXT-18) are gratefully acknowledged.