Amine-Reactive BODIPY Dye: Spectral Properties and Application for Protein Labeling

A boron-dipyrromethene (BODIPY) derivative reactive towards amino groups of proteins (NHS-Ph-BODIPY) was synthesized. Spectroscopic and photophysical properties of amine-reactive NHS-Ph-BODIPY and its non-reactive precursor (COOH-Ph-BODIPY) in a number of organic solvents were investigated. Both fluorescent dyes were characterized by green absorption (521–532 nm) and fluorescence (538–552 nm) and medium molar absorption coefficients (46,500–118,500 M−1·cm−1) and fluorescence quantum yields (0.32 – 0.73). Solvent polarizability and dipolarity were found to play a crucial role in solvent effects on COOH-Ph-BODIPY and NHS-Ph-BODIPY absorption and emission bands maxima. Quantum-chemical calculations were used to show why solvent polarizability and dipolarity are important as well as to understand how the nature of the substituent affects spectroscopic properties of the fluorescent dyes. NHS-Ph-BODIPY was used for fluorescent labeling of a number of proteins. Conjugation of NHS-Ph-BODIPY with bovine serum albumin (BSA) resulted in bathochromic shifts of absorption and emission bands and noticeable fluorescence quenching (about 1.5 times). It was demonstrated that the sensitivity of BSA detection with NHS-Ph-BODIPY was up to eight times higher than with Coomassie brilliant blue while the sensitivity of PII-like protein PotN (PotN) detection with NHS-Ph-BODIPY and Coomassie brilliant blue was almost the same. On the basis of the molecular docking results, the most probable binding sites of NHS-Ph-BODIPY in BSA and PotN and the corresponding binding free energies were estimated.


Introduction
Detection, tracking, and imaging of proteins are considerable challenges for modern science. One of the ways of tackling this challenge is to label proteins with fluorescent dyes. Fluorescent labels allow imaging various biochemical processes including protein interactions with other biomolecules, protein localization, protein dynamics, enzyme activity, etc. [1,2].
Currently, three main classes of fluorescent labels are widely used [3]: organic dyes, fluorescent proteins, and quantum dots. Each of them has its own advantages and disadvantages. In this regard, the development of new fluorophores that meet all the requirements of bioimaging continues unabated.
Bright members of organic dyes are boron-dipyrromethene (BODIPY) fluorophores and their derivatives [13]. This family of dyes possesses a number of striking properties, such as great stability, high molar absorption coefficients, excellent fluorescence quantum yields, relatively small Stokes shifts, and sharp excitation and emission peaks [14]. Moreover, its ease of chemical modification opens up great possibilities for a synthesis of new fluorescent compounds with predetermined properties.
A rather new promising area in BODIPYs chemistry is the development of dyes reactive towards principal functional groups of proteins [15][16][17][18][19][20][21][22][23][24]. Such fluorescent dyes possess certain reactive substituents capable of covalent binding with proteins, thereby providing their effective labeling. Unfortunately, most of existing reactive BODIPY derivatives have several issues with solubility in water and/or photophysical characteristics. Therefore, the search for reactive fluorescent dyes suitable for conjugation with proteins continues.
In this work, we synthesized a BODIPY derivative reactive towards amino groups of proteins. A comparative analysis of spectroscopic and photophysical properties as well as solvatochromic behavior of the synthesized reactive fluorescent dye and its non-reactive precursor was carried out by means of ultraviolet-visible spectroscopy and steady-state and time-resolved fluorescence spectroscopy supported by quantum chemical calculations. The investigated reactive fluorescent dye was used for labeling of a range of proteins. The sum of the experimental and theoretical results suggests that the synthesized amine-reactive fluorescent dye is a promising fluorescent label for molecular biology and biotechnology.
Bright members of organic dyes are boron-dipyrromethene (BODIPY) fluorophores and their derivatives [13]. This family of dyes possesses a number of striking properties, such as great stability, high molar absorption coefficients, excellent fluorescence quantum yields, relatively small Stokes shifts, and sharp excitation and emission peaks [14]. Moreover, its ease of chemical modification opens up great possibilities for a synthesis of new fluorescent compounds with predetermined properties.
A rather new promising area in BODIPYs chemistry is the development of dyes reactive towards principal functional groups of proteins [15][16][17][18][19][20][21][22][23][24]. Such fluorescent dyes possess certain reactive substituents capable of covalent binding with proteins, thereby providing their effective labeling. Unfortunately, most of existing reactive BODIPY derivatives have several issues with solubility in water and/or photophysical characteristics. Therefore, the search for reactive fluorescent dyes suitable for conjugation with proteins continues.
In this work, we synthesized a BODIPY derivative reactive towards amino groups of proteins. A comparative analysis of spectroscopic and photophysical properties as well as solvatochromic behavior of the synthesized reactive fluorescent dye and its nonreactive precursor was carried out by means of ultraviolet-visible spectroscopy and steady-state and time-resolved fluorescence spectroscopy supported by quantum chemical calculations. The investigated reactive fluorescent dye was used for labeling of a range of proteins. The sum of the experimental and theoretical results suggests that the synthesized amine-reactive fluorescent dye is a promising fluorescent label for molecular biology and biotechnology.

UV-Vis and Fluorescence Spectroscopy Experiments of BODIPYs
Absorption spectra, emission spectra, fluorescence decay curves as well as a number of photophysical characteristics of synthesized fluorescent dyes COOH-Ph-BODIPY and NHS-Ph-BODIPY in a range of organic solvents of various nature were obtained (Figures 1 and S1, Table 1).

