Next Article in Journal
Human Estrogen Receptor Alpha Antagonists, Part 3: 3-D Pharmacophore and 3-D QSAR Guided Brefeldin A Hit-to-Lead Optimization toward New Breast Cancer Suppressants
Next Article in Special Issue
Pd(II) Binding Strength of a Novel Ambidentate Dipeptide-Hydroxypyridinonate Ligand: A Solution Equilibrium Study
Previous Article in Journal
Research Update of Emergent Sulfur Quantum Dots in Synthesis and Sensing/Bioimaging Applications
Previous Article in Special Issue
Hydroxypyridinone-Based Metal Chelators towards Ecotoxicity: Remediation and Biological Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson’s Disease

by
Joseph C. J. Bergin
1,
Kean Kan Tan
1,
Anya K. Nelson
1,
Cristina-Andreea Amarandei
2,
Véronique Hubscher-Bruder
2,
Jérémy Brandel
2,
Varvara Voinarovska
3,
Annick Dejaegere
3,
Roland H. Stote
3 and
David Tétard
1,*
1
Department of Applied Sciences, Faculty of Health and Life Sciences, Northumbria University, Newcastle upon Tyne NE1 8ST, UK
2
Université de Strasbourg, CNRS, IPHC UMR 7178, F-67000 Strasbourg, France
3
Institut de Génétique et de Biologie Moléculaire et Cellulaire (IGBMC), Institut National de La Santé et de La Recherche Médicale (INSERM), U1258/Centre National de Recherche Scientifique (CNRS), UMR7104/Université de Strasbourg, 67404 Illkirch, France
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(9), 2816; https://doi.org/10.3390/molecules27092816
Submission received: 31 January 2022 / Revised: 14 March 2022 / Accepted: 18 March 2022 / Published: 28 April 2022
(This article belongs to the Special Issue Novel Functional Hydroxypyridinone-Based Derivatives)

Abstract

:
Two analogues of tolcapone where the nitrocatechol group has been replaced by a 1-hydroxy-2(1H)-pyridinone have been designed and synthesised. These compounds are expected to have a dual mode of action both beneficial against Parkinson’s disease: they are designed to be inhibitors of catechol O-methyl transferase, which contribute to the reduction of dopamine in the brain, and to protect neurons against oxidative damage. To assess whether these compounds are worthy of biological assessment to demonstrate these effects, measurement of their pKa and stability constants for Fe(III), in silico modelling of their potential to inhibit COMT and blood–brain barrier scoring were performed. These results demonstrate that the compounds may indeed have the desired properties, indicating they are indeed promising candidates for further evaluation.

1. Introduction

Parkinson’s disease (PD) is a devastating disease characterised by the death of dopaminergic neurons in the substantia nigra pars compacta. In this very complex disease, ferroptosis, accumulation of iron and oxidative damage are thought to be critically involved in the pathology [1,2,3]. PD has no cure, but symptoms can be managed by the administration of L-DOPA as a precursor to dopamine. L-DOPA is capable of crossing the blood–brain barrier (BBB) into the central nervous system (CNS) where it is metabolised by aromatic aminoacid decarboxylase (AADC) enzymes into dopamine. However, L-DOPA and dopamine are further metabolised by catechol O-methyl transferase (COMT, E.C. 2.1.1.6), thus reducing the amount of dopamine available as a neurotransmitter [4]. To limit these unwanted metabolic degradations, inhibitors of AADC and COMT are usually co-administered with L-DOPA. In the case of COMT, inhibitors such as tolcapone, opicapone and entacapone (Figure 1) were developed and used in the clinic [4,5]. Tolcapone is rather unique in this set as it is the only COMT inhibitor that is known to significantly cross the BBB and act both peripherally and in the CNS [6,7]. These inhibitors are based on a nitrocatechol coordinating group and act by chelating the Mg2+ co-factor present in COMT [8,9]. One effect of the nitro group is to lower the pKa of the hydroxyl groups (around 4.6–4.7 and 10.0–10.2 in entacapone, nitecapone and tolcapone [10]), helping to bind to Mg2+. However, this nitro group has been linked to hepatic toxicity in tolcapone and its use has been discontinued [11].
We have recently described the use of metal chelators based on the 1-hydroxy-2(1H)-pyridinone (1,2-HOPO) coordinating group as neuron protective agents against toxins mimicking the insult observed in PD on dopaminergic neurons [16,17]. We proposed that their ability to bind excess Fe(III) and prevent oxidative damage to neurons contribute to their efficacy. Hydroxypyridinones are known to be good coordinating groups for hard Lewis acids such as Mg2+ and require only one deprotonation to be able to chelate a metal cation. These are expected to make 1,2-HOPO good alternatives to nitrocatechols in COMT inhibitors, without the associated hepatic toxicity.
We therefore hypothesised that simple derivatives of 1,2-HOPO could serve a dual purpose of COMT inhibition to limit dopamine degradation and neuron rescue against oxidative damage. Both of these are expected to be beneficial for treatment of PD. For neuron rescue, the compound must act in the CNS, requiring the compound to be able to cross the BBB. For this reason, we chose tolcapone as the template for modification (Figure 2). As the charge state of the compound is a key factor for BBB permeation, the pKa values of the compounds were determined and the BBB score, representing a statistical chance to cross the BBB, was calculated. Moreover, Fe(III) chelation ability is thought to be determining for the efficacy of the compounds as neuro-protective agents, as they could limit radical oxidation species (ROS) production via inhibition of the Fenton reaction. As such, their affinity for this metal and stoichiometry of the species formed were also determined. Finally, an in silico study was conducted to assess their potential as COMT inhibitors. The results reported herein show that they indeed have the potential to inhibit Fe(III)-catalysed oxidative damage, can inhibit COMT and are likely to penetrate the CNS. Further study of these compounds, including biological screening is therefore warranted.

2. Results and Discussion

