Cholinesterase Inhibitory and Anti-Inflammatory Activity of the Naphtho- and Thienobenzo-Triazole Photoproducts: Experimental and Computational Study

New 1,2,3-triazolo(thieno)stilbenes were synthesized as mixtures of isomers and efficiently photochemically transformed to their corresponding substituted thienobenzo/naphtho-triazoles in high isolated yields. The resulting photoproducts were studied as acetyl- (AChE) and butyrylcholinesterase (BChE) inhibitors without or with interconnected inhibition potential of TNF-α cytokine production. The most promising anti-inflammatory activity was shown again by naphtho-triazoles, with a derivative featuring 4-pentenyl substituents exhibiting notable potential as a cholinesterase inhibitor. To identify interactions between ligands and the active site of cholinesterases, molecular docking was performed for the best potential inhibitors. Additionally, molecular dynamics simulations were employed to assess and validate the stability and flexibility of the protein–ligand complexes generated through docking.


Introduction
A small and simple triazole core is present in compounds that possess antimicrobial, antitumor, antitubercular, antidepressant, anti-inflammatory, and numerous other activities.Triazoles show a wide range of biological activities and are found in many powerful, biologically active compounds such as trazodone (antidepressant), hexaconazole (an antifungal drug), alprazolam (tranquilizer) and others [1].So far, modifications of triazoles have proven to be very effective, and some of the newly prepared derivatives have been observed to have better activity and less toxicity.Drugs are also being developed for the treatment of neurodegenerative diseases such as Alzheimer's disease (AD), with a characteristic 1,2,3-triazole core [2].Triazolostilbenes have been synthesized in the context of nonlinear optical (NLO) materials [3,4].In addition to the synthesis of new industrial materials, such as the triazole derivative of the natural pigment curcumin with significant biological and spectroscopic properties [5][6][7], functionalized trans-1,2,3-triazolostilbenes have been synthesized as suitable derivatives for further photochemical and photophysical research due to easily achievable change in their structures [8].
Non-steroidal anti-inflammatory drugs have a wide therapeutic application in the treatment of various types of inflammatory conditions.The long-term use of these drugs has harmful effects on the gastrointestinal tract and subsequently leads to problems such as kidney damage, stomach ulcers, and hepatotoxicity [9].This is exactly why there is an effort to synthesize new types of these compounds that should be highly effective and have improved safety so inflammations can be suppressed with as few unwanted effects as possible.Numerous studies related to such compounds have been conducted.Considering the biological and medical importance of ibuprofen (Figure 1) and 1,2,3-triazole, new chemical compounds were investigated based on these compounds' previously known biological significance.
Int. J. Mol.Sci.2023, 24, x FOR PEER REVIEW 2 of 22 for further photochemical and photophysical research due to easily achievable change in their structures [8].
Non-steroidal anti-inflammatory drugs have a wide therapeutic application in the treatment of various types of inflammatory conditions.The long-term use of these drugs has harmful effects on the gastrointestinal tract and subsequently leads to problems such as kidney damage, stomach ulcers, and hepatotoxicity [9].This is exactly why there is an effort to synthesize new types of these compounds that should be highly effective and have improved safety so inflammations can be suppressed with as few unwanted effects as possible.Numerous studies related to such compounds have been conducted.Considering the biological and medical importance of ibuprofen (Figure 1) and 1,2,3triazole, new chemical compounds were investigated based on these compounds' previously known biological significance.Ibuprofen-based compounds containing 1,4-disubstituted 1,2,3-triazole were synthesized and examined for in vivo anti-inflammatory activity, and one of the compounds (Figure 1) showed a strong effect as the reference anti-inflammatory drug ibuprofen at the same concentration (10 mg/kg body weight).The compounds in Figure 2 also showed significant anti-inflammatory activity.Their bactericidal profiles were also examined for all mentioned triazole analogs (including the ibuprofen derivative from Figure 1).All compounds showed significant bactericidal activity against Gram-positive and Gramnegative strains of bacteria.It is an interesting fact that the presence of an electron-withdrawing group or atom (NO2 or Cl) in the meta-or para-position (C3 or C4) of the benzyl or phenyl ring on triazole leads to a significant increase in anti-inflammatory and bactericidal activity [10].Potential anti-inflammatory activity was also shown by compounds that contain both triazole and tetrazole rings in their structure (Figure 3).Their anti-inflammatory effect was evaluated by a test with carrageenan (a polysaccharide obtained from red seaweed), and they showed superior anti-inflammatory activity compared to other synthesized compounds [11,12].Ibuprofen-based compounds containing 1,4-disubstituted 1,2,3-triazole were synthesized and examined for in vivo anti-inflammatory activity, and one of the compounds (Figure 1) showed a strong effect as the reference anti-inflammatory drug ibuprofen at the same concentration (10 mg/kg body weight).The compounds in Figure 2 also showed significant anti-inflammatory activity.Their bactericidal profiles were also examined for all mentioned triazole analogs (including the ibuprofen derivative from Figure 1).All compounds showed significant bactericidal activity against Gram-positive and Gram-negative strains of bacteria.It is an interesting fact that the presence of an electron-withdrawing group or atom (NO 2 or Cl) in the metaor para-position (C3 or C4) of the benzyl or phenyl ring on triazole leads to a significant increase in anti-inflammatory and bactericidal activity [10].
Int. J. Mol.Sci.2023, 24, x FOR PEER REVIEW 2 of 22 for further photochemical and photophysical research due to easily achievable change in their structures [8].Non-steroidal anti-inflammatory drugs have a wide therapeutic application in the treatment of various types of inflammatory conditions.The long-term use of these drugs has harmful effects on the gastrointestinal tract and subsequently leads to problems such as kidney damage, stomach ulcers, and hepatotoxicity [9].This is exactly why there is an effort to synthesize new types of these compounds that should be highly effective and have improved safety so inflammations can be suppressed with as few unwanted effects as possible.Numerous studies related to such compounds have been conducted.Considering the biological and medical importance of ibuprofen (Figure 1) and 1,2,3triazole, new chemical compounds were investigated based on these compounds' previously known biological significance.Ibuprofen-based compounds containing 1,4-disubstituted 1,2,3-triazole were synthesized and examined for in vivo anti-inflammatory activity, and one of the compounds (Figure 1) showed a strong effect as the reference anti-inflammatory drug ibuprofen at the same concentration (10 mg/kg body weight).The compounds in Figure 2 also showed significant anti-inflammatory activity.Their bactericidal profiles were also examined for all mentioned triazole analogs (including the ibuprofen derivative from Figure 1).All compounds showed significant bactericidal activity against Gram-positive and Gramnegative strains of bacteria.It is an interesting fact that the presence of an electron-withdrawing group or atom (NO2 or Cl) in the meta-or para-position (C3 or C4) of the benzyl or phenyl ring on triazole leads to a significant increase in anti-inflammatory and bactericidal activity [10].Potential anti-inflammatory activity was also shown by compounds that contain both triazole and tetrazole rings in their structure (Figure 3).Their anti-inflammatory effect was evaluated by a test with carrageenan (a polysaccharide obtained from red seaweed), and they showed superior anti-inflammatory activity compared to other synthesized compounds [11,12].Potential anti-inflammatory activity was also shown by compounds that contain both triazole and tetrazole rings in their structure (Figure 3).Their anti-inflammatory effect was evaluated by a test with carrageenan (a polysaccharide obtained from red seaweed), and they showed superior anti-inflammatory activity compared to other synthesized compounds [11,12].On the other side, benzene-triazole derivatives (BTA) are active against a wide range of target molecules.Also, BTA is of biological and industrial importance [13,14].Due to its pharmacological activities, 1H-benzo[d]-1,2,3-triazole is considered a favored structure.Used as a template for the design of new pharmacologically active compounds, BTA is On the other side, benzene-triazole derivatives (BTA) are active against a wide range of target molecules.Also, BTA is of biological and industrial importance [13,14].Due to its pharmacological activities, 1H-benzo[d]-1,2,3-triazole is considered a favored structure.Used as a template for the design of new pharmacologically active compounds, BTA is undergoing rapid development in the synthesis of heterocycles.Research has shown that BTA has various pharmacological activities (anti-inflammatory, antimicrobial, antifungal, and anticarcinogenic).It is also known that molecules containing the benzotriazole core have CNS (Central Nervous System) activity [15].
The first confirmation of the strong connection between the cholinergic and immune systems appeared at the beginning of the last century by researching the effects of CNI-1493 on laboratory mouse models.It has been observed that the anti-inflammatory effects of the compound are indirect and occur due to the compound's activity in the brain and not in the blood (immune cells) [16].Further testing of the compound in animal models showed a dose-dependent reduction in the systemic production of tumor necrosis factor, TNF-α, after endotoxin (LPS) challenge [17].It is known that the cholinergic anti-inflammatory system is a powerful instrument by which the body regulates the amount of inflammation.The key participant in this pathway is the neurotransmitter acetylcholine (ACh).Increased calcium levels in cells in that pathway lead to activation of the transcription factor NFαB and suppression of the immune response [18].This mechanism can be interrupted by cholinesterase (ChE) enzymes that break down ACh, acetyl-(AChE), and butyrylcholinesterase (BChE) enzymes.Inhibitors of these enzymes could also have antiinflammatory effects [19][20][21].AChE and BChE also represent pharmacologically suitable targets in neurodegenerative disorders (Alzheimer's disease, AD), given their physiological roles in the body.Treatment of neurodegenerative disorders currently includes common reversible cholinesterase enzyme inhibitors, such as galantamine, with proven efficacy in improving cognitive function [2,22,23].
In our previous research, we have proven that some of the thienobenzo-triazoles obtained by photochemical cyclization have a dual activity.They showed very good cholinesterase enzyme inhibition interconnected with anti-inflammatory activity (Structure A, Figure 4) [24,25].On the other hand, naphtho-triazoles (Structures B and C, Figure 4) showed a much more potent anti-inflammatory effect but not cholinesterase inhibition.This work evaluated a new series of naphtho/thienobenzo-triazoles to confirm previous experimental and computational data on biological activity.On the other side, benzene-triazole derivatives (BTA) are active against a wide range of target molecules.Also, BTA is of biological and industrial importance [13,14].Due to its pharmacological activities, 1H-benzo[d]-1,2,3-triazole is considered a favored structure.Used as a template for the design of new pharmacologically active compounds, BTA is undergoing rapid development in the synthesis of heterocycles.Research has shown that BTA has various pharmacological activities (anti-inflammatory, antimicrobial, antifungal, and anticarcinogenic).It is also known that molecules containing the benzotriazole core have CNS (Central Nervous System) activity [15].
The first confirmation of the strong connection between the cholinergic and immune systems appeared at the beginning of the last century by researching the effects of CNI-1493 on laboratory mouse models.It has been observed that the anti-inflammatory effects of the compound are indirect and occur due to the compound's activity in the brain and not in the blood (immune cells) [16].Further testing of the compound in animal models showed a dose-dependent reduction in the systemic production of tumor necrosis factor, TNF-α, after endotoxin (LPS) challenge [17].It is known that the cholinergic anti-inflammatory system is a powerful instrument by which the body regulates the amount of inflammation.The key participant in this pathway is the neurotransmitter acetylcholine (ACh).Increased calcium levels in cells in that pathway lead to activation of the transcription factor NFαB and suppression of the immune response [18].This mechanism can be interrupted by cholinesterase (ChE) enzymes that break down ACh, acetyl-(AChE), and butyrylcholinesterase (BChE) enzymes.Inhibitors of these enzymes could also have antiinflammatory effects [19][20][21].AChE and BChE also represent pharmacologically suitable targets in neurodegenerative disorders (Alzheimer's disease, AD), given their physiological roles in the body.Treatment of neurodegenerative disorders currently includes common reversible cholinesterase enzyme inhibitors, such as galantamine, with proven efficacy in improving cognitive function [2,22,23].
In our previous research, we have proven that some of the thienobenzo-triazoles obtained by photochemical cyclization have a dual activity.They showed very good cholinesterase enzyme inhibition interconnected with anti-inflammatory activity (Structure A, Figure 4) [24,25].On the other hand, naphtho-triazoles (Structures B and C, Figure 4) showed a much more potent anti-inflammatory effect but not cholinesterase inhibition.This work evaluated a new series of naphtho/thienobenzo-triazoles to confirm previous experimental and computational data on biological activity.