UV-Vis and Fluorescence Spectroscopy Experiments of BODIPYs
Absorption spectra, emission spectra, fluorescence decay curves as well as a number of photophysical characteristics of synthesized fluorescent dyes COOH-Ph-BODIPY and NHS-Ph-BODIPY in a range of organic solvents of various nature were obtained (Figures 1 and S1, Table 1).   In Table 1, λ abs max is the maximum absorption wavelength, λ em max is the maximum emission wavelength, ∆ν is the Stokes shift, and ε is the molar absorption coefficient at the maximum absorption wavelength. The asterisk (*) indicates the average fluorescence lifetime τ calculated using Equation (1) [26]: where α is the amplitude. The absorption and emission spectra of COOH-Ph-BODIPY and NHS-Ph-BODIPY are typical for this class of dyes ( Figure 1). The absorption spectra gave two bands in the ranges of 521-532 and 375-385 nm due to S 0 -S 1 and S 0 -S 2 electron transitions, respectively, while the emission spectra gave one band in the region of 538-552 nm due to S 1 -S 0 electron transitions (Section 2.3). The absorption and emission bands maxima of NHS-Ph-BODIPY are 2-7 nm red shifted compared with the spectra of COOH-Ph-BODIPY. The fluorescent dyes are characterized by medium molar absorption coefficients and fluorescence quantum yields which change non-linearly from COOH-Ph-BODIPY to NHS-Ph-BODIPY. The fluorescence lifetimes of NHS-Ph-BODIPY are shorter than those of COOH-Ph-BODIPY, while the radiative and non-radiative rate constants change non-linearly.
Solvent effects on spectral characteristics of COOH-Ph-BODIPY and NHS-Ph-BODIPY, namely the maximum absorption ν abs (max) and emission ν em (max) wavenumbers and Stokes shift ∆ν, were analyzed in terms of the Kamlet-Taft [27], Catalán [28], and Lippert-Mataga [29] equations. Tables 2 and S1 list the estimated regression coefficients y 0 , a − d and their standard errors as well as coefficients of determination R 2 for the multiple linear regression analyses of the maximum absorption ν abs (max) and emission ν em (max) wavenumbers and Stokes shifts ∆ν of COOH-Ph-BODIPY and NHS-Ph-BODIPY according to Equations (5) and (S2) for the solvents chosen. The analyses of the ν abs (max) and ν em (max) data within the Kamlet-Taft model, in which solvent polarizability and dipolarity effects are combined in the single parameter π*, show poor fits with the R 2 = 0.44-0.66 and large standard errors on the estimated regression coefficients a α , b β , and c π * as goodness-of-fit criteria (Table S1). In contrast, the analyses of the same data within the Catalán model, in which solvent polarizability and dipolarity effects are separated by the two parameters SP and SdP, give good to perfect fits with R 2 = 0.75-0.96 (Table 2). To visualize the goodness-of-fits of ν abs (max) and ν em (max) as functions of the Catalán solvent parameters {SA, SB, SP, SdP}, the plots of ν abs (max) and ν em (max) of COOH-Ph-BODIPY and NHS-Ph-BODIPY calculated according to Equation (5) using the estimated values of y 0 , a SA , b SB , c SP , and d SdP vs. the corresponding experimental ν abs (max) and ν em (max) values were drawn ( Figure 2).  Figure 2).  Table 1).
The Catalán model further reveals solvent properties that are mainly responsible for the observed spectral shifts. The very large (negative) estimated regression coefficients cSP compared to the aSA, bSB, and dSdP ones as well as the relatively large standard errors on the aSA, bSB, and dSdP compared to those on the cSP (Table 2) indicate that the changes of νabs(max) and νem(max) may primarily reflect changes in polarizability of the environment of COOH-Ph-BODIPY and NHS-Ph-BODIPY [30]. The other supporting evidence comes from the multiple linear regression analyses of the νabs(max) and νem(max) according to Equation (5) (Table S2). The analyses give good fits with an R 2 equal to 0.91, 0.89, and 0.72 and 0.75, 0.69, and 0.67 for νabs(max) and νem(max) of COOH-Ph-BODIPY, respectively, as well as 0.91, 0.92, and 0.78 and 0.85, 0.82, and 0.90 for νabs(max) and νem(max) of NHS-Ph-BODIPY, respectively. It is significant to note that the lowest coefficients of determination (indicated in bold) were obtained in case of {SA, SB, SP} as independent variables (except for νem(max) of NHS-Ph-BODIPY), indicating that solvent dipolarity should not be neglected as an influencing factor. In contrast, the multiple linear regression analyses of the νabs(max) and νem(max), according to Equation (5) with {SA, SB, SdP} as independent variables (Table S2), show poor fits with an R 2 equal to 0.54 and 0.52 for νabs(max) and νem(max) of COOH-Ph-BODIPY, respectively, as well as 0.38 and 0.39 for νabs(max) and νem(max) of NHS-Ph-BODIPY, respectively. Thus, the crucial solvent property affecting the maximum absorption νabs(max) and emission νem(max) wavenumbers of COOH-Ph-BODIPY and NHS-Ph-BODIPY is solvent polarizability with a small contribution of solvent dipolarity.  Table 1). The Catalán model further reveals solvent properties that are mainly responsible for the observed spectral shifts. The very large (negative) estimated regression coefficients c SP compared to the a SA , b SB , and d SdP ones as well as the relatively large standard errors on the a SA , b SB , and d SdP compared to those on the c SP (Table 2) indicate that the changes of ν abs (max) and ν em (max) may primarily reflect changes in polarizability of the environment of COOH-Ph-BODIPY and NHS-Ph-BODIPY [30]. The other supporting evidence comes from the multiple linear regression analyses of the ν abs (max) and ν em (max) according to Equation (5) (Table S2). The analyses give good fits with an R 2 equal to 0.91, 0.89, and 0.72 and 0.75, 0.69, and 0.67 for ν abs (max) and ν em (max) of COOH-Ph-BODIPY, respectively, as well as 0.91, 0.92, and 0.78 and 0.85, 0.82, and 0.90 for ν abs (max) and ν em (max) of NHS-Ph-BODIPY, respectively. It is significant to note that the lowest coefficients of determination (indicated in bold) were obtained in case of {SA, SB, SP} as independent variables (except for ν em (max) of NHS-Ph-BODIPY), indicating that solvent dipolarity should not be neglected as an influencing factor. In contrast, the multiple linear regression analyses of the ν abs (max) and ν em (max), according to Equation (5) with {SA, SB, SdP} as independent variables (Table S2), show poor fits with an R 2 equal to 0.54 and 0.52 for ν abs (max) and ν em (max) of COOH-Ph-BODIPY, respectively, as well as 0.38 and 0.39 for ν abs (max) and ν em (max) of NHS-Ph-BODIPY, respectively. Thus, the crucial solvent property affecting the maximum absorption ν abs (max) and emission ν em (max) wavenumbers of COOH-Ph-BODIPY and NHS-Ph-BODIPY is solvent polarizability with a small contribution of solvent dipolarity.
The Lippert plots of the Stokes shifts ∆ν of COOH-Ph-BODIPY and NHS-Ph-BODIPY vs. orientation polarizabilities ∆f of the solvents chosen are presented in Figure S2. There are extremely poor linear relationships between ∆ν and ∆f with an R 2 equal to 0.03 and 0.05, respectively. The small slopes (−68 and 56 cm −1 , respectively) imply that the dipole moments of COOH-Ph-BODIPY and NHS-Ph-BODIPY do not change noticeably between the ground and excited states [31]. The other supporting evidence comes from the multiple linear regression analyses of the ∆ν within the Kamlet-Taft and Catalán models (Tables S1 and 2), both of which show poor fits with the R 2 equal to 0.31 and 0.29 for COOH-Ph-BODIPY, respectively, as well as 0.61 and 0.70 for NHS-Ph-BODIPY, respectively. It is interesting to note that quantum chemical calculations gave another explanation of the abovementioned trends (Section 2.3).