The synthesis of compound 1a is depicted in Scheme 1. The synthetic strategy relied primarily on the formation of the 1,2-HOPO group in two steps from the 2-chloropyridine derivative, first forming the 2-chloropyridine N-oxide (4a, step ii). Subsequent substitution of the chlorine atom for a hydroxyl group (step iv) resulted in the formation of 2-hydroxypyridine N-oxide (6a) that is a tautomer of 1a.
The formation of the ketone was planned in two steps. First was the formation of an alcohol using a Grignard reaction on commercially available aldehydes (step i) and a subsequent oxidation into the carbonyl group (step iii). We speculated that the oxidation needed to be performed after the formation of the N-oxide to avoid the possible formation of the ester via a Baeyer–Villiger oxidation when using a peroxyacid. It was deduced that the oxidation of the alcohol to the ketone was to be performed either after the 1,2-HOPO moiety had been fully formed or between the two steps required to synthesise that group. In the former option, this may have required the protection then deprotection of the 1,2-HOPO’s hydroxyl group to avoid possible interference with metal-based reagents used in the oxidation (e.g., MnO2) or with organic oxidation reagents (e.g., oxalyl chloride in the Swern oxidation). To limit the need for protection and deprotection, the oxidation of the alcohol was therefore planned between the formation of the 2-chloropyridine N-oxide (step ii) and its reaction with hydroxide (step iv).
The first step of the synthesis was therefore the reaction of a Grignard reagent from 4-bromotoluene on the aldehyde to form the alcohols 3a in 93% yield. The formation of the 2-chloropyridine N-oxide 4a using mCPBA resulted in the formation of a small amount of 5a. We did not attempt to separate these compounds. Instead, the impure alcohol was oxidised into the diaryl ketone 5a using manganese dioxide. Target compound 1a was then obtained by reacting potassium hydroxide on the 2-chloropyridinine N-oxide.
The synthesis of compound 1b followed a similar strategy with one exception (Scheme 2). To our surprise, the formation of the pyridine N-oxide using mCPBA was also accompanied by the unexpected full oxidation of the alcohol into the wanted carbonyl group, thus giving 5b directly and saving one step in our synthetic route. The unambiguous oxidation of the alcohol into the carbonyl was demonstrated by the appearance of a 13C NMR signal at 190.9 ppm, the expected MS fragmentations observed at m/z = 119.04927 (C8H7O, calc. 119.04969) and 91.05464 (C7H7, calc. 91.05478) and no evidence of any CH group in the 1H NMR around 4.5 ppm. We speculate that this rather mild and efficient oxidation might have been catalysed by traces of contaminant, possibly transition metals in the glassware (although this has been a systematic observation made with different glassware and different people) and proceeded via a radical mechanism [18].
The formation of compounds 5a and 5b in our synthetic strategy allowed us to use them as control compounds in our biological evaluations as they are expected to have similar electron distribution as 1a and 1b, but without the ability to bind Mg2+ or Fe3+. An additional control compound (7bScheme 3), which is not expected to have any biological function [19], was prepared by the methylation of 1b using methyl iodide in excellent purity without any need for specific purification procedures [20].
Physicochemical measurements were carried out on ligands 1a and 1b to evaluate their acido-basic behaviour and Fe(III) complexation properties in solution with respect to pH. The pKa of these compounds will have a strong influence on their ability to penetrate the brain [21,22]. It will also, together with the Fe(III) stability constant, influence their ability to redox-inactivate Fe(III), thought to be key for neuron rescue [23]. Acid-base properties of both ligands were determined by potentiometric titrations and UV-visible spectrophotometric titrations vs. pH (Figure 3). For solubility reasons, the studies were carried out in a mixed MeOH/H2O (80/20 w/w) solvent.
The titrations of both ligands showed the presence of one protonation equilibrium in the 2 to 12 pH range for each ligand. Analysis of the spectral variations [24,25] allowed us to calculate the protonation constants of the formed species (Table 1).
From these protonation constants, the electronic spectra (1b on Figure 3C and 1a on Figure S1A) and distribution curves (1b on Figure 4C and 1a on Figure S1B) of the species vs. pH were calculated. The distribution curves suggested that, at physiological pH (pH 7.4), both ligands are predominantly negatively charged (1a 76%, 1b 100%) with the N-hydroxyl group of the 1-hydroxypyrydin-2(1H)-one ring being deprotonated.
As the ability to bind Fe(III) is a key factor for the evaluation of the potential use of these compounds against PD and oxidative damage, Fe(III) complexation studies vs. pH were carried out in the same way on compounds 1a and 1b, via spectrophotometric titrations between pH 2 and 12 (Figure 4).
The spectra of Fe(III) complexation titrations showed the appearance of a ligand to metal charge transfer (LMCT) band between 400 and 600 nm for both ligands, a sign that complexation started at very low pH. The LMCT band underwent a hyperchromic and hypsochromic shift with increasing pH suggesting an increase in the number of ligands in the coordination sphere of Fe(III), in accordance with the bidentate structure of the ligands. From a pH around 6, the LMCT band started decreasing and an increase in the baseline was observed, indicating the progressive decomposition of the complexes and precipitation of iron hydroxides (Fe(OH)3). Analysis of the spectral variations prior to precipitation suggested the successive formation of FeL2 and FeL3 complexes (L = 1a or 1b) in the pH range 2 to 12, in accordance with the bidentate character of the ligands. The calculated stability constants (log β) of the species are reported in Table 1. From the values of these stability constants, the electronic spectra of the complex species (Figure 4C and Figure S2A) and their distribution curves (Figure 4D and Figure S2B) were calculated. As the complexation of a metal by a ligand is dependent on several parameters such as metal ion concentration, ionic strength, ionic medium, pH, temperature and protonation constants of the ligands, the stability constants of the complexes cannot be used as such to compare the sequestering ability of a series of ligands for a given metal. The sequestering power of a ligand can be evaluated by the determination of the empirical and quantitative parameter pL0.5, which represents the total concentration of ligand required for the sequestration of 50% of the metal. [26]. A higher pL0.5 value will indicate a higher sequestering ability.
The pL0.5 values are presented in Table 1 (sequestering diagrams in Figure S3) and suggest that ligand 1b is a substantially better Fe(III) chelator than 1a. As a comparison, the 3,4-HOPO ligand DFP shows a pL0.5 value of 6.81. Previously studied 1-hydroxypyrazin-2(1H)-ones (1,2-HOPY) ligands 8 and 9 (structure in Figure S3), showed pL0.5 values of 7.34 and 5.98, respectively, and encouraging neuroprotection against the PD-relevant neurotoxin, 6-OHDA [27], which suggests that compounds 1a and 1b have a suitable Fe(III) chelation power and might help limit the formation of ROS in neurons.
Numerous methods have been described to predict the BBB permeation (via passive diffusion) of an organic compound [28,29,30]. Unfortunately, the molecular descriptors identified in these different methods as being key to BBB permeation are not always the same and some discrepancies exist, for example, regarding the influence of the partition coefficient log P usually taken as the n-octanol/water partition coefficient. Until a consensus is reached, we therefore have to assess the compounds using several tools to obtain a qualitative idea of their ability. The first tool we used was SwissADME [31]. Calculated parameters relevant to BBB Score of tolcapone, 1a and 1b obtained by SwissADME are given in Table 2.
Tolcapone is known to be able to cross the BBB [7] despite the negative prediction of SwissADME. Interestingly, our compounds are predicted to be BBB permeants by the same tool. In particular, 1a and 1b have far fewer numbers of hydrogen bond donors and acceptors. As the ability to cross the BBB has been proposed to be hampered by hydrogen bond formation [28], this suggests that our compounds are more likely to be BBB permeants than tolcapone, as required for neuron rescue.
Another recent method has defined the BBB Score based on the following parameters: number of aromatic rings (Aro_R), number of heavy atoms (HA), molar mass, number of hydrogen bond donors and acceptors, topological surface area (TPSA) and pKa [21]. The BBB scores of tolcapone, 1a and 1b are given in Table 3.
Tolcapone has a BBB score of 3.25 which gave it a 21.9% chance of being BBB permeant. Our compounds have a BBB Score of 4.8–5.0 that puts them in the higher category of compounds with a statistical chance of being BBB permeant of 54.5% and close to the category of compounds having a statistical chance of being CNS permeant of over 90%.
Another recent method to predict BBB permeation identified log P (partition coefficient of the neutral species, typically between n-octanol and water), distribution coefficient log D[pH 7.4] (typically the n-octanol/water distribution of the ionised species at pH 7.4), TPSA, neutrality/acidity/basicity, hydrogen bond donor (HBD), hydrogen bond acceptor (HBA), number of nitrogen atoms (NC), number of oxygen atoms (OC) and total number of nitrogen and oxygen atoms (NOC) as key parameters. [22]. The calculations for tolcapone, 1a and 1b are given in Table 4.
Again, this analysis indicates that our compounds are better than tolcapone with regard to HBA and HBD numbers. TPSA, OC and NOC are also much better although still classified in the red regions.
The pKas of tolcapone are 4.64 and 10.20. The pKa of our compounds are 6.9 and 5.69 for 1a and 1b, respectively. It has to be noted that the pKa values of compounds 1a and 1b were determined in a MeOH/H2O (80/20 w/w) solvent, which results in slightly higher values than in water. Thus, assuming that the pKa in water is not too much lower for our compounds compared to the ones measured in a mixed solvent, the percentage of non-ionised species present at physiological pH is far greater for compounds 1a (24%) and 1b (1.9%) than for tolcapone (0.2%), favouring brain penetration as required for neuron rescue, especially for 1a. These properties suggest that our compounds could reach the CNS in a higher concentration than tolcapone.
Molecular modelling studies were carried out to investigate the nature of compounds 1a and 1b binding to COMT. More specifically, focus was on whether there are any major structural or energetic factors that would suggest non-binding of the ligands. Following the protocols described in the Methods (Supplementary Data), molecular dynamics (MD) simulations of the known inhibitor tolcapone in complex with the COMT enzyme for rat (crystal structure available) and human COMT (model constructed by comparative docking) were carried out. The structural stability and free energy binding scores were compared to those obtained by simulations of both rat and human COMT complexed to molecules 1a and 1b (models constructed by comparative docking).
Molecular dynamics (MD) simulations of tolcapone were carried out in two different protonation states and in a single orientation (based on the crystallographic data, see Figure 5), while simulations of compounds 1a and 1b were performed in a single protonation state but in two different orientations with respect to the SAM cofactor (see Figure 6). Representations of the structures of the different tolcapone and 1a and 1b complexes are given in Figures S4–S6, respectively.
The root-mean-square-difference (RMSD) time series of the backbone atoms calculated from the trajectories were compared for the rat and human complexes (Figure 7). The results show that the global conformations do not deviate significantly from the initial reference structures, which demonstrates the general stability of the different complexes. Even though the RMSD time series for 1b shows slightly higher values that the other compounds in rat COMT, the results suggest that compounds 1a and 1b do not introduce any strong interactions that would perturb the protein structure. To further confirm this, average structures were calculated over the final 10 ns of the simulations and the backbone RMSDs were calculated by-residue for both rat and human COMT, Figures S7 and S9, and the average structures were superposed, Figures S8 and S10, respectively. These results show that all helices maintain their secondary structure (Figures S8A and S10A) and the main structural changes are limited to the loops in the vicinity of the ligand binding site (Figures S8B and S10B). As suggested by the RMSD times series shown in Figure 7, the comparison of the average structures further supports the conclusion that compounds 1a and 1b do not introduce any significant structural changes. A closer examination of the ligand conformation in the binding pocket shows that in some of the simulations, the ligand moved from a bidentate coordination starting configuration to a more monodentate coordination of the Mg2+ ion. An earlier study also observed a movement from bidentate to monodentate coordination [33]. In this work, we do not attempt to resolve this question of mono- vs. bidentate coordination, rather we focus on comparing relative behaviours between the complexes of the different ligands.
The root-mean-square-fluctuations (RMSF) for the crystallographic complex of rat COMT with tolcapone were calculated from the MD simulations and compared those calculated from the experimental temperature factors (B factors) determined for the crystal structure. In Figure 8, a comparison of the backbone RMSFs averaged by residue shows good agreement in the overall pattern of flexibility. Only residues 170–190 (220–240 in human COMT) in the MD simulations show motions notably greater than those determined from crystallographic B factors. Using the symexp option in the Pymol program (ref: The PyMOL Molecular Graphics System, Version 2.5 Schrödinger, LL) to visualise the position of the enzyme in the crystal lattice, we conclude that the lower values of fluctuations observed in the experimental structure are due to crystal lattice contacts. Visualisation confirmed that, in the crystal structure, these residues form an alpha helix that is in contact with other proteins in the lattice. In the simulations, this region is solvent-exposed and displays a wider range of motions.
The RMSF pattern for the docked complexes with compounds 1a and 1b is overall similar to the complex with tolcapone.
From each simulation, we calculated an average free energy of binding score by the MM/GBSA method [34] using the CHARMM program combined with in-house programs. The MM/GBSA method is a rapid calculation of the free energy of binding based on the molecular mechanics force field and an implicit representation of the solvent environment. The calculation was repeated for 1000 structures extracted from each simulation. Here, we do not consider the absolute values of the numbers, but their values relative to that of the doubly deprotonated tolcapone (rat and human). The results are presented in Table 5. The absolute values from the MM/GBSA calculations are given in Supplementary Materials.
The relative binding free energy scores with respect to doubly deprotonated tolcapone in rat and human COMT were calculated. Between singly and doubly deprotonated tolcapone, the doubly deprotonated form is the stronger binder because of favourable electrostatic interactions with the Mg2+ ion. In the simulations of rat COMT, molecules 1a and 1b showed binding free energies in the same range as singly or doubly deprotonated tolcapone, depending on the orientation. In the simulations of human COMT, the ligand showed more varied free energy scores, but one orientation of 1a (1a N-O→SAM) had a binding free energy score more favourable than the doubly deprotonated tolcapone.
Interestingly, one orientation of compound 1a (1a = O→SAM) gave a positive free energy score in both rat and human COMT (Table S1), which suggests that there may be some orientational preference in the binding of this compound to COMT. The large error bars reflect the strong electrostatic contribution to the binding free energy, which leads to large fluctuations in the binding free energy.
The results provide the first investigation into the potential binding capabilities of a new class of COMT inhibitors, compounds 1a and 1b. Besides the structural stability of complexes of COMT with these ligands, as compared to simulations of a complex for which the experimental structure is known, the energetic rankings of these ligands were in the same range as the known COMT inhibitor, tolcapone. The simulation results suggest the possibility of a preferred orientation in the ligand binding pocket for at least one of the compounds (1a). While further investigations will be necessary, these first results do not provide any evidence against the potential binding capabilities of these compounds to the COMT enzyme.