Synthesis of New Thienobenzo/Naphtho-Triazoles 1-13
New triazolo-stilbenes 1a-7a and triazolo-thienostilbenes 8a-13a (Scheme 1) were prepared by the Wittig reaction from phosphonium bromide and the 1-substituted-1,2,3triazole-4-carbaldehydes, according to the previously described procedure [26].Wittig reaction provided the new triazolo-stilbenes 1a-7a and triazolo-thienostilbenes 8a-13a as mixtures of cisand trans-isomers (isolated yields on the mixtures were 42-78%, Scheme 1) which were not separated in this research into pure geometrical isomers but immediately converted into photoproducts, naphtho-triazoles 1-7 and thienobenzo-triazoles 8-13 (Scheme 1) as new biological targets by an intramolecular electrocyclization reaction in high isolated yields (29-60%, Scheme 1).According to 1 H NMR spectroscopy, the substituent on the triazole ring directed the ratio of geometric isomers in the Wittig reaction (see Materials and Methods).In the 1 H NMR spectra of the mixture of geometrical isomers 1a-7a, it can be seen the resolved patterns for ethylenic protons with the characteristic coupling constants, the signals for the protons on various substituents, and the singlets for the protons on the triazole rings.
New triazolo-stilbenes 1a-7a and triazolo-thienostilbenes 8a-13a (Scheme 1) were prepared by the Wittig reaction from phosphonium bromide and the 1-substituted-1,2,3triazole-4-carbaldehydes, according to the previously described procedure [26].Wittig reaction provided the new triazolo-stilbenes 1a-7a and triazolo-thienostilbenes 8a-13a as mixtures of cis-and trans-isomers (isolated yields on the mixtures were 42-78%, Scheme 1) which were not separated in this research into pure geometrical isomers but immediately converted into photoproducts, naphtho-triazoles 1-7 and thienobenzo-triazoles 8-13 (Scheme 1) as new biological targets by an intramolecular electrocyclization reaction in high isolated yields (29-60%, Scheme 1).According to 1 H NMR spectroscopy, the substituent on the triazole ring directed the ratio of geometric isomers in the Wittig reaction (see Materials and Methods).In the 1 H NMR spectra of the mixture of geometrical isomers 1a-7a, it can be seen the resolved patterns for ethylenic protons with the characteristic coupling constants, the signals for the protons on various substituents, and the singlets for the protons on the triazole rings.In aerobic conditions, preparative irradiations of 1a-13a at 313 nm in toluene solutions gave the naphthotriazoles 1-7 and thienobenzo-triazoles 8-13; after that, they were successfully isolated and fully characterized by NMR spectroscopy and HRMS analyses (see Materials and Methods and Suppl.Material).The formation of the electrocyclization photoproducts 1-13 was generally accompanied by the formation of some high-molecular-weight products, which were not investigated.In their 1 H NMR spectra, it can be seen the disappearance of the ethylenic protons and the singlets on the 1,2,3-triazole rings.
As explained in previous research, the inhibitors of enzyme cholinesterases could also have anti-inflammatory effects [24].For the stated reason, thienobenzo-triazoles 14-31 (Figure 5), which stood out as very good inhibitors of these enzymes, were also tested here for their possible interconnected inhibition potential of TNF-α cytokine production.Most of these molecules showed better inhibition of the enzyme butyrylcholinesterase (BChE) [27].Also, the binding affinity of BChE for even seven new compounds was similar to that reported for common cholinesterase inhibitors.In aerobic conditions, preparative irradiations of 1a-13a at 313 nm in toluene solutions gave the naphthotriazoles 1-7 and thienobenzo-triazoles 8-13; after that, they were successfully isolated and fully characterized by NMR spectroscopy and HRMS analyses (see Materials and Methods and Suppl.Material).The formation of the electrocyclization photoproducts 1-13 was generally accompanied by the formation of some high-molecularweight products, which were not investigated.In their 1 H NMR spectra, it can be seen the disappearance of the ethylenic protons and the singlets on the 1,2,3-triazole rings.
As explained in previous research, the inhibitors of enzyme cholinesterases could also have anti-inflammatory effects [24].For the stated reason, thienobenzo-triazoles 14-31 (Figure 5), which stood out as very good inhibitors of these enzymes, were also tested here for their possible interconnected inhibition potential of TNF-α cytokine production.Most of these molecules showed better inhibition of the enzyme butyrylcholinesterase (BChE) [27].Also, the binding affinity of BChE for even seven new compounds was similar to that reported for common cholinesterase inhibitors.

