Flow-Through Macroporous Polymer Monoliths Containing Artiﬁcial Catalytic Centers Mimicking Chymotrypsin Active Site

: Synthetic catalysts that could compete with enzymes in term of the catalytic e ﬃ ciency but surpass them in stability have a great potential for the practical application. In this work, we have developed a novel kind of organic catalysts based on ﬂow-through macroporous polymer monoliths containing catalytic centers that mimic the catalytic site of natural enzyme chymotrypsin. It is known that chymotrypsin catalytic center consists of L-serine, L-histidine, and L-aspartic acid and has speciﬁcity to C-terminal residues of hydrophobic amino acids (L-phenylalanine, L-tyrosine, and L-tryptophan). In this paper, we have prepared the macroporous polymer monoliths bearing grafted polymer layer on their surface. The last one was synthesized via copolymerization of N-methacryloyl-L-serine, N-methacryloyl-L-histidine, and N-methacryloyl-L-aspartic acid. The spatial orientation of amino acids in the polymer layer, generated on the surface of monolithic framework, was achieved by coordinating amino acid-polymerizable derivatives with cobalt (II) ions without substrate-mimicking template and with its use. The conditions for the preparation of mimic materials were optimized to achieve a mechanically stable system. Catalytic properties of the developed systems were evaluated towards the hydrolysis of ester bond in a low molecular substrate and compared to the results of using chymotrypsin immobilized on the surface of a similar monolithic framework. The e ﬀ ect of ﬂow rate increase and temperature elevation on the hydrolysis e ﬃ ciency were evaluated for both mimic monolith and column with immobilized enzyme. visualization—E.K.-V.; scientiﬁc supervision—E.K.-V.; writing—original E.K.-V.; writing—review and M.S.,


Introduction
Being biocatalysts, enzymes are widely used in many technological processes including food and drug industry, biotechnology, pulp and paper industry, etc. [1,2]. The catalytic efficiency of enzymes is very high. At the same time, the protein nature of enzymes is very sensitive to some conditions (organic solvents, pH, elevated temperature, etc.) and the catalytic activity can be easily lost due to distortion of the active site conformation. In addition, the time-consuming production and isolation of enzymes cause their high cost. In this regard, the development of artificial synthetic enzymes or moieties mimicking enzyme active site is the focus of many researches [3,4].
Currently, a number of enzyme mimics based on inorganics [5] or metal-organic frameworks (MOFs) [6,7] have been reported possessing with catalytic activity similar to that of enzymes. Mainly, inorganic systems and MOFs simulate the activity of oxidoreductases like peroxidase, superoxide In order to prepare macroporous monolithic MICs, we selected the chymotrypsin-mimic system and considered two strategies for generation catalytic centers. Co(II)-coordinated monomers containing hydroxyl, carboxyl, and imidazole groups and N-nicotinoyltyrosylbenzyl ester as a template were applied to create the chymotrypsin-mimicking active site. The catalytic properties of the developed mimics were evaluated with the use of Z-L-tyrosine p-nitrophenyl ester (Z-Tyr-OPNP) as known specific low molecular substrate for chymotrypsin. In addition, the catalytic properties of mimic monoliths were compared to those determined for chymotrypsin immobilized on the surface of monolithic column.