3. Experimental

All reagents were purchased from Sigma-Aldrich (St. Louis, MO, USA) and Fluorochem (Glossop, UK) and used without further purification unless stated. Solvents were obtained from Fisher Scientific (Waltham, MA 02451 USA) and were of either reagent or HPLC Grade and column chromatography was performed using Fisher Scientific silica 60Å (35–70 μm). THF was dried over 3Å molecular sieves for 48 h. Reactions were monitored with thin layer chromatography (TLC) performed on Merck plastic foil plates pre-coated with silica gel 60 F254. Infrared-red (IR) spectra were recorded using a SensIR Technologies Durascope diamond anvil cell mounted on a Perkin-Elmer Paragon 1000 FTIR Spectrometer (Texas, TX, USA), wavenumbers are reported in cm−1 and band intensities are described using the following abbreviations: w = weak, m = medium, s = strong. NMR spectra were recorded using a JEOL ECS400 Delta spectrometer (Tokyo, Japan) at frequencies of 399.78 MHz for 1H NMR and 100.53 MHz for 13C NMR. All chemical shifts are reported in parts per million (ppm) relative to a tetramethylsilane (TMS) internal standard. Chemical shifts are reported as follows: Shift (δ) in ppm (multiplicity, coupling constant (J) in Hz, normalised integral). Multiplets are described as follows: s = singlet, d = doublet, t = triplet and dd = doublet of doublets. High-resolution mass spectrometry was performed using a Q Exactive Mass Spectrometer (New York, NY, USA) under electrospray ionisation (ESI) operated in positive and negative ion modes as specified (scan range 85–800 m/z). All observed peak values are recorded as mass–charge ratio (m/z).

3.1. 2-Chloro-5-([4-methylphenyl]hydroxymethyl)-pyridine (3a)

Under an argon atmosphere, magnesium (0.99 g, 41 mmol) was added to a rapidly stirring solution of 4-bromotoluene (7.57 g, 44.3 mmol) in dry tetrahydrofuran (95 mL). The mixture was then stirred at reflux for three hours until all of the Mg had been reacted. The Grignard mixture was then cooled to 0 °C and a solution of 2-chloropyridine-5-carboxaldehyde (5.22 g, 36.9 mmol) in dry THF (40 mL) was then added dropwise. The reaction was stirred for a further thirty minutes at 0 °C and then quenched with water. The product was then extracted with ethyl acetate and the organic layer was then washed with brine, dried with magnesium sulphate filtered and evaporated under vacuum. The solid was triturated with petroleum ether and evaporated to yield 8.05 g of yellowish crystals (34.4 mmol, 93%). 1H NMR (400 MHz, CDCl3) δ 8.37 (d, J = 2.3 Hz, 1H), 7.65 (dd, J = 8.2, 2.75 Hz, 1H), 7.27 (d, J = 8.2 Hz, 1H), 7.22 (d, J = 7.8 Hz, 2H), 7.17 (d, J = 8.2 Hz, 2H), 5.82 (d, J = 2.8 Hz, 1H), 2.56 (d, J = 3.2 Hz, 1H), 2.34 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 150.3 (quat.), 148.0 (CH), 139.9 (quat.), 138.5 (quat.), 138.2 (quat.), 137.2 (CH), 129.7 (CH), 126.6 (CH), 124.1 (CH), 73.4 (CH), 21.2 (CH3).
Molecules 27 02816 i001

3.2. 2-Chloro-5-([4-methylphenyl]hydroxymethyl)-pyridine N-oxide (4a)

