Pass and Swiss ADME collaborated in silico docking approach to the synthesis of certain pyrazoline spacer compounds for dihydro- folate reductase inhibition and antimalarial activity

New series of pyrazoline spacer compounds were prepared by the reaction between benzimidazole chalcones and (2-methyl-5-nitro-imidazole-1-yl)acetic acid hydrazide by the sensible use of Michael addition. The building blocks used for the synthesis of pyrazoline derivatives were opted by using virtual screening by molinspiration search engine. The hypothetically resulted pyrazoline spacer compounds from this list are checked for their reliability on other in silico drug designing online web services like PASS online bioactivity, Swiss ADME predictor. The docking study on final four pyrazoline compounds was carried out using Accelrys Discovery Studio 3.5. These synthesized compounds were, later, characterized with the help of UV, IR, mass and 1H NMR techniques. These compounds were further screened for their in vitro antimalarial effect. The PASS, Swiss ADME assisted docking approach and the use of combo heterocyclic ring with pyrazoline scaffold were found to be beneficial to derive and synthesize effective antimalarial agents in the present study. Article Info Received: 19 August 2017 Accepted: 12 December 2017 Available Online: 25 January 2017 DOI: 10.3329/bjp.v13i1.33625 Cite this article: Lohidakshan K, Rajan M, Ganesh A, Paul M, Jerin J. Pass and Swiss ADME collaborated in silico docking approach to the synthesis of certain pyrazoline spacer compounds for dihydrofolate reductase inhibition and antimalarial activity. Bangladesh J Pharmacol. 2018; 13: 23-29. Pass and Swiss ADME collaborated in silico docking approach to the synthesis of certain pyrazoline spacer compounds for dihydrofolate reductase inhibition and antimalarial activity Krishnakumar Lohidashan1, Manju Rajan2, Andhale Ganesh1, Mathew Paul3 and Jithu Jerin1 Department of Pharmaceutical Chemistry, Acharya and B. M. Reddy College of Pharmacy, Bengaluru 560107, Karnataka, India; Department of Pharmaceutical Chemistry, Nirmala College of Health Science, Meloor 680311, Kerala, India; Department of Chemistry, Christ College, Irinjalakuda 680125, Kerala, India. This work is licensed under a Creative Commons Attribution 4.0 License. You are free to copy, distribute and perform the work. You must attribute the work in the manner specified by the author or licensor. the effect of substitutions in the interaction of this molecule to the binding pockets of the enzyme. The title compounds subjected to PASS online bioactivity prediction (Parasuraman, 2011). The result came as synthesized compounds showed more probabilities to exhibit the antiprotozoal activity. Our recent study on dihydrofolate reductase enzyme is a good receptor for the lipophilic derivative gave a direction to do a check on these molecules (Krishnakumar et al., 2017). Malarial dihydrofolate reductase (DHFR) is the target of antifolate antimalarial drugs. Inhibition of this enzyme results in disruption of DNA, RNA, and protein synthesis, with consequent cell death. Recent studies on available antimalarials say efficacy of these drugs has decreased because of mutations in the enzyme that have led to drug resistance (Yuthavong et al., 2012). This fostered us to screen the compounds for effective dihydrofolate reductase inhibition. Further web based in silico methods like SWISS ADME calculator, PASS bioactivity tester are implemented to study the efficacy of these compounds in the due course of synthesis. Drug design involves the estimation of absorption, distribution, metabolism, and excretion (ADME) increased in the contemporary discovery process (Daina et al., 2017). In that context, computer models constitute valid alternatives to experiments. The identified molecules are further subjected to Docking study using CDocker program of Accelrys Discovery Studio 3.5 against the dihydrofolate reductase enzyme. Since the in silico drug design studies are carried out in an intention to minimize the animal sacrifice and economic burden due to clinical studies, the pharmacological screening in this work performed with nonclinical evaluation of the antimalarial activity. Materials and Methods


