Molecular Recognition of Imidazole-Based Drug Molecules by Cobalt(III)- and Zinc(II)-Coproporphyrins in Aqueous Media

The methods of 1H NMR, spectrophotometric titration, mass spectrometry and elemental analysis are applied to determine the selective binding ability of Co(III)- and Zn(II)-coproporphyrins I towards a series of imidazole-based drug molecules with a wide spectrum of pharmacological activity (metronidazole, histamine, histidine, tinidazole, mercazolil, and pilocarpine) in phosphate buffer (pH 7.4) simulating the blood plasma environment. It is shown that in aqueous buffer media, Co(III)-coproporphyrin I, unlike Zn(II)-coproporphyrin I, binds two imidazole derivatives, and the stability of mono-axial Co-coproporphyrin imidazole complexes is two to three orders of magnitude higher than that of similar complexes of Zn-coproporphyrin I. The studied porphyrinates are found to have the highest binding ability to histamine and histidine due to the formation of two additional hydrogen bonds between the carboxyl groups of the porphyrinate side chains and the binding sites of the ligands in the case of histidine and a hydrogen bond between the amino group of the ligand and the carbonyl oxygen atom of the carboxyl group of the porphyrinate in the case of histamine. The structures of the resulting complexes are optimized by DFT quantum chemical calculations. The results of these studies may be of use in the design of biosensors, including those for the detection, control and verification of various veterinary drug residues in human food samples.


Introduction
Complexes of porphyrins with coordinationally unsaturated metal cations (MP, M = Fe(II, III), Co(II, III), Mn(II, III, IV), etc.) are widespread in the living nature surrounding us. Coproporphyrins (tetrapyrrole macroheterocycles, which have propionic acid residues in four β-pyrrole positions) are intermediates in the biosynthesis of many vital substrates. In particular, the level of copro-and uroporphyrins in the urine can be used to identify heme biosynthesis disorders in patients, which can be either congenital, as, for example, in case of hereditary fermentopathies, or acquired, as, for example, with liver diseases and hemolytic anemia [1][2][3]. At the same time, coproporphyrinates have the properties of stable photosensitizers, are characterized by low toxicity and can be easily eliminated from the body. These compounds can be used to treat arthrosis and arthritis and in photodynamic therapy of cancer and surface-localized infections [4][5][6][7].
The aim of this work was to study the processes of the molecular recognition of imidazole-containing drug molecules of different natures by complexes of coproporphyrin I with coordinationally unsaturated metal cations. Due to the electronic effects of substituents and structural features, imidazole and its derivatives exhibit a wide range of pharmacological and/or biological activities and are effectively used in the development of new broad-spectrum drugs, including those used for veterinary medicine. In this work, we studied the binding ability of Co(III)-and Zn(II)-coproporphyrinates I (CoCP and ZnCP, respectively) to antimicrobial and antiprotozoal metronidazole and tinidazole, antithyroid mercazolil, m-cholinomimetic pilocarpine, proteinogenic amino acid histidine, work, we studied the binding ability of Co(III)-and Zn(II)-coproporphyrinates I (CoCP and ZnCP, respectively) to antimicrobial and antiprotozoal metronidazole and tinidazole, antithyroid mercazolil, m-cholinomimetic pilocarpine, proteinogenic amino acid histidine, and histamine, one of the mediators that regulate the vital functions of the human body ( Figure 1) [8]. The ability of metalloporphyrin molecules to selectively interact with substrates of different natures depends on the ability of the macrocycle and the substrate to form several intermolecular bonds with each other. For example, in addition to the donoracceptor bond between the coordinationally unsaturated metal cation of the macrocycle and the electron donor atom of the ligand, there may be additional binding sites formed due to van der Waals and/or dispersion interactions, hydrogen bonds, π-π stacking, etc. [9][10][11]. For example, as [12,13] show, π-π interactions between the porphyrin receptor and the aromatic part of an amino acid significantly stabilize the Zn-porphyrin complex with tryptophan, in comparison with a similar complex with glycine (the stability constant of the complex increases from 60 L/mol in glycine to 1200 L/mol in tryptophan). In our previous work [14], we studied the recognition ability of two hydrophilic Co-tetraarylporphyrins with hydroxyl groups in the phenyl rings of the macrocycle towards imidazole and its derivatives (metronidazole and histidine) to find out whether additional hydrogen bonds could be formed between the hydroxyl groups of the receptor (porphyrinate) and the oxygen atoms of the nitro and carboxyl groups of the substrates. In the coproporphyrinates, the groups that are responsible for their selectivity in organic substrate binding are conformationally flexible propionic acid residues in four β-pyrrole The ability of metalloporphyrin molecules to selectively interact with substrates of different natures depends on the ability of the macrocycle and the substrate to form several intermolecular bonds with each other. For example, in addition to the donor-acceptor bond between the coordinationally unsaturated metal cation of the macrocycle and the electron donor atom of the ligand, there may be additional binding sites formed due to van der Waals and/or dispersion interactions, hydrogen bonds, π-π stacking, etc. [9][10][11]. For example, as [12,13] show, π-π interactions between the porphyrin receptor and the aromatic part of an amino acid significantly stabilize the Zn-porphyrin complex with tryptophan, in comparison with a similar complex with glycine (the stability constant of the complex increases from 60 L/mol in glycine to 1200 L/mol in tryptophan). In our previous work [14], we studied the recognition ability of two hydrophilic Co-tetraarylporphyrins with hydroxyl groups in the phenyl rings of the macrocycle towards imidazole and its derivatives (metronidazole and histidine) to find out whether additional hydrogen bonds could be formed between the hydroxyl groups of the receptor (porphyrinate) and the oxygen atoms of the nitro and carboxyl groups of the substrates. In the coproporphyrinates, the groups that are responsible for their selectivity in organic substrate binding are conformationally flexible propionic acid residues in four β-pyrrole positions of the macrocycle. The structures of the resulting complexes were optimized by DFT quantum chemical calculations. The calculated and experimental data were in good agreement with each other, confirming the ability of the coproporphyrinates to molecularly recognize imidazole-based drug molecules.
In this work, zinc and cobalt were taken as metal ions coordinating the porphyrin ligand, since in buffer media both of these metals have the ability to additionally coordinate nitrogen-containing electron-donor molecules. Moreover, in the case of cobalt, it is possible to attach two ligand molecules with very high stability constants. The Zn(II)-porphyrin was used for comparison with the Co(III)-porphyrin to further demonstrate the differences in the behavior of these metals in aqueous and non-aqueous media.

