Optimizing the optical and biological properties of 6-(1H-benzimidazole)-2-naphthalenol as a fluorescent probe for the detection of thiophenols: a theoretical study

The study presents the influence of structure modulation by introduction of selected donor and acceptor substituents on the properties of 6-(1H-benzimidazole)-2(2,4-dinitrobenzenesulfonate)naphthalene used in thiophenol identification. The presence of –OH and –OR groups enhances the non-linear optics (NLO) response of the marker. The –NO2 substituent maximizes the non-linear response and increases the amount of transferred charge and the charge-transfer distance. The introduction of the –OH, –NO2 and –CN groups into the marker structure significantly improves the solubility and optical availability. The –NO2 group however contributes to mutagenicity and carcinogenicity. The –OH and –OR groups can be successfully used in bioimaging to detect specific molecules containing the –SH group in their structure. At the same time, the –OR group minimizes the energy barrier necessary to break the bond between the chromophore and the linker. The paper also includes a comparison of optical and biological properties of structures before and after identification of thiophenols.


Introduction
Thiophenols are widely used in the synthesis of many industrial products, primarily polymers, pesticides and pharmaceutical products. 1,2 Despite their widespread use, they are highly toxic and the average lethal dose (LC 50 ) for sh is only 0.01-0.04 mM. 3 Also, the human body is very exposed to the toxic effects of these substances due to inhalation and rapid absorption through the skin. Exposure of the human body to these compounds results in many serious illnesses, including polyposis, nausea, vomiting, muscle weakness and in special cases even death. 4,5 Against this background, it is particularly important to search for new solutions for the easy and quick detection of thiophenols, especially in pharmaceuticals and in biological material. In this context, uorescent probes, which are dyes relatively easy to manufacture with a wide spectrum of Stokes shis, have become of practical importance. It is easy to modulate the structure of such compounds in order to optimize desirable photophysical features and they can be used in in vivo and in vitro tests. For detecting thiophenol derivatives, uorescent probes are designed in such a way that the uorophore is linked to an appropriate quencher unit to transfer the transformed molecule to another group with no quenching capability or cleaved from the uorophore permanently. In this way, uorescence is recovered and the last derivative is responsible for marker uorescence. An example is 2,4-dinitrophenylsulphonyl (DNS). [6][7][8][9] Many probes tailored for the detection of thiophenols are designed on the basis of the reaction between thiophenols and the sulfamide unit (RNH-SO 2 R 0 ). However, thiophenol detection via such probes is relatively long and takes on average 20-30 min. [10][11][12] The search for more effective solutions leads to the use of a sulfonyl ester (RO-SO 2 R 0 ), which signicantly reduces the detection time. 13 Based on this second strategy of uorophore construction, the W. Lv 14 group has synthesized 6-(1H-benzimidazole)-2naphtalenol (BIN) as a satisfactory probe for the detection of thiophenols. It has been shown to have excellent photostability, good biocompatibility and a large Stokes' shi in aqueous solutions. By connecting the BIN to the sulfonyl ester group, the 6-(1H-benzimidazole)-2(2,4-dinitrobenzenesulfonate)naphtalene (BIN-T) probe was synthesized. In this construction, BIN is a uorophore, while the 2,4-dinitrobenzenesulfonyl group is both a thiophenol recognition group and a uorescence quenching unit (Scheme 1). For the marker constructed in this way, the average detection time of thiophenols was about 10 min.
Due to the possibility of using BIN for the detection of thiophenols in pharmaceutical products and in biological material, it was decided to check which factors enhance the linear and non-linear optical response and improve biological properties. For this purpose, electron-donating and electronwithdrawing substituents have been attached to the naphthalene part. In this way, six derivatives (Scheme 1) with different photochemical and biological properties were obtained. The presented research results allow to answer the question how the probe structure should be modulated to obtain a marker with the desired photochemical properties required for molecular imaging. The study of marker structure modulation by changing substituents has a very important cognitive purpose. On the one hand, it allows achieving the desired optical and biological characteristics. On the other hand, it allows the creation of a uorescent probe operating in a specic absorption and uorescence range, depending on the research requirements.