Quantum Chemical Calculations of BODIPYs
In order to gain more insight into the abovementioned relations, a series of quantum chemical calculations was performed. First, a semiempirical metadynamic screening revealed the existence of four different conformers for both COOH-Ph-BODIPY and NHS-Ph-BODIPY, of which only two were found to be non-redundant according to further DFT optimization. Those conformations differ only in the relative orientation of ethyl groups, which have limited impact on spectroscopic properties. As was later confirmed, all conformers have approximately the same excitation energies ( Figure S3). It is, thus, reasonable to expect that existence of those conformers does not lead to any complications and further analysis can be performed only on the most thermodynamically stable conformer ( Figure 3). The Lippert plots of the Stokes shifts Δν of COOH-Ph-BODIPY and NHS-Ph-BODIPY vs. orientation polarizabilities Δf of the solvents chosen are presented in Figure  S2. There are extremely poor linear relationships between Δν and Δf with an R 2 equal to 0.03 and 0.05, respectively. The small slopes (−68 and 56 cm −1 , respectively) imply that the dipole moments of COOH-Ph-BODIPY and NHS-Ph-BODIPY do not change noticeably between the ground and excited states [31]. The other supporting evidence comes from the multiple linear regression analyses of the Δν within the Kamlet-Taft and Catalán models (Tables S1 and 2), both of which show poor fits with the R 2 equal to 0.31 and 0.29 for COOH-Ph-BODIPY, respectively, as well as 0.61 and 0.70 for NHS-Ph-BODIPY, respectively. It is interesting to note that quantum chemical calculations gave another explanation of the abovementioned trends (Section 2.3).