Results and Discussion
The binding ability of ZnCP and CoCP(H 2 O) 2 to imidazole (Im) and a number of imidazole-containing drug molecules (L1-L6) in phosphate buffer (pH 7.4) at 20 • C was studied by the methods of 1 H NMR, spectrophotometric titration, mass spectrometry and elemental analyses.

Formation of Bis-Axial CoCP(H 2 O) 2 Complexes with Imidazole Derivatives
Co(III) complexes of water-soluble tetraarylporphyrins (Co(III)P or CoP) are characterized by bis-axial coordination of both nitrogen-and oxygen-containing small organic molecules [15,16].
The Co 2p spectrum recorded from CoCP(Im) shows the main Co 2p3/2 photoemission line at 788.9 eV and its satellite at 793.8 eV. According to [17], this feature allows us to identify the Co ion in the compound as an intermediate spin (S = 1) Co(III) species. The X-ray photoelectron spectra of the CoCP(L) samples were studied using an X-ray photoelectron spectrometer Kratos Analytical AXIS Supra, with a high-power monochromatic Al Kα radiation source with a photon energy of 1486.6 eV ( Figure S1).
Earlier, we studied the processes of complexation of Co(III)-tetraarylporphyrins (CoTCPP and CoTSPP) with various drugs based on nitrogen-containing heterocycles [18,19]. Table 1 and Figure 2a show data on the stability constants of CoTCPP axial complexes, taken in this work as the object of comparison, with a number of imidazole derivatives. The location of the carboxyl groups in this macrocycle (4-positions of the phenyl rings) makes it impossible for additional hydrogen bonds to be formed with these ligands. The stability constants of the mono-and diaxial CoTCPP complexes with different L differ insignificantly.
Like CoTCPP and CoTSPP, CoCP forms bis-aqua complexes in aqueous media at pH 7.4, and the replacement of the aqua-groups with organic bases represents a twostage process: Table 1. Total energies of formation of ZnCP(L), CoCP(L) and CoCP(L) 2 (−∆G g , kJ·mol −1 ), M-N(Im) bond energies (−∆G M-L , kJ·mol −1 ) and hydrogen bond energies (−∆G HB , kJ·mol −1 ).  plexes, taken in this work as the object of comparison, with a number of imidazole derivatives. The location of the carboxyl groups in this macrocycle (4-positions of the phenyl rings) makes it impossible for additional hydrogen bonds to be formed with these ligands. The stability constants of the mono-and diaxial CoTCPP complexes with different L differ insignificantly. Like CoTCPP and CoTSPP, CoCP forms bis-aqua complexes in aqueous media at pH 7.4, and the replacement of the aqua-groups with organic bases represents a two-stage process: Each of the process stages is accompanied by its own characteristic spectral changes in the UV-Vis spectra of the reaction system. In particular, the replacement of the first water molecule with the ligand is characterized by a decrease in the intensity of the Soret band and its red shift by 3.5 nm (Figure 3a). The replacement of the second water molecule is characterized by less noticeable spectral changes ( Figure 3b).
Theoretical UV-Vis spectra were calculated using the ZINDO/S semi-empirical method. This method is adapted for transition element compounds and also helps in predicting and describing electronic transitions of metalloporphyrin molecules and related compounds ( Figure S2).
It is known that the electronic transitions of porphyrins corresponding to the Soret band (B-band) in the visible region are carried out mainly due to the π→π* transitions between HOMO-1 and LUMO/LUMO + 1 orbitals. Each of the process stages is accompanied by its own characteristic spectral changes in the UV-Vis spectra of the reaction system. In particular, the replacement of the first water molecule with the ligand is characterized by a decrease in the intensity of the Soret band and its red shift by 3.5 nm (Figure 3a). The replacement of the second water molecule is characterized by less noticeable spectral changes ( Figure 3b).
Theoretical UV-Vis spectra were calculated using the ZINDO/S semi-empirical method. This method is adapted for transition element compounds and also helps in predicting and describing electronic transitions of metalloporphyrin molecules and related compounds ( Figure S2).
It is known that the electronic transitions of porphyrins corresponding to the Soret band (B-band) in the visible region are carried out mainly due to the π→π* transitions between HOMO-1 and LUMO/LUMO + 1 orbitals.
x FOR PEER REVIEW 5 of 18 The calculated constants of the water molecules' replacement for the imidazole compounds (K 1 and K 2 ) are listed in Table 1 and in Figure 2c,d. As the presented data show, the CoCP binding ability towards imidazole in aqueous media is lower than that of hydrophilic CoTCPP. Moreover, the stability of the mono-and diaxial Co(III)-coproporphyrin complexes with the imidazole derivatives (L1-L6) strongly depends on the presence of additional binding sites in the latter ( Figure 2).
The obtained dependence of the stability constants on the ligand nature and the calculated thermodynamic parameters of the axial coordination process, which will be discussed below, suggest that at least four of the six studied drug molecules, when coordinated on the macrocycle zinc cation, form additional hydrogen bonds with the binding sites of the conformationally mobile propionic acid branches of CoCP. The most stable complexes are formed with histidine and histamine, probably due to the formation The results of quantum chemical calculations are in good agreement with the experimental red shift of the Soret band during complexation. Figure S2 shows the calculated electronic spectra for the complexes, as well as the scheme of electronic transitions.
The calculated constants of the water molecules' replacement for the imidazole compounds (K 1 and K 2 ) are listed in Table 1 and in Figure 2c,d. As the presented data show, the CoCP binding ability towards imidazole in aqueous media is lower than that of hydrophilic CoTCPP. Moreover, the stability of the mono-and diaxial Co(III)-coproporphyrin complexes with the imidazole derivatives (L1-L6) strongly depends on the presence of additional binding sites in the latter ( Figure 2). The obtained dependence of the stability constants on the ligand nature and the calculated thermodynamic parameters of the axial coordination process, which will be discussed below, suggest that at least four of the six studied drug molecules, when coordinated on the macrocycle zinc cation, form additional hydrogen bonds with the binding sites of the conformationally mobile propionic acid branches of CoCP. The most stable complexes are formed with histidine and histamine, probably due to the formation of fairly strong hydrogen bonds between the porphyrinate binding sites and amine and/or carboxyl groups of the ligand (6.5 and 5.8 kJ/mol, respectively). The stability constant of the CoCP(L2) and CoCP(L3) complexes is higher than that of the CoCP(Im) complex by an order of magnitude.
In the Co(III)-5,15-bis-(4-sylfophenyl)-10,20-bis-(3-hydroxyphenyl)-porphyrin and Co(III)-5,15-bis-(4-sylfophenyl)-10,20-bis-(2-hydroxyphenyl)-porphyrin described by us in [14], the selectivity of binding imidazole derivatives is only observed when the first ligand is bound. When the second axial ligand is attached, the binding process is not selective. This is probably explained by the fact that in the case of these porphyrins, the second complexation step is only due to the single-center binding of the ligands through the metal cation of the macrocycle (i.e., the sulfo-and hydroxy-binding sites rigidly fixed on the macrocycle periphery are not capable of forming additional bonds with the ligand molecule).
At the same time, it was established that the CoCP's ability to molecularly recognize imidazoles is also retained upon coordination of the second ligand, although the energy of additional binding for the second ligand is somewhat lower. This, apparently, means that the coordination of the second imidazole derivative is also accompanied by the formation of additional binding centers.