To a solution of 3a (3.09 g, 12.4 mmol) in dichloromethane (50 mL), 10.32 g of meta-chloroperoxybenzoic acid (77 wt.%, 46 mmol) was added and stirred at reflux overnight. The organic phase was then cooled to 0 °C and then washed with a solution of sodium thiosulphate (10%, 3 × 50 mL) followed by a saturated solution of sodium bicarbonate and brine. The organic layer was then dried over magnesium sulphate and evaporated under reduced pressure to yield 0.75 g of a solid contaminated with 5a which was used without further purification. 1H NMR (400 MHz, CDCl3) δ 8.35 (d, J = 1.3 Hz, 1H), 7.31 (d, J = 8.2 Hz, 1H), 7.12–7.21 (m, 5H), 5.70 (s, 1H), 4.70 (br, 1H), 2.32 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 143.0 (quat.), 140.0 (quat.), 138.9 (CH), 138.9 × 2 (CH + quat.), 138.1 (CH), 129.5 (CH), 126.5 (CH), 126.3 (CH), 72.3 (CH), 21.2 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C13H12ClNO2 + H+) calculated: 250.06348; observed: 250.06314.
Molecules 27 02816 i002

3.3. 2-Chloro-5-([4-methylphenyl]carbonyl)-pyridine N-oxide (5a)

To a solution of impure 4a (ca 750 mg, ca 3.2 mmol) in dichloromethane (30 mL), manganese dioxide (5.2 g, 60 mmol) was added at 0 °C. The reaction mixture was stirred at room temperature overnight. The MnO2 was removed by filtration and the solvent was evaporated under reduced pressure to yield 670 mg of off-white crystals (2.71 mmol, 84%). 1H NMR (400 MHz, CDCl3) δ 8.68 (d, J = 8.2 Hz, 1H), 7.71 (d, J = 7.8 Hz, 2H), 7.64 (d, J = 8.2 Hz, 1H), 7.58 (dd, J = 8.7, 1.83 Hz, 1H), 7.33 (d, 2H, J = 8.2 Hz), 2.47 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 190.8 (quat.), 145.3 (quat.), 145.0 (quat.), 141.7 (CH), 134.8 (quat.), 132.9 (quat.), 130.2 (CH), 129.7 (CH), 127.0 (CH), 126.1 (CH), 21.9 (CH3). IR: 1650 (C=O, s); 1292 (N-O, m). HRMS (Positive ESI), m/z: [MH]+ (C13H10ClNO2 + H+) calculated: 248.04783; observed: 248.04715.
Molecules 27 02816 i003

3.4. 1-Hydroxy-5-([4-methylphenyl]carbonyl)-2(1H)-pyridinone (1a)

A solution of aqueous potassium hydroxide (10% w/v, 7.5 mL) was added to a round bottom flask along with 5a (0.670 g, 2.7 mmol) and stirred at 70 °C for 72 h. The reaction mixture was then cooled to 0 °C and concentrated hydrochloric acid was added dropwise until pH 1. The resulting precipitate was collected by Büchner filtration, washed with water and finally the crude product was recrystallised in ethanol to afford 232 mg of fine yellow needles (1.01 mmol, 37%). 1H NMR (400 MHz, CDCl3) δ 9.42 (s, 1H), 8.06 (d, J = 9.4 Hz, 1H), 7.38–7.43 (m, 4H), 6.80 (d, J = 9.4 Hz, 1H), 2.48 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 187.1 (CH), 158.1 (quat.), 149.9 (quat.), 141.6 (quat.), 135.4 (CH), 130.1 (CH), 129.4 (CH), 122.8 (quat.), 116.0 (quat.), 114.8 (CH), 21.6 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C13H11NO3 + H+) calculated: 230.08172; observed: 230.08101.
Molecules 27 02816 i004

3.5. 2-Chloro-4-([4-methylphenyl]hydroxymethyl)-pyridine (3b)

To 95 mL of dry THF under an argon atmosphere, 0.99 g of magnesium turnings (41.2 mmol) and 7.61 g of 4-bromotoluene (44.5 mmol) were mixed and the mixture was heated to reflux for 2 h. The solution was then cooled to 0 °C and 40 mL of dry THF containing 5.0 g of 2-chloropyridine-4-carboxaldehyde (35.3 mmol) was added dropwise. The reaction was stirred at 0 °C for 30 min then quenched with 75 mL of water. The organic phase was separated and the aqueous phase was extracted 3 times with 50 mL of ethyl acetate. The combined organic phases were washed twice with 75 mL of water, dried over magnesium sulphate and evaporated. The solid residue was triturated with petroleum ether (60–80) then recrystallised in acetonitrile to give 4.69 g of the pure product (20.1 mmol, 57%). 1H NMR (400 MHz, CDCl3) δ 8.20 (d, J = 5.5 Hz, 1H), 7.40 (s, 1H), 7.13–7.21 (m, 5H), 5.68 (d, J = 3.7 Hz, 1H), 4.47 (d, J = 3.7 Hz, 1H). 13C NMR (100 MHz, CDCl3) δ 156.2 (quat.), 151.7 (quat.), 149.4 (CH), 139.2 (quat.), 138.5 (quat.), 129.7 (CH), 126.8 (CH), 121.6 (CH), 120.1 (CH), 74.4 (CH), 21.2 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C13H12ClNO + H+) calculated: 234.06857 and 236.06562; observed: 234.06761 and 236.06429.
Molecules 27 02816 i005

3.6. 2-Chloro-4-([4-methylphenyl]carbonyl)-pyridine N-oxide (5b)

In total, 3.0 g of 3b (12.8 mmol) and 10.0 g of meta-chloroperoxybenzoic acid (77 wt.%, 44.6 mmol) were added to 50 mL of CH2Cl2 and the mixture was heated at 60 °C overnight. The solid that formed upon cooling of the solution was filtered and the organic phase was washed with three times 50 mL of saturated sodium hydrogencarbonate. The organic phase was then dried with magnesium sulphate and evaporated. The compound was purified on silica gel using 2.5% methanol in dichloromethane to yield 1.97 g of product (7.95 mmol, 62%). 1H NMR (400 MHz, CDCl3) δ 8.40 (d, J = 6.9 Hz, 1H), 7.93 (d, J = 2.3 Hz, 1H), 7.69 (d, J = 8.4 Hz, 2H), 7.63 (dd, J = 6.9, 2.3 Hz, 1H), 7.35 (d, J = 8.4 Hz, 2H), 2.47 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 190.9 (quat.), 144.7 (quat.), 142.3 (quat.), 140.4 (CH), 133.8 (quat.), 133.0 (quat.), 129.9 (CH), 129.6 (CH), 127.7 (CH), 124.3 (CH), 21.8 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C13H10ClNO2 + H+) calculated: 248.04783 and 250.04488; observed: 248.04700 and 250.04359.
Molecules 27 02816 i006

3.7. 1-Hydroxy-4-([4-methylphenyl]carbonyl)-2(1H)-pyridinone (1b)

A mixture of 5b (0.32 g, 1.29 mmol) in 4 mL of a 10% potassium hydroxide solution was heated to 70 °C for 72 h. Upon cooling, a brown precipitate formed that was filtered. This solid was then redissolved in 5 mL of warm water and carefully acidified with concentrated HCl to pH 1. The yellow precipitate was filtered, rinsed with cold water to give 0.25 g of pure product (1.09 mmol, 84%). 1H NMR (400 MHz, CDCl3) δ 7.92 (d, J = 7.0 Hz, 1h), 7.74 (d, J = 8.0 Hz, 2H), 7.31 (d, J = 8.0 Hz, 2H), 6.98 (d, J = 2.0 Hz, 1H), 6.67 (dd, J = 7.0, 2.0 Hz, 1H), 2.45 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 193.5 (quat.), 157.5 (quat.), 146.7 (quat.), 145.2 (quat.), 132.9 (quat.), 132.3 (CH), 130.5 (CH), 129.6 (CH), 118.8 (CH), 105.8 (CH), 22.0 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C13H11NO3 + H+) calculated: 230.08172; observed: 230.08063.
Molecules 27 02816 i007

3.8. 1-Methoxy-5-([4-methylphenyl]carbonyl)-2(1H)-pyridinone (7b)