Quantum Chemical Calculations of BODIPYs
In order to gain more insight into the abovementioned relations, a series of quantum chemical calculations was performed. First, a semiempirical metadynamic screening revealed the existence of four different conformers for both COOH-Ph-BODIPY and NHS-Ph-BODIPY, of which only two were found to be non-redundant according to further DFT optimization. Those conformations differ only in the relative orientation of ethyl groups, which have limited impact on spectroscopic properties. As was later confirmed, all conformers have approximately the same excitation energies ( Figure S3). It is, thus, reasonable to expect that existence of those conformers does not lead to any complications and further analysis can be performed only on the most thermodynamically stable conformer ( Figure 3). Firstly, we have obtained vibrationally resolved absorption and emission spectra of the studied molecules. Calculations were performed using the VG-FC model. As shown in the literature [32], neglect of the excited state geometry relaxation is not only a reasonable, but sometimes crucially important (for calculation convergence) assumption. The comparison of the ground and excited state geometries ( Figure 3) shows that although there is a noticeable change in structure between S0 and S1 (reduction of C7-C8-C1Ph-C2Ph dihedral angle leading to improved conjugation between BODIPY and C8substituent), this change is consistent between the two studied dyes, so the general trends for relevant properties should be preserved.
As can be seen from the VG-FC absorption and emission spectra (Figure 4), one of the preserved trends is the presence of a reasonably slight bathochromic shift resulting from the addition of succinimidyl group to BODIPY core. The presented theoretical spectra also give the same overall shape and Stokes shifts seen in experimental data ( Figure 1). Firstly, we have obtained vibrationally resolved absorption and emission spectra of the studied molecules. Calculations were performed using the VG-FC model. As shown in the literature [32], neglect of the excited state geometry relaxation is not only a reasonable, but sometimes crucially important (for calculation convergence) assumption. The comparison of the ground and excited state geometries ( Figure 3) shows that although there is a noticeable change in structure between S 0 and S 1 (reduction of C 7 -C 8 -C 1Ph -C 2Ph dihedral angle leading to improved conjugation between BODIPY and C 8 -substituent), this change is consistent between the two studied dyes, so the general trends for relevant properties should be preserved.
As can be seen from the VG-FC absorption and emission spectra (Figure 4), one of the preserved trends is the presence of a reasonably slight bathochromic shift resulting from the addition of succinimidyl group to BODIPY core. The presented theoretical spectra also give the same overall shape and Stokes shifts seen in experimental data ( Figure 1).  Next, we have compared the results of the multiple linear regression analysis (Section 2.2) with the data obtained from the quantum chemical calculations (Table 3). In particular, relatively large (~ 3.6 D) ground state dipole moments of both COOH-Ph-BODIPY and NHS-Ph-BODIPY as well as large (~ 1.5 D) dipole moment changes upon excitation show why such Catalán solvent parameters as polarizability SP and dipolarity SdP dominate the abovementioned regressions. Moreover, the poor applicability of the Lippert relations may be explained not by a small dipole moment change upon excitation but by an indifference of the change towards solvent nature. Indeed, for both COOH-Ph-BODIPY and NHS-Ph-BODIPY, the dipole moment change is about 0.71 a.u. and for both non-polar n-hexane and slightly polar DMSO, the dipole moment magnitude change is about 0.68 a.u. Table 3. Ground S0 and first excited S1 total dipole moments ||µGS|| and ||µTD|| of COOH-Ph-BODIPY and NHS-Ph-BODIPY, difference between their absolute values, and absolute value of their difference in n-hexane and DMSO. To gain a better understanding of the nature of this excitation, we have opted to compute the partitioned transition density matrix [33]. This technique has been used to great effect to explain the similarities and differences between compounds [34,35]. Moreover, recent papers [36] have even compared different DFT functionals in terms of the quality of their transition density matrices.

Compound
We have defined five common units for COOH-Ph-BODIPY and NHS-Ph-BODIPY-two pyrrole rings (L-pyr and R-pyr), coordination center (CC), C8-carbon at- Next, we have compared the results of the multiple linear regression analysis (Section 2.2) with the data obtained from the quantum chemical calculations (Table 3). In particular, relatively large (~3.6 D) ground state dipole moments of both COOH-Ph-BODIPY and NHS-Ph-BODIPY as well as large (~1.5 D) dipole moment changes upon excitation show why such Catalán solvent parameters as polarizability SP and dipolarity SdP dominate the abovementioned regressions. Moreover, the poor applicability of the Lippert relations may be explained not by a small dipole moment change upon excitation but by an indifference of the change towards solvent nature. Indeed, for both COOH-Ph-BODIPY and NHS-Ph-BODIPY, the dipole moment change is about 0.71 a.u. and for both non-polar n-hexane and slightly polar DMSO, the dipole moment magnitude change is about 0.68 a.u. Table 3. Ground S 0 and first excited S 1 total dipole moments ||µ GS || and ||µ TD || of COOH-Ph-BODIPY and NHS-Ph-BODIPY, difference between their absolute values, and absolute value of their difference in n-hexane and DMSO.

Compound
Solvent ||µ GS ||, a.u. ||µ TD ||, a.u. ||µ GS || − ||µ TD ||, a.u. ||µ GS − µ TD ||, a.u. To gain a better understanding of the nature of this excitation, we have opted to compute the partitioned transition density matrix [33]. This technique has been used to great effect to explain the similarities and differences between compounds [34,35]. Moreover, recent papers [36] have even compared different DFT functionals in terms of the quality of their transition density matrices.
As can be seen from the visual representation of the matrices (Figure 5b,c), the first singlet excitation is localized mostly on the pyrrole rings and has a local character, although there is also a significant electron transfer between the rings, which, nevertheless, is symmetric. The change in the dipole moment arises due to the density transfer between the pyrrole rings as well as C8-carbon atom and its substituent. However, the position of the phenyl ring relative to dipyrrin is almost perpendicular due to the existence of the methyl groups in 1-and 7-positions of BODIPY core. As such, steric factors around this position should prevent any solvent molecules from getting close enough to provide a significant local impact on electron density via polarization. Comparing two matrices, it becomes clear that the nature of the phenyl substituent does not play a strong role when it comes to the first excited state-transition density on the succinimidyl group, as well as density change between it and other units, are small. As such, as will be shown below (Sections 2.4 and 2.5), any alterations in the absorption and emission bands positions of COOH-Ph-BODIPY and NHS-Ph-BODIPY upon conjugation with proteins will come as a result of intermolecular interactions, such as polarization due to nearby amino acids or loss of polarization induced by solvent.

