Francisco Mainardi
Martins
a,
Daniele Cocco
Durigon
bc,
Otávio Augusto
Chaves
d,
Rosely Aparecida
Peralta
c,
Davi Fernando
Back
*a and
Hernán
Terenzi
*b
aLaboratory of Inorganic Materials, Department of Chemistry, CCNE, Federal University of Santa Maria, Santa Maria, RS 97105-900, Brazil. E-mail: davi.f.back@ufsm.br
bLaboratory of Structural Molecular Biology, Department of Biochemistry, Federal University of Santa Catarina, Florianópolis, SC 88040-900, Brazil. E-mail: hernan.terenzi@ufsc.br
cLaboratory of Bioinorganic and Crystallography, Department of Chemistry, Federal University of Santa Catarina, Florianópolis, SC 88040-900, Brazil
dCQC-IMS, Department of Chemistry, University of Coimbra, Coimbra, 3004-535, Portugal
First published on 26th September 2024
Research on metallodrugs with biological properties remains at the forefront, focusing on the development of compounds that interact non-covalently with deoxyribonucleic acid (DNA) and possess the ability to cleave the double helix strands of this biomacromolecule. In this context, vanadium(V) imine complexes featuring a fluorinated triphenylphosphonium group, (3-formyl-4-hydroxybenzyl)tris(4-fluorophenyl)phosphonium chloride ([AF]Cl), were studied for their targeting and accumulation in mitochondria, in addition to their interactions and ability to cleave DNA. The solid-state structures of complexes C1–C3 were elucidated using single-crystal X-ray diffraction and were characterized using vibrational techniques and elemental analysis, along with extensive characterization in solution. These studies revealed that the complexes contain cis-dioxidovanadium(V) species and are zwitterionic species. It was shown that C1–C3 can interact with and are capable of cleaving plasmid DNA through oxidative mechanisms without the need for photoinduction. When potential interactions with bovine serum albumin were analyzed, it was revealed that interactions in the order of 104 M−1 (Stern–Volmer quenching constant, KSV) were observed. Additionally, in silico molecular docking studies showed that C1–C3 can preferentially interact with the minor grooves of DNA and with domain IB (site III) of bovine and human serum albumins.
The cleavage of both strands of the DNA double helix is medically significant because it causes damage to cells, resulting in cell death.9 Therefore, the study of complexes capable of cleaving the DNA chain at specific sites has garnered interest due to their potential contributions to genomics and their applications in photodynamic therapy and chemotherapy against neoplasms.10,11
Mechanistically, two major types of cleavage can occur, namely hydrolytic or induced by reactive oxygen species (ROS) and the last one being categorized into oxidative or photocleavage. Generally speaking, synthetic metallonucleases that cleave DNA through hydrolytic mechanisms typically feature metal ions with inert redox properties and strong Lewis acidity via a reversible reaction, whereas those with ROS-induced mechanisms usually occur via irreversible reactions involving redox-active metal ions or photosensitizers.9,12,13 The interaction and cleavage of DNA by redox-active metal complexes are crucial for developing antineoplastic therapeutic compounds.14 In this context, certain cations have garnered increased attention, such as manganese(III),15 iron(III),15 cobalt(II),16 ruthenium(II),17,18 and vanadium(IV/V).19,20
Vanadium is a redox-active transition metal abundantly dissolved in marine waters,21 exhibiting ions in variable oxidation states,22 with III, IV, and V being the most common and biologically relevant.23,24 When tetra- or pentavalent, vanadium metal centers are often stabilized in their oxido ([VO]2+ or [VO]3+) or cis-dioxidovanadium (cis-[VO2]+) forms due to their high oxophilicities and Lewis acidity in these oxidation states.25,26 Vanadium ions perform various biological functions,27–29 such as constituting nitrogenase enzymes30,31 and haloperoxidases.32,33
Several types of vanadium-containing compounds have been studied for their potential therapeutic properties,34–36 including insulin-enhancing properties37 (notably, complexes bis(maltolate)oxidovanadium(IV), BMOV, and bis(ethylmaltolate)oxidovanadium(IV), BEOV),38,39 antibacterial,40,41 antiviral,42–44 antiparasitic,45,46 and antineoplastic activities47–49 (e.g., bis(4,7-dimethyl-1,10-phenanthroline)sulfateoxidovanadium(IV), metvan).50 In recent years, several vanadium(V) complexes containing different imine ligands have been synthesized and evaluated for their interactions with biomacromolecules, including serum albumins and DNA,51–53 cytotoxic potential,35,54–56 and/or photoinduced DNA cleavage.57–59
Recently, the focus on vanadium complexes has extended to specifically target mitochondria,60 offering a potential strategy for effectively inducing tumor cell death. Mitochondria play key roles in cellular energy production and in the execution of cell death through apoptotic mechanisms. Therefore, the production of reactive oxygen species, mitochondrial changes, upregulation of anti-apoptotic factors, and the acquired inability of the tumor suppressor protein p53 to activate the pro-apoptotic proteins Bak and Bax connect this organelle to the formation and growth of neoplasms.61
For mitochondrial targeting of a compound of interest, derivatization with triphenylphosphonium groups (+TPP) has been a strategy employed by the pharmaceutical industry for years,62 as these groups preferentially accumulate in the anionic mitochondrial matrix due to their cationic charge.61 These groups, widely regarded as the “gold standard” for mitochondrial vectors, also exhibit high lipophilicity, exceptional stability under physiological conditions, and facile conjugation with the compound of interest.63 Consequently, the potential for derivatizing +TPP groups with 4-fluoro substituents is assessed.
Fluorination of molecules is a prevalent step in modern medicinal chemistry aimed at developing compounds with enhanced metabolic stability, bioavailability, and protein interactions. Fluorine atoms can modulate lipophilicity and facilitate weak hydrogen bonds and electrostatic interactions in organofluorine compounds.64 The unique physical and chemical characteristics of fluorinated compounds arise from the small atomic radius and high electronegativity of fluorine.65 Significant advances in fluorinated drug development with varying properties have occurred in recent years, including antifungals and antivirals, as well as in antineoplastic therapies.66
This study reports the synthesis of (3-formyl-4-hydroxybenzyl)tris(4-fluorophenyl)phosphonium chloride ([AF]Cl), along with its use as a building block in the synthesis, characterization, and biological applications of three new zwitterionic complexes of cis-dioxidovanadium(v) (C1–C3) containing imine hydrazone ligands ([H2L1]Cl–[H3L3]Cl). The structures of complexes C1–C3 were elucidated by single-crystal X-ray diffraction and characterized by various techniques in solution and in the solid state. Concerning biologically relevant applications, these complexes were evaluated for their interactions and cleavage potentials with DNA and bovine serum albumin through gel electrophoresis assays, fluorescence and UV-vis spectroscopies, and also molecular docking studies.
The nuclear magnetic resonance (NMR) spectra of nucleus 1H (600 MHz), 13C (151 MHz), 19F (565 MHz), 31P (243 MHz), and 51V (158 MHz) were recorded using a Bruker AVANCE III 600 MHz spectrometer at 25 °C. Dimethyl sulfoxide (DMSO-d6) and chloroform (CDCl3) were used as deuterated solvents. Tetramethylsilane (TMS, Me4Si) was used as the internal reference, while trichlorideoxidovanadium(V) in deuterated benzene (90% VOCl3/C6D6 solution), triphenylphosphane (PPh3), and trifluorotoluene (CF3Ph) were used as external references. Chemical shifts were reported in parts per million (δ, ppm) and were referenced to residual DMSO-d6 or TMS (CDCl3) peaks. The multiplicities were expressed as follows: s, singlet; d, doublet; t, triplet, and m, multiplet.
The Fourier-transform infrared (FT-IR) spectra in transmission mode were recorded using a Bruker Vertex 70 spectrometer equipped with an ATR diamond accessory in the 4000–30 cm−1 region with 64 scans and 4 cm−1 resolution. On the other hand, the Fourier-transform Raman (FT-Raman) spectra of the complexes (crystalline materials) in absorbance mode were recorded using a Bruker SENTERRA confocal Raman microscope. The scattered beam was collected using an Olympus 20× objective and the spectra were recorded in the 3500–50 cm−1 range with a laser line of 785 nm. The spectra were recorded at a power of 25 mW with 6 coadditions of 10 seconds each.
The mass spectra of the inorganic complexes were obtained using an Amazon ion trap mass spectrometer via electrospray ionization (ESI-MS). The analysis was carried out in an ultrapure acetonitrile solution at a concentration of 500 ppb at 180 μL min−1. The capillary temperature was maintained at 180 to 200 °C and a voltage of −400 to −500 V. The simulated spectra were calculated using the mMass software.67
Finally, the electronic spectra in the UV-Vis region were obtained using a Varian Cary 50 BIO UV-vis spectrometer coupled to a thermostatic bath at 25 ± 1 °C, using N,N-dimethylformamide (DMF) as the solvent and a cuvette with 1.0 cm optical path.
The gels were imaged and quantitatively analysed by densitometry using the ImageJ software. Given that ethidium bromide intercalates to a lesser extent with the supercoiled form of the plasmid DNA, a correction factor of 1.47 was applied when assessing this specific conformation of plasmid DNA.71
To determine the DNA interaction and cleavage preferences by the complexes, external agents were incorporated into the reaction mixtures. This included lithium perchlorate (LiClO4, to augment the ionic strength of the reaction media) and minor-groove binders (4′,6-diamidino-2-phenylindole, DAPI, and netropsin, NET; 50 μM) along with a DNA major-groove binder (methyl green; 50 μM).
Studies on the interactions between the complexes and calf thymus DNA (ct-DNA, stock solution 1.0 mM) were conducted using a Varian Cary 50 Bio UV-vis spectrophotometer connected to a thermostatic bath maintained at 25 ± 1 °C. Complexes C1–C3 were examined in a DMF/buffer (TRIS 0.5 mM) solution (5% v/v) at pH 7.40 and an ionic strength (I) of 0.5 mM (NaCl).
For the DNA structure, an 8 Å radius around the two main possible binding sites (major and minor grooves) was explored,75 while for BSA and HSA, an 8 Å radius around the three main binding pockets (subdomains IIA, IIIA and IB – known as sites I, II, and III, respectively)76,77 was defined and explored. For all biomacromolecules, the standard ChemPLP was used as the scoring function. Protein–ligand interaction profiler (PLIP) webserver (https://plip-tool.biotec.tu-dresden.de/plip-web/plip/index)78 was used for the identification of protein–ligand interactions and the 3D-figures were generated using the PyMOL Molecular Graphics System 1.0 level software (Delano Scientific LLC software, Schrödinger, New York, NY, USA).79
The complexes feature an iminic hydrazone ligand unity arising from the condensation of (3-formyl-4-hydroxybenzyl)tris(4-fluorophenyl) chloride and the respective hydrazide, resulting in a tridentate chelating cavity that comprises phenol and carbonyl oxygens, and an iminic nitrogen as potential donor atoms for coordination.80 However, the vanadium metal centers in C1–C3 (mononuclear complexes) are coordinated to the respective ligands through the anionic forms of such oxygen atoms. One is in the phenolate form (ranging from 1.898(5) to 1.915(2) Å) and the other is in the enolate form (ranging from 1.973(5) to 1.982(4) Å).
The enolate-form oxygen atom is generated by deprotonation of the amide nitrogen of the hydrazone moiety, as indicated by the bond lengths of C–N(amide) (ranging from 1.292(7) to 1.311(3) Å) and C–O(enolate) (ranging from 1.288(6) to 1.305(3) Å), and the C–N(amide)–N(imine) bond angles81 (ranging from 108.1(2)° to 109.0(5)°) in C1–C3. Moreover, a neutral iminic nitrogen (ranging from 2.138(5) to 2.156(2) Å) acts as a third site in ligands characterized by an O,N,O tridentate dianionic motif. These values for C1–C3, along with the bond length values of their coordination spheres, bear similarities to those of selected complexes in the literature (Tables S2 and S3 in the ESI†).82,83
In addition to the ligand donor atoms, the vanadium metal centers in C1–C3 coordinate with two oxido ligands (ranging from 1.588(6) to 1.646(2) Å), adopting the form of cis-dioxidovanadium (cis-[VO2]+), a distinctive species of vanadium(V). Consequently, C1–C3 exhibit pentacoordinate coordination spheres. The geometries of C1–C3 have been determined based on the values of the structural index parameter τ, calculated as (β − α)/60, with β and α being, respectively, the largest and second largest bond angle values of the coordination spheres. The values of this parameter range from 0, representing a perfect square-based pyramid, to 1, denoting a perfect trigonal bipyramid. Non-integer values within this range indicate distorted geometries.84 The β and α angles (°), calculated τ values (dimensionless), and coordination geometries of C1–C3 are summarized in Table 1, with comparisons to some dioxidovanadium(V) complexes from the literature presented in Table S4 in the ESI.†82,83
Complex | β anglea (°) | α angleb (°) | τ | Coordination geometriesd |
---|---|---|---|---|
a Highest bond angle value. b Second highest bond angle value. c Structural index parameter τ (dimensionless) calculated as τ = (β − α)/60. d SBP = square based pyramidal geometry and TB = trigonal bipyramidal geometry. e Two independent units of the complex in the asymmetric unit. | ||||
C1 | 147.59(8) | 145.77(11) | 0.03 | Slightly distorted SBP |
C2 | 155.46(17)/153.50(18) | 125.6(2)/136.2(2) | 0.50/0.29 | Distorted between SBP and TB/distorted SBP |
C3 | 154.9(2) | 129.8(3) | 0.42 | Distorted between SBP and TB |
Furthermore, given the sum of cationic and anionic charges, alongside the absence of counterions, these three complexes are classified as neutral zwitterionic complexes,85 and this is due to the presence of the delocalized cationic charge intrinsic to the ligand of the tris(4-fluorophenyl)phosphonium group. This characteristic distinguishes them from the complexes reported by Maurya et al.86,87 and Dash et al.,57 wherein protonation and coordination of the pyridinic nitrogen atoms occur.
The stability of C1–C3 in DMF/buffer – (TRIS 0.5 mM) 5% v/v; pH 7.40, I = 0.5 mM (NaCl) at 25 °C (Fig. S46–S48 of the ESI†) – was examined over 24 h using UV-vis absorption spectroscopy. The absorption spectra of all three compounds remained unchanged, indicating their stability and structural integrity during the investigated period.
Initially, the impact of pH on the reaction medium was assessed (Fig. 2). For all three complexes, a notable increase in DNA cleavage was observed at pH 6.0 and 7.0, after 8 hours of reaction at a complex concentration of 500 μM. Additionally, we noted that C1–C3 exhibited similar activity. Consequently, pH 7.0 was selected for further assays.
Notably, unlike similar compounds described by Dash and coworkers,58,59 which involve mononuclear dioxovanadium(V) complexes containing aroylhydrazone derivative ligands and have shown DNA cleavage activity through photoinduced reactions, the compounds described herein, C1–C3, can cleave DNA in the absence of light and under mild pH conditions, specifically at pH 7.0.
The effect of complex concentrations was also investigated. For C1–C3, as shown in Fig. S52–S54 in the ESI,† it was possible to observe that the complexes exhibit a capacity to cleave DNA in a concentration-dependent manner at 500–30 μM. This demonstrates that increasing the concentration of the compound also increases its DNase activity.
The cleavage reaction was carried out in both aerobic and anaerobic environments, in a glovebox using argon as an inert gas. As shown in Fig. S55 in ESI,† the absence of oxygen significantly decreased the DNase activity of compounds C1–C3. Therefore, it can be inferred that the cleavage activity is dependent on oxidative processes. To obtain information on the preference for cleavage sites, the cleavage was investigated using DNA groove blockers, and the results of which are shown in Fig. S56 in the ESI.† Generally, complexes C1–C3 showed an increase in DNA cleavage activity in the presence of minor groove binders, DAPI and NET, suggesting higher affinities of the vanadium complexes with the minor groove. Conversely, methyl green, a major groove binder, did not affect the DNase activity.
To comprehend the forms of interaction between the complexes and DNA, the impact of the ionic strength within the reaction medium was initially explored. We found that as the concentration of lithium perchlorate increased from 0 to 50 mM for all three complexes (Fig. 3 for C1 and Fig. S57 and S58 in the ESI† for C2 and C3, respectively); the cleavage of DNA by the complexes also exhibited a gradual decrease for all three complexes. Considering the negative charge of the DNA molecule at pH 7.0, this finding suggests that electrostatic interactions, probably the positive charge centered on the tris(4-fluorophenyl)phosphonium group, between the complex and DNA might play a role in their binding and subsequent cleavage.
The nature of the binding interaction between complexes C1–C3 and ct-DNA was also examined using the UV-Vis absorption titration method, observing changes in absorbance and shifts in the wavelength of the complex spectra with increasing amounts of ct-DNA. The results are presented in Fig. 4 for C3 and Fig. S59 and S60 in the ESI† for C1 and C2, respectively. As a result, significant hypochromic changes were observed for all three compounds, along with slight bathochromic shifts.
In the absorption spectrum of a compound in the presence of DNA, changes such as hypochromism or hyperchromism and a wavelength shift (red or blue) frequently occur. Hypochromism and bathochromism (red shift) in UV-vis spectra are typically linked to the binding of a metal complex to DNA through intercalation involving π-stacking, van der Waals, and hydrophobic interactions between the metal complex and DNA.91
BSA fluorescence is primarily governed by the aromatic amino acid residues, namely tyrosine, phenylalanine, and tryptophan (W134 and W213), with excitation at 280 nm. It is well-known that changes in the fluorescence behavior of albumin result from the interactions and/or changes in the microenvironment of these fluorophores caused by the quencher.93
Fig. 5 depicts the protein's maximum emission band at approximately 340 nm when excited at 280 nm. Observations from the steady-state fluorescence spectra revealed that an increase in the complexes’ concentration (Fig. 5 for C1 and C3, and Fig. S61 in the ESI† for C2) led to a significant decrease in the albumin's fluorescence intensity. This evidence suggests the binding between the complexes and the fluorophores of BSA.
Fig. 5 Fluorescence spectra of BSA (2 μM) quenched by complexes C1 (up) and C3 (down) in the concentration range of 0–20 μM. Inset: The Stern–Volmer plot. |
An isosbestic point was also observed along with a new fluorescence band, indicating an increased fluorescence intensity in the 425–525 nm range. This phenomenon suggests the formation of an exciplex with BSA, which can be considered further evidence of more complex interactions between albumin and vanadium complexes.94,95 This phenomenon was more evident for C1 and C2 than for C3, indicating the impact of the nicotinic and isonicotinic moieties on the spectroscopic signature of the binding of the inorganic complexes under study. Rio et al.95 previously reported the capacity of albumin in forming exciplex in the presence of heterocyclic compounds (e.g., 1,2,4-triazolo-3-thiones), highlighting the impact of the sulfide and sulfoxide groups.
To further explore these interactions, the Stern–Volmer quenching (KSV) and binding (Kb) constants of the complexes were calculated and are summarized in Table 2. The Stern–Volmer equation is F0/F = 1 + KSV[Q] = 1 + kqτ0[Q], where F0 and F represent the fluorescence intensity without and with the complexes, respectively, KSV is the Stern–Volmer quenching constant, [Q] is the concentration of the complexes, kq is the bimolecular quenching rate, and τ0 is the fluorescence lifetime of BSA (6.06 × 10−9 s).96 This calculation used corrected fluorescence data to account for the effect of dilution. The Kb values and the number of binding sites (n) were obtained from the Scatchard equation using double logarithm plots through the equation log(F0 − F/F) = nlog[Q] + logKb.
Complex | K SV × 104 (M−1) | k q × 1012 (s−1 M−1) | K b × 104 (M−1) | n | ΔGe (kJ mol−1) |
---|---|---|---|---|---|
a Stern–Volmer quenching constant. b Bimolecular quenching rate constant with BSA (τ = 6.06 ns). c Double logarithm binding constant. d Number of binding sites. e Gibbs free-energy values (R = 8.314 J mol−1 and T = 298.15 K). | |||||
C1 | 4.48 ± 0.09 | 7.34 | 6.03 ± 0.13 | 1.03 ± 0.03 | −27.28 |
C2 | 3.88 ± 0.07 | 6.36 | 7.76 ± 0.12 | 1.06 ± 0.02 | −27.91 |
C3 | 7.62 ± 0.26 | 12.49 | 36.31 ± 0.18 | 1.22 ± 0.04 | −31.73 |
In Table 2, one can observe that the number of interaction sites for C1–C3 is preferably one. The results (both KSV and Kb) demonstrated that the binding constants of C1 and C2 are very similar, which is reasonable given their structural similarities. Conversely, compound C3 exhibited the highest binding constant among the three complexes, reflecting the structural effects of incorporating naphthalenol as a pendant arm on the ligand. The values observed are in the same range or higher than those for similar [VVO2]-complexes reported in the literature, indicating that C1–C3 might be distributed in the human bloodstream by the interaction with albumin.26,52,59,97
The bimolecular quenching rate constant (kq) of all the complexes are larger than the diffusion rate constant in water (kdiff ≈ 7.40 × 109 M−1 s−1 at 298 K, following the Smoluchowski–Stokes–Einstein theory98), indicating a ground-state association.99 The negative sign for the Gibb's free energy (ΔG), calculated based on the Van’t Hoff equation (ΔG = −RTlnKb, where R is the gas constant that is equal to 8.314 J mol−1 and T is the temperature in Kelvin), indicates that the binding process is spontaneous at 298.15 K.100
Complex | DNA | Bovine serum albumin | Human serum albumin | |||||
---|---|---|---|---|---|---|---|---|
Major groove | Minor groove | Site I | Site II | Site III | Site I | Site II | Site III | |
C1 | 26.1 | 63.1 | 44.9 | 40.5 | 88.6 | 60.4 | 11.9 | 84.2 |
C2 | 32.1 | 65.1 | 41.4 | 25.2 | 84.9 | 59.0 | 7.76 | 87.8 |
C3 | 25.6 | 92.0 | 47.8 | 16.8 | 91.1 | 57.0 | 10.1 | 98.7 |
Thus, according to the computational results, for DNA which has two main binding sites (minor and major grooves), the [VVO2]-complexes fit preferentially into minor grooves than into major grooves (Table 3), agreeing with the experimental competitive binding assays. The same trend was previously in silico reported for the interaction between DNA and VV-complexes derived from pyridoxal/salicylaldehyde.52 On the other hand, BSA has three main binding pockets for both endogenous and exogenous compounds: subdomain IIA (site I) located in a hydrophobic binding pocket, subdomain IIIA (site II) also located in a hydrophobic binding pocket, and subdomain IB (site III) located on the surface of the protein.101,102
All binding sites were able to accommodate the [VVO2]-complexes; however, site III is the main pocket (Table 3) with docking score values close to 2-fold higher than for sites I and II. Interestingly, subdomain IB was not reported as the main binding site to VV-complexes derived from pyridoxal or salicylaldehyde,52,103 indicating dependence on the chemical nature of the ligands complexed with vanadium(V) species.
To verify if C1–C3 will also have the same binding trend obtained from bovine to the human analogue, in silico calculations for HSA were also carried out and the results are summarized in Table 3. In this case, subdomain IB was also identified as the main binding site, reinforcing site III as the feasible cavity to the interaction of the [VVO2]-complexes. However, after a deep analysis of the docking score values, it is possible to note that BSA and HSA had practically the same score in site III, with differences in the other binding sites. This was expected due to the highest superposition of the 3D structure of the two proteins into subdomain IB compared to the other subdomains IIA and IIIA (Fig. 7A) with a global root mean square deviation (RMSD) value of 1.634 dimensionless. Interestingly, albumin-C3 had the highest docking score value than albumin-C1 or albumin-C2, following the same experimental binding trend.
In the case of DNA, molecular docking results suggested van der Waals interactions as the main intermolecular force responsible for the interaction with [VVO2]-complexes C1–C3 (Table S5 in the ESI†), and despite the high steric volume of the inorganic complexes under study, they had good fit capacity in the double helix of DNA (Fig. 6) with the positive charge centered in the tris(4-fluorophenyl)phosphonium group of C1–C3 close to the negative charged phosphate moieties of the DNA structure, agreeing with the experimental trend.
On the other hand, for BSA, interactions via hydrophobic, π-cation, and hydrogen bonding were suggested as the main binding forces, while for HSA hydrophobic interaction was detected and hydrogen bonding with π-stacking interaction was detected only for HSA-C3 (Tables S6, S7 in the ESI† and Fig. 7). Interestingly, Dias et al.104 reported that the interaction of the mixed species cis-[VO(carrier)2-(HSA)] (carrier is picolinato, maltolato, or 3-hydroxy-1,2-dimethyl-4-pyridinone) involves hydrogen bonding and/or hydrophobic interactions with the protein surface, corroborating with the data described above. Finally, in silico calculations suggested that the naphthalenol moiety of C3 is more buried into the protein cavity than the nicotinic and isonicotinic moieties of C1 and C2, respectively, reinforcing the obtained highest experimental binding capacity for C3 to albumin.
Fig. 7 (A) Superposition of the 3D structures of BSA and HSA (PDB codes 4F5S and 3JRY, respectively). The best docking pose for the interaction: (B) BSA:C1, (C) BSA:C2, (D) BSA:C3, (E) HSA:C1, (F) HSA:C2, and (G) HSA:C3 in the three main binding sites with a zoom representation for subdomain IB. Selected amino acid residues for BSA and HSA, C1, C2, and C3, are in stick representation in green, duck blue, beige, pink, and brown, respectively. Elements’ color: oxygen, nitrogen, fluorine, phosphorus, and vanadium(V) are in red, dark blue, cyan, orange, and gray, respectively. For better interpretation, all hydrogen atoms were omitted. |
Both experimental and in silico results indicated that the chemical nature of the pendant arms impacts not only the binding capacity to the biomacromolecules but also the steady-state fluorescence signal. For example, for C1, the isosbestic point detected an exciplex formation in the presence of albumin (Fig. 5). Molecular docking calculations suggested that it occurs due to the capacity of the nicotinic moiety from C1 to interact with the fluorophore Tyr-161 via π-stacking force (Fig. 7B). In other words, since albumin was excited at 280 nm, the amino acid residue Tyr in its excited state might transfer energy to the nicotinic moiety of C1. In the case of C2 and C3, a weak and non-exciplex formation was detected, respectively (Fig. S61 and Fig. 5B, respectively). The in silico results suggested that the isonicotinic and naphthalenol moieties from C2 and C3, respectively, interact differently in the subdomain IB of albumin, i.e., the isonicotinic moiety interacts with Tyr-161 via hydrophobic forces (Fig. 7C). In contrast, naphthalenol does not interact with aromatic amino acid residues (Fig. 7D).
Finally, despite the number of connecting points for BSA-C3 is lower than that for BSA-C1 and BSA-C2 (Fig. 7), the inorganic complex C3 presented the highest binding affinity (around 6- and 4.7-fold higher than C1 and C2, respectively, Table 2) probably due to the capacity of the naphthalenol moiety in interaction via π-cation forces with the charge center of Arg-186 driven by the hydrogen bond between the hydroxyl group of the naphthalenol moiety and Glu-183.
In vitro tests with DNA revealed that complexes C1–C3 are capable of interacting with and cleaving DNA. Competitive binding studies combined with in silico modeling indicated a preference for all the complexes to bind to the minor groove of DNA, significant for the specificity of DNA cleavage. All the new vanadium(V) compounds exhibited significant interactions with the protein, as demonstrated by the binding constant values from emission spectroscopy, indicating a substantial interaction between BSA and the complexes.
Experimental and molecular docking studies showed that compound C3, which incorporates naphthalenol as a pendant arm on the ligand, interacts more strongly with the protein than C1 and C2. This heightened interaction is attributed to π-stacking, enhancing the attraction between the amino acid (the naphthalenol moiety of C3 is more buried into the protein cavity than the nicotinic and isonicotinic moieties of C1 and C2, respectively) residues and the compound.
Footnote |
† Electronic supplementary information (ESI) available. CCDC 2367124–2367126 for C1–C3. For ESI and crystallographic data in CIF or other electronic format see DOI: https://doi.org/10.1039/d4nj03087g |
This journal is © The Royal Society of Chemistry and the Centre National de la Recherche Scientifique 2024 |