Design, synthesis, and in silico-in vitro antimalarial evaluation of 1,2,3-triazole-linked dihydropyrimidinone quinoline hybrids

In response to the malaria parasite’s resistance towards quinoline-based antimalarial drugs, we have employed quinoline-containing compounds in combination with dihydropyrimidinone (DHPM) analogues as resistance reversal agents (RAs) and investigated their antimalarial activities based on DHPM’s resistance reversal abilities. The present study employed click chemistry to link DHPM and quinoline compounds which offered several synthetic advantages over the previously used amide coupling for the same hybrids. Among the synthesised compounds, 4 hybrids with the 7-chloroquinoline moiety showed antimalarial activity below 1 µM while compounds with the mefloquine moiety showed lower antimalarial activity than chloroquine (CQ) and the 7-chloroquinoline hybrids. Among the tested hybrids for the IC50 determination, four compounds displayed good antimalarial activity with increased sensitivity against the CQ-resistant K1 strain between 421 and 567 nM and showed higher activity between 138 and 245 nM against the NF54 CQ-sensitive strain, while three compounds have IC50 values greater than 5 µM. Additionally, in silico molecular docking and molecular dynamics studies were conducted to investigate the binding affinities of all the synthesised compounds as glutathione reductase protein competitive inhibitors. Further optimisation of the compound with the highest binding affinity generated 16 compounds with higher binding affinities than the flavine adenine dinucleotide (FAD) cofactor.