Anti-Inflammatory Activity of 1,2,3-Triazole Derivatives 1-31
The potential to act anti-inflammatory was evaluated for 1,2,3-triazole derivatives 1-31 in vitro in an assay where PBMCs from 2 healthy donors were stimulated with LPS to induce the immune response.As a measurement of response, the pro-inflammatory cytokine TNF-α was determined.The highest potency in TNF-α production inhibition was observed for naphtho-triazoles 1-4 and 6 (Figure 6) and thienobenzo-triazoles 18-22 and 26-27 (Table 1).However, cell viability was determined following incubation with compounds by measuring ATP levels inside the cells to ensure that the lower TNF-α production did not result from cytotoxicity.While the naphthotriazoles did not affect cell viability, for the most potent thienobenzo-triazoles 19, 20, and 26, cell viability reduction was observed at the highest tested concentration.Since previously studied thienobenzo/naphtho-triazoles displayed activity toward cholinesterase inhibition, further functionalization and examination of these core units were imposed.The inhibition of AChE and BChE was assessed for seven naphtho-tria-   Since previously studied thienobenzo/naphtho-triazoles displayed activity toward cholinesterase inhibition, further functionalization and examination of these core units were imposed.The inhibition of AChE and BChE was assessed for seven naphtho-triazoles and six thienobenzo-triazole derivatives in a wide range of concentrations.If more than 50% inhibition is achieved, the IC 50 value is determined (Table 2).Results were compared to a reference compound galantamine.
Naphtho-triazoles with methoxy substituent at aryl moiety and various substituents at triazole ring were studied previously [25] and showed no significant activity.Change in the type of substituent improved activity since almost all of the tested compounds achieved IC 50 values.Derivatives 1, 2, 3, and 7 showed activity with IC 50 values in similar ranges of concentrations toward both enzymes.Among naphtho-triazoles, structure 3 with pentenyl substituent showed the most potent activity with IC 50 values in a very good range (IC 50 51.3µM for AChE and 53.5 µM for BChE), but still less than galantamine (IC 50 7.9 µM for BChE and 0.15 µM for AChE).Selective inhibition toward AChE was demonstrated by derivatives 5 and 6, with meta substituents at phenyl moiety.The one with m-fluorophenyl substituent 5 achieved a very good IC 50 value, while change of the fluorine atom to para position in derivative 4 diminished inhibitory activity.A similar structure to 1, with a methyl group at aryl moiety, was studied earlier [24], and it can be concluded that the change to methoxy derivative improved activity toward AChE.At the same time, affinity toward BChE remained the same.Previous results show that introducing a thiophene ring in the core of these structures improves their biological activity [25].Numerous thienobenzo-triazoles with different substituents at the triazole ring were examined.They showed higher affinity toward BChE, where some of them, particularly derivative A (Figure 4), achieved better IC 50 values than reference galantamine.Structure 13, examined here, represents one in a series of such derivatives.The results are consistent with the previous ones, i.e., it showed selective inhibition for BChE with a very good IC 50 value.The subsequent modification in the structure of the tested compounds refers to introducing a methyl substituent on the thiophene ring (compounds 8-12).Derivatives with p-fluorophenyl substituent showed the best results among this group, while the same substituent had the opposite effect on naphtho-triazoles (derivative 4).Generally, none of the compounds in this group displayed remarkable results since IC 50 values, where achieved, are in the good to moderate range.A comparison of structure 8 with its analogs without methyl group (structure A) shows that introducing methyl at thiophene decreases the activity for BChE five times and completely reduces activity toward AChE.Reduction of activity toward BChE inhibition by introducing the methyl group was also observed for 9, 10, and 12 compared to their non-methylated analogs [27].The presence of a furan ring at the triazole unit diminished enzyme activity, and 50% of inhibition was not achieved.
Among the tested compounds in this research, the two with pentenyl substituent, 3 and 13, showed the best results, with one significant difference: naphtho derivative is active toward both enzymes while thienobenzo derivative is completely selective toward BChE.The introduction of methyl substituent at the thiophene ring is not a desirable structure modification since no improvement was observed for the inhibition of enzymes.
Concerning the evidence of a strong link between the cholinergic and immune systems [28], it is known that the cholinergic anti-inflammatory system, a part of the vagus nerve system, is a powerful tool by which the organism regulates the inflammation extent to ensure the elimination of the threat and prevent any potential tissue disruption or damage.The key player in this pathway is acetylcholine (ACh).As cholinesterases degrade ACh to acetic acid and choline, the inhibitors of these enzymes could also have anti-inflammatory effects.
As a result of experimental data in this research obtained by testing inhibition of cholinesterases and inhibition of TNF-α cytokine production, it can be concluded that naphtho-triazoles show better anti-inflammatory activity than thienobenzo-triazoles, especially compounds 1, 2, 3 and 6.The best candidate showing the associated inhibition of cholinesterase and anti-inflammatory effect is naphtho-triazole 3, for which the possibility of the described mechanism in which ACh is the link can be assumed.As a result of experimental data in this research obtained by testing inhibition of cholinesterases and inhibition of TNF-α cytokine production, it can be concluded that naphtho-triazoles show better anti-inflammatory activity than thienobenzo-triazoles, especially compounds 1, 2, 3 and 6.The best candidate showing the associated inhibition of cholinesterase and anti-inflammatory effect is naphtho-triazole 3, for which the possibility of the described mechanism in which ACh is the link can be assumed.