To 3 mL of dry DMF, 80 mg of 1b (0.35 mmol), 170 mg of methyl iodide (1.20 mmol) and 150 mg of anhydrous potassium carbonate (1.08 mmol) were added and the mixture was stirred at room temperature overnight. The solid was then filtered, the solvent evaporated and water and ethyl acetate were added to the residue (10 mL each). The organic phase was separated. The aqueous phase was extracted twice with 10 mL of ethyl acetate and the combined organic phases were dried over magnesium sulphate and evaporated to give 54.4 mg of pure 7b (0.22 mmol, 63%). 1H NMR (400 MHz, CDCl3) δ 7.75 (d, J = 8.0 Hz, 2H), 7.66 (d, J = 7.3 Hz, 1H), 7.30 (d, J = 8.0s Hz, 2H), 6.91 (d, J = 2.3 Hz, 1H), 6.48 (dd, J = 7.3, 2.3 Hz, 1H), 4.14 (s, 3H), 2.45 (s, 3H). 13C NMR (100 MHz, CDCl3) δ 193.4 (Cq), 157.8 (Cq), 147.2 (Cq), 145.0 (Cq), 135.8 (CH), 132.5 (Cq), 130.2 (CH), 129.4 (CH), 124.0 (CH), 104.5 (CH), 64.9 (CH3), 21.7 (CH3). HRMS (Positive ESI), m/z: [MH]+ (C14H13NO3 + H+) calculated: 244.09736; observed: 244.09732.
Molecules 27 02816 i008

3.9. Physicochemical Studies

Water was purified via reverse osmosis, followed by passing through a mixed bed of ion exchange resin (Bioblock Scientific R3-83002, M3-83006) and activated carbon (Bioblock Scientific ORC-83005). All of the stock solutions were prepared by weighing solid products using an AG 245 Mettler Toledo analytical balance (precision 0.01 mg). Fe3+ cation stock solution was prepared from its perchlorate salt and its concentration was determined spectrophotometrically [35]. All the experiments were carried out at least in duplicate.
The acid-base properties (log K) of ligands 1a and 1b and their affinity for Fe3+ were determined via spectrophotometric titrations versus pH between pH 2 and 12. The titrations were carried out in a mixed MeOH/H2O 80/20 w/w (I = 0.1 M NaClO4) solvent due to their very low solubility in water. Sodium hydroxide (NaOH) and perchloric acid (HclO4) were used to adjust pH during titrations. The ionic strength of all the solutions was fixed to 0.1 M with sodium perchlorate (NaClO4). The measurement of pH was achieved by the use of combined glass electrodes (Metrohm, 6.0234.100, Long Life) filled with 0.1 M NaCl in MeOH/H2O (80/20 w/w). The electrodes were calibrated daily as hydrogen concentration probes by titrating known amounts of hydrochloric acid with CO32− free sodium hydroxide solutions. The GLEE program [36] was used for the glass electrode calibration with a pKw of 14.42 for studies in MeOH/H2O (80/20 w/w).
Typically, an aliquot of 40 mL of ligand solution was introduced into a thermostated jacketed titration vessel (25.0(2) °C), with an additional 0.3 equiv of Fe3+ in the case of complexation titrations. A known volume of perchloric acid solution was added to adjust the pH value to around 2, and the titrations were carried out between pH 2 and 12 by addition of known volumes of sodium hydroxide solution with a Metrohm 904 DMS Titrino automatic titrator (Methrom AG; Herisau, Switzerland) equipped with a 2mL Dosino 800 burette. After each addition, the pH was allowed to equilibrate and a UV-visible spectrum was recorded automatically with an Agilent Cary 60 UV-visible spectrophotometer.
The spectrophotometric data were fitted with the HypSpec software [24,25], to calculate the protonation constants of the ligands (log K), the stability constants (log β) of the formed complex species and the coordination model of the studied systems. The data for Fe3+ hydrated species and their solubility products were taken into account in the equilibrium model.

3.10. Modelling Methods

Structures used in this study were Protein Data Bank entries 3s68 [9] (rat COMT + Tolcapone as inhibitor) and 5lsa [37] (human COMT + dinitrocatechol, DNC, as inhibitor). Rat and human S-COMT enzymes show a high degree of similarity; they differ by 40 amino acids, of which 3 are in the active site; these three are MET89(ILE139), MET91(ILE141) and TYR95(CYS145) in rat (human COMT numbering in parenthesis as given in the PDB files). All structures used for the simulations presented here were derived from these two structures. Protonation states of the proteins were determined using the PROPKA program [38].
In preparation of the structures for molecular dynamics simulations, crystallographic water molecules, including the one water molecule that coordinates the Mg2+ ion, were kept in the initial structure. Initial preparations were carried out using the CHARMM program [39].
The CHARMM36 all-atom force field was used for the protein [40] and the CHARMM generalized force field [41] was used to represent the small molecule ligands. The parameters for the ligands were obtained from the PARAMCHEM website [42] which generates parameters for novel ligands based on the CHARMM generalized force field. The partial atomic charges were then refined through quantum mechanical calculations using the Gaussian program [43]. To be consistent with the level used in the development of the CHARMM36 all-atom force field, the quantum calculations were carried out at the MP2/6-31G(d) level.
The system was energy minimised for 700 steps using the steepest descent algorithm with a non-bonded cutoff distance of 12 Å and the dielectric constant equal to 4. This minimisation had as its purpose to reduce any strong van der Waals clashes that may be present in the crystal structure. Following the initial minimisation, the system was solvated in an 85 × 85 × 85 Å side box of explicit TIP3P water molecules. Sodium and chlorine ions were added to the system to neutralise the charge and introduce a physiological ionic strength to the solution.
The NAMD program version 2.13 was used for the simulations which were carried out using the NPT isothermal–isobaric ensemble using the standard options in the NAMD program [42]. Periodic boundary conditions were used for all simulations. Long-range electrostatic interactions were accounted for the particle mesh Ewald (PME) method. The water molecules were first relaxed around the fixed protein by 1000 steps of Conjugate Gradient (CG) energy minimisation using a constant dielectric coefficient with ε = 1. Under the constraint of the fixed protein, the water and ions were heated to 600 K over 23 ps followed by 250 steps of CG minimisation, and finally heating over 25 ps to reach the temperature of 300 K. Next, the positional constraints were removed and the entire system was subject to 2000 steps of CG, followed by heating to 300 K over 15 ps. An integration time step of 2 fs was used for all simulations. Each simulation was run for 200 ns.
Calculations were carried out for rat and human COMT in complex with tolcapone, a known inhibitor of COMT. Simulations for rat and human COMT in complex with compounds 1a and 1b were also carried out. No experimental structures exist for these complexes, so they were built by comparative docking, i.e., using tolcapone or DNC as structural guides for placing 1a and 1b in the COMT binding sites.

4. Conclusions

In conclusion, we have designed and synthesised analogues of tolcapone that could constitute the first generation of a multi-functional treatment against the symptoms and progress of the disease. Modelling of the binding of Mg2+ in COMT indicated that they might indeed be inhibitors of the enzyme. These compounds are expected to act as COMT inhibitors in the CNS without its known side effects caused by the nitro group. Calculation of the BBB Score indicated that they are likely to penetrate the CNS. Calculation of their pKa and stability constant for Fe(III) suggests that they are suitable Fe(III) chelators and thus might limit iron-mediated formation of ROS. The biological evaluation of these molecules as COMT inhibitors and as neuro-rescue agents will be reported in due course.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27092816/s1, Figure S1: Electronic spectra and distribution curves of the protonated species of ligand 1a; Figure S2: Electronic spectra and distribution curves of Fe3+ complexes of 1a; Figure S3: Sequestration diagrams; Figure S4: Starting conformation; Figure S5. Starting conformations; Figure S6: Starting conformations; Figure S7: Backbone RMSD by-residue; Figure S8: Superposition of average structures; Figure S9. Backbone RMSD by-residue; Figure S10: Superposition of average structures; Table S1: Binding Free energy scores with absolute scores with standard deviations and differences between doubly protonated tolcapone score and the scores of the other compounds and orientations.

Author Contributions

D.T.: Conceptualization of the molecules, writing—review and editing; J.C.J.B., K.K.T., A.K.N.: Synthesis of the molecules; V.V.; Investigation, Formal analysis, writing—review and editing, visualization, A.D.; Investigation, writing—review and editing, resources, R.H.S.; Conceptualization, Formal analysis, Investigation, resources, supervision, visualization, writing—original draft preparation, J.B.; Investigation, Formal analysis, supervision, writing—review and editing, C.-A.A. Investigation, Formal analysis, writing—review and editing; V.H.-B. writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