Preparation and Characterization of Macroporous Monolithic MICs
In contrast to the preparation of coatings and particles, the synthesis of macroporous monoliths is more complicated task. The polymerization mixture for the preparation of continuous media must satisfy to several requirements. First of all, besides monomers, template, and initiator, the polymerization mixture must contain porogenic solvents to generate a macroporous structure. To be porogenic, the solvents should be poor for the final polymer. At the same time, all components of the polymerization mixture must form a homogenous system. However, some components (monomers/templates/initiator) are solid compounds, which must be dissolved in the appropriate solvents. In turn, the introduction of a good solvent into the polymerization mixture can negatively affect the formation of a macroporous structure due to the increase in system compatibility. Finally, macroporous monoliths should be mechanically stable to resist the flow of the mobile phase.
It is well known that the catalytic site of chymotrypsin is formed by three key amino acids: Ser 195 , Asp 102 , and His 57 . Therefore, the catalytic function of the enzyme is provided by three main functional groups: hydroxyl, carboxyl, and imidazole ( Figure 1). The mechanism of enzymatically catalyzed hydrolysis of low-molecular-amino acid-based substrates is well known [32,33]. A schematic representation of the hydrolysis of N-modified tyrosine p-nitrophenyl ester catalyzed by chymotrypsin is shown in Figure 1.
Initially, attempts were made to synthesize nonimprinted macroporous monoliths using N-methacryloyl-L-histidine, methacrylic acid, 2-hydroxyethyl methacrylate, and EDMA/PEGDA, self-assembled into pre-polymerization complex with Co(II) ions. The tested compositions can be found in Table S1 of the Supplementary materials. Unfortunately, some of the tested polymerization mixtures were unstable and quickly delaminated, while others did not provide the necessary macroporous structure. The main problem when preparing a target catalytic system using a one-step approach was related to poor compatibility of components. The selected MAA and EDMA are moderately hydrophobic liquids, while CoCl 2 and MA-His are solid compounds, which are soluble in water and methanol, respectively. EDMA and MAA are poorly compatible with water so its amount should not be very high to avoid liquid-liquid separation. In turn, the introduction of hydrophobic solvents like toluene or 1-dodecanol as porogens favored to the rapid delamination of the system. Among the tested polymerization systems, only one consisting of MA-His/MAA/HEMA/EDMA = 10/10/10/70 (% mol/mol/mol/mol) (monomers), 10 mol% of CoCl 2 and 1,4-butanediol/methanol/water = 35/43/22 (% v/v/v) (porogens) allowed the formation of monoliths with necessary characteristics: the average pore size and porosity were equal to 1090 ± 85 nm and 90%, respectively. However, the mechanical properties of these monoliths were found to be inappropriate.
Catalysts 2019, 9, x FOR PEER REVIEW 4 of 20 monoliths with necessary characteristics: the average pore size and porosity were equal to 1090 ± 85 nm and 90%, respectively. However, the mechanical properties of these monoliths were found to be inappropriate. In order to overcome that problem, another approach to obtain macroporous monolithic mimics was proposed. The second way was based on a two-step approach, which included the preparation of macroporous monolithic framework with subsequent grafting of the functional layer on the surface of monolithic skeleton ( Figure 3). In this case, the macroporous monolithic framework must be rigid and have a suitable pore size. All of the above advantages of monoliths can be realized with an average pore size of at least 1 m. Since grafting the polymer layer on the walls of macropores can reduce the pore size, the initial pore size in monolithic framework should be not less than 1.5 m. Given that framework performs only a supporting function, to simplify the preparation of monoliths, we used EDMA, DEGDMA, or PEGDA as a single monomer. Detailed information about the composition of the polymerization mixture and the characteristics of monoliths is given in Table In order to overcome that problem, another approach to obtain macroporous monolithic mimics was proposed. The second way was based on a two-step approach, which included the preparation of macroporous monolithic framework with subsequent grafting of the functional layer on the surface of monolithic skeleton (Figure 3). In this case, the macroporous monolithic framework must be rigid and have a suitable pore size. All of the above advantages of monoliths can be realized with an average pore size of at least 1 µm. Since grafting the polymer layer on the walls of macropores can reduce the pore size, the initial pore size in monolithic framework should be not less than 1.5 µm. Given that framework performs only a supporting function, to simplify the preparation of monoliths, we used EDMA, DEGDMA, or PEGDA as a single monomer. Detailed information about the composition of the polymerization mixture and the characteristics of monoliths is given in Table  S2 of the Supplementary materials. The monolith based on poly(ethylene glycol dimethacrylate) (poly(EDMA)) synthesized using 1-dodecanol/toluene = 70/30 (% v/v) as a porogenic mixture (sample 4,  Grafting a polymer layer onto a macroporous polymer framework requires a certain amount of unreacted double bonds on the monolith surface. In this regard, it was necessary to choose the time of polymerization of the framework to achieve the required mechanical stability at dynamic flow of the mobile phase, on the one hand, and to preserve some double bonds for grafting, on the another one. In order to select the optimal time, the materials were prepared within 3, 4, 5, 6, and 9 h of polymerization. Residual surface double bonds were determined by Raman spectroscopy (Figure 4 and Figure S1 of the Supplementary Materials), where the characteristic band at 1639 cm −1 corresponded to the valent oscillation of the C=C group. As seen from the data presented in Figure 4, the amount of residual surface double bonds varied from 2.9% to 2.2% at the polymerization time of 3 and 9 h, respectively. Polymerization within 5 h was considered as optimal, since the shorter polymerization time did not provide the required rigidity (the monolith was compressed by 15-20% after a day of washing at flow rate of 1.0 mL·min −1 ). At the same time, after 5 h polymerization, about 2.5% of double bonds remained at monolith surface. The developed macroporous framework was characterized by an average pore size and total porosity of 1590 ± 80 nm and 82 ± 5%, respectively. Grafting a polymer layer onto a macroporous polymer framework requires a certain amount of unreacted double bonds on the monolith surface. In this regard, it was necessary to choose the time of polymerization of the framework to achieve the required mechanical stability at dynamic flow of the mobile phase, on the one hand, and to preserve some double bonds for grafting, on the another one. In order to select the optimal time, the materials were prepared within 3, 4, 5, 6, and 9 h of polymerization. Residual surface double bonds were determined by Raman spectroscopy (Figure 4 and Figure S1 of the Supplementary Materials), where the characteristic band at 1639 cm −1 corresponded to the valent oscillation of the C=C group. As seen from the data presented in Figure 4, the amount of residual surface double bonds varied from 2.9% to 2.2% at the polymerization time of 3 and 9 h, respectively. Polymerization within 5 h was considered as optimal, since the shorter polymerization time did not provide the required rigidity (the monolith was compressed by 15-20% after a day of washing at flow rate of 1.0 mL·min −1 ). At the same time, after 5 h polymerization, about 2.5% of double bonds remained at monolith surface. The developed macroporous framework was characterized by an average pore size and total porosity of 1590 ± 80 nm and 82 ± 5%, respectively.   Optimization of functional layer grafting was performed using MA-His, MAA, and HEMA as sources of imidazole, carboxyl, and hydroxyl functionalities and CoCl2 as source of cobalt ions for coordination of the monomers. EDMA was used as a cross-linking monomer. At this step, we did not use the template. The formation of a Co 2+ -coordinated assembly of stoichiometric amounts of monomers containing hydroxyl, carboxyl, and imidazole groups was previously testified by Lele at al. with the use of electron spin resonance (ESR) method [19]. As a result, formed assembly undergoes polymerization as a monomeric Co 2+ -coordinated complex but not a random mixture of the three monomers.
The composition of the polymerization mixtures and the polymerization time used for grafting the functional layer are presented in Table 1. Methanol was used as a solvent for monomers and initiator. In all cases, we used equal amounts of functional monomers, whereas the amount of cross-linking monomer was different. Initially, we used equimolar quantities of components to prepare the grafted polymer layer. In this case, more or less appropriate grafting was achieved after 15 h of polymerization (sample 1). Two-three times better grafting of the functional layer was achieved with a threefold increase in the amount of cross-linking agent and a reduction in the polymerization time (sample 3). As expected, the grafting of the polymer layer on the surface of    Optimization of functional layer grafting was performed using MA-His, MAA, and HEMA as sources of imidazole, carboxyl, and hydroxyl functionalities and CoCl2 as source of cobalt ions for coordination of the monomers. EDMA was used as a cross-linking monomer. At this step, we did not use the template. The formation of a Co(II)-coordinated assembly of stoichiometric amounts of monomers containing hydroxyl, carboxyl, and imidazole groups was previously testified by Lele at al. with the use of electron spin resonance (ESR) method [19]. As a result, formed assembly undergoes polymerization as a monomeric Co(II)-coordinated complex but not a random mixture of the three monomers.
The composition of the polymerization mixtures and the polymerization time used for grafting the functional layer are presented in Table 1. Methanol was used as a solvent for monomers and initiator. In all cases, we used equal amounts of functional monomers, whereas the amount of cross-linking monomer was different. Initially, we used equimolar quantities of components to prepare the grafted polymer layer. In this case, more or less appropriate grafting was achieved after 15 h of polymerization (sample 1). Two-three times better grafting of the functional layer was achieved with a threefold increase in the amount of cross-linking agent and a reduction in the polymerization time (sample 3). As expected, the grafting of the polymer layer on the surface of monoliths was accompanied by the slight decrease in their average pore size while the total porosity Optimization of functional layer grafting was performed using MA-His, MAA, and HEMA as sources of imidazole, carboxyl, and hydroxyl functionalities and CoCl 2 as source of cobalt ions for coordination of the monomers. EDMA was used as a cross-linking monomer. At this step, we did not use the template. The formation of a Co(II)-coordinated assembly of stoichiometric amounts of monomers containing hydroxyl, carboxyl, and imidazole groups was previously testified by Lele at al. with the use of electron spin resonance (ESR) method [19]. As a result, formed assembly undergoes polymerization as a monomeric Co(II)-coordinated complex but not a random mixture of the three monomers.
The composition of the polymerization mixtures and the polymerization time used for grafting the functional layer are presented in Table 1. Methanol was used as a solvent for monomers and initiator. In all cases, we used equal amounts of functional monomers, whereas the amount of cross-linking monomer was different. Initially, we used equimolar quantities of components to prepare the grafted polymer layer. In this case, more or less appropriate grafting was achieved after 15 h of polymerization (sample 1). Two-three times better grafting of the functional layer was achieved with a threefold increase in the amount of cross-linking agent and a reduction in the polymerization time (sample 3). As expected, the grafting of the polymer layer on the surface of monoliths was accompanied by the slight decrease in their average pore size while the total porosity kept at the same level. In particular, the average pore size and porosity for samples 1, 2, and 3 were found to be 1480 ± 20, 1430 ± 50, and 1410 ± 20 nm and 81 ± 1, 75 ± 4 and 78 ± 4%, respectively. * the maximal quantities of components, which were managed to dissolve in 1.5 mL of methanol; ** determined by quantitative analysis of His released after acidic hydrolysis of amide bond between MA and His moieties (for details see Section 3.2.4).
The developed approach was applied to the synthesis of mimic layers on the surface of a macroporous monolithic material using MA-Asp and MA-Ser as functional monomers. Additionally, N-nicotinoyltyrosylbenzyl ester was introduced into polymerization mixture as a template. Functional monomers, template, and CoCl 2 were taken in equimolar quantities. Since the part of solid compounds, which must be dissolved for introduction into polymerization, was increased, the complete dissolution of the monomers and template in methanol was achieved at about 4 times reduced concentration. The polymerization was carried out for 8 h as previously optimized. The compositions of polymerization mixtures, as well as the characteristics of the prepared materials are given in Table 2. One can see that the introduction of Co(II) ions (nonimprinted catalyst (NIC)) as well as the template (molecularly imprinted catalyst (MIC)) favored to the higher inclusion of amino acid-based monomers in the grafted polymer layer. This result can be explained by the copolymerization of monomers in the Co(II)-assembled complex in the case of mimic systems (MIC and NIC) in contrast to the polymer formed from a disordered monomer mixture (control sample). Reducing the concentration of monomers provided the formation of a fairly thin grafted layer, which practically did not affect the average pore size. The average pore size was found to be 1560-1570 nm for all systems.
* polymerization time was 8 h; ** determined by quantitative analysis of amino acids released after acidic hydrolysis of amide bond between MA and amino acid moieties (for details see Section 3.2.4); the results of the analysis are given in terms of the number of catalytic sites where 3 equivalents of amino acids correspond to 1 equivalents of catalytic centers.
In addition, the morphology of the samples obtained was studied by the scanning electron microscopy (SEM). Figure 5 shows images of the surfaces of a macroporous framework and MIC monolithic material. The obtained images testify that all the monoliths were characterized by a homogeneous and fairly similar structure with a large number of interconnected macropores.
The mechanical stability of materials towards the flow of the mobile phase is one of the key properties of the flow-through systems. For monoliths, it is known that increasing the flow rate of the mobile phase increases the efficacy of dynamic interfacial processes, such as separation performance and the efficiency of biocatalytic processes based on immobilized enzymes [35,36]. The range of appropriate operating flow rates is determined by the linear interval of backpressure dependence on the flow rate of a mobile phase. Significant deviations from linearity, caused by backpressure growth,   Figure 6 illustrates the dependence of the backpressure on the flow rate for the developed MIC, NIC, control, and immobilized enzyme (IME) columns. A mixture of acetonitrile/water = 60/40 (% v/v) was used as the mobile phase. All tested materials showed a linear dependence of backpressure on flow rate up to 8.5 mL·min −1 . At this flow rate, the backpressure did not exceed 7.0 MPa. The result obtained clearly indicated that all materials were characterized by high mechanical stability towards the flow of mobile phase that agrees with tendency known for macroporous monoliths [35]. Moreover, this indirectly confirms that the grafted functional polymer is covalently fixed on the walls of the framework, but is not located inside the pores. Otherwise, the deviation from linearity would have occurred much earlier, and a significant increase in backpressure would have been observed due to the loss of material permeability caused by clogging the macropores by a polymer not attached to the framework.
A comparison of the tested columns discovers that the MIC column demonstrated an earlier deviation from linearity than other columns. This fact can be explained by the distortion of grafted layer whose quantity in MIC was higher than in NIC monolith. A comparison of the tested columns discovers that the MIC column demonstrated an earlier deviation from linearity than other columns. This fact can be explained by the distortion of grafted layer whose quantity in MIC was higher than in NIC monolith.

Study of Catalytical Properties
It is well known that chymotrypsin has specificity to ester and amide bonds formed with L-isomers of tyrosine, phenylalanine, and tryptophan [37]. In this study, hydrolytic activity of the prepared monoliths was studied towards chymotrypsin's low molecular substrate Z-Tyr-OPNP (Figure 7). In addition, Fmoc-Ala-OPNP (Figure 7) was used as a nonspecific substrate to evaluate the specificity of the developed chymotrypsin mimic systems towards aromatic amino acid-based substrates. A comparison of the tested columns discovers that the MIC column demonstrated an earlier deviation from linearity than other columns. This fact can be explained by the distortion of grafted layer whose quantity in MIC was higher than in NIC monolith.

Study of Catalytical Properties
It is well known that chymotrypsin has specificity to ester and amide bonds formed with L-isomers of tyrosine, phenylalanine, and tryptophan [37]. In this study, hydrolytic activity of the prepared monoliths was studied towards chymotrypsin's low molecular substrate Z-Tyr-OPNP ( Figure 7). In addition, Fmoc-Ala-OPNP (Figure 7) was used as a nonspecific substate to evaluate the specificity of the developed chymotrypsin mimic systems towards aromatic amino acid-based substrates.

Z-Tyr-OPNP
Fmoc-Ala-OPNP Comparison of the substrate hydrolysis rate for MIC and NIC performed at 37 C for 240 min to noncatalyzed system at the same conditions testified that both mimic monoliths surpassed the noncatalyzed hydrolysis ( Figure 8A). Previously it was established that for monoliths an increase in the flow rate often leads to the improvement of efficiency of interphase mass transfer [36,38]. Figure  8B illustrates the efficiency of Z-Tyr-OPNP hydrolysis for reactions catalyzed with the use of MIC monolith performed at 37 C for 60 min at two flow rates and fixed ratio of molar amounts of catalytic centers to substrate ([Cat]/[S]). As expected, the efficiency of the MIC system was strongly affected by the flow rate: an increase in the flow rate shifted the reaction balance towards product formation by enhancing the number of effective contacts with an increase in the number of recirculation cycles. Catalysis using MIC at a flow rate of 0.5 mL·min −1 was less effective than that measured for free chymotrypsin in solution ( Figure 8B). However, at a flow rate of 1.5 mL·min −1 , the hydrolysis efficacy detected for the MIC monolith was higher than for the native enzyme. Comparison of the substrate hydrolysis rate for MIC and NIC performed at 37 • C for 240 min to noncatalyzed system at the same conditions testified that both mimic monoliths surpassed the noncatalyzed hydrolysis ( Figure 8A). Previously it was established that for monoliths an increase in the flow rate often leads to the improvement of efficiency of interphase mass transfer [36,38]. Figure 8B illustrates the efficiency of Z-Tyr-OPNP hydrolysis for reactions catalyzed with the use of MIC monolith performed at 37 • C for 60 min at two flow rates and fixed ratio of molar amounts of catalytic centers to substrate ([Cat]/[S]). As expected, the efficiency of the MIC system was strongly affected by the flow rate: an increase in the flow rate shifted the reaction balance towards product formation by enhancing the number of effective contacts with an increase in the number of recirculation cycles. Catalysis using MIC at a flow rate of 0.5 mL·min −1 was less effective than that measured for free chymotrypsin in solution ( Figure 8B). However, at a flow rate of 1.5 mL·min −1 , the hydrolysis efficacy detected for the MIC monolith was higher than for the native enzyme. Considering that heterogeneous catalysis may differ from homogenous one, we prepared a column with immobilized chymotrypsin (IME column) and applied it in flow-through mode to compare the catalysis with mimics under the same conditions. Figure 9 shows the dependences of the Z-Tyr-OPNP hydrolysis rates on time plotted for IME, MIC, and NIC monolithic columns. Both the MIC and NIC monoliths showed the comparable efficiency, which was twice lower than that for IME column. Considering that heterogeneous catalysis may differ from homogenous one, we prepared a column with immobilized chymotrypsin (IME column) and applied it in flow-through mode to compare the catalysis with mimics under the same conditions. Figure 9 shows the dependences of the Z-Tyr-OPNP hydrolysis rates on time plotted for IME, MIC, and NIC monolithic columns. Both the MIC and NIC monoliths showed the comparable efficiency, which was twice lower than that for IME column.
Considering that heterogeneous catalysis may differ from homogenous one, we prepared a column with immobilized chymotrypsin (IME column) and applied it in flow-through mode to compare the catalysis with mimics under the same conditions. Figure 9 shows the dependences of the Z-Tyr-OPNP hydrolysis rates on time plotted for IME, MIC, and NIC monolithic columns. Both the MIC and NIC monoliths showed the comparable efficiency, which was twice lower than that for IME column.  Figure 10 illustrates the dependence of the product amount on time at different catalyst-to-substrate ratios. The enzyme-like behavior was observed for both MIC and NIC systems, that allow us to apply the Michaelis-Menten theory to determine the maximum reaction velocity (Vmax) and the Michaelis constant (KM). According to the theory for enzymatically catalyzed processes, a rapid increase in the rate of hydrolysis is observed at the first step of the reaction [37]. After that, the active sites of the enzyme become occupied by substrate molecules, and as a result, the hydrolysis rate reaches a plateau.  Figure 10 illustrates the dependence of the product amount on time at different catalyst-to-substrate ratios. The enzyme-like behavior was observed for both MIC and NIC systems, that allow us to apply the Michaelis-Menten theory to determine the maximum reaction velocity (V max ) and the Michaelis constant (K M ). According to the theory for enzymatically catalyzed processes, a rapid increase in the rate of hydrolysis is observed at the first step of the reaction [37]. After that, the active sites of the enzyme become occupied by substrate molecules, and as a result, the hydrolysis rate reaches a plateau.
In order to determine the kinetic parameters for Z-Tyr-OPNP hydrolysis using MIC and NIC systems, the Michaelis-Menten plots were built. The reaction of Z-Tyr-OPNP hydrolysis was carried out with the use of a 20-100-fold molar excesses of the substrate relative to the catalytic centers in the recirculation mode ( Figure S2, the Supplementary Materials). The Michaelis constant (K M ) and the maximum reaction velocity (V max ) were determined using a graphical approach based on linearization by Lineweaver-Burk (double reciprocal coordinates, Figure S3, the Supplementary Materials) and Hanes methods ([S]/V on [S] plot). A monolithic column with immobilized chymotrypsin was used as a control. The summarized data on K M and V max are presented in Table 3. Both mimic systems demonstrated comparable catalytic properties but were inferior in those to natural biocatalyst. The results for the control system were random and poorly reproducible and did not give a Michaelis-Menten plot. An attempt to linearize the data set obtained for the control system gave the positive intercept of x axis. Thus, the Michaelis-Menten theory was not applicable for the disordered polymer system prepared without pre-coordination of functional monomers to mimic the active center of the enzyme. In order to determine the kinetic parameters for Z-Tyr-OPNP hydrolysis using MIC and NIC systems, the Michaelis-Menten plots were built. The reaction of Z-Tyr-OPNP hydrolysis was carried out with the use of a 20-100-fold molar excesses of the substrate relative to the catalytic centers in the recirculation mode ( Figure S2, the Supplementary Materials). The Michaelis constant (KM) and the  It is known that in contrast to the Lineweaver-Burk method, which suffers from a great magnification of an error in determining the linear reaction velocity at a low substrate concentration, the Hanes method is less prone to this obstacle [39]. The results obtained by Hanes method seem to be more trustworthy. Indeed, in other experiments, IME column demonstrated higher hydrolysis rate than MIC. Obviously, it should have had a higher maximal reaction velocity. In turn, it seems reasonable that MIC, prepared with the use of the template, had K M close to the enzyme and lower than for the NIC system. Increasing the flow rate to 1.5 mL·min −1 favored to the enhancement of catalytic activity ( Figure 11) that is in agreement with the experiment performed at fixed [Cat]/[S] ratio ( Figure 8B).
To evaluate the ability of mimics to hydrolyze p-nitrophenyl esters of other amino acid derivatives, Fmoc-Ala-OPNP was used as a substrate. The comparison of the hydrolysis efficiency of p-nitrophenyl esters of Tyr-and Ala-derivatives by MIC and NIC monoliths allowed the conclusion that for both mimics, the hydrolysis of the Tyr-based ester was more effective than Ala-based one ( Figure 12A). Furthermore, the application of the molecular imprinting technique, involving the use of a hydrophobic template contributed to the formation of a kind of "hydrophobic pocket" (Figure 3), resulted in the more pronounced hydrolysis efficiency by MIC towards Fmoc-Ala-OPNP. However, even for the MIC column, the determined maximum reaction velocity of the hydrolysis reaction of the Ala-based substrate by MIC was twice lower than for Tyr-based one ( Figure 12B) that indicates a higher selectivity of MIC to Tyr-based substrates similar to the applied template. To evaluate the ability of mimics to hydrolyze p-nitrophenyl esters of other amino acid derivatives, Fmoc-Ala-OPNP was used as a substrate. The comparison of the hydrolysis efficiency of p-nitrophenyl esters of Tyr-and Ala-derivatives by MIC and NIC monoliths allowed the conclusion that for both mimics, the hydrolysis of the Tyr-based ester was more effective than Ala-based one ( Figure 12A). Furthermore, the application of the molecular imprinting technique, involving the use of a hydrophobic template contributed to the formation of a kind of "hydrophobic pocket" (Figure  3), resulted in the more pronounced hydrolysis efficiency by MIC towards Fmoc-Ala-OPNP. However, even for the MIC column, the determined maximum reaction velocity of the hydrolysis that for both mimics, the hydrolysis of the Tyr-based ester was more effective than Ala-based one ( Figure 12A). Furthermore, the application of the molecular imprinting technique, involving the use of a hydrophobic template contributed to the formation of a kind of "hydrophobic pocket" (Figure  3), resulted in the more pronounced hydrolysis efficiency by MIC towards Fmoc-Ala-OPNP. However, even for the MIC column, the determined maximum reaction velocity of the hydrolysis reaction of the Ala-based substrate by MIC was twice lower than for Tyr-based one ( Figure 12B) that indicates a higher selectivity of MIC to Tyr-based substrates similar to the applied template. It is known that natural biocatalysts are highly effective but only at certain pH and temperature optimum. Deviation from the optimal values can lead to deterioration of the catalytic properties or even their total disappearance due to the distortion of the spatial structure of the enzyme [40]. As seen from Figure 13, an elevation of temperature from 37 to 50 °C increased the catalytic efficiency of the immobilized enzyme during 40 min of the reaction. However, after 60 min, the rate of hydrolysis It is known that natural biocatalysts are highly effective but only at certain pH and temperature optimum. Deviation from the optimal values can lead to deterioration of the catalytic properties or even their total disappearance due to the distortion of the spatial structure of the enzyme [40]. As seen from Figure 13, an elevation of temperature from 37 to 50 • C increased the catalytic efficiency of the immobilized enzyme during 40 min of the reaction. However, after 60 min, the rate of hydrolysis slowed down and finally decreased from 55% to 45% after 240 min. In turn, as expected, the artificial mimic system demonstrated better stability at elevated temperature. In the case of the MIC monolith, the rate of hydrolysis was doubled with increasing temperature. After 240 min, the hydrolysis rate reached 36% at 50 • C instead of 18% at 37 • C. This result may be related to the enhanced diffusivity of the substrate molecules at elevated temperatures. Contrary to the enzyme active site, the artificial sites are stable and, thus, do not depend on temperature.  Evaluation of the stability of the MIC column over time showed that the catalytic activity remained constant even after 40 catalytic cycles that indicates high regeneration of the catalyst after use and good stability of artificial catalytic sites. Evaluation of the stability of the MIC column over time showed that the catalytic activity remained constant even after 40 catalytic cycles that indicates high regeneration of the catalyst after use and good stability of artificial catalytic sites.
Macroporous monoliths were synthesized inside the stainless-steel housings of 50 mm × 4.6 mm produced by Supelco (Bellefonte, PA, USA). All mobile phases were purified by ultrafiltration through nylon filters with pore size of 0.45 µm produced by Merck Millipore (Moscow, Russia). All sample solutions loaded into the monolithic columns were filtered with the use of MCE or PTFE syringe filters with pore size of 0.22 µm produced by Jet Biofil (Guangzhou, China).

Synthesis of Amino Acid-based Monomers (MA-His, MA-Ser, and MA-Asp)
The synthesis of MA-His was carried out according to the protocol published elsewhere [41] with some changes. Briefly, 5.00 g of L-histidine (0.032 mol), 0.20 g of sodium nitrite (0.003 mol), and 30 mL of potassium carbonate solution (water/saturated solution of K 2 CO 3 = 95/5 (vol%/vol%)) adjusted to a value of 9.8 (control by a pH meter) with 0.83 g of sodium hydroxide (0.021 mol) were placed into a three-necked flask with a volume of 100 mL. To remove oxygen from the reaction mixture, the synthesis was preceded by passing an inert gas (argon) through the reaction mixture for 10-15 min. Then, the reaction flask was placed in an ice bath (0-2 • C) and 4 mL (0.041 mol) of methacryloyl chloride was added dropwise under stirring of reaction mixture. After that, the reaction was carried out with constant stirring for 2 h at room temperature. At the end of the reaction, the pH of the reaction mixture was adjusted to 7 with 4 M aqueous sodium hydroxide solution. Thereafter, the reaction mixture was extracted to remove unreacted methacryloyl chloride in portions (1/2 of the resulting mixture) with 20 mL of ethyl acetate for each portion. As a result of three repetitions of the extraction, a yellow aqueous solution was obtained, and then, it was evaporated to dryness on a rotary evaporator. Furthermore, 10 mL of methanol was added to the flask to dissolve the product and separate it from salts by filtration. The product dissolved in methanol was precipitated in 750 mL of fry acetone. After two reprecipitations, the isolated yellow hygroscopic solid was placed in a vacuum dryer without heating for 12-14 h. The product yield was 49%. IR (KBr, cm MA-Asp was synthesized and purified as described for MA-His using 0.032 mol of L-aspartic acid. In this case, a twofold increase in the amount of sodium hydroxide was required to neutralize the β-carboxyl group of aspartic acid and achieve the required pH value. The product yield was 59%. IR (KBr, cm

Synthesis of Template and Substrates
Template (N-nicotinoyltyrosylbenzyl ester) and substrate (Z-Tyr-OPNP) were synthesized without changes as described elsewhere [19].
The template yield was 35%. IR (KBr, cm Fmoc-Ala-OPNP was synthesized using the following procedure. For this, 1.5567 g (0.005 mol) of Fmoc-Ala and 0.6951 g (0.005 mol) of p-nitrophenol were placed and dissolved in 25 mL of ethyl acetate freshly distilled over sodium sulfate. The mixture obtained was added to the solution of 1.0599 g (0.005 mol) DCC in 20 mL of freshly distilled ethyl acetate and the reaction medium was left overnight at room temperature under stirring. After that, the formed dicyclohexylurea precipitate was filtered off with a Schott filter and the resulting yellow solution was evaporated on a rotor until a precipitate began to form. Then, the mixture was left at room temperature until the precipitation stopped and the resulting precipitate was separated using a Schott filter. Then, the product on the Schott filter was washed with 50 mL of cold methanol (4-5 • C) and dried in a vacuum for 24 h. The washing and drying procedures were repeated 3 times. Fmoc-Ala-OPNP yield was 51%.

Synthesis of Macroporous Mimic Monoliths
Precisely, 3 mL of a homogeneous polymerization mixture, consisting of functional and cross-linking monomers, cobalt (II) chloride, porogens, and an initiator (AIBN, 1 wt% from total mass of monomers) was purged for 3-5 min with argon and loaded into a housing for in situ polymerization. The filled housing was placed in a water thermostat and polymerized at 70 • C. For details of one-step mimic synthesis and framework synthesis optimization, see Tables S1 and S2, respectively (the Supplementary Materials). After the polymerization, the column with the obtained polymer was connected to a peristaltic pump and washed with a mixture of methanol/hydrochloric acid = 95/5 (% v/v) until the total removal of cobalt ions (approx. 100 mL) for one-step mimic system. The degree of washing of the material from cobalt ions was monitored spectrophotometrically at 673 nm. The polymer was washed until there was no change in the absorption of the wash solution. After that, the column was washed sequentially with 15 mL of ethanol, 20 mL of ethanol/water = 50/50 (% v/v), and then 20 mL of water. The last procedure was used for washing of monolithic frameworks in the two-step mimic preparation. The flow rate in all cases was 0.5 mL·min −1 . All washing solutions were filtered through syringe filters before use.
Before grafting the mimic polymer layer, a weighed portion of cobalt (II) chloride was dissolved in methanol. After that, the weighed portions of monomers, a template, an aliquot of a cross-linking monomer, and a weighed portion of an initiator (AIBN, 1 wt% from the total mass of monomers) were sequentially introduced into cobalt chloride solution. Homogenization of the mixture after the addition of each component was carried out using an ultrasonic bath. The resulting polymerization mixture was filtered through a Schott filter and purged with argon for 3-5 min to remove oxygen. The prepared mixture was left for 30 min and after that injected with a syringe into macroporous monolithic framework, in situ prepared in the stainless-steal column. The grafting was carried out at 70 • C during different time ( Table 1). The resulting columns were washed using a peristaltic pump in the same way as described above. The control monolith was prepared in a similar manner but without the addition of cobalt (II) chloride.

Characterization of Macroporous Monoliths
The pore characteristics (average pore size and porosity) of all macroporous monoliths were determined based on the data of hydrodynamic permeability calculated according to the Poiseuille-Darcy law as described elsewhere [42].
The surface morphology of the obtained macroporous monoliths was investigated using scanning electron microscopy (SEM) with secondary electron detection. Images were recorded at a magnification 1000× and 5000× using a scanning electron microscope Zeiss AURIGA Laser (Oberkochen, Germany).
The mechanical stability of macroporous monoliths to the mobile phase flow was investigated by the installation of columns into a Shimadzu LC-20 Prominence HPLC system (Kyoto, Japan) and measuring the backpressure at different flow rates of the mobile phase (acetonitrile/water=60/40, % v/v). The flow rate was varied from 0.5 to 10 mL·min −1 .
The quantity of double bonds was determined by Raman spectroscopy with the use of a Nicolet 6700 Fourier spectrometer (Thermo Scientific, Waltham, MA, USA) with an NXR Raman attachment at room temperature. Experimental parameters-spectral resolution 4 cm −1 , number of scans-800, exciting laser-Nd: YVO4, 1064 nm, and laser power-1 W. The spectra were taken from the rounded surfaces of the cylinders of monoliths placed on the bases. Dimensions of cylinders were height-3-5 mm and diameter-2-3 mm. The relative quantity (%) of double bonds was calculated by referring the areas of the characteristic bands for each sample, corresponding to the vibrations of the C=C-groups, to the areas of the characteristic bands of the CH-groups for the corresponding samples.
The quantitative determination of catalytic centers was carried after hydrolysis of grafted functional polymer. For this, 0.1 g, of grafted monolith, prewashed with ethyl alcohol and dried, was ground, placed in a 25 mL round-bottom flask, and refluxed for 24 h with 10 mL of 0.1 M aqueous hydrochloric acid under constant stirring. Then, the solution was cooled to room temperature, the pH was adjusted to 6-7 using an aqueous solution of sodium hydroxide (0.4 M) and evaporated to dryness on a rotary evaporator. Furthermore, 5 mL of water was added to the obtained dry residue and stirred for 1 h at room temperature. The resulting aqueous solution of amino acid(s) was filtered from the polymer particles with a glass fiber syringe filter. The qualitative analysis of amino acid(s) was carried out with use of 2,4,6-trinitrobenzenesulfonic acid (TNBS) specific to primary amino groups. Analysis was performed as follows. Precisely, 1.5 mL of 0.0125 M borate buffer solution (pH = 9.4) and 0.1 mL of 8% aqueous TNBS solution were added to 0.5 mL of the analyzed solution, and kept for 50 min at room temperature. After that, 1.5 mL of 0.01 M citrate buffer solution (pH = 4.8) was added and mixed, and the optical density of the solution was determined at λ = 425 nm against the control sample (prepared in the same way, but with the use of 0.5 mL of water instead of the analyzed solution).
To determine the amount of amino acid(s) in the solution, the corresponding calibration curves were preliminarily built ( Figure S4, the Supplementary Materials).
3.2.5. Chymotrypsin Immobilization 25 mg (10 −6 mol) of the α-chymotrypsin was dissolved in 1.5 mL of 0.0125 M borate buffer solution, pH = 8.5. The resulted solution was loaded into macroporous monolithic column based on poly(glycidylmethacrylate-co-ethylene glycol dimethacrylate). The monolithic column had the following characteristics: length 50 mm, diameter 4.6 mm, average pore size 1590 nm, and void volume 0.67 mL. The immobilization was carried out under shaking at 37 • C for 16 h. After that, the columns with the immobilized enzyme were washed by the same buffer solution, an aqueous solution of hydrochloric acid (pH 2), and again with buffer solution to remove unreacted enzyme. The volume of each washing portion was 3 mL. The amount of bound enzyme was determined by the Lowry-Folin assay. The amount of immobilized enzyme was calculated as difference between quantity of enzyme in initial solution and enzyme quantity determined in solutions after immobilization and washing. The amount of immobilized enzyme was equal to 3.4 mg per column. The average pore size after chymotrypsin immobilization practically did not change.

Evaluation of Catalytic Properties
After preparation of monolithic column, the test materials were washed using a peristaltic pump with 10 mL of acetonitrile/water = 50/50 (% v/v), and then, the liquid was removed from the column and blown with air for 2 min.
A catalytic reaction was carried out by pumping a substrate solution through the column in a recirculation mode using a peristaltic pump. The molar ratios of substrate/catalytic centers ranged from 20 to 100. The volume of reaction medium was equal to 5.6 mL for NIC and IME and 10 mL for NIC. The reaction mixtures for the catalytic reaction represented the substrate solution in acetonitrile/0.005 M sodium phosphate buffer (pH 7.8) = 60/40 (% v/v). The reaction was performed at 37 • C at different flow rates depending on experiment (0.5 or 1.5 mL·min −1 ). Time of catalytic reaction ranged from 10 to 240 min. The amount of the resulting reaction product (p-nitrophenol) was determined using the previously plotted calibration curve ( Figure S5, the Supplementary Materials). The absorbance of p-nitrophenol in solution was measured at a wavelength of 405 nm. The determination of Michaelis constant (K M ) and maximum reaction velocity (V max ) was performed with use of the initial velocity of the catalytic reaction calculated at 40 min. Lineweaver-Burk and Hanes coordinates were used for linearization.

Conclusions
Flow-through macroporous polymer monoliths containing catalytic centers mimicking chymotrypsin active sites were obtained by surface grafting of Co(II)-coordinated assembly of N-methacryloyl-L-serine, N-methacryloyl-L-histidine, and N-methacryloyl-L-aspartic acid as monomers, with or without template. Both kinds of mimics prepared in the presence or absence of the template demonstrated the catalytic enzyme-like properties towards Z-Tyr-OPNP that followed Michaelis-Menten kinetics. Polymerization of the same monomers in the absence of Co-ions provided a disordered polymer layer, which demonstrated poor catalytic properties different from the Michaelis-Menten kinetics. The efficiency of the developed macroporous mimic monoliths depended on the flow rate of recirculation of the substrate solution that is typical for immobilized enzyme reactors based on monoliths. The hydrolysis efficiency for the mimic system operating at a flow rate of 1.5 mL·min −1 was higher than that for soluble enzyme under the same conditions but lower than for the immobilized chymotrypsin operating in flow-through mode. At the same time, an elevation of the temperature to 50 • C led to a decrease in the hydrolysis efficiency for the immobilized enzyme, while it increased for the mimic system. The mimics obtained were characterized by excellent stability of artificial catalytic sites. The catalytic activity kept constant even after 40 catalytic cycles of hydrolysis. Despite the fact that mimics were less efficient than immobilized enzymes in catalysis at standard conditions for the enzyme, they can surpass natural biocatalysts in long-term stability and applicability to harsher reaction conditions. Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4344/10/12/1395/ s1, Table S1: Compositions of polymerization mixtures tested for the one-step synthesis of macroporous polymer monolith containing catalytic centers mimicking chymotrypsin active site. Table S2: Compositions of polymerization mixtures tested for the synthesis of macroporous monolithic framework. Figure S1: Raman spectra of macroporous polyEDMA synthesized during different polymerization time. Figure S2: Scheme of catalysis performance in the recirculation mode. Figure S3: Lineweaver-Burk plots for hydrolysis of Z-L-tyrosine p-nitrophenyl ester (Z-Tyr-OPNP) by macroporous nonimprinted mimic (NIC) (A) and immobilized enzyme (IME) (B) monoliths, and hydrolysis of Fmoc-Ala-OPNP by macroporous imprinted mimic (MIC) monolith (C). Figure S4: Calibration plot for L-His and mixture of amino acids (L-Ala, L-His, and L-Ser) determination with the use of 2,4,6-trinitrobenzenesulfonic acid (TNBS) solution. Figure S5: Calibration plot for p-nitrophenol.

Author Contributions:
Conceptualization