Docking and Molecular Dynamics Study of 1,2,3-Triazole Derivatives
To elucidate the interactions underpinning the stability of complexes formed between the triazole derivatives and cholinesterases, we conducted molecular docking on the compounds that exhibited the most potent inhibitory activity in experimental assays.Notably, for AChE, the most promising outcomes were observed with compounds 3 and 5.The structures of complexes obtained by docking 3 and 5 into the active site of AChE are shown in Figure 7a,b, respectively.These structural representations unveil the presence of non-covalent interactions between each ligand and residues within the active site, which are responsible for the complex's stability.The placement of the most stable ligand pose into the active site of AChE for both molecules exhibits a striking similarity, which is to be expected given that their only structural divergence lies in the substituent at the triazole subunit.Both compounds 3 and 5 occupy the peripheral anionic site (PAS), thereby enabling the naphthalene subunit to engage in - stacking interactions with residues Trp286 and Tyr72, and the triazole ring to interact with Tyr341.The substituent at the triazole subunit in compound 3, pentenyl, engages in a hydrophobic interaction with Phe338, while its counterpart in compound 5, m-fluorobenzyl, interacts with the same residue through orthogonal - stacking.No hydrogen bonds are observed; this orientation places the triazole ring in a position where there are no residues that could act as hydrogen bond donors to nitrogen atoms of triazole.The lowest free energies of binding (ΔGbind) estimated by molecular docking of compounds 3 and 5 into AChE are −8.91 and −9.53 kcal mol −1 , respectively (Table S1).Expectedly, these values are not as good as those estimated for the reference ligand galantamine, whose ΔGbind of −10.11 kcal mol −1 was obtained using the same docking procedure.These structural representations unveil the presence of non-covalent interactions between each ligand and residues within the active site, which are responsible for the complex's stability.The placement of the most stable ligand pose into the active site of AChE for both molecules exhibits a striking similarity, which is to be expected given that their only structural divergence lies in the substituent at the triazole subunit.Both compounds 3 and 5 occupy the peripheral anionic site (PAS), thereby enabling the naphthalene subunit to engage in π-π stacking interactions with residues Trp286 and Tyr72, and the triazole ring to interact with Tyr341.The substituent at the triazole subunit in compound 3, pentenyl, engages in a hydrophobic interaction with Phe338, while its counterpart in compound 5, m-fluorobenzyl, interacts with the same residue through orthogonal π-π stacking.No hydrogen bonds are observed; this orientation places the triazole ring in a position where there are no residues that could act as hydrogen bond donors to nitrogen atoms of triazole.The lowest free energies of binding (∆G bind ) estimated by molecular docking of compounds 3 and 5 into AChE are −8.91 and −9.53 kcal mol −1 , respectively (Table S1).Expectedly, these values are not as good as those estimated for the reference ligand galantamine, whose ∆G bind of −10.11 kcal mol −1 was obtained using the same docking procedure.
For BChE, the most potent inhibitory activity was observed with compounds 3 and 13.The structures of the complexes formed by docking molecules 3 and 13 into the active site of BChE are depicted in Figure 8a,b.Despite the difference in the scaffold, where one of the phenyl rings in compound 3 is replaced by thiophene in compound 13, both compounds exhibit similar placement within the active site.
For BChE, the most potent inhibitory activity was observed with compounds 3 and 13.The structures of the complexes formed by docking molecules 3 and 13 into the active site of BChE are depicted in Figure 8a,b.Despite the difference in the scaffold, where one of the phenyl rings in compound 3 is replaced by thiophene in compound 13, both compounds exhibit similar placement within the active site.This alignment enables - stacking interactions between the aromatic scaffold of both compounds and Trp82.Furthermore, on the opposing side of the ligand scaffold, - interaction occurs with His 438, residue belonging to the catalytic triad.In both of these ligands, the substituent at the triazole subunit is the same, 4-pentenyl, and its conformation allows for hydrophobic alkenyl- interaction with Tyr332, a residue located within the PAS of BChE.Docking scores for these systems are presented in Table S2: the lowest free energies of binding estimated by molecular docking of compounds 3 and 13 into BChE are −7.56 and −6.81 kcal mol −1 , respectively, while, for reference, the ligand galantamine ΔGbind was estimated to be −7.49kcal mol −1 .
Additionally, we conducted molecular dynamics (MD) simulations on studied systems to assess the stability of protein-ligand complexes suggested by molecular docking.For each investigated protein-ligand complex, the structure with the lowest estimated free energy of binding was chosen, and a 30 ns MD simulation was performed.Root-meansquare deviation values (RMSD) were computed for each system to gauge the structural alterations of the protein-ligand complex throughout the simulation time.RMSD analysis was conducted on all atoms within each protein-ligand complex, excluding hydrogen atoms.Additionally, RMS fluctuation values (RMSF) are calculated as average quadratic fluctuations of positions of  carbons over the trajectory.The RMSD quantifies how much a structure diverges from the initial one, and the RSMF values here show which parts of the protein backbone are the most mobile.Finally, the gyration radius (Rg) was computed to estimate the proteins' compactness.Rg is calculated as a root-mean-square average of the distance of all atoms from the center of mass of the protein-ligand system.Systems where Rg does not increase significantly during the simulation can be described as compact.The results of MD simulations analysis for complexes of AChE with compounds 3 and 5, respectively, are shown in Figure 9.
After an initial 5-10 ns period, both systems achieve convergence, with no significant fluctuations in subsequent RMSD values.Notably, the complex containing compound 3 This alignment enables π-π stacking interactions between the aromatic scaffold of both compounds and Trp82.Furthermore, on the opposing side of the ligand scaffold, π-π interaction occurs with His 438, residue belonging to the catalytic triad.In both of these ligands, the substituent at the triazole subunit is the same, 4-pentenyl, and its conformation allows for hydrophobic alkenyl-π interaction with Tyr332, a residue located within the PAS of BChE.Docking scores for these systems are presented in Table S2: the lowest free energies of binding estimated by molecular docking of compounds 3 and 13 into BChE are −7.56 and −6.81 kcal mol −1 , respectively, while, for reference, the ligand galantamine ∆G bind was estimated to be −7.49kcal mol −1 .
Additionally, we conducted molecular dynamics (MD) simulations on studied systems to assess the stability of protein-ligand complexes suggested by molecular docking.For each investigated protein-ligand complex, the structure with the lowest estimated free energy of binding was chosen, and a 30 ns MD simulation was performed.Root-meansquare deviation values (RMSD) were computed for each system to gauge the structural alterations of the protein-ligand complex throughout the simulation time.RMSD analysis was conducted on all atoms within each protein-ligand complex, excluding hydrogen atoms.Additionally, RMS fluctuation values (RMSF) are calculated as average quadratic fluctuations of positions of α carbons over the trajectory.The RMSD quantifies how much a structure diverges from the initial one, and the RSMF values here show which parts of the protein backbone are the most mobile.Finally, the gyration radius (Rg) was computed to estimate the proteins' compactness.Rg is calculated as a root-mean-square average of the distance of all atoms from the center of mass of the protein-ligand system.Systems where Rg does not increase significantly during the simulation can be described as compact.The results of MD simulations analysis for complexes of AChE with compounds 3 and 5, respectively, are shown in Figure 9.
an average RMSD of 1.34 Å and a maximum value of 1.58 Å , reducing the difference between them to 0.24 Å .For AChE with ligand 5 (illustrated by the red line in Figure 9a), the average RMSD value is 1.74 Å , accompanied by a maximum RMSD of 2.19 Å , resulting in a 0.45 Å difference.This difference decreases to 0.32 Å when the analysis is limited to the last 20 nanoseconds of the trajectory, where the average RMSD is 1.46 Å with a maximum value of 1.78 Å .In both complexes, the RMSF values obtained during the entire simulation indicate that -carbon positions in the protein backbone exhibit small fluctuations compared to their reference positions.Specifically, in the system containing ligand 3, these fluctuations range from 0.34 to 1.75 Å , while in AChE complexed with ligand 5, -carbon positions vary between 0.32 and 2.13 Å .Similarly, the radius of gyration does not change significantly throughout the simulation; for a complex of AChE with 3, this parameter varies within the range of 22.57 to 22.93 Å , while in the AChE with ligand 5, the values span from 22.60 to 23.08 Å .
To gain a visual understanding of the complexes under study in molecular dynamics (MD), we analyzed the final 20 ns of each simulation to derive their average structure.After an initial 5-10 ns period, both systems achieve convergence, with no significant fluctuations in subsequent RMSD values.Notably, the complex containing compound 3 (represented by the blue line in Figure 9a) exhibits slightly better performance.Across the entire blue trajectory, the average RMSD stands at 1.69 Å, with a maximum value of 2.01 Å, resulting in a 0.32 Å difference between the maximum and average RMSD values.An examination of the last 20 ns of this trajectory, where the system had converged, reveals an average RMSD of 1.34 Å and a maximum value of 1.58 Å, reducing the difference between them to 0.24 Å.For AChE with ligand 5 (illustrated by the red line in Figure 9a), the average RMSD value is 1.74 Å, accompanied by a maximum RMSD of 2.19 Å, resulting in a 0.45 Å difference.This difference decreases to 0.32 Å when the analysis is limited to the last 20 nanoseconds of the trajectory, where the average RMSD is 1.46 Å with a maximum value of 1.78 Å.
In both complexes, the RMSF values obtained during the entire simulation indicate that α-carbon positions in the protein backbone exhibit small fluctuations compared to their reference positions.Specifically, in the system containing ligand 3, these fluctuations range from 0.34 to 1.75 Å, while in AChE complexed with ligand 5, α-carbon positions vary between 0.32 and 2.13 Å.Similarly, the radius of gyration does not change significantly throughout the simulation; for a complex of AChE with 3, this parameter varies within the range of 22.57 to 22.93 Å, while in the AChE with ligand 5, the values span from 22.60 to 23.08 Å.
To gain a visual understanding of the complexes under study in molecular dynamics (MD), we analyzed the final 20 ns of each simulation to derive their average structure.
Then, the structure with the least RMSD compared to the calculated average structure was selected for each of the two complexes.Figure 10 showcases detailed representations of the AChE active site with ligands 3 and 5 taken from these selected protein-ligand complex structures.
Int. J. Mol.Sci.2023, 24, x FOR PEER REVIEW 11 of 22 Then, the structure with the least RMSD compared to the calculated average structure was selected for each of the two complexes.Figure 10 showcases detailed representations of the AChE active site with ligands 3 and 5 taken from these selected protein-ligand complex structures.Similar to the complexes generated through docking, structures closest to the average ones obtained through MD resemble ligand placement within the peripheral anionic site of the AChE active site.The predominant stabilizing interaction observed is again - stacking.Ligand 3 maintains contact with Trp286 and Tyr341, whereas ligand 5 shifts away from Trp286 towards the His447 of the esteratic site and establishes contact with Trp86.
The analysis of MD simulation results for BChE complexes with compounds 3 and 13 is shown in Figure 11.These complexes exhibit structural stability throughout the simulations, with average RMSD values of 1.97 Å and 1.74 Å , and maximum values of 2.43 Å and 2.17 Å , respectively.During the last 20 nanoseconds of the trajectories, the average RMSD decreases to 1.52 Å and 1.64 Å , with maxima at 1.90 Å and 2.16 Å .Ligand 3 induces more pronounced structural changes in the protein than molecule 13, as evident from the RMSF values: RMS fluctuation reaches 2.82 Å in the complex with 3 and 2.47 Å for compound 13.The compactness of the protein remains stable during the simulations, with gyration radius values similar to those observed in the complexes of AChE.Similar to the complexes generated through docking, structures closest to the average ones obtained through MD resemble ligand placement within the peripheral anionic site of the AChE active site.The predominant stabilizing interaction observed is again π-π stacking.Ligand 3 maintains contact with Trp286 and Tyr341, whereas ligand 5 shifts away from Trp286 towards the His447 of the esteratic site and establishes contact with Trp86.
The analysis of MD simulation results for BChE complexes with compounds 3 and 13 is shown in Figure 11.These complexes exhibit structural stability throughout the simulations, with average RMSD values of 1.97 Å and 1.74 Å, and maximum values of 2.43 Å and 2.17 Å, respectively.During the last 20 nanoseconds of the trajectories, the average RMSD decreases to 1.52 Å and 1.64 Å, with maxima at 1.90 Å and 2.16 Å. Ligand 3 induces more pronounced structural changes in the protein than molecule 13, as evident from the RMSF values: RMS fluctuation reaches 2.82 Å in the complex with 3 and 2.47 Å for compound 13.The compactness of the protein remains stable during the simulations, with gyration radius values similar to those observed in the complexes of AChE.
The active site of BChE with ligands 3 and 13 that are structurally closest to the calculated average structures obtained by MD is presented in Figure 12.
Despite the displacements and structural fluctuations during the simulation, the dominant structure in both BChE-ligand complexes derived by the docking study remained similar to the initial one, thus preserving all relevant stabilizing interactions observed earlier: the π-π stacking between residue Trp82 and ligand, the distance from His438 in the esteratic subsite, and, finally, the contact between the 4-pentenyl, substituent attached to the triazole ring, andTyr332 of PAS.The active site of BChE with ligands 3 and 13 that are structurally closest to the calculated average structures obtained by MD is presented in Figure 12.Despite the displacements and structural fluctuations during the simulation, the dominant structure in both BChE-ligand complexes derived by the docking study remained similar to the initial one, thus preserving all relevant stabilizing interactions observed earlier: the - stacking between residue Trp82 and ligand, the distance from