Fluorescent Labeling of Proteins
The synthesized amine-reactive fluorescent dye NHS-Ph-BODIPY was used for labeling of a range of proteins, such as bovine serum albumin (BSA), PII-like protein PotN (PotN), PotA subunit of polyamine ABC transporter (PotAc), and glutamine synthetase (GS).
The fluorescent labeling is based on a conjugation reaction between NHS-Ph-BODIPY possess an active succinimidyl ester moiety and a protein possess aminecontaining residues (i.e., ε-amines of lysine side chains and α-amines at N-termini) to form a stable amide linkage (Scheme 2). The conjugation reaction conditions, namely slightly basic pH and low temperature, were chosen to ensure amine groups in proteins were unprotonated and to prevent hydrolysis of NHS-Ph-BODIPY to remove a reactive group [37]. Labeling details are given in Section 3.4. As can be seen from the visual representation of the matrices (Figure 5b,c), the first singlet excitation is localized mostly on the pyrrole rings and has a local character, although there is also a significant electron transfer between the rings, which, nevertheless, is symmetric. The change in the dipole moment arises due to the density transfer between the pyrrole rings as well as C 8 -carbon atom and its substituent. However, the position of the phenyl ring relative to dipyrrin is almost perpendicular due to the existence of the methyl groups in 1-and 7-positions of BODIPY core. As such, steric factors around this position should prevent any solvent molecules from getting close enough to provide a significant local impact on electron density via polarization.
Comparing two matrices, it becomes clear that the nature of the phenyl substituent does not play a strong role when it comes to the first excited state-transition density on the succinimidyl group, as well as density change between it and other units, are small. As such, as will be shown below (Sections 2.4 and 2.5), any alterations in the absorption and emission bands positions of COOH-Ph-BODIPY and NHS-Ph-BODIPY upon conjugation with proteins will come as a result of intermolecular interactions, such as polarization due to nearby amino acids or loss of polarization induced by solvent.

Fluorescent Labeling of Proteins
The synthesized amine-reactive fluorescent dye NHS-Ph-BODIPY was used for labeling of a range of proteins, such as bovine serum albumin (BSA), PII-like protein PotN (PotN), PotA subunit of polyamine ABC transporter (PotAc), and glutamine synthetase (GS).
The fluorescent labeling is based on a conjugation reaction between NHS-Ph-BODIPY possess an active succinimidyl ester moiety and a protein possess amine-containing residues (i.e., ε-amines of lysine side chains and α-amines at N-termini) to form a stable amide linkage (Scheme 2). The conjugation reaction conditions, namely slightly basic pH and low temperature, were chosen to ensure amine groups in proteins were unprotonated and to prevent hydrolysis of NHS-Ph-BODIPY to remove a reactive group [37]. Labeling details are given in Section 3.4.
At the first stage, the ability of NHS-Ph-BODIPY to label biomolecules was tested using a globular protein BSA as a model compound. At the first stage, the ability of NHS-Ph-BODIPY to label biomolecules was tested using a globular protein BSA as a model compound.
Absorption spectra, emission spectra, fluorescence decay curves, as well as a number of photophysical characteristics of NHS-Ph-BODIPY and its conjugate with BSA (BSA-NHS-Ph-BODIPY) were obtained in a mixture of DMSO and bicarbonate buffer with pH 8.3 (1:9) (Figures 6 and S4, Table 4).  Absorption spectra, emission spectra, fluorescence decay curves, as well as a number of photophysical characteristics of NHS-Ph-BODIPY and its conjugate with BSA (BSA-NHS-Ph-BODIPY) were obtained in a mixture of DMSO and bicarbonate buffer with pH 8.3 (1:9) (Figures 6 and S4, Table 4). At the first stage, the ability of NHS-Ph-BODIPY to label biomolecules was tested using a globular protein BSA as a model compound.
Absorption spectra, emission spectra, fluorescence decay curves, as well as a number of photophysical characteristics of NHS-Ph-BODIPY and its conjugate with BSA (BSA-NHS-Ph-BODIPY) were obtained in a mixture of DMSO and bicarbonate buffer with pH 8.3 (1:9) (Figures 6 and S4, Table 4).    The absorption and emission spectra of BSA-NHS-Ph-BODIPY have several distinctive features ( Figure 6). The absorption spectrum gave three bands with the maxima at 280, 380, and 530 nm, while the emission spectrum gave two bands with the maxima at 350 (not shown) and 554 nm. The new bands in the UV region are certain to relate with the absorption and emission of the protein moiety of the conjugate. The characteristic absorption and emission bands maxima of the BODIPY moiety of the conjugate are red shifted compared with the spectra of the free dye. Along with this, it is observed a decrease (about 1.5 times) of the fluorescence quantum yield of BSA-NHS-Ph-BODIPY in contrast to that of NHS-Ph-BODIPY as well as a change of the intensity decay from single exponential in case of NHS-Ph-BODIPY to double exponential in case of BSA-NHS-Ph-BODIPY. These phenomena can be explained by a dynamic quenching by aromatic amino acids located at the immediate environment of the dye in the conjugate [38,39], which is confirmed via synchronous fluorescence spectroscopy ( Figure S5) and molecular docking (Section 2.5). As regards the synchronous fluorescence spectroscopy, it was found that the change from a pure BSA to a labeled BSA results in a decrease of fluorescence intensities (from 3 times in case of ∆λ = 15 nm to 10 times in case of ∆λ = 60 nm) and a blue shift (from 2 nm in case of ∆λ = 15 nm to 10 nm in case of ∆λ = 60 nm) of the bands maxima at both wavelength shifts. This suggests that the addition of NHS-Ph-BODIPY causes changes in the environment of tyrosine and tryptophan residues of BSA, namely an increase of its hydrophobicity [40].
Furthermore, on the basis of the absorption spectra of NHS-Ph-BODIPY and BSA-NHS-Ph-BODIPY, the degree of labeling of BSA was estimated (Section 3.7). The value amounts to 1 which means that each protein molecule possesses one fluorescent label.
At the next stage, the possibility of NHS-Ph-BODIPY utilization in protein analysis was examined. To