ZnCP Mono-Axial Complexes with Imidazole Derivatives
In both non-aqueous and aqueous media, the zinc cation in the composition of porphyrin macrocycles is in the oxidation state (II) and forms mono-axial ZnCP(L) complexes with strong organic bases, such as imidazole and its derivatives [20][21][22]. The data presented in Figure 4 indicate that complexation of Zn-porphyrin with the organic ligand results in the formation of 1:1 complexes (there is one family of spectral curves in the UV-Vis spectrum of the reaction system, Figure 4a, and there is one step on the respective curves of spectrophotometric titration, Figure 4b). The shift in the Soret band that accompanies the formation of the ZnCP mono-ligand complex is~7 nm. , x FOR PEER REVIEW 6 of 18 In the Co(III)-5,15-bis-(4-sylfophenyl)-10,20-bis-(3-hydroxyphenyl)-porphyrin and Co(III)-5,15-bis-(4-sylfophenyl)-10,20-bis-(2-hydroxyphenyl)-porphyrin described by us in [14], the selectivity of binding imidazole derivatives is only observed when the first ligand is bound. When the second axial ligand is attached, the binding process is not selective. This is probably explained by the fact that in the case of these porphyrins, the second complexation step is only due to the single-center binding of the ligands through the metal cation of the macrocycle (i.e., the sulfo-and hydroxy-binding sites rigidly fixed on the macrocycle periphery are not capable of forming additional bonds with the ligand molecule).
At the same time, it was established that the CoCP's ability to molecularly recognize imidazoles is also retained upon coordination of the second ligand, although the energy of additional binding for the second ligand is somewhat lower. This, apparently, means that the coordination of the second imidazole derivative is also accompanied by the formation of additional binding centers.