Genotoxicity Testing on 1,2,3-Triazole Derivatives 1-31
Pharmaceutical development of new drugs requires a thorough investigation into all impurities present in the active pharmaceutical compound (API) and the finished drug product.Impurities that are or can be theoretically present in the API, as well as in each of the intermediates in every step of the manufacturing process of the API, as well as the finished drug, have to be evaluated for their genotoxic potential.The primary evaluations are always done with the in silico models.So, within the scope of impurities of drug substances and drug products, this special subcategory of mutagenic/cancerogenic impurities is extremely important, and all possible risks must be evaluated.These compounds are more strictly regulated and controlled at much lower levels than other impurities.The regulation that is followed for them is the ICH M7 Guideline, and the levels that can be present in the drug substance or the drug product have to be calculated individually for each identified PMI (possible mutagenic compound) based on their determined acceptable daily intake (AI) and the maximum daily dose (MDD) of the final dosage form in question.Suppose toxicological studies on animals have not determined the AI.In that case, the most conservative approach must be taken with the most strict presumed AI described in the guideline itself.
When developing new API and finished drugs, it is expected that the impurities will also be new compounds and that no experimental data will be available for them in the genotoxicity databases.In these cases, the Q(SAR) approach is of vital importance.(Q)SAR models predict biological activity based on structural components [29].This approach used to determine the mutagenic potential of impurities can also be used during the early stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.