Introduction
The rule of drug design is more focused into a receptor based design.The structure of the target is usually known to high resolution and the binding site is well defined.Binding of the drug to these hypothetical interaction sites is further influenced by hydrogen bonds, electrostatic and other non-covalent interactions (Wade and Goodford, 1989).The utilization of molecular docking studies made its intervention at this point of drug synthesis.Here comes the scope of utilization of other web based in silico screening as a method of the multilevel screening process.The results obtained from Molinspiration and checking for its match on the principles of Lipinski rule are the one of first evolved of this kind (Mabkhot et al., 2014) used as the choice of present study.The resulted compounds were a fusion of two major heterocyclic entities benzimidazole and imidazole connected with a pyrazoline spacer.The nitrogen heterocycles are well-known about their ability to participate in the diverse biological activities (Akhtar et al., 2017).Considerable attention was given here on the synthesizing molecules, which have three nitrogen heterocyclic rings.Chalcones required for synthesis were prepared from Benzimidazole and aromatic aldehydes.The acid hydrazide was prepared from imidazole.
The pictorial representation of this heterocyclic combination, spacer and the way bulky groups can be linked (Figure 1).The bulky group with substitutions is planned on this pyrazoline spacer in an intention of studying the effect of substitutions in the interaction of this molecule to the binding pockets of the enzyme.The title compounds subjected to PASS online bioactivity prediction (Parasuraman, 2011).The result came as synthesized compounds showed more probabilities to exhibit the antiprotozoal activity.
Our recent study on dihydrofolate reductase enzyme is a good receptor for the lipophilic derivative gave a direction to do a check on these molecules (Krishnakumar et al., 2017).Malarial dihydrofolate reductase (DHFR) is the target of antifolate antimalarial drugs.Inhibition of this enzyme results in disruption of DNA, RNA, and protein synthesis, with consequent cell death.Recent studies on available antimalarials say efficacy of these drugs has decreased because of mutations in the enzyme that have led to drug resistance (Yuthavong et al., 2012).This fostered us to screen the compounds for effective dihydrofolate reductase inhibition.
Further web based in silico methods like SWISS ADME calculator, PASS bioactivity tester are implemented to study the efficacy of these compounds in the due course of synthesis.Drug design involves the estimation of absorption, distribution, metabolism, and excretion (ADME) increased in the contemporary discovery process (Daina et al., 2017).In that context, computer models constitute valid alternatives to experiments.The identified molecules are further subjected to Docking study using CDocker program of Accelrys Discovery Studio 3.5 against the dihydrofolate reductase enzyme.Since the in silico drug design studies are carried out in an intention to minimize the animal sacrifice and economic burden due to clinical studies, the pharmacological screening in this work performed with nonclinical evaluation of the antimalarial activity.

Molinspiration
Various heterocyclic chalcones and acid hydrazide were considered for virtual screening method to check the bioactivity score.When these chalcones were compared, benzimidazole chalcones and imidazole hydrazide were found to have interesting molecular properties, which match with Lipinski prediction.Once a bioactivity model is generated, the actual virtual screening was performed for the predicted pyrazoline molecules.

ADME studies
The ADME study was carried out using SWISS ADME predictor.This is a free web tool to evaluate the pharmacokinetics, drug likeness and medicinal chemistry friendliness small molecules.As mentioned earlier, the attention was given to design the molecules which fit into the rule of drug likeness.The properties like molecular weight less than 500 g/mol, less than 5 numbers of hydrogen bond donors, less than 10 numbers of hydrogen bond acceptors and less than 10 rotatable bonds were chosen as criteria, while the selection of molecules to be synthesized (Doak et al., 2014).The search engine further gave a compiled result on lipophilicity and hydrophilicity of these molecules by integrating results obtained from various Log P and S prediction programs called ILOGP, XLOGP3, WLOGP, ESOL, and SILICOS-IT.Log P, a measure of lipophilicity of the molecule is the logarithm of the ratio of the concentration of drug substance between two solvents in an unionized form.Lipinski rule prescribes an upper limit of 5 for druggable compounds.The lower the log P values the stronger the lipophilicity for the chemical substance.The aqueous solubility of a compound significantly affects its absorption and distribution characteristics.On the other side, low water solubility goes along with a bad absorption, and therefore, the general aim is to avoid poorly soluble compounds.Log S is a unit of expressing solubility in itself, which is the 10-based logarithm of the solubility measured in mol/L.Distribution of Log S in traded drugs reveals a value somewhere between -1 to -4, will be optimum for better absorption and distribution of drugs in the body.