ZnCP Mono-Axial Complexes with Imidazole Derivatives
In both non-aqueous and aqueous media, the zinc cation in the composition of porphyrin macrocycles is in the oxidation state (II) and forms mono-axial ZnCP(L) complexes with strong organic bases, such as imidazole and its derivatives [20][21][22]. The data presented in Figure 4 indicate that complexation of Zn-porphyrin with the organic ligand results in the formation of 1:1 complexes (there is one family of spectral curves in the UV-Vis spectrum of the reaction system, Figure 4a, and there is one step on the respective curves of spectrophotometric titration, Figure 4b). The shift in the Soret band that accompanies the formation of the ZnCP mono-ligand complex is ~7 nm.  The stability constants of the mono-axial ZnCP complexes with various imidazole-containing drugs and the enthalpies of their formation are shown in the diagram ( Figure 5) and in Table 1.  The stability constants of the mono-axial ZnCP complexes with various imidazolecontaining drugs and the enthalpies of their formation are shown in the diagram ( Figure 5) and in Table 1.   As the presented data show, ZnCP's binding ability to imidazole in aqueous medi is three to four times lower than that of hydrophilic Zn-tetraarylporphyrins. Moreove the stability of the mono-axial ZnCP complexes with the imidazole derivatives (L1-L6 also strongly depends on the presence of additional binding sites in the latter ( Figure 5 The ZnCP(L) stability constants for the imidazole derivatives have a similar relationshi to those of CoCP(L) and CoCP(L)2. This means that upon coordination of the ligands o ZnCP, the studied substrates form similar additional bonds with the propionic aci fragments of the porphyrinate. The most stable complexes are formed with histidine an histamine, probably due to the formation of fairly strong hydrogen bonds between th porphyrinate binding sites and amine and/or carboxyl groups of the ligand (6.5 and 5. kJ·mol −1 , respectively). The stability constant of the ZnCP(L2) and ZnCP(L3) complexes i higher than that of the ZnCP(Im) complex by more than an order of magnitude. Table 2 shows the calculated energy and geometric parameters of ZnCP complexe with axial ligands (L1-L6). An NBO analysis was performed to predict the natural mo lecular charge distribution as well as to calculate the stabilization energies between th As the presented data show, ZnCP's binding ability to imidazole in aqueous media is three to four times lower than that of hydrophilic Zn-tetraarylporphyrins. Moreover, the stability of the mono-axial ZnCP complexes with the imidazole derivatives (L1-L6) also strongly depends on the presence of additional binding sites in the latter ( Figure 5). The ZnCP(L) stability constants for the imidazole derivatives have a similar relationship to those of CoCP(L) and CoCP(L) 2 . This means that upon coordination of the ligands on ZnCP, the studied substrates form similar additional bonds with the propionic acid fragments of the porphyrinate. The most stable complexes are formed with histidine and histamine, probably due to the formation of fairly strong hydrogen bonds between the porphyrinate binding sites and amine and/or carboxyl groups of the ligand (6.5 and 5.8 kJ·mol −1 , respectively). The stability constant of the ZnCP(L2) and ZnCP(L3) complexes is higher than that of the ZnCP(Im) complex by more than an order of magnitude. Table 2 shows the calculated energy and geometric parameters of ZnCP complexes with axial ligands (L1-L6). An NBO analysis was performed to predict the natural molecular charge distribution as well as to calculate the stabilization energies between the ligand and the porphyrinate zinc cation. The interaction between the donor and acceptor orbitals was analyzed based on the second-order perturbation stabilization energy (E st ) and the charge transfer value (q CT ) presented in Table 2. The intermolecular interaction energy of the studied complexes with the axial ligands was also calculated by taking into account the BSSE superposition error (−E int ). The NBO analysis shows that the Zn-N bond in the coordination center was identical for all of the calculated structures and was formed by the partial transfer of the nitrogen atom electron pair to the free orbitals of the zinc cation (LP N → LP * Zn). The interaction between the donor and acceptor orbitals was analyzed based on the second-order perturbation stabilization energy. The E st value increased as the interaction energy between the host and guest molecules became higher. The stabilization energy between the central zinc cation and the nitrogen atom of the macrocycle pyrrole ring, according to the NBO analysis data, was 43-49 kcal·mol −1 , which is higher than the stabilization energy of Zn-porphyrinate with axial ligands.