Table 3.
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Structure ICH M7 Class Derek Prediction Sarah Prediction
Experimental stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.x FOR PEER REVIEW 14 of 22 stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.stages of searching for potentially active drug substances.The elimination of all compounds that, even if they have biological activity, can also have mutagenic potential can be made fast and easy.The most commonly used tool is the Lhasa Nexus v.2.5.2 software because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.

General Remarks
Nuclear magnetic resonance (NMR) spectroscopic data for 1 H and 13 C nuclei were recorded at room temperature on Bruker Avance 300 and 600 MHz (Coventry, UK) spectrometers.Deuterated chloroform, CDCl3, with TMS as a standard, was used to record NMR spectra.All solvents used are commercially available and purified by distillation.Anhydrous MgSO4 was used to dry the organic layers after extraction.Column chromatography was performed on silica gel columns (60 Å , technical grade).Abbreviations used in this experimental procedure were NMR-nuclear magnetic resonance; EtOAc-ethyl acetate; PE-petroleum ether; E-diethyl ether; EtOH-ethanol; MeOH-methanol; and DCM-dichloromethane. Preparative photochemical reactions were carried out in a closed quartz cuvette in a Rayonet photochemical reactor equipped with 313 nm UV lamps.HRMS analyses were performed with a mass spectrometer (MALDI TOF/TOF analyzer) (Agilent Technologies, Santa Clara, CA, USA) equipped with a Nd:YAG laser operating at 355 nm with a firing rate of 200 Hz in a positive (H + ) or negative (H -) ion reflector.All solvents were removed from the solutions using a rotary evaporator under reduced pressure.

Synthesis of Starting Thieno-Triazole Stilbenes 1a-13a
Starting compounds 1a-13a were obtained as mixtures of cis-and trans-isomers of heterostilbene synthesized by the Wittig reaction.The reaction apparatus was purged with nitrogen for 15 min before adding the reagents.Solutions of 2-thienyl-phosphonium salt (11 mmol) were dissolved in 50 mL of absolute EtOH (dried on a 3 Å sieve) in threenecked round-bottomed flasks (100 mL).Solutions of sodium ethoxide (11 mmol, 1.1 equiv Na dissolved in 10 mL absolute ethanol) were added dropwise under strictly anhydrous conditions under a nitrogen atmosphere.Various triazole aldehydes (11 mmol) were added directly to the mixed solutions.The reaction mixtures were allowed to stir for 24 h at room temperature under a nitrogen bubble.After removing the solvent with a rotary evaporator under reduced pressure, the solid reaction mixtures were extracted with toluene p.a. (3 × 25 mL).The organic layers were dried over anhydrous MgSO4.The final products, a mixture of cis-and trans-isomers 1a-13a, were isolated by silica gel column chromatography using PE/E as eluent and confirmed via 1 H NMR and 13 C NMR spectroscopy and HRMS analyses.

.Synthesis of Cyclization Photoproducts 1-13
Mixtures of previously synthesized compounds 1a-13a were dissolved in toluene p.a. (~ 2.5 × 10 −3 M) and transferred to a quartz cuvette (50 mL) with the addition of a catalytic amount of iodine and illuminated with 10 UV lamps at 313 nm in a Rayonet photochemical reactor (The Southern New England Ultraviolet Co., Branford, USA) for 1-3 h to achieve almost complete conversion.After removing the solvent with a rotary because it uses two complementary models, and their predictions are then reviewed one more time by an expert.
In the case of compounds 1-31 investigated in this paper, several compounds had the Sarah Nexus v.3.2.1 software predictions positive (Table 3).Still, only one structure was found to have a very high risk of being genotoxic.Structures 5 and 6 had very low probability, according to Sarah's prediction, and all of the others (7, 9-11, 15, 18, 21, 22, 24, 27, 29, and 30) had the model rely on the data for the thiophene ring, with also a relatively low percentage of certainty.
As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.As for structure 31, here the probability is very high, and this structure also has a Derek (Derek Nexus v.6.2.1 software); positive prediction, with a much higher certainty.Here, the Q(SAR) found a very similar compound (condensed aromatic system with the amide functionality linked directly) with a positive AMES test for mutagenicity in vitro.This can be taken as a compound that can be eliminated from further research because the probability of its genotoxic potential is high.