Activity prediction and toxicity
PASS online predicts over 4000 kinds of biological activity, including pharmacological effects been utilized for this study to check the bioactivity score of synthesized compounds as antimalarial drugs.Further preADMET ver.2.0 toxicity predictor predicts test values for Ames test and rodent carcinogenicity assay from the available database of compounds of similar structure.This web service is used for finding mutagenic and carcinogenicity of synthesized compounds.

Docking studies
Molecular docking study for the inhibition of dihydrofolate reductase was studied.A number of inhibitions made by the different derivatives were studied by evaluating the result obtained after docking the com- pounds using CDocker program.The scoring functions and hydrogen bonds formed with the surrounding amino acids were used to predict their level of inhibition on binding sites of dihydrofolate reductase receptor.The structure of 1-[3-(1H-benzoimidazol-2-yl)-4, 5dihydro-pyrazol-1-yl]-2-(2-methyl-5-nitro-imidazol-1yl)-ethanone derivatives were docked using Accelrys Drug Discovery Studio 3.5.The structure of the enzyme dihydrofolate reductase was obtained from Protein Data Bank (PDB code: 7 dfr) and was used for docking.The pose with the most negative CDocker interaction value of the compound was considered to be the most favorable pose having highest inhibition of DHFR.
ethanone derivatives are described here as three schemes.
A mixture of ortho-phenylenediamine (0.25M) and lactic acid (0.35M) was refluxed for 2.5 hours.The reaction mixture was cooled and made alkaline by gradual addition of 10% NaOH.The residue was collected.The crude product obtained was dissolved in 400 mL of boiling water.To this, 2 g of decolorizing carbon was added and heated for 15 min.The mixture was filtered rapidly at the pump through preheated Buchner funnel.The filtrate was cooled to about 10°C.The product obtained was filtered and washed with 25 mL of cold water and dried at 100°C.The purity of the compound was confirmed by getting a single spot in thin layer chromatography (TLC).

Step 2: Synthesis of 2-acetyl benzimidazole
To the solution of 2-(hydroxy ethyl) benzimidazole (50 mM) in diluted H2SO4 (5%, 40 mL), a solution of K2Cr2O7 (150 mM) in aqueous H2SO4 (25%, 80 mL) was added over a period of 20 min.The mixture was then stirred further for 2 hours.The separated orange solid was washed with water, then suspended in water (50 mL) and treated with aqueous NH3 (1:1) to a pH 6.0-6.5.The separated solid was washed with water, dried and recrystallized from boiling ethyl acetate.A single spot in the TLC plate confirmed the purity.

Step 3: Synthesis of benzimidazole chalcones
To a solution of step 2 (0.01 M) and aromatic aldehyde (0.01 M) in ethanol (20 mL) was added into a solution of sodium hydroxide (8 mL, 10% solution).The reaction mixture was stirred at room temperature over a period of 24 hours, diluted with water (100 mL) and acidified with concentrated HCl.The product obtained was filtered, washed with water and recrystallized from ethanol (Scheme 1).

Antimalarial screening
The in vitro antimalarial assay was carried out in 96well microliter plates (Lee et al., 2002).The cultures of Plasmodium falciparum RKL-2 strain were maintained in medium RPMI 1640 supplemented with 25 mM HEPES, 1% D-glucose, 0.2% sodium bicarbonate and 10% heat inactivated human serum.180 mL of 1% D-sorbitol synchronized ring stage parasitemia in 3% hematocrit in a total volume of 200 µL of medium RPMI-1640 was assessed after converting to final dose.A stock solution of 5 mg/mL of each of the test samples was prepared in DMSO, and subsequent dilutions were prepared with culture medium to obtain 5 and 50 mg/mL.The test plates, as well as the positive and negative control (chloroquine and without sample respectively), were incubated at 37°C for 40 hours after treating with 5% CO2.

Pass and Swiss ADME
The results obtained from in silico studies clearly indicate compound code-3 with hydroxy substitution is the most druggable substance with a minimum violation from any of drug likeness rules discussed above.It was interesting to note that the results from the SWISS ADME predictor values of Log P, molar refractivity, and the total polar surface area in these molecules were in excellent agreement with the most important rules of drug likeness.eg.Lipinski, Ghose, Veber, Egan.Though these compounds were exhibiting good hydrophiliclipophilic balance and same predicted bioavailability, the hydroxy derivative with high lipophilicity was expecting to show decent GI absorption.This hydroxy derivative with a higher value of probability of antimalarial activity and non-carcinogenic and mutagenic properties were predicted as the lead in the study.The results obtained from the Swiss ADME search engine are listed in Table I.The PASS activity prediction for pyrazoline compounds given a higher value of probability for antiprotozoal action (Figure 2).The results obtained from activity prediction are given in Table I.