Molecular Docking of Protein-BODIPY Conjugates
To gain insight into protein-BODIPY interactions, a two-stage molecular docking of BSA and PotN with NHS-Ph-BODIPY was carried out.
At the first stage, a blind docking was performed to determine the most probable binding site of NHS-Ph-BODIPY in BSA. It was shown that the fluorescent dye is localized in a cavity between the IB, IIA, and IIB subdomains of the protein (Figure 9a) with a binding free energy of −26.1 kJ/mol. It is noteworthy that these subdomains are known to contain major drug binding sites [41][42][43]. The amino acid composition of the binding site of NHS-Ph-BODIPY in BSA is shown in Figure 9b. The binding site was found to contain three lysine residues (Lys187, Lys221, Lys294), each of which was selected for a subsequent covalent docking. In addition, it contains tyrosine (Tyr149, Tyr156, Tyr451)

Molecular Docking of Protein-BODIPY Conjugates
To gain insight into protein-BODIPY interactions, a two-stage molecular docking of BSA and PotN with NHS-Ph-BODIPY was carried out.
At the first stage, a blind docking was performed to determine the most probable binding site of NHS-Ph-BODIPY in BSA. It was shown that the fluorescent dye is localized in a cavity between the IB, IIA, and IIB subdomains of the protein (Figure 9a) with a binding free energy of −26.1 kJ/mol. It is noteworthy that these subdomains are known to contain major drug binding sites [41][42][43]. The amino acid composition of the binding site of NHS-Ph-BODIPY in BSA is shown in Figure 9b. The binding site was found to contain three lysine residues (Lys187, Lys221, Lys294), each of which was selected for a subsequent covalent docking. In addition, it contains tyrosine (Tyr149, Tyr156, Tyr451) and tryptophan (Trp213) residues, which is consistent with the results of synchronous fluorescence spectroscopy (Section 2.4). and tryptophan (Trp213) residues, which is consistent with the results of synchronous fluorescence spectroscopy (Section 2.4). At the second stage, a covalent docking was performed to estimate a binding free energy of NHS-Ph-BODIPY and BSA bound together via an amide linkage through one of the lysine residues (Lys187, Lys221, Lys294) found in the binding site at the previous stage (Figures 10 and S6). The binding free energies values of Lys187-BSA-NHS-Ph- At the second stage, a covalent docking was performed to estimate a binding free energy of NHS-Ph-BODIPY and BSA bound together via an amide linkage through one of the lysine residues (Lys187, Lys221, Lys294) found in the binding site at the previous stage (Figures 10 and S6). The binding free energies values of Lys187-BSA-NHS-Ph-BODIPY, Lys221-BSA-NHS-Ph-BODIPY, and Lys294-BSA-NHS-Ph-BODIPY systems being −127, −129, and −108 kJ/mol, respectively, indicate an exergonic type of the conjugation reaction. The first two systems are characterized by almost equal binding free energies, which does not allow one to get an unambiguous answer about an anchor residue involved in the conjugation. Therefore, it was estimated a distance between the tyrosine and tryptophan residues (Tyr149, Tyr156, Tyr451, Trp213) and the lysine residues bound to NHS-Ph-BODIPY. The distance values are 5.669 (average) and 12.201 Å for Tyr Lys187-BSA-NHS-Ph-BODIPY and Trp···Lys187-BSA-NHS-Ph-BODIPY, respectively, and 3.524 (average) and 3.357 Å for Tyr···Lys221-BSA-NHS-Ph-BODIPY and Trp···Lys221-BSA-NHS-Ph-BODIPY, respectively. Thus, a synergy of synchronous fluorescence spectroscopy and molecular docking made it possible to establish that it is Lys221 that participate in the conjugation reaction between NHS-Ph-BODIPY and BSA. A similar two-stage approach was used for molecular docking of PotN with NHS-Ph-BODIPY. A blind docking results shown that the fluorescent dye is localized near the B-loop of the protein ( Figure 11) with a binding free energy of −31.9 kJ/mol. The binding site was found to contain one lysine residue (Lys92). A covalent docking results ( Figure S7) gave a binding free energy value of −122 kJ/mol for the Lys92-PotN-NHS-Ph-BODIPY system. A similar two-stage approach was used for molecular docking of PotN with NHS-Ph-BODIPY. A blind docking results shown that the fluorescent dye is localized near the B-loop of the protein ( Figure 11) with a binding free energy of −31.9 kJ/mol. The binding site was found to contain one lysine residue (Lys92). A covalent docking results ( Figure S7) gave a binding free energy value of −122 kJ/mol for the Lys92-PotN-NHS-Ph-BODIPY system.
It is noteworthy that the molecular docking results give explanation for a difference in sensitivity of BSA and PotN detection with NHS-Ph-BODIPY (Section 2.4). As stated above, the BSA binding site contains three lysine residues, while the PotN binding site contains only one lysine residue. Thus, lower detection limit for BSA than for PotN may be due to a higher probability of NHS-Ph-BODIPY to bind lysine residues at the binding site of the first protein than the second one.
A similar two-stage approach was used for molecular docking of PotN with NHS-Ph-BODIPY. A blind docking results shown that the fluorescent dye is localized near the B-loop of the protein ( Figure 11) with a binding free energy of −31.9 kJ/mol. The binding site was found to contain one lysine residue (Lys92). A covalent docking results ( Figure S7) gave a binding free energy value of −122 kJ/mol for the Lys92-PotN-NHS-Ph-BODIPY system. It is noteworthy that the molecular docking results give explanation for a difference in sensitivity of BSA and PotN detection with NHS-Ph-BODIPY (Section 2.4). As stated above, the BSA binding site contains three lysine residues, while the PotN binding site contains only one lysine residue. Thus, lower detection limit for BSA than for PotN may be due to a higher probability of NHS-Ph-BODIPY to bind lysine residues at the binding site of the first protein than the second one.
Dialysis tubing MD25 with 8.0-14.0 kDa molecular weight cut-off was obtained from a commercial source and prepared for use according to the procedure [46].