The synthesis was supported by internal funding at Northumbria University. VV, AD and RHS acknowledge support by the Centre National de la Recherche Scientifique, the Institut National de la Santé et de la Recherche Médicale and the University of Strasbourg. We acknowledge particularly the support of the Strasbourg University High Performance Computing Center). This work was granted access to the HPC resources of TGCC under the allocation 2021-A0070705109 made by GENCI. J.B, C.-A.A. and V.H.-B. acknowledge support by the Centre National de la Recherche Scientifique and the University of Strasbourg.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article or Supplementary Material.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

No samples are available from the authors as they will be tested in the context of developing collaborations.

References

  1. Barnham, K.J.; Bush, A.I. Biological metals and metal-targeting compounds in major neurodegenerative diseases. Chem. Soc. Rev. 2014, 43, 6727–6749. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Guiney, S.J.; Adlard, P.A.; Bush, A.I.; Finkelstein, D.I.; Ayton, S. Ferroptosis and cell death mechanisms in Parkinson’s disease. Neurochem. Int. 2017, 104, 34–48. [Google Scholar] [CrossRef] [PubMed]
  3. Masaldan, S.; Bush, A.I.; Devos, D.; Rolland, A.S.; Moreau, C. Striking while the iron is hot: Iron metabolism and ferroptosis in neurodegeneration. Free Rad. Biol. Med. 2019, 133, 221–233. [Google Scholar] [CrossRef] [PubMed]
  4. Kiss, L.E.; Soares-da-Silva, P. Medicinal chemistry of catechol O-methyltransferase (COMT) inhibitors and their therapeutic utility. J. Med. Chem. 2014, 57, 8692–8717. [Google Scholar] [CrossRef]
  5. Akhtar, M.J.; Yar, M.S.; Grover, G.; Nath, R. Neurological and psychiatric management using COMT inhibitors: A review. Bioorg. Chem. 2020, 94, 103418. [Google Scholar] [CrossRef]
  6. Forsberg, M.M.; Huotari, M.; Savolainen, J.; Mannisto, P.T. The role of physicochemical properties of entacapone and tolcapone on their efficacy during local intrastriatal administration. Eur. J. Pharm. Sci. 2005, 24, 503–511. [Google Scholar] [CrossRef]
  7. Forsberg, M.; Lehtonen, M.; Heikkinen, M.; Savolainen, J.; Järvinen, T.; Männistö, P.T. Pharmacokinetics and Pharmacodynamics of Entacapone and Tolcapone after Acute and Repeated Administration: A Comparative Study in the Rat. J. Pharmacol. Exp. Ther. 2003, 304, 498–506. [Google Scholar] [CrossRef]
  8. Ehler, A.; Benz, J.; Schlatter, D.; Rudolph, M.G. Mapping the conformational space accessible to catechol-O-methyltransferase. Acta Cryst. D 2014, 70 Pt 8, 2163–2174. [Google Scholar] [CrossRef]
  9. Ellermann, M.; Lerner, C.; Burgy, G.; Ehler, A.; Bissantz, C.; Jakob-Roetne, R.; Paulini, R.; Allemann, O.; Tissot, H.; Grünstein, D.; et al. Catechol-O-methyltransferase in complex with substituted 3′-deoxyribose bisubstrate inhibitors. Acta Cryst. D 2012, D68, 253–260. [Google Scholar] [CrossRef]
  10. Novaroli, L.; Bouchard Doulakas, G.; Reist, M.; Rolando, B.; Fruttero, R.; Gasco, A.; Carrupt, P.-A. The Lipophilicity Behavior of Three Catechol-O-Methyltransferase Inhibitors and Simple Analogues. Helv. Chim. Acta 2006, 89, 144–152. [Google Scholar] [CrossRef]
  11. Smith, K.S.; Smith, P.L.; Heady, T.N.; Trugman, J.M.; Harman, W.D.; Macdonald, T.L. In Vitro Metabolism of Tolcapone to Reactive Intermediates: Relevance to Tolcapone Liver Toxicity. Chem. Res. Toxicol. 2003, 16, 123–128. [Google Scholar] [CrossRef]
  12. Ernst, G.; Akuma, D.; Au, V.; Buchler, I.P.; Byers, S.; Carr, G.V.; Defays, S.; de Leon, P.; Demaude, T.; DePasquale, M.; et al. Synthesis and Evaluation of Bicyclic Hydroxypyridones as Inhibitors of Catechol O-Methyltransferase. ACS Med. Chem. Lett. 2019, 10, 1573–1578. [Google Scholar] [CrossRef] [PubMed]
  13. Zhao, Z.; Harrison, S.T.; Schubert, J.W.; Sanders, J.M.; Polsky-Fisher, S.; Zhang, N.R.; McLoughlin, D.; Gibson, C.R.; Robinson, R.G.; Sachs, N.A.; et al. Synthesis and optimization of N-heterocyclic pyridinones as catechol-O-methyltransferase (COMT) inhibitors. Bioorg. Med. Chem. Lett. 2016, 26, 2952–2956. [Google Scholar] [CrossRef] [PubMed]
  14. Harrison, S.T.; Poslusney, M.S.; Mulhearn, J.J.; Zhao, Z.; Kett, N.R.; Schubert, J.W.; Melamed, J.Y.; Allison, T.J.; Patel, S.B.; Sanders, J.M.; et al. Synthesis and Evaluation of Heterocyclic Catechol Mimics as Inhibitors of Catechol-O-methyltransferase (COMT). ACS Med. Chem. Lett. 2015, 6, 318–323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Robinson, R.G.; Smith, S.M.; Wolkenberg, S.E.; Kandebo, M.; Yao, L.; Gibson, C.R.; Harrison, S.T.; Polsky-Fisher, S.; Barrow, J.C.; Manley, P.J.; et al. Characterization of non-nitrocatechol pan and isoform specific catechol-O-methyltransferase inhibitors and substrates. ACS Chem. Neurosci. 2012, 3, 129–140. [Google Scholar] [CrossRef] [Green Version]
  16. Workman, D.G.; Tsatsanis, A.; Lewis, F.W.; Boyle, J.P.; Mousadoust, M.; Hettiarachchi, N.T.; Hunter, M.; Peers, C.; Tétard, D.; Duce, J.A. Protection from Neurodegeneration in the 6-Hydroxydopamine (6-OHDA) Model of Parkinson’s with Novel 1-Hydroxypyridin-2-one Metal Chelators. Metallomics 2015, 7, 867–876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Lewis, F.W.; Fairooz, S.; Elson, J.L.; Hubscher-Bruder, V.; Brandel, J.; Soundararajan, M.; Smith, D.; Dexter, D.T.; Tetard, D.; Pienaar, I.S. Novel 1-hydroxypyridin-2-one metal chelators prevent and rescue ubiquitin proteasomal-related neuronal injury in an in vitro model of Parkinson’s disease. Arch. Toxicol. 2020, 94, 813–831. [Google Scholar] [CrossRef] [Green Version]
  18. Bravo, A.; Bjorsvik, H.-R.; Fontana, F.; Minisci, F.; Serri, A. Radical versus “Oxenoid” Oxygen Insertion Mechanism in the Oxidation of Alkanes and Alcohols by Aromatic Peracids. New Synthetic Developments. J. Org. Chem. 1996, 61, 9409–9416. [Google Scholar] [CrossRef]
  19. Gutbier, S.; Kyriakou, S.; Schildknecht, S.; Uckert, A.K.; Brull, M.; Lewis, F.; Dickens, D.; Pearson, L.; Elson, J.L.; Michel, S.; et al. Design and evaluation of bi-functional iron chelators for protection of dopaminergic neurons from toxicants. Arch. Toxicol. 2020, 94, 3105–3123. [Google Scholar] [CrossRef]
  20. Ando, M.; Sato, N.; Nagase, T.; Nagai, K.; Ishikawa, S.; Takahashi, H.; Ohtake, N.; Ito, J.; Hirayama, M.; Mitobe, Y.; et al. Discovery of pyridone-containing imidazolines as potent and selective inhibitors of neuropeptide Y Y5 receptor. Bioorg. Med. Chem. 2009, 17, 6106–6122. [Google Scholar] [CrossRef]
  21. Gupta, M.; Lee, H.J.; Barden, C.J.; Weaver, D.F. The Blood-Brain Barrier (BBB) Score. J. Med. Chem. 2019, 62, 9824–9836. [Google Scholar] [CrossRef] [PubMed]
  22. Dichiara, M.; Amata, B.; Turnaturi, R.; Marrazzo, A.; Amata, E. Tuning Properties for Blood-Brain Barrier Permeation: A Statistics-Based Analysis. ACS Chem. Neurosci. 2020, 11, 34–44. [Google Scholar] [CrossRef] [PubMed]
  23. Crumbliss, A.L.; Harrington, J.M. Iron Sequestration by Small Molecules: Thermodynamic and Kinetic Studies of Natural Siderophores and Synthetic Model Compounds. Adv. Inorg. Chem. 2009, 61, 179–250. [Google Scholar]
  24. Gans, P.; Sabatini, A.; Vacca, A. Investigation of equilibria in solution. Determination of equilibrium constants with the HYPERQUAD suite of programs. Talanta 1996, 43, 1739–1753. [Google Scholar] [CrossRef]
  25. Gans, P.; Sabatini, A.; Vacca, A. Stability Constants. Available online: http://www.hyperquad.co.uk (accessed on 19 January 2022).
  26. Crea, F.; De Stefano, C.; Foti, C.; Milea, D.; Sammartano, S. Chelating Agents for the Sequestration of Mercury(II) and Monomethyl Mercury(II). Curr. Med. Chem. 2014, 21, 3819–3836. [Google Scholar] [CrossRef]
  27. Lewis, F.W.; Bird, K.; Navarro, J.-P.; El Fallah, R.; Brandel, J.; Hubscher-Bruder, V.; Tsatsanis, A.; Duce, J.A.; Tétard, D.; Bourne, S.; et al. Synthesis, Physicochemical Characterization and Neuroprotective Evaluation of Novel 1-hydroxypyrazin-2(1H)-one Iron Chelators in an In Vitro Cell Model of Parkinson’s Disease. Dalton Trans. 2022, 59, 3590–3603. [Google Scholar] [CrossRef] [PubMed]
  28. Patel, N.C. Methods to optimize CNS exposure of drug candidates. Bioorg. Med. Chem. Lett. 2020, 30, 127503. [Google Scholar] [CrossRef]
  29. Vilar, S.; Sobarzo-Sanchez, E.; Santana, L.; Uriarte, E. Ligand and Structure-based Modeling of Passive Diffusion through the Blood-Brain Barrier. Curr. Med. Chem. 2018, 25, 1073–1089. [Google Scholar] [CrossRef]
  30. Morales, J.F.; Scioli Montoto, S.; Fagiolino, P.; Ruiz, M.E. Current State and Future Perspectives in QSAR Models to Predict Blood-Brain Barrier Penetration in Central Nervous System Drug R&D. Mini Rev. Med. Chem. 2017, 17, 247–257. [Google Scholar]
  31. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
  32. Xing, L.; Glen, R.C. Novel Methods for the Prediction of log P, pKa, and logD. J. Chem. Inf. Comput. Sci. 2002, 42, 796–805. [Google Scholar] [CrossRef] [PubMed]
  33. Patra, N.; Ioannidis, E.I.; Kulik, H.J. Computational Investigation of the Interplay of Substrate Positioning and Reactivity in Catechol O-Methyltransferase. PLoS ONE 2016, 11, e0161868. [Google Scholar] [CrossRef] [PubMed]
  34. Kollman, P.A.; Massova, I.; Reyes, C.; Kuhn, B.; Huo, S.; Chong, L.; Lee, M.; Lee, T.; Duan, Y.; Wang, W.; et al. Calculating Structures and Free Energies of Complex Molecules:  Combining Molecular Mechanics and Continuum Models. Acc. Chem. Res. 2000, 33, 889–897. [Google Scholar] [CrossRef]
  35. Fish, W.W. Rapid colorimetric micromethod for the quantitation of complexed iron in biological samples. Methods Enzymol. 1988, 158, 357–364. [Google Scholar] [PubMed]
  36. Gans, P.; O’Sullivan, B. GLEE, a new computer program for glass electrode calibration. Talanta 2000, 51, 33–37. [Google Scholar] [CrossRef]
  37. Ehler, A.; Lerner, C.; Rudolph, M.G. Human Catechol O-Methyltransferase in Complex with SAM and DNC at 1.50A. Available online: https://www.rcsb.org/structure/5LSA (accessed on 19 January 2022).
  38. Dolinsky, T.J.; Nielsen, J.E.; McCammon, J.A.; Baker, N.A. PDB2PQR: An automated pipeline for the setup of Poisson–Boltzmann electrostatics calculations. Nucleic Acids Res. 2004, 32, W665–W667. [Google Scholar] [CrossRef] [PubMed]
  39. Brooks, B.R.; Brooks III, C.L.; Mackerell, A.D., Jr.; Nilsson, L.; Petrella, R.J.; Roux, B.; Won, Y.; Archontis, G.; Bartels, C.; Boresch, S.; et al. CHARMM: The biomolecular simulation program. J. Comput. Chem. 2009, 30, 1545–1614. [Google Scholar] [CrossRef]
  40. MacKerell, A.D., Jr.; Bashford, D.; Bellott, M.; Dunbrack, R.L., Jr.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef]
  41. Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; et al. CHARMM General Force Field (CGenFF): A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields. J. Comput. Chem. 2010, 31, 671–690. [Google Scholar]
  42. Vanommeslaeghe, K.; Mackerell, A.D., Jr. The CHARMM General Force Field. Available online: https://cgenff.umaryland.edu/ (accessed on 19 January 2022).
  43. Frisch, M.J. Gaussian 09, Revision A.02; Gaussian Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