General Remarks
Nuclear magnetic resonance (NMR) spectroscopic data for 1 H and 13 C nuclei were recorded at room temperature on Bruker Avance 300 and 600 MHz (Coventry, UK) spectrometers.Deuterated chloroform, CDCl 3 , with TMS as a standard, was used to record NMR spectra.All solvents used are commercially available and purified by distillation.Anhydrous MgSO 4 was used to dry the organic layers after extraction.Column chromatography was performed on silica gel columns (60 Å, technical grade).Abbreviations used in this experimental procedure were NMR-nuclear magnetic resonance; EtOAc-ethyl acetate; PE-petroleum ether; E-diethyl ether; EtOH-ethanol; MeOH-methanol; and DCM-dichloromethane. Preparative photochemical reactions were carried out in a closed quartz cuvette in a Rayonet photochemical reactor equipped with 313 nm UV lamps.HRMS analyses were performed with a mass spectrometer (MALDI TOF/TOF analyzer) (Agilent Technologies, Santa Clara, CA, USA) equipped with a Nd:YAG laser operating at 355 nm with a firing rate of 200 Hz in a positive (H + ) or negative (H -) ion reflector.All solvents were removed from the solutions using a rotary evaporator under reduced pressure.

Synthesis of Starting Thieno-Triazole Stilbenes 1a-13a
Starting compounds 1a-13a were obtained as mixtures of cisand trans-isomers of heterostilbene synthesized by the Wittig reaction.The reaction apparatus was purged with nitrogen for 15 min before adding the reagents.Solutions of 2-thienyl-phosphonium salt (11 mmol) were dissolved in 50 mL of absolute EtOH (dried on a 3 Å sieve) in threenecked round-bottomed flasks (100 mL).Solutions of sodium ethoxide (11 mmol, 1.1 equiv Na dissolved in 10 mL absolute ethanol) were added dropwise under strictly anhydrous conditions under a nitrogen atmosphere.Various triazole aldehydes (11 mmol) were added directly to the mixed solutions.The reaction mixtures were allowed to stir for 24 h at room temperature under a nitrogen bubble.After removing the solvent with a rotary evaporator under reduced pressure, the solid reaction mixtures were extracted with toluene p.a. (3 × 25 mL).The organic layers were dried over anhydrous MgSO 4 .The final products, a mixture of cisand trans-isomers 1a-13a, were isolated by silica gel column chromatography using PE/E as eluent and confirmed via 1 H NMR and 13 C NMR spectroscopy and HRMS analyses.

In Vitro Biological Activity
PBMC isolation.Human peripheral blood mononuclear cells (PBMCs) were isolated from buffy coats obtained from healthy adult volunteers.The buffy coats were diluted with sterile PBS and layered over Lymphoprep (Axis-Shield Diagnostics, Dundee, UK).Tubes were centrifuged for 35 min at 400×g (RT; acceleration and break turned off).After centrifugation, the mononuclear ring was collected, transferred into a new 50 mL tube, and washed thrice with sterile PBS.The remaining erythrocytes were lysed with an isotonic solution of ammonium chloride (150 mM; Kemika, Zagreb, Croatia).The cell pellet was resuspended in Roswell Park Memorial Institute (RPMI) 1640 medium (Lonza, Basel, Switzerland) supplemented with 10% FBS (Biowest, Nuaillé France).
Treatment and stimulation of human PBMCs.Immediately after isolation, 200,000 PBMCs were seeded per well of a 96-well plate.Compounds were dissolved in 100% dimethyl sulfoxide (DMSO, Sigma, Darmstadt, Germany), and serial dilutions in DMSO were prepared and added to cells, with starting concentration for test compounds 100 µ M or 1 µ M for reference dexamethasone.After pre-incubation for 1h with the compounds, lipopolysaccharide (LPS) from E. coli 0111:B4 (Sigma, Darmstadt, Germany) was added to cells at 1 ng/mL final concentration.Plates were incubated for 24 h at 37°C, 5% CO2, followed by the collection of supernatants for measurement of TNF-α and cell viability assessment.
Enzyme-linked immunosorbent assay (ELISA) and cytotoxicity analysis.TNFα concentration in PBMC supernatants was measured using ELISA assay, using antibodies obtained from BD Pharmingen (San Diego, CA, USA) and recombinant human TNFα protein (standard) from R&D Systems (Minneapolis, MN, USA).Lumitrac 600 384-well plates

In Vitro Biological Activity
PBMC isolation.Human peripheral blood mononuclear cells (PBMCs) were isolated from buffy coats obtained from healthy adult volunteers.The buffy coats were diluted with sterile PBS and layered over Lymphoprep (Axis-Shield Diagnostics, Dundee, UK).Tubes were centrifuged for 35 min at 400× g (RT; acceleration and break turned off).After centrifugation, the mononuclear ring was collected, transferred into a new 50 mL tube, and washed thrice with sterile PBS.The remaining erythrocytes were lysed with an isotonic solution of ammonium chloride (150 mM; Kemika, Zagreb, Croatia).The cell pellet was resuspended in Roswell Park Memorial Institute (RPMI) 1640 medium (Lonza, Basel, Switzerland) supplemented with 10% FBS (Biowest, Nuaillé France).
Treatment and stimulation of human PBMCs.Immediately after isolation, 200,000 PBMCs were seeded per well of a 96-well plate.Compounds were dissolved in 100% dimethyl sulfoxide (DMSO, Sigma, Darmstadt, Germany), and serial dilutions in DMSO were prepared and added to cells, with starting concentration for test compounds 100 µM or 1 µM for reference dexamethasone.After pre-incubation for 1h with the compounds, lipopolysaccharide (LPS) from E. coli 0111:B4 (Sigma, Darmstadt, Germany) was added to cells at 1 ng/mL final concentration.Plates were incubated for 24 h at 37 • C, 5% CO 2 , followed by the collection of supernatants for measurement of TNF-α and cell viability assessment.
Enzyme-linked immunosorbent assay (ELISA) and cytotoxicity analysis.TNFα concentration in PBMC supernatants was measured using ELISA assay, using antibodies obtained from BD Pharmingen (San Diego, CA, USA) and recombinant human TNFα protein (standard) from R&D Systems (Minneapolis, MN, USA).Lumitrac 600 384-well plates (Greiner Bio-One, Kremsmünster, Austria) were coated overnight with 1 µg/mL of TNFα capture antibody diluted in phosphate-buffered saline (PBS, Thermo Fisher Scientific, Waltham, MA, USA), followed by a washing step and blocking with 5% sucrose in assay diluent (1% bovine serum albumin (BSA; Sigma, Darmstadt, Germany) in PBS) for 4 h at RT.After the blocking step, the plates were washed, and samples/standards were added to the plates.After overnight incubation at 4 • C, detection antibodies diluted in assay diluent to a final concentration of 250 ng/mL were added.Following incubation for 2 h at RT, plates were washed, and streptavidin-HRP (Thermo Fisher Scientific, Waltham, MA, USA) solution was added (final concentration 500 ng/mL, prepared in assay diluent).Plates were incubated for 20 min at RT in the dark and washed, followed by adding Luminol substrate (Sigma, Darmstadt, Germany).Luminescence was measured using EnVision 2104 Multilabel Plate Reader (PerkinElmer, MA, USA), with an exposition time of 0.1 s.To evaluate cell viability, Cell Titer Glo reagent was used (Promega, Madison, WI, USA) according to the manufacturer's instructions.
Data analysis.For ELISA, the blank was subtracted from measured RLU values, and TNFα concentrations in samples were interpolated from the standard curve.QC parameters for each plate were calculated from raw data.Percentages of inhibition of each compound were calculated from obtained cytokine concentrations using the formula: PIN = 100 − (((compound − no trigger)/(trigger − no trigger)) * 100).For viability assessment, the average luminescence value was calculated from all trigger vehicle samples, and the percentage of vehicle value was calculated for all samples.The compound effect was considered cytotoxic if reduction from a vehicle was ≥20%, and these concentrations were excluded from further analysis.IC 50 values of tested compounds were determined in GraphPad Prism 9 software using nonlinear regression curve fit (four parameters, with variable slope).

