Polytherapy with a combination of three repurposed drugs (PXT3003) down-regulates Pmp22 over-expression and improves myelination, axonal and functional parameters in models of CMT1A neuropathy
- Ilya Chumakov1,
- Aude Milet1,
- Nathalie Cholet1,
- Gwenaël Primas1,
- Aurélie Boucard1,
- Yannick Pereira1,
- Esther Graudens1,
- Jonas Mandel1,
- Julien Laffaire1,
- Julie Foucquier1,
- Fabrice Glibert1,
- Viviane Bertrand1,
- Klaus-Armin Nave2,
- Michael W Sereda2, 3,
- Emmanuel Vial1,
- Mickaël Guedj1,
- Rodolphe Hajj1,
- Serguei Nabirotchkin1 and
- Daniel Cohen1Email author
© Chumakov et al.; licensee BioMed Central. 2014
Received: 9 September 2014
Accepted: 24 November 2014
Published: 10 December 2014
Charcot-Marie-Tooth disease type 1A (CMT1A) is the most common inherited sensory and motor peripheral neuropathy. It is caused by PMP22 overexpression which leads to defects of peripheral myelination, loss of long axons, and progressive impairment then disability. There is no treatment available despite observations that monotherapeutic interventions slow progression in rodent models. We thus hypothesized that a polytherapeutic approach using several drugs, previously approved for other diseases, could be beneficial by simultaneously targeting PMP22 and pathways important for myelination and axonal integrity. A combination of drugs for CMT1A polytherapy was chosen from a group of authorised drugs for unrelated diseases using a systems biology approach, followed by pharmacological safety considerations. Testing and proof of synergism of these drugs were performed in a co-culture model of DRG neurons and Schwann cells derived from a Pmp22 transgenic rat model of CMT1A. Their ability to lower Pmp22 mRNA in Schwann cells relative to house-keeping genes or to a second myelin transcript (Mpz) was assessed in a clonal cell line expressing these genes. Finally in vivo efficacy of the combination was tested in two models: CMT1A transgenic rats, and mice that recover from a nerve crush injury, a model to assess neuroprotection and regeneration. Combination of (RS)-baclofen, naltrexone hydrochloride and D-sorbitol, termed PXT3003, improved myelination in the Pmp22 transgenic co-culture cellular model, and moderately down-regulated Pmp22 mRNA expression in Schwannoma cells. In both in vitro systems, the combination of drugs was revealed to possess synergistic effects, which provided the rationale for in vivo clinical testing of rodent models. In Pmp22 transgenic CMT1A rats, PXT3003 down-regulated the Pmp22 to Mpz mRNA ratio, improved myelination of small fibres, increased nerve conduction and ameliorated the clinical phenotype. PXT3003 also improved axonal regeneration and remyelination in the murine nerve crush model. Based on these observations in preclinical models, a clinical trial of PTX3003 in CMT1A, a neglected orphan disease, is warranted. If the efficacy of PTX3003 is confirmed, rational polytherapy based on novel combinations of existing non-toxic drugs with pleiotropic effects may represent a promising approach for rapid drug development.
KeywordsSystems Biology Repurposing Combination therapy Baclofen Naltrexone Sorbitol Synergy CMT1A Low dose
Charcot-Marie-Tooth disease (ORPHA166), although rare, is the most common hereditary peripheral neuropathy with an estimated prevalence of 1 in 2,500 . Seventy disease genes have now been identified to cause defects primarily in long axons or myelinating Schwann cells . Approximately half of the patients belong to the CMT1A subtype which is caused in the majority of cases by the duplication of a 1.5 megabase region of human chromosome 17 containing the myelin gene PMP22 encoding the peripheral myelin protein of 22 kDa. Overexpression of this gene by 50% results in abnormal Schwann cell differentiation and dysmyelination , eventually leading to axonal loss and muscle wasting. This striking dosage sensitivity points to the important regulatory function of PMP22 protein both for integrity as well as the function of peripheral nerves. Tight reciprocal interactions of neurons and Schwann cells are essential for peripheral nervous system (PNS). Axons provide signalling clues necessary for normal differentiation of Schwann cells, while Schwann cells not only insulate growing axons, but also deliver trophic factors supporting neuronal functions . We hypothesized that polytherapeutic intervention to treat CMT1A should include agents able both to normalise PMP22 gene expression, and to improve axonal dysfunction. As a first step we identified possible targets for intervention, looking at signalling pathways. Then we selected compounds known as to interfere with these pathways, to finally test these compounds in a rat model of the disease.
Systems biology analysis
Genes associated with development of Charcot-Marie-Tooth disease were retrieved from OMIM database (Online Mendelian Inheritance in Man database: http://www.ncbi.nlm.nih.gov/omim).
The initial identification and analysis of the potentially affected CMT1A-relevant signalling pathways, integrating these genes, was made with the help of Ingenuity (https://analysis.ingenuity.com) and GeneGO (https://portal.genego.com) databases/bioinformatics tools and proprietary bioinformatics programs, with subsequent filtering and manual curation of the data. This initial analysis was complemented by manual updating of the most recent literature findings. The choice of the drugs was made through matching the emerging relevant targets using with the help of DrugBank, Reaxys and Integrity databases (http://www.drugbank.ca, http://www.elsevier.com/online-tools/reaxys, https://integrity.thomson-pharma.com). Final selection of candidate drugs for functional validation studies included also analysis of their pharmacological efficacy and safety profiles (https://www.medicinescomplete.com).
(RS)-baclofen (B5399), naltrexone hydrochloride (N3136) and D-sorbitol (S3889) were all obtained from Sigma-Aldrich. PXT3003 consisted in the combination of these 3 drugs.
Co-cultures of sensory neurons and Schwann cells
15 days gestation pregnant Wild Type (WT) Sprague Dawley female rats, bred with CMT1A heterozygous rats , were killed by cervical dislocation and embryos (E15) were removed from the uterus. Rat Dorsal Root Ganglia (DRG) cultures were obtained as previously described - and performed at Neuronexperts laboratories (Marseille, France). The cultures were maintained in standard Neurobasal medium for 7 days to allow Schwann cells to populate sensory neurites. On day 7, the cultures were incubated with standard neuronal medium containing 50 μg/mL ascorbic acid (in order to initiate basal lamina formation and myelination) and drugs until 19 days. Our analyses of myelination were performed after 10–11 days of incubation. Three separate and independent cultures of DRG (from Transgenic (TG) embryos male rats) were performed with 6 replicates per condition.
Cells were fixed by a cold solution of ethanol (95%) and acetic acid (5%) for 10 min, permeabilised, and blocked with PBS containing 0.1% saponin (Sigma) and 10% FBS (Invitrogen) for 15 min. Then, cells were incubated with a rabbit polyclonal antibody specific to Myelin Basic Protein (MBP) (Sigma-Aldrich; M3821). MBP staining was revealed with Alexa Fluor 568 goat anti-rabbit IgG (Molecular Probes). Nuclei were counter-stained with Hoechst (Sigma-Aldrich). 20 pictures per well were taken in the same conditions using the InCell Analyzer™ 1000 (GE Healthcare) with 20 × magnification. Analysis of total length of myelinated axonal segments and total number of nuclei was automatically performed with the InCell Analyzer Software. Total nuclei count showed consistency in all wells and at each time point of the experiment.
RT4 schwannoma cultures
RT4-D6P2T cell line was provided by ATCC (CRL-2768, batch 3993689, mycoplasma-free). Cells were thawed and cultured at Neurofit facilities (Illkirch, France) in DMEM medium (ATCC) containing 10% Foetal Bovine Serum (FBS; ATCC) and 1% antibiotic-antimicotic mixture (Gibco), and maintained in a humidified incubator at 37°C in 5% CO2-95% air atmosphere. Two days later, cells were trypsinised and transferred to 12-well plates (Nunc) at 30,000 cells per well. After 48 h, culture medium was replaced by DMEM containing drugs without FBS. After 8 h of treatment, cells were washed and harvested in 100 μL Phosphate Buffer Saline (PBS), then centrifuged 10 min at 13,000 rpm (4°C). The supernatant was discarded and the pellet was stored at −80°C until use. Three experiments with two independent cultures of RT4 were performed and each condition was done in triplicate.
Animals and housing
CMT1A transgenic rats initial colony was established by a couple of heterozygous transgenic rats belonging to the original source . We used only male animals to reduce potential variability in our experiments . Nevertheless, male rats are clinically indistinguishable from female CMT1A rats . Animals were housed and maintained at Key-Obs (Orléans, France). Two to three animals were housed per type III polycarbonate cages (Techniplast, Italy) under standard conditions at constant temperature (22 ± 1.5°C), hygrometry (50 ± 25%), and lighting conditions (50 Lux in housing room and 13 Lux in experimental room) with a 12/12 h daylight cycle. CD-1 mice were provided by Janvier Labs (Saint-Berthevin, France). Mice were housed at the animal facility of Neurofit (Illkirch, France). They were group-housed (5–10 mice per cage) and maintained in a room with controlled temperature (21–22°C) and a reversed light–dark cycle (12h/12h). Animals had free access to food and tap water. Animal procedures were conducted in strict adherence to the EU Directive of September 22, 2010 (2010/63/UE). Only male animals were included in the experiments.
Rat (E15) embryos
A piece of each E15 embryo head (3 mm3) was used for genotyping. DNA was extracted with the SYBR Green Extract-N-Amp tissue PCR kit (Sigma), then, quantitative PCR (qPCR) was performed using the 7500 fast RT-PCR system and analysed following the manufacturer’s instructions (Applied Biosystems). The gender of each embryo was determined using the male Sry [NM_012772] primers  (F: 5′-GAGAGAGGCACAAGTTGGC-3′; R: 5′-GCCTCCTGGAAAAAGGGCC-3′). The genotyping of each embryo was determined using the Pmp22 [NM_008885] primers . PCR was performed (20 s at 95°C, 45 cycles of 10 s at 95°C, 10 s at 65°C then 30 s at 72°C) with the 7500 fast RT-PCR system (Applied Biosystems).
CMT1A rats (aged 3 weeks)
Total DNA was isolated from ear biopsies by the single-step purification method with DNeasy blood & Tissue kit (Qiagen Gmbh) according to the manufacturer’s protocol. PCR was performed (10 min at 95°C, 45 cycles of 10 s at 95°C, 20 s at 64°C then 20 s at 72°C) with a rapid thermal cycler system (LightCycler® 480 II, 96 wells, Roche) using LightCycler® 480 SYBR Green I Master (Roche). Rat Pmp22 gene DNA [NM_17037.1] was used as a quality control of DNA extraction (F: 5′-GACAAACCCCAGATGGCC-3′, R: 5′-CCGCAGCCACCAGCTATTGGT-3′). Mouse-specific Pmp22 gene [NM_008885.2] PCR primers (F: 5′-GACAAACCCCAGACAGTTGA-3′; R: 5′-CAGGAGCCACCAGCTATTACT-3′) was used to identify transgenic rats. Primers were synthetised by Eurofins (MWG Operon, Germany).
Sciatic nerve crush
Experiments were performed at Neurofit facilities. Sample sizes used for analyses were determined based on previous experience on assay variability. 4–5 week-old CD-1 male mice were anaesthetised using isoflurane (2.5–3% in air). The right thigh was shaved and the sciatic nerve was exposed at mid-thigh level (5 mm proximal to the bifurcation of the sciatic nerve) and crushed for 10 s twice with a microforceps (Holtex) with a 90° rotation between each crush. For sham operated animals, sciatic nerves were exposed but not crushed. Finally, the skin incision was secured with wound clips. Forty two days after the crush, the tibial nerve was taken out from 6 mice per group to perform morphometric analyses.
Stratification of the animals was based on the weight, the parent origin and the muscular performances obtained in behavioural tests (bar test and inclined plane test) performed one week before drug treatment. Randomization of animals was done to avoid bias in animal studies ,. PXT3003 contained 30 μg/kg (RS)-baclofen, 3.5 μg/kg naltrexone hydrochloride and 1.05 mg/kg D-sorbitol. Drugs were dissolved in distilled water and were freshly prepared daily before each per os administration. PXT3003 was administered once daily by gavage with inox steel canula (Dominique Dutscher, 075486) in a volume of 1 mL/kg. At the dates of behavioural testing, gavages were performed after the test. Ages and treatment durations of rats are stated in the behaviour section.
Drugs were dissolved in distilled water and were freshly prepared daily before each per os administration. First administration of treatment (PXT3003: mix of (RS)-baclofen (60 or 600 μg/kg), naltrexone (7 or 70 μg/kg) and D-sorbitol (2.1 or 21 mg/kg)) was performed 30 min after the crush. From day 1 to day 42, administration was performed once daily. On test days, mice were treated 1.5 h before Compound Muscle Action Potential (CMAP) recording.
Pmp22mRNA expression in RT4 cells and in sciatic nerves of CMT1A rats (9 weeks of treatment, 17 weeks of age)
Total RNA was isolated from RT4 cells using the RNeasy Micro Kit (Qiagen Gmbh) as described by the manufacturer’s protocol (Qiagen-RNeasy Micro Handbook), while total RNA was isolated from sciatic nerves using Qiazol (Qiagen Gmbh) followed by the single-step purification method with RNeasy Mini Kit (Qiagen Gmbh) as described by the manufacturer’s protocol (Qiagen-RNeasy Fibrous tissue Handbook). DNA contamination was removed by RNase-free DNase I (Qiagen-RNase-free DNase set 1,500 Kunits). RNA concentrations were estimated by NanoDrop ND-1000. RNA quality control was performed by Agilent RNA 6000 nano chips on Agilent 2100 Bioanalyzer. 80 ng of total RNA was reverse transcribed using SuperScript™ II Reverse Transcriptase with Oligo(dT)12–18 (Invitrogen), then, qPCR was performed with a rapid thermal cycler system (LightCycler® 480 II, Roche). Amplifications were performed (10 min at 95°C, 45 cycles of 10 s at 95°C, 40 s at 60°C then 10 s at 72°C) using LightCycler® 480 SYBR Green I Master (Roche). The sequences of the primers (synthesised by Eurofins) used for the Reverse Transcription quantitative PCR (RTqPCR) analyses were: Rps9 [NM_031108.2]: F: 5′-ATCCGCCAACGTCACATC-3′ and R: 5′-CCGCCACCATAAGGAGAAC-3′, Pmp22 rat and mouse [NM_17037.1] and [NM_008885.2]: F: 5′-TGTACCACATCCGCCTTGG-3′ and R: 5′-GAGCTGGCAGAAGAACAGGAAC-3′, Mpz [NM_017027.1]: F: 5′-TGTTGCTGCTGTTGCTCTTC-3′ and R: 5′-TTGTGAAATTTCCCCTTCTCC-3′, Actb [NM_031144.2]: F: 5′-CACCATGTACCCAGGCATT-3′ and R: 5′-ACTTGCGCTCAGGAGGAG-3′.
Animals were tested in a random and blind manner for treatment and outcome measurements. Behavioural experiments and readouts (bar, inclined plane, hot plate, weight bearing and electrophysiological tests) were performed and validated at Key-Obs and Neurofit facilities by the examiners who were blinded for the treatment. Sample sizes used for behavioural analyses were also determined based on previous experience on assay variability ,.
Bar and inclined plane tests were performed on CMT1A rats after 9 weeks of treatment, 13–17 weeks of age.
Bar test: This test evaluated the muscular strength of the four paws and the equilibrium performances on a fixed rod. The rat was placed on its four paws on the middle of the wooden rod (diameter: 2.5 cm; length: 50 cm). The time spent on the bar (fall latency) in each trial and the number of falls were recorded. Five successive trials were performed (180 s max) .
Inclined plane test: The sliding apparatus  of 30 × 50 cm Plexiglas plane that could be inclined from an angle of 0° (horizontal) up to 60° was used. For each angle, two trials separated by 1 min were performed. Each rat was initially placed on the 25° inclined plane in the up-headed position (head-up orientation). From 30 min up to 1 h later, the same experiment was performed on a 35° inclined plane. Between the tests, rats were returned to their cages. The Plexiglas plane was cleaned after each trial. The performances of rats were evaluated by a 3-level score from 0 (no slide) to 3 (the rat slides to the bottom of the plane). The analysis was performed on the mean score at 25° and 35°.
Hot plate: The animal (after 4 months of treatment, 5 months of age) was placed into a glass cylinder on a hot plate adjusted to 48°C. The latency of the first reaction of the hind paw either left or right was recorded (paw lifting or licking, leaps or a jump to escape the heat). The cut-off time was set to 30 s .
Adverse events were monitored along animal experiments by assessment of the following parameters: loss of weight, lack of appetite, diarrhoea, sneezing, sniffling, laboured breathing, rough hair coat, hair loss, inactivity, abnormal grooming behaviour, aggressiveness.
CD-1 mice, hind limb weight bearing
Animals’ weight distribution on the four limbs was assessed on day 13 using the dynamic weight bearing test (BioSeb, France) . This disability test consists in a continuous measurement of all pressure points applied by a freely moving animal, allowing a quantitative evaluation of the weight imbalance caused by pain. Each mouse was placed in an 11 × 11 × 22 cm cage with a 44 × 44 sensor cells grid on the floor. The pressure applied on sensor cells by animal’s paws was recorded at a 10 Hz frequency over a 5 min period. Pressure and surface detection thresholds were determined automatically for each animal by the Dynamic weight bearing 1.3.2 h software (Bioseb, France). After the manual attribution of each pressure point to the corresponding paw, the mean weight applied on each paw was calculated for surface ratio. Unilateral gait dysfunction was finally evaluated through the ipsi/contralateral hind paw weight and surface ratio.
Electrophysiological studies and readouts were done and validated at Neurofit laboratories. Rats or mice were anaesthetised by 2.5–3% isoflurane-air mixture. CMAP was recorded by needle electrodes placed into the intrinsic foot muscles of the plantar surface (Keypoint electromyography, Medtronic, France). Subcutaneous monopolar needle electrodes (Medtronic, 9013R0312) were used for both stimulation and recording. 12.8 mA square wave pulses of 0.2 ms duration were used to stimulate sciatic nerve. For the CMT1A rats (after 8 months of treatment, 9 months of age), the latencies of CMAPs elicited by stimulation at both proximal (hip) and distal (hock) sites were measured. Motor Nerve Conduction Velocity (MNCV) was estimated by the time interval between two stimulation sites latencies relative to distance between these sites (with leg fully extended) ,. For the CD-1 crushed mice (42 days of treatment), the right sciatic nerve (ipsilateral) was stimulated with single pulse applied at the sciatic notch. CMAP was recorded by needle electrodes placed at the gastrocnemius muscle. The amplitude of the action potential was determined 30 days after the crush.
Mice (42 days of treatment) or CMT1A rats (after 4 months of treatment, 5 months of age) were anaesthetised with intraperitoneal mixture of ketamine and xylazine (80 mg/kg, and 15 mg/kg, i.p. respectively). 5 mm segment of the right sciatic nerve of rats (just above the bifurcation of sciatic nerve into two branches, in contrast to other studies collecting sciatic nerves at their distal (1/3) ends ), or tibial nerve of mice were excised. Tissues were fixed in 4% glutaraldehyde in PBS (Sigma-Aldrich) at 4°C overnight. Then, nerve samples were rinsed with PBS at room temperature and post-fixed in 1% osmium tetroxide (EMS, France) for 1 h, rinsed in PBS and dehydrated in serial alcohol solutions (VWR Prolabo), and embedded in Epon. Embedded tissues were kept at 70°C during the 3 days of polymerisation. The tissue block was cut with a microtome (Microm, HM360). Transverse sections 5 μm thick were performed. Then, 1.5 μm cross section was performed and stained with 1% of toluidine blue (EMS) for 2 min, rinsed in PBS, dehydrated in serial alcohol solutions (95% and 100%) followed by xylene solution substitute (EMS) for 10 min at room temperature, and mounted in Eukitt (EMS).
Composite image of the entire nerve section of each sample was obtained using an optical microscope (Nikon optiphot-2) equipped with a digital camera (Nikon DS-Fi1). Morphometric analyses were performed with digital images using Image-Pro Plus software (Media Cybernetics Inc., USA), which computes the axonal and myelin pixel sizes after individualisation of each myelinated fibre via the grey level of toluidine blue-stained myelin sheath. These two parameters were used to calculate the number of myelinated axons and the g-ratio (ratio of the inner axonal diameter to the total outer).
Statistical tests were two-tailed and conducted at a 5% significance level. Sample sizes used for behavioural analyses were defined on the basis of previous experience on assay variability ,. Data distribution and within-group variation were preliminary assessed in order to guide our methodological choices. Statistical analyses were performed with Prism (http://www.graphpad.com/scientific-software/prism) and R (http://cran.r-project.org). We applied an Analysis of Variance (ANOVA) with Dunnett’s test for comparison of more than one experimental group to a reference and a Welch’s t-test for comparison of two experimental groups. These tests are described to be robust enough to overcome changes in distribution and variation ,. We performed a combined analysis of all cultures and replicates of in vitro experimental data to assess myelination in DRG co-cultures after 10–11 days of treatment, and to study gene expression in RT4 schwannoma cells. In each in vitro experimental system (DRG and RT4), data were normalised to the non-treated vehicle control. In DRG co-cultures, synergy was assessed by drug combination analysis with isobologram and calculation of Combination Index (CI) that compares dosages for a given combination to those expected to obtain the same combination effect under a simple additive assumption, and offers a quantitative definition for additive effect (CI = 1), synergism (CI < 1), or antagonism (CI > 1) -. To study the superiority of PXT3003 to its single drugs in RT4 cells, an ANOVA model was fitted and the difference between the effect of PXT3003 and the effect of the most efficient drug (i.e. the drug having the largest effect when compared to control) was tested using a post-hoc test of contrast. Kaplan-Meier curves were analysed with a Cox model. Analyses were adjusted to the age of starting treatment when appropriate (bar test and inclined plane test). Statistical analysis of linear regression of g-ratio with respect to the axon diameter was performed on slope.
Systems biology analysis and drug choice
The first drug, baclofen (BCL), is a specific agonist of GABAB receptors. These receptors decrease the activity of adenylate cyclases and therefore reduce the levels of intracellular cAMP that positively regulate PMP22 expression ,,-. GABAB receptors are also abundant in neurons . BCL is a safe drug currently used to treat spasticity. The second selected drug is the opioid receptor antagonist naltrexone (NTX). It is approved for treating drug addiction, and was shown at low non-toxic doses, to potentiate, rather than to block, cell signalling through opioid receptors that are coupled to inhibitory G alpha protein subunit (Gαi), thereby reducing intracellular levels of cAMP. NTX could mediate this action either through interaction with accessory protein Filamin A (FLNA)  or by induction of endogenous opioid agonists . The third compound D-sorbitol (SRB), a natural metabolite playing an important role in the energy production/storage (polyol pathway), was chosen as another safe drug involved in processes that are deregulated in CMT1A. No specific receptor for SRB is currently known although muscarinic acetylcholine GPCRs, present both on neurons and Schwann cells and affecting cAMP levels, were reported to bind SRB with unexpected high affinity . Moreover, SRB acting as chaperone might also improve PMP22 protein folding , that is impaired when overexpressed in CMT1A Schwann cells ,. Importantly, we hypothesised that a combination of these drugs (PXT3003), by acting on different receptor systems, could decrease the toxic effects of the overexpression of PMP22 gene and also improve downstream consequences on myelination and nerve function. This hypothesis was then tested in a battery of ex vivo and in vivo tests modelling different aspects of CMT1A disease.
PXT3003 improves myelination ex vivo
PXT3003 modulates Pmp22 mRNA levels in vitro
PXT3003 down-regulates Pmp22 expression and is active in the in vivotransgenic rat model of CMT1A disease
Efficacy of 9 weeks of oral daily PXT3003 treatment given to CMT1A TG rats was assessed by (i) measuring muscle strength in the bar test, (ii) performance on an inclined plane and (iii) thermal sensitivity in the hot plate test. We observed a substantial heterogeneity of TG rats with respect to muscular strength, as previously described ,. After 9 weeks of treatment, the overall mean latency to fall (bar test) was significantly prolonged (Figure 4B). The same data set represented as Kaplan-Meier curves allowed us to confirm the significant improvement of PXT3003-treated animals (Figure 4C). It is important to note that the CMT1A-like phenotype was already manifested before the start of the treatment (Additional file 2A). We also assessed motor performance with the inclined plane test  that measures the ability of animals to rest on the plane fixed at two different angles. Before treatment, there was a substantial and significant difference between TG vehicle and WT groups in this test (Additional file 2B). After 9 weeks of treatment with PXT3003, CMT1A rats showed significantly improved performance on the inclined plane that even reached the level of performance of non-affected WT rats (Figure 4D).
CMT1A neuropathy is also characterised by sensory impairment , which we measured in the CMT1A rats using the hot plate . We confirmed the reduced heat sensitivity reported previously in this model . The daily treatment with PXT3003 over 4 months restored sensitivity of CMT1A rats almost to WT levels (Figure 4E).
PXT3003 improves histological and electrophysiological parameters in CMT1A rats
One consequence of dysmyelination in CMT1A disease is an impairment of electrophysiological parameters such as decrease of MNCV or of CMAP amplitude ,. After 8 months of PXT3003 treatment in CMT1A rats, MNCV was improved (Figure 5C). However, CMAP amplitudes had only a slight statistical trend for improvement (Additional file 3B). All above treatment effects and durations in CMT1A rats are summarised in Additional file 4.
PXT3003 improves axonal myelination and regeneration in a nerve crush model
Next, we tested the potency of the PXT3003 combination compared to single drugs in the functional behavioural assay of paw surface bearing, where the surface generated by paw pressure against the floor was measured . In operated mice, unable to set down properly their affected paw to perform normal movements, a substantial decrease in the surface bearing of the affected hind paw was confirmed (Additional file 5E). After treatment of mice with PXT3003, using doses used in the experiments described above, no functional improvement was observed. However, when we increased drug doses to a level that was still less than the range of clinical use for their principal indications, i.e. 1/27 of the approved dosage for BCL (80 mg/day), 1/147 of approved NTX (50 mg/day) and 1/147 of approved SRB (15 g/day), PXT3003 was able to normalise completely this functional defect (Figure 6E) in contrast to single drugs that had no detectable effect (Figure 6F to H). These results highlight a neuro-regenerative and promyelinating potential of our drug combination.
We found no evidence of side-effects in any of the animal experiments described above (see Methods).
It is striking that a small change in the copy number of PMP22 (three compared to two) is responsible for axonal loss in peripheral nerves, and consecutive muscle atrophy and impaired sensitivity in affected individuals. Functional integrity of myelin sheaths is highly sensitive to the stoichiometry of its components, and any dysfunction in the transcriptional network regulating expression of myelin genes is potentially deleterious. Systems biology analysis of available experimental data allowed us to prioritise an integral network of cellular signalling pathways implicated in transcriptional regulation of the PMP22 gene. The converging point of this network is cAMP-dependent silencer upstream of the PMP22 promoter. Three drugs: baclofen, naltrexone and sorbitol, all binding to GPCR receptors and affecting intracellular cAMP levels, were chosen among available authorised drugs with additional criteria of safety and ability to affect other aspects of nerve function. The underlying idea was that their combined action will be more potent which could permit their use at lower doses further minimising potential unwanted secondary effects. Other myelin genes can be potentially down-regulated by these drugs, which necessitated an experimental verification of our underlying hypothesis.
We performed this primary analysis by testing the three molecules in a disease-relevant ex vivo myelination model derived from Pmp22 transgenic rats. While separately the three molecules were able to increase myelination, their common effect demonstrated a positive synergistic interaction that provided a rationale for further combination assays. In pilot experiments using cultured Schwannoma cells, we were able to detect a limited but significant down-regulation of Pmp22 gene expression. The RT4-D6P2T Schwannoma cell line that we used displayed features that resemble those of normal myelinating Schwann cells with the expression of myelin genes at relatively high levels ,. The extent of the down-regulation observed was significantly higher for the combination compared to the action of individual drugs. Expression of the major myelin gene Mpz was not affected by these treatments, while Pmp22 gene expression (when normalised to that of Mpz) was decreased by the presence of PXT3003 but not by its constitutive drugs. Interestingly, this specific effect on Pmp22 but not on Mpz expression is similar to the effect observed for the progesterone antagonist onapristone . We note that CMT1A is a gene dosage disease where translation, processing and transport of PMP22 protein in Schwann cells is not yet fully elucidated and only a small fraction of the protein is incorporated into mature myelin . Moreover, PMP22 protein and mRNA expression in CMT1A patients is variable and no correlation was found between PMP22 overexpression and the Charcot-Marie-Tooth disease neuropathy score (CMTNS) clinical score ,,. Future studies may shed light on the effect of PXT3003 on the regulation of PMP22 protein and other myelin proteins in peripheral tissue of CMT1A patients.
Since PXT3003 proved efficacy over its single compounds in the cellular models, we further tested the combination drug in the transgenic rat model of CMT1A. As patients included in clinical trials could have manifested CMT1A for years before the beginning of the trial, we used CMT1A adult rats for PXT3003 treatment. We also asked whether electrophysiological amelioration could be expected after a long time of treatment, i.e. 8 months, when the disease has reached an advanced stage in the animal model. Here, we demonstrated the beneficial effect of PXT3003 on the clinically relevant muscle motor performance and sensory nerve function. Several biochemical, electrophysiological and morphological characteristics such as the ratio of Pmp22/Mpz gene expression in sciatic nerves, small fibre myelination and MNCV were significantly improved by PXT3003. Of note, the mild reduction in Pmp22 expression after PXT3003 treatment (15%) is reminiscent of the therapeutic effect of onapristone in CMT1A rats . We here argue that the reduction of Pmp22 by PXT3003 improves peripheral nerve myelination and/or trophic support to axons leading to the increase of the number of myelinated axons in the sciatic nerve by 31%. At the same time, PXT3003 combination had limited effect on CMAP and myelin g-ratio in this CMT1A TG rat model in contrast to the acute in vivo nerve crush model. Further work will be needed to see if this increase in NCV without changes in myelin g-ratio may be due to the changes in internodal lengths. We note that the speed of nerve impulse propagation is also determined by the internodal length, which could explain ameliorated NCV and unaltered CMAP. In addition, we hypothesise that the absence of the effect on CMAP at this moment may be due to the fact that the treatment of TG rats had been started late, in early adulthood when the disease was already at an advanced, chronic stage. We note that in the progesterone antagonist therapy trial, a stronger therapeutic effect was observed when treatment was initiated early postnatally ,. Also, early postnatal treatment of CMT1A rats with soluble NRG corrected the dysmyelinating and dedifferentiated Schwann cell phenotype . Experiments with an early onset treatment paradigm with PXT3003 in CMT1A rats are in progress and will also focus on its effects on the differentiation state of Pmp22 transgenic Schwann cells. Finally, we note that in the nerve crush axonal regeneration model, PXT3003 administrated very early, i.e. 30 minutes after nerve crush, accelerated the process of axonal regeneration thus leading to the significantly decreased proportion of small myelinated axons and causing improvement in the amplitude of CMAP. These observations might also raise the possibility that PXT3003 could modify directly axons activity enabling their protection or regeneration, in agreement with our initial hypothesis (Figure 1B).
The synergy observed in models permits the use of individual drugs at lower doses up to several hundred times less than what is being used in their primary indications, diminishing potential undesirable secondary effects.
Throughout this study, we focused on the validation of the efficacy of PXT3003 in cellular and animal models by assessing several CMT1A disease-relevant endpoints. Our goal was more on the translational side to provide the rationale for further testing in an affected human population. The major point of this work was to provide a proof of principle for the repurposing of approved and clinically safe drugs, using reproducible preclinical systems with the main accent to genetic disease model. Many features of CMT1A disease are far from being understood, neither in rats nor in human patients. It is important to keep in mind that, although for some endpoints the effect in some models was significant but rather limited, it can be considered as successful test of the hypothesis that approved drugs can be newly formulated to ameliorate a disease that they were not designed for.
Although designed from a series of mechanistic hypotheses, our preclinical data remains for the moment mostly observational. More efforts must be performed for a better understanding of exact mechanism of action and the combinational synergy observed for PXT3003. It will be particularly relevant to verify experimentally the involvement of GABAB, opioid or muscarinic receptors in the therapeutic effect reported here. Investigating the mechanism of synergy could be even more challenging. Launching these sophisticated and laborious mechanistic experiments was out of scope of this study also because of limited correlations generally observed for efficacy when translating data in animal models to the human clinical studies . Finally, as the profile of PXT3003 seemed rather safe, taking into account the nature and low doses of the 3 compounds, we decided to concentrate our efforts on its testing in an exploratory phase 2. This has recently provided first indications of therapeutic activity in mild to moderately affected adult CMT1A patients . Altogether these data indicate that PXT3003 deserves further clinical investigation as well as experimental elucidation of its detailed mechanism of action.
In conclusion, the novel combination of three well-known and approved drugs (baclofen, naltrexone and sorbitol), identified by our Systems Biology approach, is able to improve myelination in the ex vivo myelination model of CMT1A, while single drugs displayed a lower efficacy. The combination is also able to down-regulate the expression of Pmp22 in cultured Schwannoma cells more efficiently than single drugs. The tests performed in vivo to assess the efficacy of the combination in the nerve crush mouse model demonstrate the acute neuro-regenerative and promyelinating potential of PXT3003 with the functional test demonstrating superiority of the combined action. Moreover, PXT3003 combination improves relevant parameters in the CMT1A rat model (muscular performances, heat sensitivity, histology and electrophysiology). No toxic effect of PXT3003 has been detected in animal studies.
Rational polytherapy based on a combinational repositioning of existing non-toxic drugs that act on pleiotropic pathways, may represent an important novel approach for rapid drug development in a variety of disorders. Future work will show if a similar approach could be useful in the treatment of other neglected orphan diseases where therapies are urgently needed.
We would like to thank Prof Eric R. Kandel for his critical reviewing and editorial suggestions on the manuscript. We thank Dr. Oxana Guerassimenko for her input in the preliminary experiments. We thank Dr. Noëlle Callizot and Dr Rémy Steinschneider for the scientific advice on in vitro studies. We would also like to thank Dr. Fabrice Trovero, Dr. Jean-Charles Bizot, Dr. Fabienne Massé, Sabrina David and Nicolas Cogrel for the scientific and technical support on the behavioural tests performed in the CMT1A rat model. We thank Dr. Stéphanie Wagner and Dr. Emile Andriambeloson for their scientific advice on electrophysiology and on the nerve crush model. We are also grateful to Ms Linda Belo for reviewing the English of the manuscript.
- Patzkó A, Shy ME: Update on Charcot-Marie-Tooth disease. Curr Neurol Neurosci Rep. 2011, 11: 78-88. 10.1007/s11910-010-0158-7.View ArticlePubMedPubMed Central
- Rossor AM, Polke JM, Houlden H, Reilly MM: Clinical implications of genetic advances in Charcot-Marie-Tooth disease. Nat Rev Neurol. 2013, 9: 562-571. 10.1038/nrneurol.2013.179.View ArticlePubMed
- Fledrich R, Stassart RM, Klink A, Rasch LM, Prukop T, Haag L, Czesnik D, Kungl T, Abdelaal TAM, Keric N, Stadelmann C, Brück W, Nave K-A, Sereda MW: Soluble neuregulin-1 modulates disease pathogenesis in rodent models of Charcot-Marie-Tooth disease 1A. Nat Med. 2014, 20: 1055-1061. 10.1038/nm.3664.View ArticlePubMed
- Nave K-A, Trapp BD: Axon-glial signaling and the glial support of axon function. Annu Rev Neurosci. 2008, 31: 535-561. 10.1146/annurev.neuro.30.051606.094309.View ArticlePubMed
- Sereda M, Griffiths I, Pühlhofer A, Stewart H, Rossner MJ, Zimmerman F, Magyar JP, Schneider A, Hund E, Meinck HM, Suter U, Nave KA: A transgenic rat model of Charcot-Marie-Tooth disease. Neuron. 1996, 16: 1049-1060. 10.1016/S0896-6273(00)80128-2.View ArticlePubMed
- Cosgaya JM, Chan JR, Shooter EM: The neurotrophin receptor p75NTR as a positive modulator of myelination. Science. 2002, 298: 1245-1248. 10.1126/science.1076595.View ArticlePubMed
- Rangaraju S, Madorsky I, Pileggi JG, Kamal A, Notterpek L: Pharmacological induction of the heat shock response improves myelination in a neuropathic model. Neurobiol Dis. 2008, 32: 105-115. 10.1016/j.nbd.2008.06.015.View ArticlePubMedPubMed Central
- Callizot N, Combes M, Steinschneider R, Poindron P: A new long term in vitro model of myelination. Exp Cell Res. 2011, 317: 2374-2383. 10.1016/j.yexcr.2011.07.002.View ArticlePubMed
- Sereda MW, Meyer zu Hörste G, Suter U, Uzma N, Nave K-A: Therapeutic administration of progesterone antagonist in a model of Charcot-Marie-Tooth disease (CMT-1A). Nat Med. 2003, 9: 1533-1537. 10.1038/nm957.View ArticlePubMed
- Tashiro H, Fukuda Y, Kimura A, Hoshino S, Ito H, Dohi K: Assessment of microchimerism in rat liver transplantation by polymerase chain reaction. Hepatology. 1996, 23: 828-834. 10.1002/hep.510230425.View ArticlePubMed
- Hirst JA, Howick J, Aronson JK, Roberts N, Perera R, Koshiaris C, Heneghan C: The need for randomization in animal trials: an overview of systematic reviews. PLoS One. 2014, 9: e98856-10.1371/journal.pone.0098856.View ArticlePubMedPubMed Central
- Kilkenny C, Browne WJ, Cuthill IC, Emerson M, Altman DG: Improving bioscience research reporting: the ARRIVE guidelines for reporting animal research. J Pharmacol Pharmacother. 2010, 1: 94-99. 10.4103/0976-500X.72351.View ArticlePubMedPubMed Central
- Meyer Zu Horste G, Prukop T, Liebetanz D, Mobius W, Nave K-A, Sereda MW: Antiprogesterone therapy uncouples axonal loss from demyelination in a transgenic rat model of CMT1A neuropathy. Ann Neurol. 2007, 61: 61-72. 10.1002/ana.21026.View ArticlePubMed
- Rivlin AS, Tator CH: Objective clinical assessment of motor function after experimental spinal cord injury in the rat. J Neurosurg. 1977, 47: 577-581. 10.3171/jns.1977.47.4.0577.View ArticlePubMed
- Beyreuther BK, Callizot N, Brot MD, Feldman R, Bain SC, Stöhr T: Antinociceptive efficacy of lacosamide in rat models for tumor- and chemotherapy-induced cancer pain. Eur J Pharmacol. 2007, 565: 98-104. 10.1016/j.ejphar.2007.02.041.View ArticlePubMed
- Tétreault P, Dansereau M-A, Doré-Savard L, Beaudet N, Sarret P: Weight bearing evaluation in inflammatory, neuropathic and cancer chronic pain in freely moving rats. Physiol Behav. 2011, 104: 495-502. 10.1016/j.physbeh.2011.05.015.View ArticlePubMed
- Bordet T, Buisson B, Michaud M, Drouot C, Gale P, Delaage P, Akentieva NP, Evers AS, Covey DF, Ostuni MA, Lacape J, Massaad C, Schumacher M, Steidl E, Maux D, Delaage M, Henderson CE, Pruss RM: Identification and characterization of cholest-4-en-3-one, oxime (TRO19622), a novel drug candidate for amyotrophic lateral sclerosis. J Pharmacol Exp Ther. 2007, 322: 709-720. 10.1124/jpet.107.123000.View ArticlePubMed
- Markowski C, Markowski E: Conditions for the effectiveness of a preliminary test of variance. Am Stat. 1990, 44: 322-326.
- Sawilowsky SS, Blair RC: A more realistic look at the robustness and Type II error properties of the t test to departures from population normality. Psychol Bull. 1992, 111: 352-360. 10.1037/0033-2909.111.2.352.View Article
- Grabovsky Y, Tallarida RJ: Isobolographic analysis for combinations of a full and partial agonist: curved isoboles. J Pharmacol Exp Ther. 2004, 310: 981-986. 10.1124/jpet.104.067264.View ArticlePubMed
- Chou T-C: Theoretical basis, experimental design, and computerized simulation of synergism and antagonism in drug combination studies. Pharmacol Rev. 2006, 58: 621-681. 10.1124/pr.58.3.10.View ArticlePubMed
- Geary N: Understanding synergy. Am J Physiol Metab. 2013, 304: E237-E253.
- Tallarida RJ: An overview of drug combination analysis with isobolograms. J Pharmacol Exp Ther. 2006, 319: 1-7. 10.1124/jpet.106.104117.View ArticlePubMed
- Suter U, Snipes GJ, Schoener-Scott R, Welcher AA, Pareek S, Lupski JR, Murphy RA, Shooter EM, Patel PI: Regulation of tissue-specific expression of alternative peripheral myelin protein-22 (PMP22) gene transcripts by two promoters. J Biol Chem. 1994, 269: 25795-25808.PubMed
- Sabéran-Djoneidi D, Sanguedolce V, Assouline Z, Lévy N, Passage E, Fontés M: Molecular dissection of the Schwann cell specific promoter of the PMP22 gene. Gene. 2000, 248: 223-231. 10.1016/S0378-1119(00)00116-5.View ArticlePubMed
- Ogata T, Iijima S, Hoshikawa S, Miura T, Yamamoto S, Oda H, Nakamura K, Tanaka S: Opposing extracellular signal-regulated kinase and Akt pathways control Schwann cell myelination. J Neurosci. 2004, 24: 6724-6732. 10.1523/JNEUROSCI.5520-03.2004.View ArticlePubMed
- Arthur-Farraj P, Wanek K, Hantke J, Davis CM, Jayakar A, Parkinson DB, Mirsky R, Jessen KR: Mouse schwann cells need both NRG1 and cyclic AMP to myelinate. Glia. 2011, 59: 720-733. 10.1002/glia.21144.View ArticlePubMed
- Pereira JA, Lebrun-Julien F, Suter U: Molecular mechanisms regulating myelination in the peripheral nervous system. Trends Neurosci. 2012, 35: 123-134. 10.1016/j.tins.2011.11.006.View ArticlePubMed
- Faroni A, Magnaghi V: The neurosteroid allopregnanolone modulates specific functions in central and peripheral glial cells. Front Endocrinol (Lausanne). 2011, 2: 103-
- Procacci P, Ballabio M, Castelnovo LF, Mantovani C, Magnaghi V: GABA-B receptors in the PNS have a role in Schwann cells differentiation?. Front Cell Neurosci. 2012, 6: 68-PubMedPubMed Central
- Glenn TD, Talbot WS: Signals regulating myelination in peripheral nerves and the Schwann cell response to injury. Curr Opin Neurobiol. 2013, 23: 1041-1048. 10.1016/j.conb.2013.06.010.View ArticlePubMed
- Towers S, Princivalle A, Billinton A, Edmunds M, Bettler B, Urban L, Bowery NG: GABA B receptor protein and mRNA distribution in rat spinal cord and dorsal root ganglia. Eur J Neurosci. 2000, 12: 3201-3210. 10.1046/j.1460-9568.2000.00237.x.View ArticlePubMed
- Wang H-Y, Frankfurt M, Burns LH: High-affinity naloxone binding to filamin a prevents mu opioid receptor-Gs coupling underlying opioid tolerance and dependence. PLoS One. 2008, 3: e1554-10.1371/journal.pone.0001554.View ArticlePubMedPubMed Central
- Hytrek SD, McLaughlin PJ, Lang CM, Zagon IS: Inhibition of human colon cancer by intermittent opioid receptor blockade with naltrexone. Cancer Lett. 1996, 101: 159-164. 10.1016/0304-3835(96)04119-5.View ArticlePubMed
- Zhu M, Li RC: Receptor binding activities of Schefflera triterpenoids and oligosaccharides. Planta Med. 1999, 65: 99-103. 10.1055/s-1999-13967.View ArticlePubMed
- Singer MA, Lindquist S: Multiple effects of trehalose on protein folding in vitro and in vivo. Mol Cell. 1998, 1: 639-648. 10.1016/S1097-2765(00)80064-7.View ArticlePubMed
- Kumar R: Role of naturally occurring osmolytes in protein folding and stability. Arch Biochem Biophys. 2009, 491: 1-6. 10.1016/j.abb.2009.09.007.View ArticlePubMed
- Notterpek L, Ryan MC, Tobler AR, Shooter EM: PMP22 accumulation in aggresomes: implications for CMT1A pathology. Neurobiol Dis. 1999, 6: 450-460. 10.1006/nbdi.1999.0274.View ArticlePubMed
- Schlebach JP, Peng D, Kroncke BM, Mittendorf KF, Narayan M, Carter BD, Sanders CR: Reversible folding of human peripheral myelin protein 22, a tetraspan membrane protein. Biochemistry. 2013, 52: 3229-3241. 10.1021/bi301635f.View ArticlePubMed
- Nobbio L, Mancardi G, Grandis M, Levi G, Suter U, Nave KA, Windebank AJ, Abbruzzese M, Schenone A: PMP22 transgenic dorsal root ganglia cultures show myelin abnormalities similar to those of human CMT1A. Ann Neurol. 2001, 50: 47-55. 10.1002/ana.1034.View ArticlePubMed
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: RESEARCH0034-10.1186/gb-2002-3-7-research0034.View ArticlePubMedPubMed Central
- Fledrich R, Schlotter-Weigel B, Schnizer TJ, Wichert SP, Stassart RM, Meyer zu Hörste G, Klink A, Weiss BG, Haag U, Walter MC, Rautenstrauss B, Paulus W, Rossner MJ, Sereda MW: A rat model of Charcot-Marie-Tooth disease 1A recapitulates disease variability and supplies biomarkers of axonal loss in patients. Brain. 2012, 135: 72-87. 10.1093/brain/awr322.View ArticlePubMed
- Shy ME, Chen L, Swan ER, Taube R, Krajewski KM, Herrmann D, Lewis RA, McDermott MP: Neuropathy progression in Charcot-Marie-Tooth disease type 1A. Neurology. 2008, 70: 378-383. 10.1212/01.wnl.0000297553.36441.ce.View ArticlePubMed
- Pareyson D, Scaioli V, Laurà M: Clinical and electrophysiological aspects of Charcot-Marie-Tooth disease. Neuromolecular Med. 2006, 8: 3-22. 10.1385/NMM:8:1-2:3.View ArticlePubMed
- Martini R, Klein D, Groh J: Similarities between inherited demyelinating neuropathies and Wallerian degeneration: an old repair program may cause myelin and axon perturbation under nonlesion conditions. Am J Pathol. 2013, 183: 655-660. 10.1016/j.ajpath.2013.06.002.View ArticlePubMed
- Castorina A, Scuderi S, D’Amico AG, Drago F, D’Agata V: PACAP and VIP increase the expression of myelin-related proteins in rat schwannoma cells: involvement of PAC1/VPAC2 receptor-mediated activation of PI3K/Akt signaling pathways. Exp Cell Res. 2014, 322: 108-121. 10.1016/j.yexcr.2013.11.003.View ArticlePubMed
- Hai M, Muja N, DeVries GH, Quarles RH, Patel PI: Comparative analysis of Schwann cell lines as model systems for myelin gene transcription studies. J Neurosci Res. 2002, 69: 497-508. 10.1002/jnr.10327.View ArticlePubMed
- Nobbio L, Visigalli D, Radice D, Fiorina E, Solari A, Lauria G, Reilly MM, Santoro L, Schenone A, Pareyson D: PMP22 messenger RNA levels in skin biopsies: testing the effectiveness of a Charcot-Marie-Tooth 1A biomarker. Brain. 2014, 137: 1614-1620. 10.1093/brain/awu071.View ArticlePubMed
- Katona I, Wu X, Feely SME, Sottile S, Siskind CE, Miller LJ, Shy ME, Li J: PMP22 expression in dermal nerve myelin from patients with CMT1A. Brain. 2009, 132: 1734-1740. 10.1093/brain/awp113.View ArticlePubMedPubMed Central
- Hay M, Thomas DW, Craighead JL, Economides C, Rosenthal J: Clinical development success rates for investigational drugs. Nat Biotechnol. 2014, 32: 40-51. 10.1038/nbt.2786.View ArticlePubMed
- Attarian S, Vallat J-M, Magy L, Funalot B, Gonnaud P-M, Lacour A, Péréon Y, Dubourg O, Pouget J, Micallef J, Franques J, Lefèvre M-N, Ghorab K, Al-Moussawi M, Tiffreau V, Magot A, Leclair-Visonneau L, Stojkovic T, Bossi L, Lehert P, Walter G, Bertrand V, Mandel J, Milet A, Hajj R, Boudiaf L, Scart-Grès C, Nabirotchkin S, Guedj M, Chumakov I, et al: An Exploratory randomised double-blind and placebo-controlled phase 2 study of a combination of baclofen, naltrexone and sorbitol (PXT3003) in patients with Charcot-Marie-Tooth disease type 1A. Orphanet J Rare Dis 2014, 9:199.,
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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.