Docking
The extent of inhibition of all the synthesized ligands was compared with inhibition produced by standard drug chloroquine.Efforts were made to study and compare the amino acid residues involved in enzyme inhibition.The results obtained were given in Table II.The four successful synthetic derivatives screened for in silico ADME and toxicity studies when subjected to docking studies the most affluent outcome was again found to match with the results obtained from ADME analysis.The compound code-3 exhibits highest -CDocker energy and -CDocker interaction values comparable with standard drug chloroquine.Unlike all other ligands, the hydroxy terminal substitution exhibited a high number of stronger electrostatic interactions.
The three hydrogen bond interactions associated with Pi bond between benzimidazole ring and amino acid residue in chain A of dihydrofolate seems to be the reason for the significance of compound-3.The compound code-3 is in good correlation with standard drug chloroquine in binding energy and number of hydrogen bond interaction at the drug receptor.The types and number of interactions between the synthesized ligands and enzyme are shown in Figure 2, and the data are tabulated in Table II.Figure 2 exhibits the prominent lipophilicity of compound 3 by the way it fit into the hydrophobic pocket of the receptor.This clearly represents hydrophobic amino acid terminals is located within that region of receptor and facilitates stronger binding with a lipophilic ligand.

Antimalarial activity
The in vitro antimalarial assay results showed that the representative compounds exhibited antimalarial  activity at 50 mg/mL (Table III).The smears when accessed for a number of dead cells, as expected the compounds showed good correlation with docking results.The smears were assessed for live and dead parasites.The significance of computer in the prediction of activity can be substantiated by correlating the activity profile with the in silico data (Table II).The antimalarial activity is compared with the dock scores (binding energy) of the synthesized compounds using GraphPad Prism 7. The correlation coefficient value "r" 0.9929 found too significant for correlating antimalarial activity and predicted docking score.

Discussion
The web-based ADME studies, bioactivity prediction, docking and in vitro studies are carried out for those molecules, which filtered off from drug rules.The authors are successful in deriving a correlation for synthesized compounds among the bioactive parameters.The correlation coefficient value obtained after comparing the docking score and antimalarial activity clearly reflect a positive equation exist between all tested parameters.It was interesting to note the Swiss ADME prediction and PASS bioactivity prediction results also in coherence with final results.This combination of Pass-Swiss-Docking in silico prediction method found to more sensitive and can work with other molecules too.
When it comes to Scaffold chemistry, the docking results are in agreement with the selection of pyrazoline part as a scaffold in the entire molecule.The scaffold is an integral part of the molecule, where a medicinal chemist can try for structural alternation to extract the bioactivity to the molecule in fullest (Mathew et al., 2010).The docking interaction describes none of the cases pyrazoline has much interaction with the amino acid residues in enzyme binding with it.Whereas substitution in pyrazoline with chloro, fluoro, hydroxy and cinnamoyl group found to bring a lot of variations in binding energies as discussed earlier.Assessment of the predictive accuracy of in silico prediction tools, alone or in combination was tried by many researchers (Leong et al., 2015).The uniqueness of this work lies with the prediction of a scaffold in a big molecule and its establishment by using binding configuration obtained after docking study.The lead optimization techniques showed that the compound code-3 with hydroxyl substitution showed good drug likeness and bioactivity score in the initial level of in silico screening.Further, the same molecule showed good binding to dihydrofolate reductase enzyme with stronger electrostatic interaction.In the later stage of work, this in silico results found matching with the results obtained from in vitro antimalarial screening.
It was interesting to note that hydroxy derivative showed good binding to dihydrofolate reductase enzyme with prominent electrostatic interactions.Whereas the other molecule, including the standard drug, showed moderate to weak Van der Walls type of interactions.Further, the optimum hydrophilic-lipophilic balance found to exist in molecule contributes to strong receptor binding.