Geometry Optimization of ZnCP(L) and CoCP(L)
The calculation results show that the interaction of the macrocycle zinc cation with the nitrogen atom of the amino acid led to the elongation of the Zn(II) ion from the macrocycle plane. This can be seen from the values of the N-Zn-N angle, which decreased proportionally to the increase in the energy of the axial ligand binding to the Zn-porphyrin ( Table 2). At the same time, the stabilization energy of the Zn-N bond decreased in inverse proportion to the energy of the porphyrinate interaction with amino acids, which was in the range of 11-20 kcal/mol. Thus, the highest binding energy of Zn-porphyrin was observed in the case of histamine (L2) and decreased in the series L2 > L3 > L6 > L4 > L5 > L1.
As the analysis of the optimized structures shows, the higher binding energy of L2 and L3 is probably associated with the formation of additional hydrogen bonds between the porphyrinate and the ligand. The lengths of these bonds are 1.9 and 2.0 Å, respectively, which allows them to be classified as medium-strength hydrogen bonds.
Similar geometric, electronic and energy parameters calculated for the CoCP(L) monoligand complexes are presented in Table 3. The geometric optimization of the extracomplexes was carried out by taking into account the above-described features of the geometric structure of the initial molecules. During the optimization, many different configurations were taken into account, differing in the presence and type of hydrogen bonds between the periphery of the macrocycle and extra-ligand. Figure 6 shows the optimized structure of the complexes corresponding to configurations with the minimal energy.
The calculated Zn-N and Co-N bond lengths and the N-Zn-N and N-Co-N angles given in Tables 2 and 3 are in good agreement with the literature data [23][24][25][26].
It is obvious that a significant contribution to the energy of the intermolecular interaction between the extra-ligand and porphyrinate is provided by hydrogen bond formation, which is observed in all of the considered extra-complexes.
According to Table 4, the strongest hydrogen bonds are observed in the CoCP(L2) and CoCP(L3) extra-complexes. In the first case, the significant energy gain is due to the formation of a cyclic hydrogen bond between the carboxyl groups, as it occurs in carboxylic acid dimers, and in the second, an abnormally strong hydrogen bond is formed between the amino group of the extra-ligand and the carboxyl group on the macrocycle periphery. The high stabilization energy can also indicate proton transfer processes in the CoCP(L3) complex, similar to the proton transfer in amino acids (Figure 7).      It is obvious that a significant contribution to the energy of the intermolecular in teraction between the extra-ligand and porphyrinate is provided by hydrogen bond formation, which is observed in all of the considered extra-complexes.
According to Table 4, the strongest hydrogen bonds are observed in the CoCP(L2 and CoCP(L3) extra-complexes. In the first case, the significant energy gain is due to the formation of a cyclic hydrogen bond between the carboxyl groups, as it occurs in car boxylic acid dimers, and in the second, an abnormally strong hydrogen bond is formed between the amino group of the extra-ligand and the carboxyl group on the macrocycle periphery. The high stabilization energy can also indicate proton transfer processes in the CoCP(L3) complex, similar to the proton transfer in amino acids (Figure 7).          In conclusion, it should be said that the calculated values of the intramolecular interaction energy for the extra-complexes in free state are consistent with the thermodynamic stability constant of the complexes calculated from the experimental data ( Figures  8 and 9). However, it should be taken into account that the real systems have been studied in buffer solutions, the components of which, together with the solvent, can also participate in the interactions. In conclusion, it should be said that the calculated values of the intramolecular interaction energy for the extra-complexes in free state are consistent with the thermodynamic stability constant of the complexes calculated from the experimental data (Figures 8 and 9). However, it should be taken into account that the real systems have been studied in buffer solutions, the components of which, together with the solvent, can also participate in the interactions.
In conclusion, it should be said that the calculated values of the intramolecular interaction energy for the extra-complexes in free state are consistent with the thermodynamic stability constant of the complexes calculated from the experimental data ( Figures  8 and 9). However, it should be taken into account that the real systems have been studied in buffer solutions, the components of which, together with the solvent, can also participate in the interactions.

1 H NMR Studies
The formation of the MCP(L) metal complexes was readily observed by 1 H NMR spectroscopy. The characteristics of the 1 H NMR spectra of a series of CoCP(L), which can be formed without an excess of a substrate, are presented in Table 5 and Figure 10.

1 H NMR Studies
The formation of the MCP(L) metal complexes was readily observed by 1 H NMR spectroscopy. The characteristics of the 1 H NMR spectra of a series of CoCP(L), which can be formed without an excess of a substrate, are presented in Table 5 and Figure 10.
The axial coordination of the L1-L6 ligands on the Co(III) cation of the macrocycle is evidenced by the characteristic high-field shifts of the signals of the ligand imidazole fragment, which are in the immediate vicinity of the ring current of the aromatic tetrapyrrole macrocycle.
The additional hydrogen bonds in the MCP(L) formed between the O-H, N-H, and S-H groups of the macrocycle and ligand are predominantly of an electrostatic nature.
Their formation is always accompanied by a strong downfield shift in the resonant signal involved in proton binding [27][28][29]. D6-DMSO is a solvent capable of forming very strong hydrogen bonds with the above protons ( Table 5). The signals of such protons in d6-DMSO are narrower and shifted downfield compared to the other solvents.
During the formation of axial CoCP(L) complexes with additional hydrogen bonds (Table 4), which are stronger than H···d6-DMSO, the signals are also shifted downfield compared to the signal positions in the L1-L6 spectra. In the 1 H NMR spectra of CoCP(L1-L4, L6), the proton signal of the porphyrinate carboxyl groups (broadened singlet) is split and partially shifted downfield ( Figure 10). Similar downfield shifts are observed in the proton signals of the CoCP(L1) hydroxy group, the CoCP(L3) amino group, and the CoP(L5) HS group ( Table 5). The fact that the -COOH resonances of the 1:1 CoCP and L2 mixture appear as four sharp signals, which are significantly shifted downfield from their original position in the individual components (CoCP and L2), confirms that two OH···O hydrogen bonds are formed in the CoP(L5) complex. 1.07 (s, 3H, -CH2-CH3) 1.12(s, 3H, -CH2-CH3) The axial coordination of the L1-L6 ligands on the Co(III) cation of the macrocycle is evidenced by the characteristic high-field shifts of the signals of the ligand imidazole fragment, which are in the immediate vicinity of the ring current of the aromatic tetrapyrrole macrocycle. Their formation is always accompanied by a strong downfield shift in the resonant signal involved in proton binding [27][28][29]. D6-DMSO is a solvent capable of forming very strong hydrogen bonds with the above protons ( Table 5). The signals of such protons in d6-DMSO are narrower and shifted downfield compared to the other solvents.
During the formation of axial CoCP(L) complexes with additional hydrogen bonds (Table 4), which are stronger than H···d6-DMSO, the signals are also shifted downfield compared to the signal positions in the L1-L6 spectra. In the 1 H NMR spectra of CoCP(L1-L4, L6), the proton signal of the porphyrinate carboxyl groups (broadened singlet) is split and partially shifted downfield ( Figure 10). Similar downfield shifts are observed in the proton signals of the CoCP(L1) hydroxy group, the CoCP(L3) amino group, and the CoP(L5) HS group ( Table 5). The fact that the -COOH resonances of the 1:1 CоCP and L2 mixture appear as four sharp signals, which are significantly shifted The existence of the most durable mono-axial complexes was also recorded by mass spectrometry (Figure 11). In the mass spectra of the CoCP(L2) and CoCP(L3) complexes obtained in the negative shooting mode, the peaks of the target ions [M] − (m/z 835.7) and (m/z 821.8), respectively, were recorded. The highest intensity peak in the mass spectra of both complexes was formed as a result of the elimination of the axial ligands and corresponds to the CoCP ion (m/z 710.7).   Figure S3 shows the powder diffraction patterns of CoCP, one of the studied ligands (L2), and the CoCP(L2)2 complex. The diffraction pattern of the CoCP(L2)2 complex differs from the diffraction pattern of the starting compounds. This sample is less crystalline (more amorphous), which does not allow any correspondence with our simulated samples.   Figure S3 shows the powder diffraction patterns of CoCP, one of the studied ligands (L2), and the CoCP(L2)2 complex. The diffraction pattern of the CoCP(L2)2 complex differs from the diffraction pattern of the starting compounds. This sample is less crystalline (more amorphous), which does not allow any correspondence with our simulated samples.

Synthesis
Zn(II)-coproporphyrin I and Co(III)-coproporphyrin I were obtained from commercially available coproporphyrin I by means of the method described in [30]. The product of the interaction of coproporphyrin I with cobalt chloride in ethanol in a nitrogen atmosphere is Co(II)-coproporphyrin. When planted in water in an air atmosphere, Co(II)coproporphyrin is oxidized to diaqua-Co(III)-coproporphyrin.

Characterization
Spectrophotometric studies. The UV-Vis spectra were recorded at 25 • C on a Jasco V-770 spectrometer. The thermodynamic constants for the complexation of CoCP and ZnCP with a ligand (L) were calculated by the formula based on the spectrophotometric titration experiment results: where λ 1 is the descending wavelength; λ 2 is the ascending wavelength; [MP-L] is the concentration of the porphyrinate with one axial ligand; and [L] is the ligand concentration.
∆A o is the maximal change in the solution optical density at the given wavelength, and ∆A i is the change in the solution optical density at the given wavelength at the given concentration [31]. NMR studies. All the NMR experiments were performed on a Bruker Avance III-500 NMR spectrometer equipped with a 5 mm probe using standard Bruker TOPSPIN Software. TMS signals were used as the internal standards to determine the chemical shifts. The temperature control was performed using a Bruker variable temperature unit (BVT-2000) in combination with a Bruker cooling unit (BCU-05) to provide chilled air. The experiments were performed at 298 K without sample spinning.
The elemental analyses were performed on a Flash EA 1112 C, H, N analyzer. The mass-spectra were obtained on a Shimadzu Biotech Axima Confidence MaldiTof mass-spectrometer (with dihydroxybenzoic acid as the matrix).
The X-ray photoelectron spectra were studied on an X-ray photoelectron spectrometer Kratos Analytical AXIS Supra, with a high-power monochromatic Al Kα radiation source with a photon energy of 1486.6 eV.

Details of Quantum Chemical Calculations
The density functional theory (DFT) calculations were carried out using the GAUS-SIAN 16 quantum chemical program package [32]. The geometry optimization, vibrational frequency calculations and conformational analysis of the complexes and individual molecules in the ground state were performed using the B97D3 functional [33,34] in combination with the 6-31 + G* [35] basis set. The analysis of the electron density distribution was carried out using the NBO program [36] implemented in GAUSSIAN 16. All of the calculations were performed for isolated molecules in vacuum, without taking into account the effect of the solvent. The results of the quantum chemical calculations were visualized with the ChemCraft program [37].

Conclusions
Thus, this study has shown that Co(III) complexes of coproporphyrin I have the ability to bind one or two molecules of imidazole-containing drugs with a wide range of pharmacological activities (metronidazole, histamine, histidine, tinidazole, mercazolil, and pilocarpine) in phosphate buffer (pH 7.4), simulating the blood plasma environment. The stability of the resulting mono-and diaxial porphyrin complexes with the ligands depends on whether the latter has a nitrogen-containing heterocycle (imidazole) and various binding sites capable of forming additional hydrogen bonds with conformationally flexible propionic acid side substituents of the macrocycle (there are two hydrogen bonds between the carboxyl groups in the case of histidine, and a hydrogen bond between the amino group of the ligand and the carbonyl oxygen atom of the carboxyl group of the porphyrinate in the case of histamine). The spectral responses to the substrate binding processes were identified and characterized, and the stability constants of the resulting complexes and concentration ranges of their existence were determined. The results of these studies may be of use in the design of biosensors, including those for the detection, control and verification of various veterinary drug residues in human food samples.