Instruments
1 H and 11 B nuclear magnetic resonance (NMR) spectra were recorded on a Bruker Avance III 500 NMR spectrometer (Billerica, MA, USA) with operating frequencies of 500.17 for 1 H and 160.48 MHz for 11 B. Deuterated chloroform was chosen as a solvent for samples preparation. Tetramethylsilane was used as an internal reference (δ = 0.00 ppm) for 1 H NMR studies, while BF 3 ·OEt 2 was used as an external reference (δ = 0.00 ppm) for 11 B NMR studies. The following abbreviations are used to designate peak multiplicities and descriptors: s for singlet, d for doublet, t for triplet, and q for quadruplet. Infrared (IR) spectra were obtained by means of a Shimadzu IRAffinity-1 Fourier transform IR spectrophotometer (Billerica, MA, USA) equipped with a Specac Quest ATR Diamond GS10800-B accessory in the mid-infrared (400-4000 cm −1 ) region. The following abbreviations are used to designate signal intensities: w for weak, m for medium, s for strong, br for broad, and sh for sharp. Matrix-assisted laser desorption/ionization (MALDI) time of flight (TOF) mass spectra (MS) were recorded on a Shimadzu AXIMA Confidence MALDI TOF-TOF mass spectrometer in positive ion reflectron mode.
Ultraviolet-visible (UV-Vis) spectroscopy experiments were performed on an Aquilon SF-104 spectrophotometer (Waltham, MA, USA). Absorption spectra were measured in the range of 190-800 nm. Fluorescence spectroscopy experiments were performed on an Agilent Cary Eclipse fluorescence spectrophotometer (Santa Clara, CA, USA). Steady-state emission spectra were measured in the ranges of 270-400 and 490-800 nm with excitation wavelengths of 260 and 480 nm, respectively. Synchronous emission spectra were measured in the range of 250-350 nm with wavelength shifts ∆λ of 15 and 60 nm. Excitation and emission slit widths values were 5.0 nm. Time-resolved fluorescence spectroscopy experiments were carried out by means of a PicoQuant FluoTime 300 high performance fluorescence lifetime and steady state spectrometer with a PicoQuant PLS 450 sub-nanosecond pulsed light-emitting diode as an excitation source. An instrument response function of a system was measured with a stray light signal of a dilute colloidal silica suspension (LUDOX ® ). Fluorescence decay curves were measured at the maximum of the emission peaks and fluorescence lifetimes were obtained by reconvolution of the decay curves using a PicoQuant EasyTau 2 software package. Standard quartz cuvettes with 10 mm light path were used for all experiments.
Labeled proteins were separated with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) and basic native-polyacrylamide gel electrophoresis (basic native-PAGE) [47], visualized with Bio-Rad ChemiDoc XRS+ System in trans-UV mode and then stained with Coomassie brilliant blue.

Fluorescent Labeling of Proteins
Fluorescent labeling of proteins was carried out on the basis of the standard aminereactive probe labeling protocol [49] (Scheme 2). A solution of NHS-Ph-BODIPY (8 equiv.) in dimethyl sulfoxide (DMSO) was slowly added to a solution of protein (1 equiv.) in bicarbonate or phosphate buffer with pH 8.3. After a thorough stirring, the mixture was incubated at 4 • C for 24 h. After completion of the reaction, the mixture was dialyzed against 70% ethanol (three times) and then distilled water (three times) to give a pure labeled protein solution.

Determination of Photophysical Characteristics
The fluorescence quantum yields Φ of the compounds investigated were estimated by comparison with Rhodamine 6G as a standard of known fluorescence quantum yield (Φ = 0.91 in ethanol [50]) using Equation (2) [50]: where S is the integrated area under the emission spectrum, A is the absorbance at the excitation wavelength, and n is the refractive index of the solvent. The subscripts x and st refer to the unknown and reference solutions, respectively. The radiative k r and non-radiative k nr rate constants were calculated using Equations (3) and (4) [26], respectively: where τ is the fluorescence lifetime.