Synthesis path
The method of synthesis of the tested derivatives is presented in Scheme 2.

Computational details
All geometrical parameters of investigated molecules in their ground (S GS ) and excited states (S CT ) were calculated using density functional theory (DFT) approach implemented in Gaussian 09 program package 15 with TIGHT threshold option and PBE0/6-311++G(d,p) basis set. In order to verify that all the structures correspond to the minima on the potential energy surface, an analysis of Hessians was performed. The electronic properties were characterized by computations of the vertical absorption and emission spectra, which were obtained using the time-dependent density functional theory (TDDFT/PBE0) 16 and by including the state-specic (SS) corrected linear response (cLR) approach. 17 Due to the high compatibility of theoretical and experimental data, [18][19][20][21] all spectroscopic calculations were performed using PBE0 functional.
For the best consideration of the solvent impact on the uorescence spectra, the ground state should be calculated with non-equilibrium solvation. 22,23 This was taken into account by including the state-specic (SS) corrected linear response (cLR) approach 24 to the theoretical calculations. In the SS approach the solvent dynamic polarizations are determined by the difference of the electron densities of the initial and nal states. [25][26][27] The dipole moments and polarities of the charge-transfer state (S CT ) were evaluated by numerical differentiation of the excitation energies (E) in the presence of an electric eld F of 0.001 a.u. strength: where i stands for the Cartesian component of the dipole moment difference and g is ground state (S g ). The isotropic average polarizability (hai), polarizability anisotropy (Da) and rst-order hyperpolarizability (b vec ) were determined based on the Gaussian 09 program and dened as: The density differences were obtained at the PBE0/6-311++G(d,p) level and are represented with a contour threshold of 0.02 a.u. In these graphs, the blue (purple) zones indicate density decrease (increase) upon electronic transition. The charge transfer parameters, namely the charge-transfer distance (D CT ) and the amount of transferred charge (q CT ), have been determined following a Le Bahers' procedure. 28 The solvent effect on the linear and nonlinear optical properties has been taken into account using the Integral Equation Formalism for the Polarizable Continuum Model (IEF-PCM). 29,30 Scheme 1 Detection pathway of thiophenol derivatives. 14 Experimentally, the two-photon absorption (TPA) can be obtained by the dissipation of the incident light, which for a single beam 2PA experiment is twice the transition rate. In this case, the two-photon cross-section of the degenerate process is written as: [31][32][33] where a is a ne structure constant, u is the frequency of absorbed photons (assuming one source of photons), G F is the broadening of the nal state (F) due to its nite lifetime and g(u) provides the spectral line prole, which oen is assumed to be a d-function and hd OF i is the two-photon transition probability for the transition from the ground state to a nal state.
In the case of a molecule absorbing two photons of the same energy in isotropic media, the degenerate hd OF i in an isotropic medium using a linearly polarized laser beam given by: 34 d OF ¼ 1 15 In this equation, S OF ij is the second-order transition moment given by: where ħu 1 + ħu 2 should satisfy the resonance condition and h0|z 1 $m _ i |Ki stands for the transition moment between electronic states |0i and |Ki, respectively. z is the vector dening polarization of photons. To describe the two-photon allowed states, the quadratic response functions formalism 35,36 within the DFT framework was used as implemented in the DALTON 2011 program. 37,38 Solvent effects were taken into account with the self-consistent reaction eld (SCRF) model. All the 2PA calculations were carried out employing the CAM-B3LYP functional and the 6-311++G(d,p) basis set.
The biological activities were simulated using a combination of the 3D/4D QSAR BiS/MC and CoCon algorithms. [39][40][41][42][43] The binding properties of considered dyes were studied by performing series of AutoDock 4.2 (ref. [44][45][46] and AutoDock Vina 47 simulations. For each active complex comprising Human Serum Albuminum (HSA), their crystal structures were taken from PDB ID: 1AO6. 48 The cubic grid box with the dimensions of 16Å and a grid spacing of 1Å was set up in such way that the reactive -NH 2 groups of lysine and -SH of cysteine were at its centre. In order to identify appropriate binding energy and conformation of compounds, the Lamarckian genetic algorithm was employed. For each atom of the receptor molecule Gasteiger charges were calculated. The investigation of the binding site was performed using a united-atom scoring function. For each amino acid the docking simulations were performed tenfold.

Physicochemical properties
For all investigated molecules, the charge-transfer (CT) excitation corresponds to the HOMO / LUMO transition ( Fig. 1 and 2). In the case of the BIN-T probe, HOMO electrons are located on the 6-(1H-benzimidazole)-2-naphtalenol part, while LUMO on the 2,4-dinitrobenzenesulfonyl group. The attachment of substituents, as well as their change, does not affect the position of these frontier orbitals. Reaction with thiophenol and obtaining BIN uorophore signicantly changes this distribution. HOMO electrons accumulate throughout the entire molecule, while LUMO mainly accumulates on the naphthalene part. This is particularly apparent for the BIN-NO 2 derivative. For BIN-T derivatives, there are slight differences in the case of energy separation between HOMO-LUMO orbital (E GAP ) - Table  SI1. † Attachment of the -NH 2 and -NHR substituents reduces E GAP by 0.55 eV on average. In turn, the presence of -OH reduces it by 0.2 eV and the presence of -OR brings this value closer to the BIN-T level. The presence of electron-acceptor substituents increases the E GAP value by more than 0.3 eV. Aer thiophenol detection (Table SI2 †), the E GAP value increases signicantly for water DE BIN-BIN-T GAP ¼ 1.5569 eV. In this case, the inuence of substituents is slightly different. The presence of -NH 2 reduces it by 0.2 eV while -NHR does not change this value. The -OH and -OR substituents slightly increase E GAP and the electronacceptor ones decrease it. For both BIN and BIN-T derivatives, no signicant effect of solvent presence on the EGAP value is observed. In general, for BIN-T derivatives with increasing polarity of the medium E GAP increases (DE H2Oe1;4Dx GAP ¼ 0.1970 eV) and for BIN decreases (DE H2Oe1;4Dx GAP ¼ 0.0160 eV). In addition, a signicant decrease (by more than 1 eV) in chemical hardness (h) is observed for BIN-T derivatives. High electronegativity (c) value suggests an easy formation of covalent bonds during various chemical processes.
In order to assess the sites of potential electrophilic and nucleophilic attack, the molecular electrostatic potential (MEP) analysis was performed ( Fig. 1 and 2). In the case of BIN-T, the place with the highest electronegativity and thus the most exposed to nucleophilic attack (positive, blue zones) is the nitrogen atom with hydrogen attached in the imidazole part. Benzene in the 2,4-dinitrobenzenesulfonyl part is also a region that may undergo nucleophilic reactions. Substitution with subsequent substituents does not change the location of these places. Only for BIN-T-NH 2 and BIN-T-NHR positive zones are also observed on the nitrogen atoms of substituents. Sites for electrophilic attack (red and yellow, negative zones) are mainly the second of the nitrogen atoms of the imidazole part and the oxygen atom in the sulfonylester group. The presence of substituents forces the negative zones to shi into the regions of the connected groups. Among the considered derivatives, the highest risk of electrophilic attack is experienced by BIN-T-CN, where the charge on the nitrogen atom of the -CN group is as much as À0.09737 a.u. Aer the detection of thiophenols and obtaining the BIN structure, the nitrogen atom of the imidazole group is still the place with the highest risk of nucleophilic attack. In addition, an important place for such a reaction is the hydrogen atom of the hydroxyl group. The place for the electrophilic attack remains the second nitrogen atom of the imidazole ring and substitution with subsequent substituents does not change these regions. Similar to BIN-T, negative zones also appear on substituents. It can also be seen that the attachment of -NH 2 slightly reduces the charge on the nitrogen atom, reducing the risk of an electrophilic reaction.
The studies have been devoted to many reections on spectral properties, corresponding to the HOMO / LUMO photoexcitation (p-p* transitions). In order to estimate contributions from other orbitals and determine the nature of electronic states, the density variation upon photoexcitation (Dr(r)) was computed for the rst electronic transitions, which is graphically depicted in Fig. 1 and 2. For BIN-T, the gure indicates that the depletion zones (blue) are located on 2-(naphthalen-2-yl)-1H-benzo[d]imidazole, while in contrast, growth zones (purple) on 2,4-dinitrobenzenosulfonyl part. The presence and change of donor-acceptor substituents does not change the position of these zones. In the case of BIN derivatives, the zones of electron density decay are visible throughout the entire molecule, and the increase zones mainly on naphthalene. Also in this case, the substituent change does not move these zones. On the other hand, the polarity of the environment affects the parameters describing Dr(r) ( Table SI3 †). For both BIN-T and BIN derivatives, the amount of the transferred charge drops as a function of solvent polarity. In addition, the value of q CT is higher for BIN-T, and the Dq BIN-T-BIN CT difference is 0.9 a.u. The same analogy applies to other derivatives. In addition, the introduction of the -NO 2 substituent increases the value of q CT and Dq BINeNO2eBIN CT is 0.5 a.u. The polarity of the environment also signicantly affects the charge-transfer distance and decreases it from 1.682Å in 1,4-Dx to 0.920Å in water for BIN. In turn, for BIN-T, this quantity increases from 7.503Å to 8.157Å. At the same time, BIN-T derivatives have a higher D CT value and DD BIN- The presented analysis indicates the chargetransfer nature of BIN and BIN-T derivatives. It also conrms the contributions from the HOMO-LUMO transition. At the same time, it indicates the possibility of the presence of contributions from other orbitals. Therefore, additional low intensity peaks can be seen in absorption and uorescence spectra.
The BIN-T molecule is characterized by good solubility in all discussed media. DG solv increases with increasing solvent polarity, going from À17.62 kcal mol À1 in 1.4Dx to À20.37 kcal mol À1 in water, with the highest value being reached in MeCN (À29.06 kcal mol À1 ). While similar behavior is observed for other molecules, the introduction of -NHR, -OH, -NO 2 and -CN signicantly improves the solubility of the compound. The presence of the remaining substituents slightly increases the DG solv value by reducing solubility. Similar conclusions can be drawn aer analyzing BIN derivatives. More importantly, the change in structure aer thiophenol identication results in a signicant reduction in solubility and for BIN DG BIN-T-BIN solv in 1,4Dx is 5.59 kcal mol À1 and in water 1.85 kcal mol À1 . It should also be clearly stated that in terms of DG solv the desired substituents will be -OH, -NO 2 and -CN, the presence of which signicantly improves the solubility of the compound and thus its optical availability.   (À23.86 nm).

Linear optical properties
Considering the effect of the solvent, negative solvatochromism is observed for BIN-T derivatives and the maximum absorption is shied in the direction of shorter wavelengths with increasing solvent polarity. This effect is observed regardless of the nature of the attached substituent, both for vertical values and the ones determined based on the cLR model. Also, for all BIN derivatives non-monotonic solvatochromism is observed. The excitation energy values E EX increase as a function of the polarity of the medium, however the transition from MeCN to DMSO is followed by a slight decrease in E EX . The only exception from the above results is the BIN-NO 2 derivative, for which positive solvatochromism is observed.
The interconnectedness of the analysis of spectroscopic parameters with MEP indicates the possibility of specic interactions in the solute-solvent system. A non-monotonic increase in excitation energy indicates on larger polarization and better S g stabilization. However, this is not consistent with the polarity of the excited state (Dm CT-g , eqn (1)). Hypsochromic shi, as the effect of the increase in environmental polarity, should result in the m g > m CT relationship. In contrast, for both BIN-T and BIN derivatives there is an inverse relationship (Table SI6 †), which is characteristic of positive solvatochromism. In any case, the polarity of the excited state also changes non-monotonously. In water, for BIN-T derivatives, the highest Dm CT-g value is characteristic for markers with electron-withdrawing substituents: Dm  (18.20 nm). In addition, in line with experimental relationships, for BIN-T uorescence intensity is negligible (Fig. 3). Identication of thiophenols and detachment of 2,4-dinitrobenzenesulfonyl group signicantly increases the intensity of l FL . At the same time, it can be observed that with the exception of BIN and BIN-NHR, a second, less intense uorescence band appears shied towards shorter wavelengths. In general, the used substituents signicantly shi the maximum l FL towards longer wavelengths, and only the -NHR substituent results in hypochromic shi relative to BIN

Nonlinear opical properties
The polarizability (Da, eqn (3)) and rst hyperpolarizability (b vec , eqn (4)) of molecules irradiated with intense laser light giving the electric eld is the subject of many research in terms of understanding various nonlinear optical properties (NLO). In particular, these studies include the interrelationship of NLO with the electronic structure to design new multifunctional uorescence probes. The calculated values of Da and b vec are collected in Tables 1 and SI8 Since the non-linear response will be derived from the structure obtained aer thiophenol identi-cation, the factor enhancing the NLO response will be the presence of the -NH 2 group and can be efficient in Second Harmonic Generation (SHG).
Tables S1 and SI9 † show the value of two-photon absorption cross section (TPA, hd OF i eqn (6)). In general, it is difficult to decisively conclude about the effect of the solvent on the hd OF i value. Non-monotonic solvatochromism is observed for BIN-T, whereas negative solvatochromism is observed for BIN.
Positive solvatochromism is revealed in the values describing BIN-NH 2 and BIN-T-OR. In addition, the values for BIN-T derivatives are on average 3.000 a.u. higher than for BIN. However, in order to compare the calculated values of the TPA with those determined experimentally (s (2) OF , eqn (5)), the relation (5) was used. In this equation, the broadening of the nal state due to its nite lifetime 0.25 eV was assumed. The effect of the solvent on the s (2) OF value remains the same as for hd OF i. More importantly, although BIN-T derivatives have a higher s (2) OF value, none of the analyzed markers meet the requirements for use in two-photon imaging. BIN-T-CN (25.33 GM) has the highest value, while BIN (0.13 GM) has the lowest one. Thus, these values indicate a very poor NLO response signal. In the case of BIN-T derivatives, only the presence of the -CN substituent is a factor increasing the non-linear response. The remaining substituents reduce the s (2) OF value, with the highest minimization observed for -NHR. For BIN derivatives, each substituent increases the two-photon absorption cross section value. The largest increase is observed for BIN-NH 2 , -NHR and -CN. Based on this analysis, to obtain a probe with a high TPA value, used as tools in real-time dynamic in vivo and in vitro research, an additional chromophore group intensifying the NLO response should be added to the structure of the analyzed dyes.

AutoDock simulations
BIN-T is a probe for detecting thiophenols in solutions. As mentioned earlier (Scheme 1), the identication reaction occurs by 6-(1H-benzimidazole) À2 (2,4-dinitrobenzenesulfonate) naphthalene reacting with thiophenol sulfur. Therefore, this marker can also be an alternative as a uorescent probe used in medical imaging. In this case, protein conjugation occurs by reaction with cysteine. Human Serum Albuminum (HSA) was selected for this study. For BIN-T derivatives, the presence of a substituent or its change does not alter the site of active binding to the protein. In each case, conjugation to HSA occurs via CYS448 (Fig. 4). The BIN-T-CN marker has the highest affinity for this active center, for which binding energy (DG b ) is À9.9 kcal mol À1 (Table SI10 †). This active binding site has an inhibition constant (K i ) of 1.00 mM. The BIN-T-CN is inserted into the aromatic cage formed by CYS448, ASN295, ARG218, LYS444, PRO447, TYR452 and VAL455. In this active site, no pp* interactions occur and the system is not stabilized by the presence of a H-bond. BIN-T-OH and BIN-T have slightly lower affinity for HSA and DG b ¼ À9.2 kcal mol À1 . Inhibition constant is 1.57 mM and 2.53 mM, respectively. The active impact cavity here is: LYS199, LYS195, LEU198, VAL455, ASN295, ARG218 and LYS444. In the case of BIN-T-OH, TRP214 also interacts here and in BIN-T -ASP451. Also in this case no hydrogen bonds are observed and no p-p* interactions. For other probes DG b is À9.1 kcal mol À1 . Only for BIN-T-NHR a hydrogen bond is formed between the oxygen of the sulfone group and the hydrogen of LYS444. In addition, the lowest K i value of 0.63 mM was recorded for BIN-T-NO 2 .
During conjugation, the O-S bond raptures and then a new one forms with the protein. The speed of biocomplex creation will therefore also be affected by the energy barrier (DE) necessary to overcome. According to Fig. SI1, † the lowest DE value when moving the oxygen atom away from the sulfur atom at a distance of 2.5Å is observed for BIN-T-OR and is 14.5 kcal mol À1 . A slightly higher value is observed for BIN-T-OH, for which DE ¼ 16.87 kcal mol À1 . The strongest bond, and thus hindering the conjugation, is the O-S in BIN-CN, where DE ¼ 26.63 kcal mol À1 .

Biological activities
The BIN-T derivatives are described by relatively good bioavailability (log P > 5). Only for BIN-T-NHR this value drops slightly and amounts to 4.70 AE 0.25. This suggests good permeability through cell membranes and achieving adequate concentration at the site of interaction with thiophenols. In turn, BIN derivatives are described with a log P value in the range of 3-4. The exceptions are BIN-NH 2 and BIN-NHR, for which log P < 3. The log BCF value calculated for BIN-T in the range from À4.2 to À4.5 indicates the lack of bioaccumulation in the tissues of living organisms and the ease of excretion with urine. Aer identifying thiophenols, log BCF decreases by one unit (in the range from À3.0 to À3.3). Despite this, BIN derivatives should not bioaccumulate aer fullling their optical role. In addition, both BIN-T and BIN derivatives are characterized by high metabolism by CYP450-2D6 and CYP450-3A4 (probability above 80%) (Tables SI11 and SI12 †). This indicates that both forms of markers will be rapidly excreted from tissues without interacting with other biomolecules and drugs.
The calculated oral toxicity value is LD 50 > 1500 mg kg À1 . Therefore BIN-T and BIN derivatives should be classied in class 4 in terms of the degree of toxicity and can be considered as practically nontoxic for humans. While the presence and change of a substituent does not affect the toxicity of presented molecules, it does, however, affect other toxicological parameters (Fig. 5). None of the derivatives shows hepatotoxicity and cytotoxicity. The presence of the NO 2 substituent causes carcinogenicity (probability P ¼ 58%) and mutagenicity (P ¼ 72%). The -OR group may affect the occurrence of immunotoxicity (P ¼ 66%). Immunotoxicity can elicit the reference BIN marker (P ¼ 57%). Other probes should not cause any toxic effects. Thus, the introduction of a substituent into the BIN structure eliminates the possible toxic effects of the molecule. In addition, both BIN-T and BIN derivatives have other biological activities suggesting their potential use in other areas of medicine (Tables SI11 and  SI12 †). First of all, BIN-T derivatives are characterized by high antioxidant activity (P> 48%), in particular BIN-T-CN (P ¼ 90%). For BIN derivatives probability of this activity drops to zero, except for BIN-NHR (P ¼ 27%) and BIN-OR (P ¼ 39%). In turn, for BIN-T derivatives no antibacterial activity is observed, while it can be attributed to BIN derivatives, in particular BIN-OH (P ¼ 79%), BIN-NO 2 (P ¼ 79%) and BIN-CN (P ¼ 80%). The following tables indicate the occurrence of many other activities, such as: alpha-radioprotector, analgetic, anti-psychotic activity diazepine site, anti-tumor alkylic, anti-tumor antimiotic, HIV1-proteaze inhibitory activities and many others. At the same time, the inuence of subsequent substituents on the maximization or reduction of these activities is demonstrated there.

Conclusion
The manuscript presents previously unavailable observations on the inuence of donor-acceptor substituents (D/A) on optical parameters and biological activity of the uorescent probe used to identify thiophenols. This type of research is the rst for the derivatives in question and is valuable information for researchers synthesizing new compounds on the benzimidazole backbone. It has been demonstrated that the used D/A substituents do not signicantly affect the position of HOMO-LUMO orbitals and only slightly affect the EGAP value. They also do not change the sites for nucleophilic and electrophilic attack and density depletion/increment zones upon photoexcitation. However, it has been shown that the -NO 2 substituent increases the amount of transferred charge and the charge-transfer distance. Studies have shown that the introduction of the -OH, -NO 2 and -CN groups into the marker structure signicantly improves the solubility and thus the optical availability of the probe. The used substituents cause the solvatochromic shi of the absorption and uorescence maxima with respect to the experimentally tested marker, and the solvatochromic relationships depend on the structure before and aer identication of thiophenols. These observations is extremely important for scientists synthesizing and studying benzimidazole derivatives. Taking into account the Stokes shi value, it can be concluded that the presence of the -OH and -OR substituents maximizes the linear response and increasing detection of the marker. However, the reducing factor will be the presence of -NO 2 . On the other hand, the -NO 2 substituent maximizes the non-linear response. Nevertheless, the presence of -OR and -OH also enhances NLO values relative to pure dye. The tested derivatives are not used in two-photon imaging. However, this observation may be a clue to search for other substituents enhancing the twophoton story. The analyzed structures show high affinity for HSA, and a change in the D/A substituent slightly reduces binding energy. Therefore, the presented compounds can be used with great success in in vivo and in vitro tests as uorescent probes for detecting specic biomolecules containing the -SH group. Therefore, the considered derivatives constitute valuable compounds for the pharmaceutical and medical market. At the same time, the -OR group minimizes the energy barrier necessary to break the bond between the chromophore and the linker, which speeds up and facilitates accessibility during conjugation This journal is © The Royal Society of Chemistry 2020 RSC Adv., 2020, 10, 24374-24385 | 24383 with protein. All derivatives are characterized by good permeability through cell membranes and no bioaccumulation. A very valuable observation is the fact that, toxicity studies excluded the attractiveness of using the -NO 2 substituent in in vivo and in vitro studies due to the potential for carcinogenicity and mutagenicity. By modulating the appropriate D/A substituents on the BIN structure, many valuable drugs can be obtained in the treatment of cancer, HIV, oxidative stress regulation, etc. Therefore, the conducted research indicates the high value of these derivatives in medicine. Considering all the above analyzes, it should be clearly stated that the presence of the -OR and -OH substituents signicantly maximizes the linear and non-linear optical response of the molecule used as a uorescent probe in medical imaging. Lethal dose, determination of the toxicity of the substance

Conflicts of interest
The authors declare that there are no potential nancial or non-nancial conicts of interest.