key: cord-0973420-h4zg0q8x authors: Singh, Kamaljit; Kaur, Hardeep; Chibale, Kelly; Balzarini, Jan title: Synthesis of 4-aminoquinoline–pyrimidine hybrids as potent antimalarials and their mode of action studies date: 2013-06-10 journal: Eur J Med Chem DOI: 10.1016/j.ejmech.2013.05.046 sha: 3a7090a9f803865dbe07b3bbcb0243741e7781ba doc_id: 973420 cord_uid: h4zg0q8x One of the most viable options to tackle the growing resistance to the antimalarial drugs such as artemisinin is to resort to synthetic drugs. The multi-target strategy involving the use of hybrid drugs has shown promise. In line with this, new hybrids of quinoline with pyrimidine have been synthesized and evaluated for their antiplasmodial activity against both CQ(S) and CQ(R) strains of Plasmodium falciparum. These depicted activity in nanomolar range and were found to bind to heme as well as AT rich pUC18 DNA. Malaria is one of the most widespread diseases besides tuberculosis and AIDS which affects more than 500 million people worldwide and results in around 1e3 million causalities every year [1] . In Africa alone, around 20% childhood deaths are due to malaria and a child dies every 30 s [2] and it is estimated that an African child has on an average 1.6e5.4 episodes of malaria fever each year. Of the four typically recognized Plasmodium species causing disease in humans, Plasmodium falciparum is most deadly to children below the age of five leading to mortality while Plasmodium vivax is most morbidity prone, and is responsible for latent infection that hampers current control and future elimination efforts [3] . The development of drug resistance for the common antimalarials such as 4-/ 8-aminoquinolines, 4-methanol quinolines, antifolate drugs, sesquiterpene lactones etc. (Fig. 1) is a rather serious issue which has stimulated considerable research efforts in the development of new drugs using different approaches [4, 5] of which the molecular hybridization approach [6, 7] is quite an attractive strategy which involves design of new chemical entities by covalent fusion of two drugs, both active compounds and/or pharmacophoric units derived from known bioactive molecules with complimentary activities and multiple pharmacological targets. The multiple target strategy led to the design of hybrid of 4-aminoquinoline with species such as triazine [8, 9] , ferrocene [10] , rhodanine [11] , thiazolidine-4-one [12] , chalcone [13] , trioxane [14] , isatin [15] and recently, pyrimidines [16e18] (Fig. 2) . Quinoline containing drugs (chloroquine and primaquine, Fig. 1 ) are known to affect parasite metabolism and cause parasite death by blocking the polymerization of the toxic heme, into an insoluble and non-toxic pigment, hemozoin, resulting in cell lysis and parasite cell auto digestion [19e21] . On the other hand, pyrimidine-based compounds are well known for their wide range of promising antiviral [22] , antitubercular [23] , anti-AIDS [24] , antinociceptive [25] , antifungal [26] , antitumor [27] and antimalarial activities [28] apart from their role in the nucleic acid synthesis. Thus, linking of the quinoline unit with pyrimidine might furnish conjugate hybrids that are capable of showing useful antimalarial activity. Recently, antimalarial activities of some quinolineepyrimidine hybrids with activities in the micromolar to nanomolar range have been reported (Fig. 3) [16e18 ,29] . In yet another report on the evaluation of quinolineepyrimidine, the activity (in micromolar range) was also reported for fixed combinations of the chloroquine and pyrimethamine. In all these reports, the pyrimidines were linked to the quinoline unit through 2-, 4-and 6-positions. We have employed rather conformationally flexible pyrimidine-5-carboxylates linked covalently to 4-aminoquinoline core. These novel pyrimidine carboxylate hybrids interact with the iron center of free heme within the physiological environment (pH 5.6), a key step in the accumulation of heme which is selectively toxic to the parasite. To enhance the possibility to accumulate within the digestive vacuole via weak-base trapping (the mechanism by which CQ and other quinoline antimalarials attain high concentrations inside this compartment), we developed a novel class of antimalarials using a pharmacophore hybridization approach in which the pyrimidine-5-carboxylate motif was hybridized with an iron-complexing, 4-aminoquinoline moiety through C-2 position [29] . To further elaborate the structureeactivity profile, here, we present additional new 4-aminoquinolineepyrimidine carboxylate hybrids. We also report on their antimalarial activity against both CQ sensitive (CQ s Dd2) as well as CQ resistant (CQ R D10) strains. Finally, the mechanism of action studies with the representative compounds has also been performed. The 4-aminoquinolineepyrimidine-5-carboxylate hybrids were synthesized in economical way using synthetic approach outlined in Scheme 1. The key starting compound, 3,4-dihydropyrimidin-2(1H)-one 1 was prepared through NH 4 Cl/TFA [30, 31] catalyzed three-component Biginelli condensation of an alkyl acetoacetate, urea and appropriate aldehyde or its formyl equivalent: 1,3oxazinane derivative, in acetonitrile or under solvent-free reaction conditions, in some cases. Subsequent oxidation of 1 using 60% nitric acid readily furnished pyrimidinones 2 which upon chlorination with POCl 3 yielded 3 [32] . The nucleophilic substitution reaction of 3 with appropriate 4-amino-7-chloroquinoline 4 which in turn was prepared from the commercially available 4,7dichloroquinoline and diaminoalkanes [33] , gave 5aeg in good yields ( Table 1) . Structures of 1e5 were unambiguously established on the basis of spectral ( 1 H NMR, 13 C NMR, MS, FT IR) as well as microanalytical analysis. We have previously established that the 4-aminoquinolinee pyrimidine hybrids 6aec intercepted by a diaminoalkyl spacer showed optimum potency (Table 1) , when the flexible spacer consisted of three or four carbon atoms [29] . Further, the introduction of nitro substituent at ortho position of the phenyl ring at the C-4 of the pyrimidine core furnished the most potent compound 6c with antimalarial activity superior to the standard CQ and close to artesunate [29] . Keeping these observations in mind, we planned to further refine the activity of these persuasive 4aminoquinolineepyrimidine by incorporating electron withdrawing substituents at C-4 phenyl of pyrimidine core, as well as by varying the C-5 ester substituent and also by altering the basicity of diaminoalkyl spacer. The in vitro antimalarial screening of the new synthesized compounds 5aeg revealed good to moderate activities in nM range against both the tested Dd2 (CQ S ) and D10 strains (CQ R ) of P. falciparum (Table 1) . Although the tested hybrids were not as active as the standard drugs viz. CQ and ASN, interesting SARs have been drawn. Analysis of the activity of the compounds recorded in Table 1 reveals that replacing C-5 ethyl ester of the previously reported [29] decrease in antimalarial activity against both the chloroquine sensitive and chloroquine resistant strains of P. falciparum. However, comparison of hybrids 5b, 5c with 6b, 6c having an identical butyl spacer showed that incorporation of isopropyl ester at C-5 of pyrimidine motif (5b and 6b) increased the antimalarial activity against the CQ S strain whereas considerable decrease in activity was observed for CQ R strain of P. falciparum. Also, the most potent compound 5b of the series displayed 2-fold increase in antimalarial activity than the standard CQ against CQ R strain of P. falciparum. When the diaminoalkyl linker of compound 6a was replaced with relatively less basic alkoxy amino linker 5d considerable decrease in antimalarial activity was observed which in turn linked to the decreased accumulation of compound via pH trapping into the digestive vacuole. It was not unexpected since the basicity of alkyl chain linker plays crucial role in determining the antimalarial activity of this class of compounds. Furthermore, the introduction of a nitro substituent on the phenyl ring at the C-4 position of the pyrimidine core to create 5c resulted in a significant increase in antiplasmodial activity in comparison to the p-chloro/p-fluoro substituents (5e and 5f). Moreover, the hybrid 5g lacking a C-4 substituent on the pyrimidine motif led to an increase in antimalarial activity against both the CQ S as well as CQ R strains of P. falciparum. However, although the antimalarial activity of 5g was superior to 5c, 5e & 5f, it was less than the corresponding C-4 phenyl substituted analogs 5b as well as 6b. Thus, the SAR study suggested that both the substitution of the C-4 phenyl group with electron withdrawing groups and alterations in basicity of linker leads to better antimalarial activity. Unfortunately, these compounds suffer from high ClogP values which are in the range 5e8 (Table 1) , which are suggestive of the fact that these possess limited aqueous solubility. However, it is not a serious limitation in view of recent advancements in formulation methods. Compounds 5aeg were evaluated for their toxicity against various (HeLa, Vero, CRFK, HEL and MDCK) cell cultures (Table 1 & SI Table S1 ). Toxicity data revealed that these compounds exhibit high toxicity (low CC 50 ) against MDCK cell cultures while CC 50 values are quiet high for other cell cultures. The CC 50 values for inhibition of MDCK cells summarized in Table 1 indicate that the strongest antimalarial compound 5b was mildly cytotoxic (Table 1) . Further, the ratio of the cytotoxicity (CC 50 in mM) and in vitro antimalarial activity (IC 50 in nM for Dd2 strain) enabled the determination of selectivity index (SI) for these compounds. Compound 5d with alkoxy amino linker and compound 5g bearing C-4 unsubstituted pyrimidine motif exhibited high CC 50 values and thus led to fairly safe selectivity index values (Table 1) . Compound 5d having less basic alkyl spacer, displayed highest SI (43.6) whereas the most potent compound 5b exhibit relatively low SI value (3.92 Table S1 ) and (vii) feline corona virus (FIPV) and feline herpes virus activity in CRFK cell cultures ( Table 2 ). The anti-viral activity of most of the compounds was not impressive except compounds 5a and 5c which exhibited relatively low EC 50 's only against the feline corona virus (FIPV) and feline herpes virus in CRFK cell cultures ( Table 2) . Quinoline antimalarials (e.g., CQ, amodiaquine and quinine) act principally by forming adducts with ferriprotoporphyrin IX, thus blocking haemozoin formation [38] . In this study, we have evaluated the mechanism of antimalarial activity of the most potent compound 5b of the series by studying its binding with heme (Fig. 4) . The titration of monomeric heme was also performed at the Plasmodial food vacuole pH 5.6 using MES buffer instead of HEPES to ensure that the compound 5b binds with heme even at acidic pH (S1). A 1:1 stoichiometry of the most stable complex of 5b with monomeric heme at pH 7.4 and 5.6 was established from the Job's plot (SI Figure S1 ). The association constants (Table 3) were calculated by analyzing the titration curves obtained at pH 7.4 using HypSpecea non-linear least square fitting programme [39] . The binding of CQ with heme under identical conditions was also determined in the similar manner and the results are presented in Table 3 for comparison. Table 3 shows that the association constants for the complexes formed between monomeric heme and 5b (log K 4.96) are comparable with those of standard antimalarial drug, CQ (log K 5.15). Furthermore, the decrease of apparent pH from 7.4 to 5.6 ( Table 3) has little effect on the binding constants indicating binding is stronger even at acidic pH. To further establish the binding of 5b with monomeric heme, 1 H NMR titrations were performed and shifts in the peaks as well as peak intensity noted. The addition of 30 mol% of heme dissolved in 40% DMSO to a solution of 5b in 40% DMSO:D 2 O\D 2 SO 4 (10 ml) caused a shift in the aromatic proton signals (Fig. 5) , indicating binding of 5b with heme but further addition of heme led to broadening of the peaks. An equimolar (3.9 mmol) solution of hemin chloride and 5b when analyzed in mass spectrometer depicted an intense molecular ion peak at 1119.3769 Da (Fig. 6a) , corresponding to the molecular formula C 62 H 62 ClFeN 9 O 6 , suggesting the formation of 1:1 complex. Thus, we propose that 5b interacts with heme by replacing chloride atom of hemin chloride and coordinating the iron atom with its endocyclic quinoline nitrogen as proposed in Fig. 6b . Similar titration of dimers of m-oxo type (10 mM) at pH 5.8 using standard procedure [29] with increasing concentration of compound 5b (0e14 mM), resulted in decrease in intensity of broad peak at 362 nm (Fig. 7a, S1) . Further, Job's plot calculations indicated a 1:1 stoichiometry for the most stable m-oxo: 5b complex (Fig. 7b) . In Table 3 the association constants of 5b (log K 5.72) are compared to that of standard CQ (log K 5.58) and also suggests that the binding of 5b is stronger with m-oxo heme (log K 5.72) than monomeric heme (log K 4.96). Thus, the compound 5b inhibits hemozoin formation by blocking the growing face of heme resulting in the observed antimalarial activity. Furthermore, the b-hematin inhibition assay (SI Table S2 ) shows that there is no correlation between antimalarial activity and b-hematin inhibition and also, all the compounds inhibit b-hematin formation although less than that of standard CQ. The mechanism of many antimalarial drugs such as CQ, quinacrine and quinine relies upon the interaction with DNA presumably through ionic interactions between phosphate groups of DNA and protonated amine in addition to the interactions between aromatic nuclei of the drug with nucleotide bases [40, 41] . Therefore, the DNA binding properties of 4-aminoquinolineepyrimidine hybrids have been evaluated using both the UVevisible spectrophotometer and fluorescence spectrophotometry in order to probe interaction of these compounds with DNA. The addition of CT-DNA (4e200 mM) to the buffered methanolic solution of 5b (30 mM) induced hyperchromic shift of 112% in absorption band at 255 nm whereas hypochromic shift of 37% in the characteristic quinoline ring absorption at 330 nm (Fig. 8) . Also, the bathochromic shift of w3 nm was observed for both the absorptions. The observed hyperchromic as well as hypochromic shifts in absorption bands of 5b upon addition of DNA results from the intercalation of 5b with CT DNA as suggested in the literature [42] . The intercalative nature of interaction of compound 5b with CT DNA was additionally supported by thermal denaturation experiment. Intercalation of molecules into the double helix is known to stabilize the DNA against thermal strand separation and thus increases thermal melting temperature (T m ) [43, 44] . The derivative melting curve presented in Fig. 9 shows an increase of 7.5 C in thermal melting temperature of CT DNA upon addition of 5b which is less than that observed for the CQ (Table S3) . Thus, both the UVevisible titrations and thermal denaturation experiment advocate partial intercalative nature of interactions between compound 5b and CT DNA. Further, to visualize the effect of DNA base composition, the fluorescence titrations of 5b were performed with both GC-rich CT DNA and AT-rich pUC18 DNA in buffered methanol. Fig. 10 shows decrease in the intensity of the emission band of 5b at 380 nm, upon addition of increasing concentration of both the DNAs. A shift of 80 nm in emission band at 380 nm was observed upon addition of CT DNA but no such shift in emission band was observed for pUC18 DNA. Comparison of binding constant of 5b with CT DNA (log K 5.76) and pUC18 DNA (log K 5.73) calculated from titration data using HypSpec [39] , suggest that 5b does not discriminate between GC rich DNA and AT rich DNA. A series of potent 4-aminoquinolineepyrimidine hybrids with antimalarial activity in nanomolar range were reported. The compound 5b exhibits lowest IC 50 value within the series against both CQ S and CQ R strains of P. falciparum. These hybrids displayed mild toxicity against MDCK cell cultures. The antiviral activity profiles of these hybrids indicate that the compound 5a and 5c effectively inhibit feline corona virus and feline herpes virus. Further, the mechanism of observed antimalarial activity was established in terms of binding with heme as well as DNA. All liquid reagents were dried/purified following recommended drying agents and/or distilled over 4 A molecular sieves. THF was dried (Na-benzophenone ketyl) under nitrogen. 1 H NMR (300 MHz) and 13 C (75 MHz) NMR spectra were recorded in CDCl 3 on a multinuclear Jeol FT-AL-300 spectrometer with chemical shifts being reported in parts per million (d) relative to internal tetramethylsilane (TMS, d 0.0, 1 H NMR) or chloroform (CDCl 3 , d77.0, 13 C NMR). Mass spectra were recorded at Department of Chemistry, Guru Nanak Dev University, Amritsar on a Bruker LC-MS MICROTOF II spectrometer. Elemental analysis was performed on FLASH EA 112 (Thermo electron Corporation) analyzer at Department of Chemistry, Guru Nanak Dev University, Amritsar and the results are quoted in %. IR spectra were recorded on Perkin Elmer FTIR-C92035 Fourier transform spectrometer in the range 400e4000 cm À1 using KBr pellets. Melting points were determined in open capillaries and are uncorrected. For monitoring the progress of a reaction and for comparison purpose, thin layer chromatography (TLC) was performed on pre-coated aluminum sheets of Merck (60F 254 , 0.2 mm) using an appropriate solvent system. The chromatograms were visualized under UV light. For column chromatography silica gel (60e120 mesh) was employed and eluents were ethyl acetate/ hexane or ethyl acetate/methanol mixtures. The steady state fluorescence experiments were carried out on Perkin Elmer LS55 fluorescence spectrometer at ambient temperature. UVevisible spectral studies were conducted on Shimadzu 1601 PC spectrophotometer with a quartz cuvette (path length, 1 cm). The absorption spectra have been recorded between 1100 and 200 nm. The cell holder of the spectrophotometer was thermostated at 25 C for consistency in the recordings. To the stirred solution of 3 (2 mmol) and potassium carbonate (5 mmol) in dry THF (30 ml), a solution of appropriate 4aminoquinoline 4 (1.0 mmol) in dry THF (50 ml) was added. The reaction mixture was stirred for 48 h at room temperature. The reaction mixture was filtered and THF was removed under vacuum. The residue was purified by column chromatography using MeOH/ ethyl acetate as eluent to obtain corresponding 5, which was recrystallized from DCM/hexane. Using this procedure the following compounds were isolated. (for optimized structure of 5b, see Figure S2 ). To the stirred solution of appropriate 4-aminoquinoline 4 in dry acetonitrile (50 ml) mixture of 3 (in a 1:2 molar ratio) and potassium carbonate in dry acetonitrile was added. The reaction mixture was refluxed for 24 h and then filtered. Acetonitrile was removed under vacuum and the residue was purified by column chromatography using MeOH/ethyl acetate as eluent to give 5 which is recrystallized from DCM/hexane. 6. Material and methods The test samples were tested in triplicate on one or two separate occasions against chloroquine sensitive (CQ S ) strain of P. falciparum (D10). Continuous in vitro cultures of asexual erythrocyte stages of P. falciparum were maintained using a modified method of Trager and Jensen [45] . Quantitative assessment of antiplasmodial activity in vitro was determined via the parasite lactate dehydrogenase assay using a modified method described by Makler et al. [46] . The test samples were prepared to a 20 mg/ml stock solution in 100% DMSO and sonicated to enhance solubility. Samples were tested as a suspension if not completely dissolved. Stock solutions were stored at À20 C. Further dilutions were prepared on the day of the experiment. Chloroquine (CQ) was used as the reference drug in all experiments. Test samples were initially tested at three concentrations (10 mg/ml, 5 mg/ml and 2.5 mg/ml) to determine the starting concentration for the full doseeresponse assay. CQ was tested at three concentrations namely 30 ng/ml, 15 ng/ml and 7.5 ng/ml. A full doseeresponse was performed for all compounds to determine the concentration inhibiting 50% of parasite growth (IC 50 -value). Test samples were tested at a starting concentration of 10 mg/ml, which was then serially diluted 2-fold in complete medium to give 10 concentrations; with the lowest concentration being 0.02 mg/ml. The same dilution technique was used for all samples. CQ was tested at a starting concentration of 1000 ng/ml. Several compounds were tested at a starting concentration of 1000 ng/ml. The highest concentration of solvent to which the parasites were exposed to had no measurable effect on the parasite viability (data not shown). The IC 50 -values were obtained using a non-linear doseeresponse curve fitting analysis via Graph Pad Prism v.4.0 software. Cytotoxicity was determined by exposing different concentrations of samples to Vero, HEL, HeLa and MDCK cells [29] . The antiviral assays were based on inhibition of virus-induced cytopathicity in HEL [herpes simplex virus type 1 (HSV-1), HSV-2 (G), vaccinia virus, and vesicular stomatitis virus], Vero (parainfluenza-3, reovirus-1, Coxsackie B4, and Punta Toro virus), HeLa (vesicular stomatitis virus, Coxsackie virus B4, and respiratory syncytial virus) and MDCK (influenza A (H1N1; H3N2) and B virus) cell cultures. Confluent cell cultures in microtiter 96-well plates were inoculated with 100 cell culture inhibitory dose-50 (CCID 50 ) of virus (1 CCID 50 being the virus dose to infect 50% of the cell cultures) in the presence of varying concentrations of the test compounds. Viral cytopathicity was recorded as soon as it reached completion in the control virus-infected cell cultures that were not treated with the test compounds [29] . Shrinking the malaria map: progress and prospects Thousands of chemical starting points for antimalarial lead identification Artemisinin resistance in Plasmodium falciparum malaria Malaria chemotherapeutics, part I: history of antimalarial drug development, currently used therapeutics, and drugs in clinical development Hybrid molecules with a dual mode of action: dream or reality? Next-generation antimalarial drugs: hybrid molecules as a new strategy in drug design Synthesis, antimalarial activity and cytotoxicity of 4-aminoquinolineetriazine conjugates Synthesis of 4-aminoquinoline-1,2,3-triazole and 4-aminoquinoline-1,2,3-triazole-1,3,5-triazine hybrids as potential antimalarial agents Synthesis and antimalarial activity in vitro and in vivo of a new ferroceneechloroquine analogue Synthesis and biological evaluation of a new class of 4-aminoquinolineerhodanine hybrid as potent anti-infective agents Synthesis and antimalarial activity of side chain modified 4-aminoquinoline derivatives Enone-and chalconeechloroquinoline hybrid analogues: in silico guided design, synthesis, antiplasmodial activity, in vitro metabolism, and mechanistic studies Preparation and antimalarial activities of trioxaquines, new modular molecules with a trioxane skeleton linked to a 4-aminoquinoline Design, synthesis and anti-plasmodial evaluation in vitro of new 4-aminoquinoline isatin derivatives Novel 4-aminoquinolineepyrimidine based hybrids with improved in vitro and in vivo antimalarial activity Synthesis, characterization and antimalarial activity of quinolineepyrimidine hybrids Substituted quinolinyl chalcones and quinolinyl pyrimidines as a new class of anti-infective agents Falcipain-2 inhibitors Quinoline antimalarial: mechanisms of action and resistance Discovery of dual function acridones as a new antimalarial chemotype Synthesis and antiviral activity of the alpha-analogues of 1,5-anhydrohexitol nucleosides Facile transformation of 3,4-dihydropyrimidin-2(1H)-ones to pyrimidines in vitro evaluation as inhibitor of mycobacterium tuberculosis and modulators of cytostatic activity Optimization of an anti-HIV hairpin ribozyme by in vitro selection Adenosine kinase inhibitors. Synthesis, water solubility, and antinociceptive activity of 5-phenyl-7-(5-deoxybeta-D-ribofuranosyl) pyrrolo[2,3-d]pyrimidines substituted at C4 with glycinamides and related compounds The synthesis of substituted pyridylpyrimidine fungicides using palladium catalysed crosscoupling reactions Nucleoside analogues and nucleobases in cancer treatment Synthesis and antiplasmodial activity of novel 2,4-diaminopyrimidines 2-Aminopyrimidine based 4-aminoquinoline anti-plasmodial agents. Synthesis, biological activity, structureeactivity relationship and mode of action studies Ammonium chloride-catalyzed one-pot synthesis of 3,4-dihydropyrimidin-2-(1H)-ones under solvent-free conditions An expedient protocol of the Biginelli dihydropyrimidine synthesis using carbonyl equivalents On the reaction of 3,4-dihydropyrimidones with nitric acid. Preparation and X-ray structure analysis of a stable nitrolic acid Chloroquine analogues: influence of side chain length and quinolyl nitrogen pK a on activity vs chloroquine resistant malaria Treatment of chronic active hepatitis B (CAH B) with chloroquine: a preliminary report Inhibition of human coronavirus 229E infection in human epithelial lung cells (L132) by chloroquine: involvement of p38 MAPK and ERK Inhibition of human immunodeficiency virus infectivity by chloroquine Mechanism of enhancement of the antiviral action of interferon against herpes simplex virus-1 by chloroquine Malarial haemozoin/betahaematin supports haem polymerization in the absence of protein Investigation of equilibria in solution. Determination of equilibrium constants with the hyperquad suite of programmes Chloroquine as intercalator: a hypothesis revived Intercalating drugs: DNA binding and molecular pharmacology Cytotoxic activity, cell imaging and photocleavage of DNA induced by a Pt(II) cyclophane bearing 1,2 diamino ethane as a terminal ligand Spectrophotometric studies of the interaction of chloroquine with deoxyribonucleic acid Spectroscopic studies on the thermodynamic and thermal denaturation of the CT-DNA binding of methylene blue Human malaria parasite in continuous culture Parasite lactate dehydrogenase as an assay for Plasmodium falciparum drug sensitivity We gratefully acknowledge financial assistance from CSIR, New Delhi (project 01(2364)/10/EMR-II) and UGC, New Delhi for Special Assistance Programme (SAP). H.K. thanks CSIR, New Delhi for senior research fellowship. JB thanks KU Leuven for financial support (GOA 10/14). Supplementary material associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.ejmech. 2013.05.046.