- Open Access
Central nervous system exposure of next generation quinoline methanols is reduced relative to mefloquine after intravenous dosing in mice
- Geoffrey S Dow1Email author,
- Erin Milner1,
- Ian Bathurst2,
- Jayendra Bhonsle1,
- Diana Caridha1,
- Sean Gardner1,
- Lucia Gerena1,
- Michael Kozar1,
- Charlotte Lanteri1,
- Anne Mannila3,
- William McCalmont1,
- Jay Moon1,
- Kevin D Read4,
- Suzanne Norval4,
- Norma Roncal1,
- David M Shackleford3,
- Jason Sousa1,
- Jessica Steuten3,
- Karen L White3,
- Qiang Zeng1 and
- Susan A Charman3
© Dow et al; licensee BioMed Central Ltd. 2011
- Received: 6 December 2010
- Accepted: 6 June 2011
- Published: 6 June 2011
The clinical use of mefloquine (MQ) has declined due to dose-related neurological events. Next generation quinoline methanols (NGQMs) that do not accumulate in the central nervous system (CNS) to the same extent may have utility. In this study, CNS levels of NGQMs relative to MQ were measured and an early lead chemotype was identified for further optimization.
The plasma and brain levels of MQ and twenty five, 4-position modified NGQMs were determined using LCMS/MS at 5 min, 1, 6 and 24 h after IV administration (5 mg/kg) to male FVB mice. Fraction unbound in brain tissue homogenate was assessed in vitro using equilibrium dialysis and this was then used to calculate brain-unbound concentration from the measured brain total concentration. A five-fold reduction CNS levels relative to mefloquine was considered acceptable. Additional pharmacological properties such as permeability and potency were determined.
The maximum brain (whole/free) concentrations of MQ were 1807/4.9 ng/g. Maximum whole brain concentrations of NGQMs were 23 - 21546 ng/g. Maximum free brain concentrations were 0.5 to 267 ng/g. Seven (28%) and two (8%) compounds exhibited acceptable whole and free brain concentrations, respectively. Optimization of maximum free brain levels, IC90s (as a measure or potency) and residual plasma concentrations at 24 h (as a surrogate for half-life) in the same molecule may be feasible since they were not correlated. Diamine quinoline methanols were the most promising lead compounds.
Reduction of CNS levels of NGQMs relative to mefloquine may be feasible. Optimization of this property together with potency and long half-life may be feasible amongst diamine quinoline methanols.
- Brain Concentration
- Equilibrium Dialysis
- Apparent Permeability
- Intermittent Preventive Treatment
The Walter Reed Army Institute of Research and collaborators are attempting to identify next generation quinoline methanols for intermittent preventive treatment (IPT) of malaria. IPT is the prevention of morbidity or mortality due to malaria through the intermittent administration of a single dose treatment of a drug at full therapeutic doses to asymptomatic, otherwise healthy infants (IPTi), pregnant women (IPTp) and travelers (IPTt) [1–3]. Drugs for IPTx indications and prophylaxis should ideally exhibit a long half-life, be very well-tolerated and safe in pregnancy. Mefloquine exhibits two of these characteristics, but will likely not find use as an IPT drug because of the adverse CNS events observed at the treatment level doses  that may be required for IPT. However, this would presumably not be an issue for next generation analogs of mefloquine without such a liability. Mefloquine accumulates in the CNS and has multiple CNS targets (see discussion in earlier papers [5, 6]). The goal is to identify a lead compound for IPT, based on a mefloquine scaffold, for which accumulation into the CNS is substantially reduced. Such a compound should have an improved CNS safety profile relative to mefloquine.
In an earlier study, it was demonstrated that non-piperidine quinoline methanols, in which the piperidine group of mefloquine was replaced with a diamine side chain, were metabolically stable, exhibited reasonable potency against Plasmodium falciparum in vitro and were less permeable across MDCK cell monolayers than their monamine counterparts [7, 8]. That study did not attempt to address whether reductions in brain concentration relative to mefloquine could be achieved in vivo. This was the goal of the present study, in which approximately 25 compounds from our original library were resynthesized, and brain and plasma concentrations were measured over 24 h in mice after i.v. dosing. Plasma concentrations were measured to generate a preliminary indication of half-life, brain concentrations to assess potential exposure relative to mefloquine, and IC90s to assess intrinsic activity against P. falciparum.
The dose-limiting CNS side effects of mefloquine at the full therapeutic doses required for IPT include dizziness, incoordination, anxiety and sleeplessness . These common side effects are largely absent at the weekly dose of mefloquine, which is five-fold lower than the treatment dose . Therefore, assuming linearity of mefloquine pharmacokinetics in humans, it makes sense that, assuming no change in affinity for the putative CNS receptors of mefloquine, a five-fold reduction in CNS total drug levels would be the minimum requirement to reasonably expect an improvement in the therapeutic index of a NGQM delivering efficacy at blood exposure equivalent to mefloquine. However, as reported elsewhere [7, 8], the lipophilicity of diamine quinoline methanols and other early lead chemotypes is lower than mefloquine. Conceivably this might alter non-specific binding in the brain leading to an increase in the free brain concentration of the drug. Since we do not know the relevant clinical CNS target(s) of mefloquine and the importance of the total and free brain concentration in relation to adverse effects, it is important that reduction relative to mefloquine be assessed with respect to both parameters.
Library synthesis and physicochemical properties
Ranges of physiochemical properties of quinoline methanols investigated in this study.
Animals, dosing and collection of samples
All in vivo studies were conducted in accordance with protocols approved by the Victorian College of Pharmacy, Animal Ethics Committee. Male, six week old FVB mice were obtained from the Animals Resources Center (Perth, Western Australia) and provided with food and water ad libitum. Mice were administered a 50 μL bolus injection of mefloquine, its enantiomers or one of twenty five quinoline methanol compounds by tail vein injection (dosing solution concentration 2.8 mg/ml) for a nominal dose of 5 mg/kg. The compounds were prepared in 5% DMSO in 5% glucose solution buffered to a pH of 3-5 using citrate or acetate as appropriate based on an assessment of the lipophilicity and pKa of the compound. At 5 min, 1 h, 6 h or 24 h post-dosing, mice were anaesthetized with inhaled isoflurane (3%) approximately 3-4 minutes prior to blood and brain harvest. Blood was collected by cardiac puncture from groups of 2-3 mice and the whole brain was removed. The whole brains were placed into preweighed polypropylene vials and snap frozen prior to storage along with the plasma samples at -20°C until analysis. On the day of analysis, the whole brain was homogenized in three parts of aqueous solution containing 0.1 M EDTA and 0.4 g/L KF (inhibitors of hydrolytic enzymes to reduce the potential for ex vivo tissue-mediated degradation of the test compounds). Calibration standards were prepared by spiking blank brain homogenate and blank plasma with the test compound. Both samples and standards were processed by adding acetonitrile (to precipitate proteins) and centrifuged to separate the supernatant for analysis. An aliquot of the supernatant from each brain homogenate and plasma samples were analysed by UPLC/MS/MS on either a Quattro Premier XE, a Quattro Ultima Pt or a Xevo XT Mass Spectrometer (Waters Corporation, USA).
Determination of unbound brain fractions by equilibrium dialysis and free brain levels
where D = dilution factor in brain homogenate and fu,apparent is the measured free fraction of diluted brain homogenate.
Determination of apparent permeability across MDCK cell layers
Permeability was determined by Absorption Systems (Exton, PA). MDR1-MDCK cell monolayers were grown to confluence on collagen-coated, microporous, polycarbonate membranes in 12-well Costar Transwell plates. Data were considered valid for a specific assay plate if TEER values were < 1400 Ω.cm2, the Papp of propanolol was between 10-30 × 10-6 cm/s and the Papp of atenolol was < 0.5 × 10-6 cm/s. The permeability assay buffer was Hanks Balanced Salt Solution containing 10 mM HEPES and 15 mM glucose at a pH of 7.4. A known p-glycoprotein inhibitor, cyclosporin A (CSA), was also added to the assay buffer at 10 μM. This was omitted for amphoteric compounds. Bovine serum albumin (1%) was added to the receiver well. The dosing solution concentrations of the test compounds were 5.0 μM in the assay buffer. All cell monolayers were first pre-incubated for 30 minutes with assay buffer. After 30 minutes, the buffer was removed, replaced with fresh buffer, and time was recorded as 0. The addition of BSA, pre-incubation, and use of a longer incubation time were employed to mitigate potential low recovery or permeability that is sometimes observed for lipophilic or 'sticky' compounds. Cell monolayers were dosed on the apical side (A-to-B) or basolateral side (B-to-A) and incubated at 37°C with 5% CO2 in a humidified incubator. After two hours, aliquots were taken from the receiver chambers. Samples were taken from the donor chamber at 0 and 2 hours. Each determination was performed in duplicate. The lucifer yellow flux was also measured for each monolayer to assess monolayer integrity during the flux period. All samples were assayed by LC/MS/MS using electrospray ionization.
Apparent permeability in the apical (A-B direction), PappA-B, and percent recovery are reported in the supplementary Information. Apparent permeability is a measure of the rate of transport across the cell monolayer. Percent recovery refers to the amount of compound recoverable at the end of the assay. Low recovery may indicate non-specific binding to assay plates, instability or accumulation in the cell pellet. In the case of mefloquine, relatively low recovery is likely a consequence of accumulation in cell membranes [15, 16] rather than non-specific binding.
d Cr/d t is the slope of the cumulative concentration in the receiver compartment versus time in μM s-1.
Vr is the volume of the receiver compartment in cm3.
Vd is the volume of the donor compartment in cm3.
A is the area of the cell monolayer (1.13 cm2 for 12-well Transwell).
C0 is the measured concentration of the donor chamber at time 0 in μM.
CN is the nominal concentration of the dosing solution in μM.
Crfinal is the cumulative receiver concentration in μM at the end of the incubation period.
Cdfinal is the concentration of the donor in μM at the end of the incubation period.
Plasmodium falciparum susceptibility assays
The in vitro activities of quinoline methanols against P. falciparum strains W2, D6, TM91C235, and TM90C2A were evaluated using the traditional labeled hypoxanthine assay of Desjardins et al as modified by Milhous et al. These four P. falciparum strains were selected since they have various levels of resistance to conventional anti-malarials. W2 is chloroquine resistant and mefloquine sensitive, D6 is chloroquine sensitive but naturally less susceptible to mefloquine, TM91C235 is resistant to mefloquine, chloroquine, and pyrimethamine as is TM90C2A, however this latter parasite is a two pfmdr1 copy strain (pfmdr1 amplification has been associated with clinical mefloquine resistance). Mefloquine is routinely screened in these assays to ensure the validity. Historical values for mefloquine against all the strains are reported elsewhere . The focus of this study was the mefloquine-sensitive strain W2, however values for the other strains for all the compounds are reported in the Supplementary Information.
Prism graphing software was used to graph key datasets in order to facilitate visual inspection for possible correlations amongst parameters of interest followed by statistical confirmation if needed. For graphical analysis of IC90 data, IC90s were assumed to be equivalent to the maximum concentrations tested where they were insufficiently potent to be determined. Our brain uptake studies yielded plasma concentration data for four time points. This was insufficient for the purposes of determining half-lives in most instances. As a surrogate for half-life, we determined the % residual plasma concentrations at 24 h versus concentrations at 5 min. Where the measured plasma concentration was below the analytical lower limit of quantitation (LLQ), an arbitrary value of 1 ng/ml was assigned for the purposes of analysis. For analysis of MDCK permeability data versus maximum brain concentration, we did not adjust Papp values on the basis of % recovery. However, percent recovery values are included in the Supplementary Information. The probability of a diamine exhibiting the desired free brain concentration, IC90s and residual plasma concentration was calculated by multiplying the marginal probabilities of each of these events since they appeared to be independent. The +/- 95% confidence limits for the probability estimate was calculated using the formula: +/- 95% CIs = 1.96*SQRT(P*(1-P)/n) where P is probability of success and SQRT is the square root of the result of the function in brackets.
Utility of MDCK permeability screen for ranking maximum brain concentrations of quinoline methanols relative to mefloquine.*
Compounds with Papp < 7.4 × 10-6 cm/s
Compounds with Papp > 7.4 × 10-6 cm/s
Compounds with lower maximum brain concentrations than mefloquine
Compounds with higher maximum brain concentrations than mefloquine
The focus of the WRAIR discovery program with next generation quinoline methanols has been to reduce partitioning into the CNS since the relevant clinical target or targets of mefloquine adverse events are not known . The incidence of the most common adverse CNS effects of mefloquine is reduced substantially when the dose is lowered five-fold from the treatment to the weekly prophylaxis level. Presumably this is also associated with a five-fold decline in CNS levels as well. Consequently, at least a five-fold reduction in whole and free brain levels relative to mefloquine at the efficacious dose is probably the minimum requirement. In this study we have shown that at least a ten-fold reduction in both parameters is feasible through replacement of the piperidine ring of mefloquine with various other four-position substituents. The prospect of such a compound being better tolerated than mefloquine is, therefore, a possibility.
While a reduction in maximum brain concentration should reduce the incidence of CNS associated adverse events, optimizing this property alone will not be sufficient for a new quinoline methanol anti-malarial. The putative late lead compound must also have a sufficiently long half-life and intrinsic efficacy to be useful if administered as a single dose. Therefore, the goal is to balance these three potentially competing pharmacological properties in the same molecule. There is an opportunity to evaluate a subset of the larger 4-position library using surrogates of these endpoints. Maximum free and whole brain concentrations represent surrogates of CNS safety if they are sufficiently reduced relative to mefloquine at the therapeutic dose. Whilst the plasma concentration-time data generated are insufficient to calculate pharmacokinetic parameters directly in most instances, the proportion of drug remaining at 24 h relative to the 5 minute time point is probably a reasonable surrogate of half-life, or at least persistence in plasma. In vitro activity (IC90) can be utilized as a surrogate for intrinsic activity.
There are no compounds amongst this sub-library that meet all the requirements in terms of potency, low brain levels and residual plasma concentrations at 24 h (see Additional file 1) suggesting that further optimization is required. The ideal chemotype around which to anchor a future lead optimization program would be one in which a balance between all these three properties is achievable. The most active compounds were generally those with aliphatic, ether or thio ether side chains (WR308245, WR308246, WR308257, WR308265, WR177000, WR308387, WR308278 and WR308412). However, these compounds exhibited much higher brain concentrations (free and whole) than mefloquine, and in most instances, 24 h residual plasma concentrations were < 1% of those at 5 min, indicating a likely short half-life relative to mefloquine. The exception to this general trend was the diamine WR308396. As might be expected, all compounds with lower whole and free brain concentrations than mefloquine had a greater number of H-bond donors and/or acceptors than mefloquine. These compounds comprised a diverse array of structures including alcohols, acids, triazines, substituted cyclic imidazoles and diamines. However, with the exceptions of the diamines WR319581 and WR318746, these compounds exhibited a lack of potency (IC90 > 300 ng/ml) or low residual plasma concentrations (< 1% at 24 h versus 5 min). Of the 25 compounds that were not mefloquine or its enantiomers, only nine had residual plasma concentrations > 1% at 24 h, and five of these were diamines. Consequently, a balance between these properties may be most likely obtained amongst diamine quinoline methanols.
Identifying a new lead compound with the requisite balance of potency, brain partitioning and half-life might not be feasible if it was true that positive trends in these characteristics were mutually exclusive, however this does not appear to be the case. Maximum free brain concentrations, IC90s, and residual plasma concentrations amongst quinoline methanols did not appear to be correlated (Figure 3A-C). Therefore, at the outset there is no reason to suspect that a balance between these properties is not achievable in a diamine quinoline methanol. There are no generic algorithms that can be used to objectively judge the probability of success of a specific lead optimization programme a priori. Thus one must extrapolate on a project-by-project basis with the available data. The lack of a general correlation between maximum free brain concentrations, residual plasma concentrations and IC90s suggests that their probability distributions may be independent. Thus, we estimated the probability of a diamine quinoline methanol possessing all these properties to be around 4%. This value was obtained by multiplying the marginal probabilities of success for the individual parameters of interest (Figure 5). The confidence intervals are wide (0-21%) which is not surprising given the small sample sizes.
The probability estimate discussed above was based on free brain concentrations being representative of CNS exposure. We have not made the same calculations for maximum whole brain concentrations, because the assumption of independence between this variable and IC90 may arguably be less justifiable. However, two lines of evidence support the notion that a compound with an acceptable maximum free brain concentrations will also have an acceptable maximum whole brain concentrations. First, whole and free brain concentrations in the broader group of quinoline methanols were strongly correlated. This suggests that they are also not independent, or put another way, that there is a strong conditional probability whole brain concentrations will be lower than mefloquine if this is the case for free brain concentrations. Second, in the same broader group of quinoline methanols it is also true that the proportion of compounds exhibiting the desired whole brain concentrations (28%) was higher than the proportion of compounds exhibiting the desired free brain concentrations (8%). Additional studies are required to determine whether these observations will hold for a larger sample of diamine quinoline methanols.
It is appropriate to execute a lead optimization programme if it can be done with a reasonable probability of success in a time frame of up to two years . For planning purposes, it has been assumed the true probability of success for an individual compound lies somewhere between the lower confidence interval and 4%. Since the synthetic route is amenable to the synthesis of 100 new compounds in 12-18 months, it is anticipated that up to four potential late lead molecules will be identified in that time frame. The assumed low probability of success with any individual compound requires an aggressive method of early triage. The permeability of a compound across MDCK cell monolayers appears to correctly categorize it relative to mefloquine in terms of its maximum brain concentration in most instances. It is, therefore, reasonable to use this assay together with more routine P. falciparum susceptibility assays to rapidly prioritize compounds for in vivo studies. It may be necessary to conduct in vivo studies on a selection of active, but permeable compounds triaged using this technique in order to control the type II error rate. A lead optimization campaign conducted broadly along these lines is now underway.
This manuscript was reviewed by the Walter Reed Army Institute of Research and the U.S. Army Medical Research and Materiel Command, and there is no objection to its publication or dissemination. The opinions expressed herein are those of the authors and do not necessarily reflect the views or opinions of the Department of the Army and the Department of Defense. All animal experiments were conducted in compliance with the Animal Welfare Act and other federal statutes and regulations relating to animals and experiments involving animals and adhere to the principles stated in the Guide for the Care and Use of Laboratory Animals (National Academy Press, 1996).
We gratefully acknowledge Medicines for Malaria Venture and the United States Military Infectious Diseases Research Program for financial support. The MDCK permeability assays were performed under contract by Absorption Systems (Exton PA). We thank Meng Shi for statistical advice.
- Dow GS, Magill AJ, Ohrt C: Clinical development of new prophylactic antimalarial drugs after the 5th Amendment to the Declaration of Helsinki. Ther Clin Risk Manag. 2008, 4: 803-819.PubMed CentralPubMedGoogle Scholar
- Schellenberg D, Cisse B, Mendendez C: The IPTi Consortium: research for policy and action. Trends Parasitol. 2006, 22: 296-300. 10.1016/j.pt.2006.05.006.View ArticlePubMedGoogle Scholar
- Shanks GD, Magill AJ, Freedman DO, Keystone JS, Bradley DJ, Steffen R: Drug-free holidays: pre-travel versus during travel malaria chemoprophylaxis. Am J Trop Med Hyg. 2007, 77: 1-2.PubMedGoogle Scholar
- Rendi-Wagner P, Noedl H, Wernsdorfer WH, Weidermann G, Mikolasek A, Kollaristch H: Unexpected frequency, duration and spectrum of adverse events after therapeutic dose of mefloquine in healthy adults. Acta Trop. 2002, 81: 167-173. 10.1016/S0001-706X(01)00210-8.View ArticlePubMedGoogle Scholar
- Caridha A, Yourick D, Cabezas M, Wolf L, Hudson TH, Dow GS: Mefloquine-induced disruption of calcium homeostasis in mammalian cells is similar to that induced by ionomycin. Antimicrob Agent Chemother. 2008, 52: 684-693. 10.1128/AAC.00874-07.View ArticleGoogle Scholar
- Dow G, Bauman R, Caridha D, Cabezas M, Du F, Gomez-Lobo R, Park M, Smith K, Cannard K: Mefloquine induces dose-related neurological effects in a rat model. Antimicrob Agent Chemother. 2006, 50: 1045-1053. 10.1128/AAC.50.3.1045-1053.2006.View ArticleGoogle Scholar
- Milner E, McCalmont W, Bhonsle J, Caridha C, Carroll D, Gardner S, Gerena L, Gettayacamin M, Lanteri C, Luong T, Melendez V, Moon J, Roncal N, Sousa J, Tangtaeng A, Wipf P, Dow GS: Structure-activity relationships amongst 4-position quinoline methanol antimalarials that inhibit the growth of drug sensitive and resistant strains of Plasmodium falciparum. Bioorg Med Chem Lett. 2010, 20: 1347-1351.View ArticlePubMedGoogle Scholar
- Milner E, McCalmont W, Bhonsle J, Caridha D, Cobar J, Gardner S, Gerena L, Goodine D, Lanteri C, Melendez V, Roncal N, Sousa J, Wipf P, Dow G: Anti-malarial activity of a non-piperidine library of next-generation quinoline methanols. Malar J. 2010, 9: 51-10.1186/1475-2875-9-51.PubMed CentralView ArticlePubMedGoogle Scholar
- Rendi-Wagner PH, Noedl H, Wernsdorfer WH, Weidermann G, Mikolasek A, Kollaritsch H: Unexpected frequency, duration and spectrum of adverse events after therapeutic dose of mefloquine in healthy adults. Acta Trop. 2002, 81: 167-173. 10.1016/S0001-706X(01)00210-8.View ArticlePubMedGoogle Scholar
- Overbosch D, Schilthuis H, Bienzle U, Behrens RH, Kain KC, Clarke PD, Toovey S, Knobloch J, Nothdurft HD, Shaw D, Roskell NS, Chulay JD: Atovaquone-proguanil versus mefloquine for malaria prophylaxis in non-immune travelers: results from a randomized, double-blind study. Clin Infect Dis. 2001, 33: 1015-1021. 10.1086/322694.View ArticlePubMedGoogle Scholar
- Lipinski CA: Drug-like properties and the causes of poor solubility and poor permeability. J Pharmacol Toxicol Methods. 2000, 44: 235-249. 10.1016/S1056-8719(00)00107-6.View ArticlePubMedGoogle Scholar
- Kerns EH, Di L: Drug-like properties: concepts, structure design and methods. 2008, Burlington, MA, Elsevier Inc, 35-100.Google Scholar
- Summerfield SG, Read K, Begley DJ, Obradovic T, Hidalgo IJ, Coggon S, Lewis AV, Porter RA, Jeffrey P: Central nervous system drug disposition: The relationship between in situ brain permeability and brain free fraction. JPET. 2007, 322: 205-213. 10.1124/jpet.107.121525.View ArticleGoogle Scholar
- Kalvass JC, Maurer TS: Influence of nonspecific brain and plasma binding on CNS exposure: implications for rational drug discovery. Biopharmaceutics Drug Disp. 2002, 23: 327-338. 10.1002/bdd.325.View ArticleGoogle Scholar
- Go ML, Ngiam TL: Thermodynamics of partitioning of the antimalarial drug mefloquine in phospholipid bilayers and bulk solvents. Chem Pharm Bull. 1997, 45: 2055-2060.View ArticlePubMedGoogle Scholar
- Zidovetzki R, Sherman IW, Attiya A, De Boeck H: A nuclear magnetic resonance study of the interactions of the antimalarials chloroquine, quinacrine, quinine and mefloquine with dipalmitoylphosphatidylcholine bilayers. Mol Biochem Parasitol. 1989, 35: 199-207. 10.1016/0166-6851(89)90206-5.View ArticlePubMedGoogle Scholar
- Desjardins RE, Canfield CJ, Haynes JD, Chulay JD: Quantitative assessment of antimalarial activity in vitro by a semiautomated microdilution technique. Antimicrob Agents Chemother. 1979, 16: 710-718.PubMed CentralView ArticlePubMedGoogle Scholar
- Milhous WK, Weatherly NF, Bowdre JH, Desjardins RE: In vitro activities of and mechanisms of resistance to antifol antimalarial drugs. Antimicrob Agents Chemother. 1985, 27: 525-530.PubMed CentralView ArticlePubMedGoogle Scholar
- Dow GS, Magill AJ: Kucer's Use of Antibiotics. Edited by: Grayson ML, Crowe SM, McCarthy JS, Mills J, Mouton JW, Norrby SR, Paterson DL and Pfaller MA. 2008, ASM Press, Washington DC, 2024-2035. "Mefloquine", SixthGoogle Scholar
- Paul MS, Mytelka DS, Dunwiddie CT, Persinger CC, Munos BH, Lindborg SR, Schacht AL: How to improve R&D productivity: the pharmaceutical industry's grand challenge. Nat Rev Drug Discovery. 2010, 9: 203-214.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.