In Vitro ChE Activity Assay
The inhibition of acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) was assessed using a modified spectrophotometric method based on Ellman's approach [30].AChE (EC from electric eel) and BChE (EC from equine serum), along with Trizma base, acetylthiocholine iodide (ATChI), S-butyrylthiocholine iodide (BTChI), and Galantamine were procured from Sigma-Aldrich (Darmstadt, Germany).Additionally, 5,50-dithiobis-(2nitrobenzoic acid) (DTNB), known as Ellman's reagent, was obtained from Zwijndrecht (Antwerpen, Belgium).Galantamine served as a reference standard in the experiment.The assessment of AChE/BChE activity was performed using a 96-well microplate reader (IRE 96, SFRI Medical Diagnostics, Saint Jean d'Illac, France) at 405 nm.The measurement lasted for 6 min at room temperature after the initiation of the enzymatic reaction.In each well of the microplate, the following components were present: 180 µL of Tris-HCl buffer (50 mM, pH 8.0), 10 µL of the respective enzyme (final concentration 0.03 U/mL), 10 µL of the tested solution with varying concentrations (final concentrations ranging from 10 to 350 µM, depending on solubility), and 10 µL of DTNB (final concentration 0.3 mM).The enzymatic reaction was initiated by adding 10 µL of ATChI/BTChI (final concentration of 0.5 mM prepared in Tris buffer) to each well.For the control measurement, a buffer was used in place of the tested compound.Additionally, for non-enzymatic hydrolysis, each measurement included a blank with the enzyme substituted by an equivalent buffer amount.To ensure accuracy, the experiment was conducted in triplicate.The solvent used in the experiment was ethanol.The percentage of enzyme inhibition was calculated using the measured data based on the equation: inhibition (%) = [(Ac − AT)/Ac] × 100, where Ac represents the enzyme activity without the test sample, and AT represents the enzyme activity with the test sample, calculated as mean values ± standard deviation.The IC 50

2. 4 .
Docking and Molecular Dynamics Study of 1,2,3-Triazole Derivatives To elucidate the interactions underpinning the stability of complexes formed between the triazole derivatives and cholinesterases, we conducted molecular docking on the compounds that exhibited the most potent inhibitory activity in experimental assays.Notably, for AChE, the most promising outcomes were observed with compounds 3 and 5.The structures of complexes obtained by docking 3 and 5 into the active site of AChE are shown in Figure 7a,b, respectively.Int.J. Mol.Sci.2023, 24, x FOR PEER REVIEW 8 of 22

Figure 7 .
Figure 7. (a) Structure of the complex between compound 3 and the active site of AChE, obtained by molecular docking; (b) compound 5 docked into the active site of AChE.Distances are given in angstroms.Hydrogen atoms of the residues were omitted for clarity.

Figure 7 .
Figure 7. (a) Structure of the complex between compound 3 and the active site of AChE, obtained by molecular docking; (b) compound 5 docked into the active site of AChE.Distances are given in angstroms.Hydrogen atoms of the residues were omitted for clarity.

Figure 8 .
Figure 8.(a) Structure of the complex between compound 3 and the active site of BChE, obtained by molecular docking; (b) compound 13 docked into the active site of BChE.All distances in angstroms.Docked molecules are presented using a ball-and-stick model, with hydrogen atoms of residues omitted for clarity.

Figure 8 .
Figure 8.(a) Structure of the complex between compound 3 and the active site of BChE, obtained by molecular docking; (b) compound 13 docked into the active site of BChE.All distances in angstroms.Docked molecules are presented using a ball-and-stick model, with hydrogen atoms of residues omitted for clarity.

Figure 9 .
Figure 9. Molecular dynamics simulations of a protein-ligand complex of AChE with compounds 3 and 5: (a) RMSD, (b) RMSF, and (c) radius of gyration plots.

Figure 9 .
Figure 9. Molecular dynamics simulations of a protein-ligand complex of AChE with compounds 3 and 5: (a) RMSD, (b) RMSF, and (c) radius of gyration plots.

Figure 10 .
Figure 10.Protein-ligand complexes derived from MD simulation as structures with the least RMSD from the calculated average structures.(a) The active site of AChE with ligand 3; (b) AChE active site with ligand 5. Ligands are presented using a ball-and-stick model, with hydrogen atoms of residues omitted for clarity.

Figure 10 .
Figure 10.Protein-ligand complexes derived from MD simulation as structures with the least RMSD from the calculated average structures.(a) The active site of AChE with ligand 3; (b) AChE active site with ligand 5. Ligands are presented using a ball-and-stick model, with hydrogen atoms of residues omitted for clarity.

Figure 12 .
Figure 12.Protein-ligand complexes derived from MD simulation as structures with the least RMSD from the calculated average structures.(a) The active site of BChE with ligand 3; (b) BChE active site with ligand 13.Ligands are presented using a ball-and-stick model, with hydrogen atoms of residues omitted for clarity.

Table 1 .
In vitro activity of novel 1,2,3-triazoles on inhibition of TNF-α production in PBMC stimulated with LPS.The IC50 values shown are average values from testing in 2 independent experiments (donors).

Table 1 .
In vitro activity of novel 1,2,3-triazoles on inhibition of TNF-α production in PBMC stimulated with LPS.The IC 50 values shown are average values from testing in 2 independent experiments (donors).

Table 2 .
Inhibition of AChE and BChE and calculated IC 50 values by compounds 1

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.

Table 3 .
Compounds with a positive mutagenic potential by Lhasa M7 evaluation.