Figure 1. Examples of COMT inhibitors based on the nitrocatechols moiety [12,13,14,15].
Figure 1. Examples of COMT inhibitors based on the nitrocatechols moiety [12,13,14,15].
Molecules 27 02816 g001
Figure 2. Target compounds 1a and 1b.
Figure 2. Target compounds 1a and 1b.
Molecules 27 02816 g002
Scheme 1. Synthesis of compounds 1a. (i) 4-bromotoluene, Mg, THF (93%), (ii) mCPBA, DCM, (iii) MnO2, DCM (84%), (iv) KOH (37%).
Scheme 1. Synthesis of compounds 1a. (i) 4-bromotoluene, Mg, THF (93%), (ii) mCPBA, DCM, (iii) MnO2, DCM (84%), (iv) KOH (37%).
Molecules 27 02816 sch001
Scheme 2. Synthesis of compound 1b. (i) 4-bromotoluene, Mg, THF (57%), (ii) mCPBA, DCM (62%), (iii) KOH (84%).
Scheme 2. Synthesis of compound 1b. (i) 4-bromotoluene, Mg, THF (57%), (ii) mCPBA, DCM (62%), (iii) KOH (84%).
Molecules 27 02816 sch002
Scheme 3. Synthesis of control compound 7b.
Scheme 3. Synthesis of control compound 7b.
Molecules 27 02816 sch003
Figure 3. Spectrophotometric titrations vs. pH of (A) Ligand 1a between 2.56 < pH < 8.68, [1a] = 5.00 × 10−5 M and (B) ligand 1b between 2.72 < pH < 8.79, [1b] = 1.00 × 10−4 M. (C) Electronic spectra of the protonated species of ligand 1b and (D) distribution curves of the protonated species of ligand 1b ([1b] = 1.00 × 10−4 M). Solvent: MeOH/H2O (80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C. N.B.: (1b)H represents the protonated ligand, i.e., compound 1b in Scheme 3 and (1b)- the deprotonated ligand.
Figure 3. Spectrophotometric titrations vs. pH of (A) Ligand 1a between 2.56 < pH < 8.68, [1a] = 5.00 × 10−5 M and (B) ligand 1b between 2.72 < pH < 8.79, [1b] = 1.00 × 10−4 M. (C) Electronic spectra of the protonated species of ligand 1b and (D) distribution curves of the protonated species of ligand 1b ([1b] = 1.00 × 10−4 M). Solvent: MeOH/H2O (80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C. N.B.: (1b)H represents the protonated ligand, i.e., compound 1b in Scheme 3 and (1b)- the deprotonated ligand.
Molecules 27 02816 g003
Figure 4. Spectrophotometric titration vs. pH of Fe(III) complexes of (A) ligand 1a (2.68 ≤ pH ≤ 11.26, [1a] = 1.01 × 10−4 M, [Fe3+] = 3.20 × 10−5 M) and (B) ligand 1b (2.95 ≤ pH ≤ 10.45, [1b] = 1.01 × 10−4 M, [Fe3+] = 4.35 × 10−5 M). (C) Electronic spectra of Fe3+ complexes of 1b and (D) distribution curves of Fe(III) complexes of 1b ([1b] = 1.01 × 10−4 M, [Fe3+] = 3.20 × 10−5 M). Solvent: MeOH/H2O (80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C.
Figure 4. Spectrophotometric titration vs. pH of Fe(III) complexes of (A) ligand 1a (2.68 ≤ pH ≤ 11.26, [1a] = 1.01 × 10−4 M, [Fe3+] = 3.20 × 10−5 M) and (B) ligand 1b (2.95 ≤ pH ≤ 10.45, [1b] = 1.01 × 10−4 M, [Fe3+] = 4.35 × 10−5 M). (C) Electronic spectra of Fe3+ complexes of 1b and (D) distribution curves of Fe(III) complexes of 1b ([1b] = 1.01 × 10−4 M, [Fe3+] = 3.20 × 10−5 M). Solvent: MeOH/H2O (80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C.
Molecules 27 02816 g004
Figure 5. Crystallographic orientation of tolcapone in the MD simulations. (A) Doubly deprotonated tolcapone and (B) singly deprotonated tolcapone.
Figure 5. Crystallographic orientation of tolcapone in the MD simulations. (A) Doubly deprotonated tolcapone and (B) singly deprotonated tolcapone.
Molecules 27 02816 g005
Figure 6. Different orientations of the compounds 1a (A,B) and 1b (C,D) used in the simulations.
Figure 6. Different orientations of the compounds 1a (A,B) and 1b (C,D) used in the simulations.
Molecules 27 02816 g006
Figure 7. RMSD time series along the simulations. Top figure, the simulations with rat COMT, bottom figure, the simulations with human COMT.
Figure 7. RMSD time series along the simulations. Top figure, the simulations with rat COMT, bottom figure, the simulations with human COMT.
Molecules 27 02816 g007
Figure 8. Backbone RMSF calculated from the MD simulations and compared to those calculated from the crystal structure of COMT in complex with tolcapone.
Figure 8. Backbone RMSF calculated from the MD simulations and compared to those calculated from the crystal structure of COMT in complex with tolcapone.
Molecules 27 02816 g008
Table 1. Protonation constants (log K) of ligands and overall stability constants (log β) of their Fe(III) complexes. Solvent: MeOH/H2O 80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C.
Table 1. Protonation constants (log K) of ligands and overall stability constants (log β) of their Fe(III) complexes. Solvent: MeOH/H2O 80/20 w/w), I = 0.1 M (NaClO4), T = 25.0 °C.
1a1bTolcapone [9]
log K LH6.9(1)5.69(6)10.20
log K LH2//4.6
log β FeL219.5(1)21.1(1)
log β FeL326.6(1)29.1(1)
pL0.5 (a)5.366.29
(a) pL0.5 represents the total concentration of ligand required for the sequestration of 50% of the metal [26].
Table 2. Moledcular descriptors relevant to BBB Score of tolcapone, 1a and 1b in their neutral form.
Table 2. Moledcular descriptors relevant to BBB Score of tolcapone, 1a and 1b in their neutral form.
CompoundMolar mass/g/molRotatable BondsHBA (a)HBD (b)TPSA (c)3Consensus log PPredicted BBB Permeant/Actual
Tolcapone273.24352103.351.91No/Yes
1a229.2423159.301.75Yes/n.d.
1b229.2423159.301.67Yes/n.d.
(a) Hydrogen bond acceptor, (b) Hydrogen bond donor, (c) Topological Polar Surface Area.
Table 3. Molecular Descriptors and BBB Score of tolcapone, 1a and 1b.
Table 3. Molecular Descriptors and BBB Score of tolcapone, 1a and 1b.
CompoundAro_RHABBB Score
Tolcapone2203.25
1a2174.97
1b2174.82
Table 4. Colour-coded descriptors linked to BBB permeation. Green indicates descriptors favourable for BBB permeation, yellow indicates borderline descriptors and red indicates descriptors having negative influence on permeation. All parameters are calculated using SwissADME [31].
Table 4. Colour-coded descriptors linked to BBB permeation. Green indicates descriptors favourable for BBB permeation, yellow indicates borderline descriptors and red indicates descriptors having negative influence on permeation. All parameters are calculated using SwissADME [31].
Tolcapone1a1b
log P (a)log D (b)TPSAHBAlog P (a)log D (b)TPSAHBAlog P (a)log D (b)TPSAHBA
1.91−0.9103.3551.750.159.3031.671.159.303
HBDNCOCNOCHBDNCOCNOCHBDNCOCNOC
215611341134
ISISIS
AcidNeutralAcid
(a) Consensus log Po/w; (b) Calculated at pH 7.4 from Ref [32].
Table 5. Relative Binding Free Energy Scores.
Table 5. Relative Binding Free Energy Scores.
SpeciesLigand OrientationRelative Scores (kcal/mole)
RatTolcapone doubly deprotonated (Figure 5a)0.0 ± 9.6
RatTolcapone singly deprotonated (Figure 5b) 17.9 ± 5.4
Rat1a = O→SAM (Figure 6a)35.6 ± 5.9
Rat1a N-O→SAM (Figure 6b)2.1 ± 5.1
Rat1b N-O→SAM (Figure 6c)10.5 ± 5.7
Rat1b = O→SAM (Figure 6d)3.3 ± 3.9
HumanTolcapone doubly deprotonated (Figure 5a)0.0 ± 5.4
HumanTolcapone singly deprotonated (Figure 5b) 5.8 ± 6.8
Human1a = O→SAM (Figure 6a) 31.4 ± 5.1
Human1a N-O→SAM (Figure 6b)−4.7 ± 5.3
Human1b N-O→SAM (Figure 6c) 23.1 ± 4.9
Human1b = O→SAM (Figure 6d)16.5 ± 4.7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bergin, J.C.J.; Tan, K.K.; Nelson, A.K.; Amarandei, C.-A.; Hubscher-Bruder, V.; Brandel, J.; Voinarovska, V.; Dejaegere, A.; Stote, R.H.; Tétard, D. 1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson’s Disease. Molecules 2022, 27, 2816. https://doi.org/10.3390/molecules27092816

AMA Style

Bergin JCJ, Tan KK, Nelson AK, Amarandei C-A, Hubscher-Bruder V, Brandel J, Voinarovska V, Dejaegere A, Stote RH, Tétard D. 1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson’s Disease. Molecules. 2022; 27(9):2816. https://doi.org/10.3390/molecules27092816

Chicago/Turabian Style

Bergin, Joseph C. J., Kean Kan Tan, Anya K. Nelson, Cristina-Andreea Amarandei, Véronique Hubscher-Bruder, Jérémy Brandel, Varvara Voinarovska, Annick Dejaegere, Roland H. Stote, and David Tétard. 2022. "1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson’s Disease" Molecules 27, no. 9: 2816. https://doi.org/10.3390/molecules27092816

APA Style

Bergin, J. C. J., Tan, K. K., Nelson, A. K., Amarandei, C. -A., Hubscher-Bruder, V., Brandel, J., Voinarovska, V., Dejaegere, A., Stote, R. H., & Tétard, D. (2022). 1-Hydroxy-2(1H)-pyridinone-Based Chelators with Potential Catechol O-Methyl Transferase Inhibition and Neurorescue Dual Action against Parkinson’s Disease. Molecules, 27(9), 2816. https://doi.org/10.3390/molecules27092816

Article Metrics

Back to TopTop