Description of Solvent Effects via Multiparameter Approach
In order to describe solvent effects on spectral characteristics of the compounds investigated, a multiparameter approach using the Catalán equation (Equation (5)) [28] was applied: y = y 0 + a SA SA + b SB SB + c SP SP + d SdP SdP where y is the value of a solvent-dependent physicochemical property in a given solvent, y 0 is the statistical quantity corresponding to the value of a solvent-dependent physicochemical property in the gas phase or in an inert solvent, SA is the empirical parameter of solvent hydrogen-bond donor acidity, SB is the empirical parameter of solvent hydrogenbond acceptor basicity, SP is the empirical parameter of solvent polarizability, SdP is the empirical parameter of solvent dipolarity, and a − d are the regression coefficients describing the sensitivity of a physicochemical property y to the different solute/solvent interaction mechanisms.

Determination of Degree of Labeling
The degrees of labeling of the proteins investigated were estimated using Equation (6) [49]: where DOL is the degree of labeling, A max conj and A max dye are the absorbances of the conjugated dye and the free dye at the maximum absorption wavelength, respectively, ε max prot and ε max dye are the molar absorption coefficients of the free protein and the free dye at the absorption wavelength, respectively, and A 280 conj and A 280 dye are the absorbances of the labeled protein and maximum the free dye at 280 nm, respectively.

Quantum Chemical Calculations Procedure
The conformational search for COOH-Ph-BODIPY and NHS-Ph-BODIPY was performed in a metadynamic approximation. The calculations were carried out by means of the xtb program [51] with the CREST add-on [52]. The more accurate GFN2-xTB method was chosen along with the iMTD-GC algorithm. The conformational screening revealed four conformers for COOH-Ph-BODIPY and four conformers for NHS-Ph-BODIPY. The lowest energy conformers were chosen for further DFT and TDDFT study. Each conformer was optimized at the CAM-B3LYP/6-31G(d,p) level, which provides a balanced description of both ground and excited states [53]. For each resulting structure, ten excitation vertical absorption spectra were computed. For the first three excitations, gradients as well as vertical gradient Franck-Condon (VG-FC) vibronic spectra within the TD formalism as implemented in the FCclasses 3.0 program [54] were also obtained. Vibronic spectra calculations used ground state vibrational frequencies, vertical excitation energies computed at ground state geometries and excited state gradients computed at the same geometries. All empirical parameters pertaining to the VG-FC algorithm, such as number of points for the fast Fourier transform (FFT), were directly estimated by the FCclasses program. All vibronic spectra are broadened with Gaussian functions with the constant width of 322 cm −1 . This value was chosen because it gives the best agreement between experimental and theoretical spectra of fully unsubstituted BODIPY. All calculations, unless otherwise mentioned, were performed within the CPCM model with DMSO parameters using the ORCA 5.0 program suite [55][56][57]. Quantities of interest, such as ground and excited state dipole moments as well as transition density matrix elements, were computed using the latest version of the Multiwfn program [58]. A visualization of the quantum chemical calculations results were made by means of both the Chemcraft 1.8 [59] and the VMD [60] programs.

Molecular Docking Procedure
The molecular docking of the proteins-BODIPY conjugates was performed in two stages. At the first stage, a blind docking was carried out by means of the Autodock 4.2 program [61]. The crystal structures of the proteins under study were taken from the Protein Data Bank: BSA [62] and PotN [63]. The structure of NHS-Ph-BODIPY was obtained by virtue of geometry optimization (Section 3.8). The calculations were performed for a 126 Å × 126 Å × 126 Å grid with BSA in its center with a step of 0.7 Å and a 90 Å × 126 Å × 126 Å grid with PotN in its center with a step of 0.4 Å. Each docking experiment included 50 runs with a maximum of 25 million energy evaluations using the Lamarckian genetic algorithm [64]. Conformations of the protein-BODIPY conjugates with minimum energies were assumed as the most stable ones. At the second stage, a covalent docking was carried out by means of the CovDock program [65]. On the basis of the blind docking results, a covalent binding of NHS-Ph-BODIPY and lysine residues (imine condensation pre-defined reaction) at the most energetically favorable binding site of the proteins was calculated. Conformations of the protein-BODIPY conjugates with minimum energies were again assumed as the most stable ones. A visualization of the molecular docking results were made by means of the UCSF Chimera program [66].

Conclusions
The development of fluorescent dyes reactive towards principal functional groups of protein is a promising area of modern chemistry. Therefore, we carried out a thorough investigation of spectroscopic and photophysical properties as well as solvatochromic behavior of the synthesized amine-reactive fluorescent dye NHS-Ph-BODIPY and its nonreactive precursor COOH-Ph-BODIPY by means of ultraviolet-visible and fluorescence spectroscopy supported by quantum chemical calculations that allowed us to better understand some specific aspects of their functioning in solutions. An inherent stable bright green fluorescence of NHS-Ph-BODIPY excited us to use this dye for fluorescent labeling of a number of proteins of various nature. NHS-Ph-BODIPY was demonstrated to act at least at the same level or even several times better than Coomassie brilliant blue when detecting proteins via SDS-PAGE. Thus, the sum of the experimental and theoretical results suggests that the synthesized amine-reactive fluorescent dye NHS-Ph-BODIPY is a promising fluorescent label for fluorescent immunoassay, cellular imaging and in vivo imaging, flow cytometry, immunohistochemical staining, etc.