Introduction
Recent data indicates that the challenge to overcome malaria still lingers [1], regardless of all the significant efforts and advances in drug developments and medical technologies.The rise of drug or multidrug resistance in Plasmodia strains has made malaria the most important tropical parasitic disease [2].Therefore, developing new antimalarial drugs is necessary to rectify this situation.Despite the development of many antimalarial leads, quinoline compounds remain important in malarial chemotherapy [3].
For many years, chloroquine (CQ) from the quinoline family was the first-line antimalarial drug.Despite the huge success recorded in using quinoline for malaria therapeutics, these drugs have been challenged by Plasmodium parasites due to their rapid development of resistance [4,5].Drug resistance in malaria occurs when the parasite strain continues to multiply or survive even though the therapeutic has been given in the recommended tolerable dose.This situation is not malaria-specific but is also observed in cancer and other parasitic diseases like the case with antibiotics.As for the case of chloroquine (CQ) resistance, there is strong evidence through different research that resistant malaria parasites accumulate less CQ in their food vacuole than sensitive parasites [6].This could result from an efflux pumping system, reduced CQ uptake by the parasite, or a combination of both processes [7].This implies that CQ concentration does not reach the parasite's lethal dose, thereby resulting in its continuous survival.
The use of resistance reversal agents (RAs) is among the strategies to reverse this drug resistance in malaria parasites.Resistance reversal refers to the use of a compound that can restore the sensitivity of resistant strains to a drug while it might have little or no therapeutic action against the parasite [8].In the context of malaria, resistance RAs ensure antimalarial drugs accumulate within the parasite's food vacuole.An interesting feature of this strategy is the ability to combine two biologically active compounds in a hybrid, maintaining the properties of both without the setback of differing elimination half-life (t 1/2 ) or pharmacokinetic properties associated with separate drugs that are used in combination therapy.Furthermore, a major advantage of this approach over the rest is the ability of the RA to competitively bind to the P. falciparum chloroquine resistance transporter (PfCRT) protein, thereby inhibiting the efflux of the CQ from the digestive food vacuole of the Plasmodium parasite.In contrast, other strategies lack this feature to bind to the mutated CQ-resistant P. falciparum [9,10].
Previously, it was discovered that concurrent exposure of resistant strains of P. falciparum to both CQ and verapamil reversed resistance in vitro [4,11].Since verapamil is a calcium channel blocker used to treat hypertension and other cardiovascular diseases, researchers have used different calcium channel blockers and other drugs in combination with quinoline-based antimalarial drugs in attempts to reverse antimalarial drug resistance.It was observed that the RA combinations were able to increase the sensitivity of cloned CQ-resistant P. falciparum strains in vitro but did not affect the sensitivity of CQ-sensitive strains [6,[11][12][13].
Before Peyton and co-workers coined the term "Reversed Chloroquine" (RCQ) to describe a hybrid compound containing CQ and an RA, the use of RAs in antimalarial chemotherapy involved the administration of the two agents as separate doses [14].They discovered that covalently linking a CQ-like moiety to imipramine was effective against both sensitive and resistant strains of the Plasmodium parasite.The hybrid molecule delivered the RA in a 1:1 ratio to the antimalarial quinoline, which would lower the dose required if the two drugs, was to be given in separate forms, with resultant lower cost and toxicity [4].This arrangement would favour the drug accumulation in the parasite's digestive vacuole and also interfere with the CQ export by the mutated CQ-resistant P. falciparum [14].
In our previous work, we covalently linked CQ to dihydropyrimidinone (DHPM) using an amide bond linker at position 3 of the DHPM.The resulting hybrids could reverse resistance in the K1 strain of P. falciparum in vitro [15].Since DHPMs present various sites of functionalisation that could be used for structure-activity relationship (SAR), several synthetic strategies could be employed to link quinolinebased antimalarials to various DHPM scaffolds.
The previous report of the synthesis of DHPM-based hybrids of chloroquine (CQ) moiety used amide coupling to link the CQ to DHPM, specifically at N3 of the reversal agent [15].The synthetic route to these hybrids involved the protection of N1 nitrogen with a methyl group.Furthermore, the synthesis of the amide carbamate ester intermediate involved overnight reflux with a large excess of phenyl chloroformate and NaH, 10.0, and 11.0 molar equivalents, respectively.This synthetic procedure generated a lot of waste.However, we report on the use of click chemistry as an alternative route to synthesise DHPM-based quinoline hybrids and their antimalarial activities.This project employs a triazole covalent linker between quinoline-based compounds and different DHPMs using the C6 of the DHPM as the point of attachment.Two series of quinoline-based compounds related to chloroquine and mefloquine antimalarial drugs were chosen (Scheme 1) It also has been demonstrated that quinoline derivatives block the P. falciparum glutathione reductase (PfGR) enzyme [16].The host erythrocytes undergo significant structural and physiological changes when invaded by P. falciparum to aid the parasites' survival and growth.The oxidative stress caused by these changes affects the malaria parasites.As a result, an effective antioxidant is necessary to inhibit reactive oxygen species from harming the parasites [17].It is thought that reduced glutathione (GSH) is a crucial antioxidant during erythrocytic infection.As an oxidoreductase, glutathione reductase in P. falciparum converts oxidised glutathione (GS-SG) to reduced GSH using hydrogen from nicotinamide adenine dinucleotide phosphate (NADPH).The transfer of hydrogen from NADPH to GS-SG is the most essential task of the FAD cofactor in PfGR because, without it, the parasite would experience oxidative stress [18,19].Thus, P. falciparum's glutathione reductase enzyme becomes a crucial target in antimalarial treatment [20].Therefore, in addition to the in vitro antimalarial testing, the binding affinities and the molecular stabilities of the synthesised compounds within the glutathione reductase enzyme were estimated using molecular docking and molecular dynamics in silico methods.

Chemistry
The synthesis of the chlorinated DHPM and the subsequent azidation to form azides 9a-d was reported earlier [21].The terminal alkynes were installed on the quinolinebased compounds using different protocols.Compound 2 was synthesised via a nucleophilic substitution reaction from 4,7-dichloroquinoline 1 according to the reported protocol (Scheme 2) [22].This involved heating a mixture of 1 and propargylamine in a sealed tube under a nitrogen atmosphere.Formation of product 2 was confirmed by 1 H NMR spectroscopy with the observance of the terminal alkyne's proton at δ H 2.66 ppm (t, 4 J HH = 2.5 Hz).The method reported by Kaval et al. was adapted to synthesise the oxyquinoline 3 (Scheme 2) [23].The formation of compound 3 was confirmed from the 1 H NMR spectrum with the terminal alkyne's proton appearing at δ H 2.62 ppm (t, 4 J HH = 2.4 Hz) and the Single crystal X-ray diffraction (SCXRD) technique (Fig. 1).The quinolinol 6 was synthesised via Conrad-Limpach synthesis of quinolines, which was readily converted to the alkynes 7 and 8 (Scheme 3).
Aniline 4 was quantitatively cyclised to the quinolinol 6 by heating it at 150 C with ethyl 4,4,4-trifluoroacetoacetate (5) for 3 h in the presence of polyphosphoric acid (PPA) as the catalyst [24].The nucleophilic alkylation of 6 to 7 was performed using the Williamson ether synthesis [25].To extend the carbon chain linker, compound 8 was synthesised from the quinolinol 6 via the Mitsunobu reaction [26] as Williamson's nucleophilic alkylation of the quinolinol 6 with excess 4-Bromo-1-butyne failed to give the desired product.From the mechanism of Conrad-Limpach synthesis of the quinolines [27,28], it is possible to obtain two tautomers for the quinolinol 6: the enol (a) and the keto (b) (Fig. 2).The simultaneous presence of both tautomers was confirmed from the SCXRD data (Fig. 2).It corroborated the earlier isolation by Sarojini et al. [29], but the subsequent deprotonations only involved the enol (a) as was observed from the SCXRD of compound 7 (Fig. 3).The involvement of the enol form for the subsequent deprotonations points to its preference over the keto form.The more stable enol tautomer is not unexpected for two reasons: enol stability by conjugation with a neighbouring pi system and aromaticity [30].The two factors are present as the quinoline moiety is highly conjugated in the two rings, and they both contain aromatic systems for additional resonance stability.The structure of compound 8 was confirmed from the 1 H NMR spectrum with the observation of all the diagnostic peaks.The terminal alkyne proton appeared at δ H 2.10 ppm (t, 4 J HH = 2.7 Hz), while the two methylene protons were assigned to the signals δ H 2.90 ppm (td, 3 J HH = 6.7 Hz, 4 J HH = 2.7 Hz) and 4.41 ppm (t, 3 J HH = 6.7 Hz).All the aromatic protons were also observed at the expected peak positions.The hybrid compounds were synthesised by click chemistry reactions (Schemes 4 and 5) as adapted from Guantai et al. [31].The reaction was first attempted in a mixture of tertbutanol and water using copper(I) iodide (CuI) as the catalyst.This gave three spots on the TLC, and no attempt was made to separate them because of their close R f values.The catalyst was changed to freshly prepared 1-M solution of copper(II) sulphate pentahydrate (CuSO 4 .5H 2 O) and sodium ascorbate which reductively produces the required Cu(I) in situ.Initially, this approach seemed to work for the reaction where only one prominent spot was observed on the TLC.However, a difficult-to-filter precipitate also formed from the solution.Another setback was the persistent appearance of the starting materials even after long hours of stirring at 25 °C.The solvent was changed to DMF with the catalyst freshly prepared in distilled water.The reaction mixture was stirred for 12 h, and the TLC analysis showed a single spot for the product with complete consumption of the starting materials.
The 1 H NMR spectra obtained for 11a-h correlated well with the proposed structures, and all the diagnostic peaks were observed (see the supplementary information for details).Additionally, the predicted connectivity of the compounds was confirmed from the SCXRD structure elucidations of 11d and 11 g (Fig. 4).

Antimalarial activity of DHPM-quinoline triazole hybrids
The percent growth inhibitory activities of DHPMs 9a-d and hybrids 10a-h and 11a-h were evaluated at 1-µM and 5-µM concentrations in an in vitro antiplasmodial assay (Fig. 5, Thermalellipsoids are drawn at a 50% probability level Table 1).An NF54 CQ-sensitive (CQS) P. falciparum asexual strain was used.The most potent compounds 10a-d were further evaluated for the IC 50 against NF54 CQS, and K1 CQ-resistant (CQR) strains, while the IC 50 of the moderately potent compounds 11c, 11e, and 11 g were determined against the K1 CQR strain only.In all the assays, CQ was used as the reference drug (Table 2).
Among the compounds tested against the NF54 strain to determine the percentage inhibition, the reversal agents (RAs) 9a-d displayed minimal inhibition (Fig. 5, Table 1).This was expected as they lack the antiplasmodial quinoline moiety, and their role in the hybrids is to inhibit the underlying process responsible for drug resistance.Within the same series, compounds 10a-d with 4-amino functionality displayed the highest inhibitions and were comparable to CQ at the two concentrations tested.Compounds 10e-f with the 4-oxyquinoline moiety exhibited very low inhibitions, even at 5 µM.The decrease in the activity of these oxygen-containing hybrids 10e-f could be attributed to the absence of the 4-amino functionality, which is present in compounds 10a-d.This agrees with previous reports that quinolines containing an amino group at the 4-position show higher antiplasmodial activities than those containing oxygen or lacking the 4-amino group [32][33][34][35].This stresses the importance of the 4-amino moiety in maintaining the basic property required for the CQ and its analogues to accumulate in the food vacuole of P. falciparum [34].Compounds 11a-h with the mefloquine moiety showed lower antiplasmodial activity than CQ and the hybrids 10a-h.Among these hybrids, compounds 11a, 11c, 11e, and 11 g, with unsubstituted acidic protons at N1 and N3, showed better activities than compounds 11b, 11d, 11f, and 11 h which instead have a methyl group at N1.The generally lower activities of hybrids 11a-h could also be ascribed to the lack of 4-amino functionality in the quinoline ring [34].
Compounds 10a-d displayed good antiplasmodial activities against the CQR K1 strain.The IC 50 s for this strain were between 421 and 567 nM, but much greater potencies of 138-245 nM were observed against the NF54 CQS strain.Compounds 11c, 11e, and 11 g were at least three orders of magnitude weaker.Hybrids 10b and 10d, containing chloroquine moiety with a methyl group at N1 of the DHPM, showed better antiplasmodial activities than hybrids 10a and 10c with unsubstituted acidic protons at N1 and N3 in both NF54 and K1 strains.This higher antiplasmodial activity could result from the decreased acidity associated with the removal of the acidic proton, which increases the basic property of the quinoline hybrid needed for its accumulation in the acidic food vacuole of the P. falciparum [36].In addition, replacing the acidic proton with a methyl group increases the lipophilicity and in turn increases the antimalarial activities [37].However, other substitutions, particularly the para-methoxy and the para-fluoro groups on the DHPM aromatic, do not seem to influence the activities of the hybrid compounds.Hybrids 10a-d also showed reduced Resistance indices (RI) factor ranging from 2.3 to 3.6 as compared to CQ with an RI factor of 8.5.The reduced RI factor suggests that the addition of the DHPM to the quinoline scaffold was able to change the sensitivity of the K1 strain observed.

In silico study
Using the Maestro tool from the Schrödinger software, we performed the in silico study utilising the crystal structure of glutathione reductase (PDB ID: 1ONF).The X-ray crystallographic structure of P. falciparum glutathione reductase was retrieved from the Protein Data Bank (PDB, https:// www.rcsb.org/).The raw protein structure was prepared using the protein preparation wizard of the Schrodinger suite.This was used to assign bond orders, remove water molecules beyond 3 Å, and add H atoms.In addition, both missing side chains and missing loops were added using the Prime.While refining, water orientations were sampled using the PROPKA at a pH of 7.0 for the energy optimisation process and the restrained minimisation was done to confine the root-mean-square deviation (RMSD) of heavy atoms to 0.30 Å using an OPLS4 force field [38].The active site of the protein residues surrounding the FAD cofactor (the bound ligand) was taken as the binding site, and these residues were used for the preparation of the docking grid using the receptor grid generation panel of the Schrodinger suite.No volume was excluded, and likewise, no constraint was defined.The input ligand structures of the synthesised compounds were built using the LigPrep module of the Schrödinger suite.For each 2D structure, the LigPrep process creates 3D structures with minimised energy, proper bond lengths, and angles.At a physiological pH of 7.0 ± 0.2, the possible ionisation states for each ligand structure were  The docking simulations of the synthesised compounds were run in the binding site of P. falciparum glutathione reductase crystal structure (PDB ID: 1ONF) co-crystallised with FAD cofactor, and the outcomes are shown in Table 3.The binding energy of the ligand-protein complex corresponds to the intensity and affinity of the ligand-protein interaction.The docked complexes were rated according to their lowest energy value (kcal/mol).The negative value of the binding energy increases with the increasing strength of the interaction and vice versa.As expected, the FAD cofactor showed the highest binding affinity at −9.644 kcal/mol, followed by compound 10a at −8.1 kcal/mol.Accidentally, no correlation was observed between the docking scores and the antiplasmodial IC 50 values of the synthesised compounds.As it is well known with competitive inhibitors, it is always a challenge to design compounds that could competitively bind at the active site with higher binding affinity than the cofactor [39]; however, rational optimisations could be a saving grace.Though all the synthesised compounds showed insignificantly lower binding affinity than the FAD cofactor, all of them displayed considerable binding energy, and five of the synthesised compounds showed higher binding affinity than the chloroquine standard.
Analysis of the binding pose of compound 10a showed different electronic interactions as shown in Fig. 6.The compound was able to make different interactions within the active site of glutathione reductase.The + NH 3 of the Lys151 side chain was able to simultaneously make two π-cation interactions with the phenyl group on the DHPM ring (at 4.69 Å) and the DHPM ring itself (at 4.94 Å), while  the benzene ring of the quinoline moiety was able to form another π-cation interaction with + NH 3 side chain of Lys32 (at 3.67 Å).In addition, several H-bond interactions were noted.These include H-bond interactions between N3H of the DHPM and the carboxylic OH of the Asp167 side chain at 1.90 Å, C = O of DHPM, and one NH of the guanidino group of Arg272 side chain at 1.89 A. Most importantly, the nitrogen of the pyridine ring of the quinoline moiety was able to form hydrophobically packed H-bond interaction with NH of Ala110 backbone at 1.92 Å (Fig. 7).This type of H-bond is very crucial in maintaining the ligand stability within a protein because of the difficulty involved in breaking the H-bond formed in a hydrophobic space [38].Another notable feature is the halogen interaction between the chlorine atom on the quinoline and the C = O backbone of Asn278 at 3.36 Å.

Optimisation of compound 10a
The synthesised compounds must exhibit stronger interactions than the FAD cofactor in order to function as competitive inhibitors at the active site of glutathione reductase.To achieve this, the protein-ligand complex of compound 10a with the highest binding affinity was loaded onto the workspace and rationally optimised using the Ligand Designer tool in the Maestro platform.Analysis of the workspace showed the 4-methoxy substituent on the phenyl ring of DHPM in a solvent-exposed region.The 4-methoxy group was therefore truncated, leading to an increase in the binding affinity from −8.1 to −8.2 kcal/mol.The growth space of compound 10a, the ligand-protein interaction, and the truncated 10a are shown in Fig. 8a, b, and c, respectively.
The light blue regions show the cavity space within the binding pocket, while the deep blue regions show the solventexposed areas.
The truncated 10a was then loaded into the workspace with the protein crystal structure.The pathfinder bonds were activated to identify possible points of R group attachments.With this, the acetate group at position 5 of the DHPM was selected and enumerated to be replaced with the default R groups in Maestro Ligand Designer.The methyl group of the acetate moiety was subsequently replaced and this generated 916 optimised derivatives of compound 10a.The top 20 ligands are shown in Table 4.
From the generated 916 new ligands, the highest scoring ligand was LD_350 with a docking score of −11.0 kcal/ mol, while ligand LD_753 had the lowest docking score     of −0.3 kcal/mol.The mean was found to be −8.2 kcal/mol, while the median was −8.4 kcal/mol.This indicates that more than 50% of the optimised ligands had higher binding affinity than the original compound 10a at −8.1 kcal/mol.Interestingly, the optimisation process generated 16 ligands with higher binding affinities than the FAD cofactor, as seen in Table 4.In comparison, 664 ligands had higher binding affinities than the original compound 10a.
The interaction diagrams of the top two ligands (LD 350 and LD 762) are shown in Fig. 9, showing the changes, especially in the electronic interactions from the added R groups.In LD 350, for example, all the original electronic interactions in compound 10a were preserved, while the added R group could generate additional H-bond interactions.These are the H-bond between 2-OH of the phenyl ring and C = O backbone of Asp167 at 2.05 Å, the H-bond between 4-OH of the phenyl ring and C = O side chain of Glu168 at 1.55 Å and H-bond between the C = O at the benzyl position and NH side chain of Asn171 at 2.01 Å.These interactions alone with the truncated 4-methoxy group on the DHPM contributed −2.9 kcal/mol of binding energy, which is a 35% increase in binding affinity compared to the original compound 10a.

Molecular dynamics (MD) simulation
Because molecular docking lacks precise information on the explicit biological water system and the inherent flexibility of the receptor, studies of molecular dynamics that imitate similar biological settings can help to understand the stability of the ligand-protein complex.Based on the docking findings, we ran molecular dynamics simulations for compounds 10a, LD_350, and LD_762 in this work.The MD simulation was carried out using the Schrödinger Desmond MD simulation software (version 2021-1).To achieve the biological solvation system, the ligand-protein complex was placed in a TIP3P solvation model and then subsequently neutralised by adding 6Cl − counter ions.The energy minimisation step was also completed to confirm that the system has no steric conflict.The simulation was run for 50 ns at 300 K and 1 bar under an "isothermalisobaric ensemble" (NPT) condition.The "Nose-Hoover chain thermostat" and "Martyna-Tobias-Klein barostat" approaches were ensembled, respectively, for isothermalisobaric conditions [40,41].At 50-ps intervals, simulation trajectories were obtained, and the resulting trajectories were evaluated.The MD simulation utilises protein equilibration, flexibility, and the average distance between backbone atoms to calculate the RMSD to estimate the fluctuation of the total protein-ligand complex.A protein-ligand complex with a lower RMSD value indicates a more stable interaction.In the case of compound 10a, both the protein and the ligand RMSD values lie at approximately 3.0 Å, though the RMSD of the protein averagely lies above 3.0.In contrast, the ligand averagely lies below 3.0 Å except at approximately 35 and 43 ns, where it went slightly above 3.0 Å (See SI Fig. S1a).In addition, the protein and the ligand attained equilibration at less than 5 ns, and no apparent separation between them could be seen throughout the 50-ns MD simulation.Indicating the stability of the original docked structures, the approximately constant RMSD of Cα atoms and minor fluctuations of the ligand atoms during the MD simulation show that the conformation of the ligand did not change or varied very little.Comparing this with the RMSD of the optimised ligands, the LD_350 complex (See SI Fig. S1b), for example, showed lower RMSD values at the beginning of the simulation for the ligand but was only able to achieve equilibration at around 37 ns.Though the RMSD of the ligand increased from ~ 1.8 Å to ~ 4.2 Å during the equilibration, no separation was observed after the equilibration till the end of the simulation period.The RMSD diagram of LD_762 revealed that the protein-ligand complex achieved equilibration at ~ 18 ns, and no separation was observed till the end of the simulation period (See SI Fig. S1c).In this case, the protein showed the slightest deviation and had a maximum RMSD value of 2.4 Å till the end of the simulation, while the RMSD value of the ligand was higher at ~ 5.6 Å at the end of the simulation.
The molecular dynamics simulation identified the major binding interactions between the simulated ligands and glutathione reductase protein.One of the most significant advantages of MD simulation is its ability to precisely identify stable binding interactions from those revealed by molecular docking studies.In most cases, the flexibility of the protein and the ligand in MD, as obtained in a physiological condition, gives rise to the observed differences.In addition, the contact time of each interaction of the ligand with different protein residues as a function of the total simulation period is shown in percentage (See SI Fig. S2).In the case of compound 10a, for example, the two π-cation interactions between + NH 3 of Lys151 and the DHPM rings, as shown in molecular docking, were not preserved.The π-cation interaction of the amino side of Lys32 was preserved at 45% of the simulation time but with the triazole ring rather than the benzene ring of the quinoline moiety.Likewise, the halogen bond interaction between Cl of the quinoline moiety and Asn278 was not preserved.
Interestingly, most of the H-bond interactions involving compound 10a were preserved, and MD was able to identify new H-bond interactions.H-bond interactions involving N3H of DHPM and Asp167 and that of the N of the pyridine ring in quinoline with Ala110 were preserved at 96 and 92%, respectively.In comparison, that of the C = O of the DHPM with Arg272 was lost, but the same C = O was able to form two new H-bond interactions with Thr38 and Asn150 simultaneously at 72 and 94%, respectively.As for the newly identified H-bond interactions, these involved N1H of the DHPM with Gly37 and N3 of the triazole ring with Gly149 at 62 and 32%, respectively.The same trend was observed with the optimised ligands.Analysing LD_350, for instance, the MD did not preserve any of the identified H-bond interactions from the added R-group except for the newly formed intramolecular H-bond interaction between the 2-OH on the benzene ring and its C = O (See SI Fig. S2b).The π-cation interactions between the quinoline rings and Lys32, as well as the DHPM ring and Lys151, are also noteworthy.These were preserved from the molecular docking interactions except for the newly identified π-cation interaction between the pyridine ring of the quinoline moiety and Lys32.Interestingly, the changes in either atomic distances or bond angles beyond the allowed limits resulted in the observed lost interactions in all the cases.
Using the root mean square fluctuation (RMSF) approach, we identified the protein regions exhibiting residue fluctuations across the simulation time.This clarifies how the flexibility of the protein is impacted by ligand binding.As expected, the graphs show the high flexibility of the protein's N and C termini [39].Although there were non-terminal residues with greater RSMF values, they were either extremely close to the terminal or very far from the protein's binding pocket, as demonstrated by Gly67 in the 10a-protein complex, which had an RMSF value of 3.29 Å.All the protein residues interacting with the ligands generally displayed variations at a mean value of 1.0 Å (See SI Fig. S3a, b, and c).

Conclusion
Click chemistry has been successfully used to link DHPM to quinoline-based compounds.The use of click chemistry offers several advantages over the previous amide coupling method.These include a better functional group tolerance, ease of synthesis, lesser waste generation, and better atom economy.From the antimalarial results obtained, we have been able to show that the hybrids containing DHPM and 4-aminoquinoline moieties with triazole linker at position 6 of the DHPM showed good antimalarial activity with IC 50 values below 1 µM and decreased resistance indices than CQ.In general, hybrids with 4-amino functionality displayed better antimalarial activities than those containing oxygen at position 4 of the quinoline ring.The hybrids with mefloquine moiety could be modified in the future to include a 4-amino group on the quinoline ring for comparative antimalarial activities.Apart from position 6 of the DHPM, other positions on the DHPM ring could also be explored to link the quinolines using triazole linkers.Using molecular docking and dynamics approaches, we further explored the synthesised compounds as glutathione reductase inhibitors.Five of the synthesised compounds showed higher binding affinity than the CQ standard, though none showed a higher binding affinity than the FAD cofactor.Rational optimisation of the compound with the highest binding affinity generated 916 ligands, with 664 having higher binding affinities than the original compound and 16 ligands higher than the FAD cofactor.The rational optimisation approach was necessary to produce ligands that could serve as competitive inhibitors of glutathione reductase.

General
A Gallenkamp melting-point apparatus was used for the melting point determination in open capillary tubes and is uncorrected. 1H and 13 C NMR spectra were recorded on either a Bruker Avance 400 (at 400.21 MHz for 1 H and 100.64 MHz for 13 C) or 300 (at 300.13 MHz for 1 H and 75.48 MHz for 13 C) spectrometers using CDCl 3 or DMSOd 6 as solvents at room temperature.Chemical shifts were recorded as part per million (ppm) using tetramethylsilane as an internal standard.2D NMR experiments were recorded on Bruker Avance 400.Mass analysis was performed on Waters® Synapt G2 High Definition Mass Spectrometry (HDMS) system with flow injection analysis (FIA) using electrospray ionisation (ESI) probe.The MS data were acquired and processed on MassLynx™ software (version 4.1).FT-IR measurements were made on a Bruker Alpha Platinum-ATR spectrometer as neat.All reagents and solvents were purchased from Sigma-Aldrich and used without further purification.

X-ray crystallography
Single-crystal diffraction experiments of compounds 3 and 6 were performed using Quazar multi-layer optics monochromated Cu Kα radiation (k = 1.54178Å) on a Bruker D8 Venture kappa geometry diffractometer with duo Iμs sources, a Photon 100 CMOS detector and APEX III control software [42].Data reduction was performed using the SAINT + [42], and the intensities were corrected for absorption using the SADABS [42].Single crystals of 7, 11d, and 11 g were analysed on a Rigaku XtaLAB Synergy R diffractometer with a rotating-anode X-ray source and a HyPix CCD detector.Data reduction and absorption were carried out using the CrysAlis Pro (version 1.171.40.23a) software package [43].All X-ray diffraction measurements were performed at 150(1) K, using an Oxford Cryogenics Cryostat.All structures were solved by direct methods with SHELXT-2016 [44] using the SHELXL-2016 algorithm [45].All H atoms were placed in geometrically idealised positions and constrained to ride on their parent atoms.For data collection and refinement parameters, see the SI (Tables S1, S2).The X-ray crystallographic coordinates for all structures have been deposited at the Cambridge Crystallographic Data Centre (CCDC), with deposition numbers CCDC 2076919-2076923.The data can be obtained free of charge from The Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/ data_ reque st/ cif.

OH
A mixture of 2-trifluoromethyl aniline (15.08 g, 120.00 mmol), ethyl 4,4,4-trifluoroacetoacetate (17.6 mL, 120.00 mmol), and polyphosphoric acid (75.40 g, 5% w/w) was heated with stirring at 150 °C for 3 h, after which the hot reaction mixture was poured into ice water (600.0mL) which precipitated the crude product.The precipitate was suction filtered and air dried.The crude product obtained was purified on column chromatography with SiO 2 gel using 2% methanol in DCM as eluent.

General procedure for the coupling of azido DHPMs to acetylenic quinolines
The general procedure for the click chemistry used for coupling azido DHPHs to acetylenic quinolines is described below using 10 h as an example.9d (0.19 g, 0.60 mmol) and 3 (0.13 g, 0.60 mmol) were dissolved in 3.0 mL DMF.To this mixture was added freshly prepared 1-M aqueous solution of Na ascorbate (240.0 µL, 0.24 mmol) and CuSO 4 (120.0µL.0.12 mmol) in sequential.The reaction mixture was stirred at ambient temperature overnight.The mixture was poured into distilled water (25.0 mL) and extracted with EtOAc (3 × 10.0 mL).The organic layers were combined and washed with brine (10.0 mL), distilled water (5 × 10.0 mL), dried with MgSO 4 , and filtered.filtrate was concentrated in vacuo to give the crude product, which was purified on column chromatography with SiO 2 using 6% MeOH in DCM as the eluent.See the supplementary information for the full characterisations of the final compounds.

Fig. 2
Fig. 2 Molecular structure of compound 6 showing the two tautomers, the keto(a) and the enol (b) forms.Thermal ellipsoids are drawn at a 50% probabilitylevel

Fig. 3
Fig. 3 Molecular structure of compound 7. Thermal ellipsoids are drawn at a50% probability level

Fig. 6
Fig. 6 Binding pose and electronic interactions of compound 10a within theactive site of glutathione reductase.3D and 2D structures are shown in a and b, respectively

Fig. 7 Fig. 8
Fig. 7 Hydrophobically packed H-bond (labelled yellow dotted lines) between Nof the pyridine ring of compound 10a and Ala110 within the hydrophobic regionof glutathione reductase (interacting residues in CPK representation)

Table 1
In vitro activity of DHPM-CQ triazole hybrids against asexual P. falciparum parasites, obtained at concentrations of 1 μM and 5 μM.(CQ was used as a positive control)

Table 2
In vitro IC 50 antiplasmodial activity of compounds 10a-d and compounds 11c, 11e, and 11 g generated using the EPIK tool.Lastly, the Glide ligand docking module and the related default force field OPLS4 were used to execute the molecular docking in an extra precision (XP) mode.

Table 3
Docking score in kcal/mol of the synthesised compounds and FAD cofactor.