PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45780.
Published online 2017 April 7. doi:  10.1038/srep45780
PMCID: PMC5384285

A staged screening of registered drugs highlights remyelinating drug candidates for clinical trials

Abstract

There is no treatment for the myelin loss in multiple sclerosis, ultimately resulting in the axonal degeneration that leads to the progressive phase of the disease. We established a multi-tiered platform for the sequential screening of drugs that could be repurposed as remyelinating agents. We screened a library of 2,000 compounds (mainly Food and Drug Administration (FDA)-approved compounds and natural products) for cellular metabolic activity on mouse oligodendrocyte precursors (OPC), identifying 42 molecules with significant stimulating effects. We then characterized the effects of these compounds on OPC proliferation and differentiation in mouse glial cultures, and on myelination and remyelination in organotypic cultures. Three molecules, edaravone, 5-methyl-7-methoxyisoflavone and lovastatin, gave positive results in all screening tiers. We validated the results by retesting independent stocks of the compounds, analyzing their purity, and performing dose-response curves. To identify the chemical features that may be modified to enhance the compounds’ activity, we tested chemical analogs and identified, for edaravone, the functional groups that may be essential for its activity. Among the selected remyelinating candidates, edaravone appears to be of strong interest, also considering that this drug has been approved as a neuroprotective agent for acute ischemic stroke and amyotrophic lateral sclerosis in Japan.

The progress in the development of immunomodulatory therapies for relapsing-remitting multiple sclerosis could not be generalized to diseases where inflammation is not prominent, including the progressive forms of multiple sclerosis1,2. This reflects the complexities often seen by those attempting to develop effective therapies3 in disorders where a neurodegenerative component is present and the experimental models are not ideal. Moreover, in progressive multiple sclerosis, several mechanisms participate in the pathophysiology of the disease, affecting different cell lineages, often with asynchronous dynamics. It is therefore possible that single treatments will bring limited benefit4,5. In such a perspective that may include the need of combination therapies, it may be advantageous to count on a substantial number of drugs targeting different pathophysiological mechanisms6.

Drug repositioning is a consolidated strategy for tackling the above problems; it cuts the timeline and costs of drug development and reduces, to some extent, the need to rely on in vivo animal models, while preserving safety7,8,9. Moreover, thanks to repurposing efforts, academic research has been able to bring, with limited resources, an increasing number of compounds to phase 2 trials in multiple sclerosis and in rare neurodegenerative diseases10,11,12,13,14,15,16. Even in the case of failure to bring these drugs to licensing, these proof-of-concept studies are providing important contributions to our understanding of “druggable” pathophysiology in neurodegenerative diseases17.

Systematic approaches are being developed to fully exploit the potential of repurposing18,19, including phenotypic screenings which may have some advantages in diseases where the pathophysiology is poorly understood. This was in fact the approach taken by three recent studies that identified pharmacological agents that were able to stimulate the differentiation of oligodendrocyte progenitor cells (OPC) into myelinating oligodendrocytes20,21,22. However, because of the above mentioned multiplicity of damage mechanisms, and because of the characteristics of the identified agents – not always compatible with prolonged treatments in chronic diseases – there is still a strong need to identify drugs that may be better tolerated, and possibly combined in polytherapies. Finally, a much-needed accomplishment is the de-risking of investments on drug development processes, including investments on repurposed drugs. Hence, the opportunity to confront data from different drug screening efforts, possibly including studies with multi-tiered evaluation systems, will be important to reduce some of the factors that currently dissuade investment in these diseases.

Here, we report a new assessment of the myelin repair capability of registered drugs through multiple phenotypic screens. Our approach was characterized by: (i) a sequential screening method that enables statistical data integration and hit picking across multiple phenotypic screens, allowing a robust ranking and prioritization of the active compounds; (ii) the validation and analysis of hits and their analogs to identify their chemical functional features, an essential step for compound optimization.

Results

Identification of compounds that stimulate oligodendrocyte development in a multiple phenotypic screen

We selected the 2,000 compounds of the Spectrum Collection library to test for the development of OPC into myelin forming oligodendrocytes. Since OPC development is a process that demands high-energy, compounds were first screened for cellular metabolic activity in mouse purified OPC cultures by the MTT assay, a rapid colorimetric test of cell reducing activity (Z-factor = 0.7). Differences in OPC ability to reduce the tetrazolium salt MTT may be due to an effect of compounds on viability, proliferation or differentiation, all of which are important components of the remyelination process. Each compound was screened once in triplicates and an MTT value was quantified by an Efficacy Ratio defined as: ER = absorbance of drug/absorbance of vehicle. Platelet derived growth factor (PDGF, 20 ng/ml) and tri-iodothyronin (T3, 30 ng/ml) and thyroxine (T4, 40 ng/ml), known inducers of OPC proliferation and differentiation, respectively, were used as positive controls (PDGF ER = 1.4 ± 0.08 and T3 + T4 ER = 1.3 ± 0.05; mean ± SEM of 5 experiments run in triplicates). On the basis of the effect of both positive controls, we defined a threshold of ER = 1.3, singling out 127 positive compounds that were further analyzed in 3–5 experiments to confirm their activity (Fig. 1a). These experiments allowed for the identification of 42 compounds presenting a mean ER value  1.3 and are suitable for human use in their present chemical form (Table 1). This set included known classes of pharmaceutical drugs such as anti-inflammatory agents, vasodilators, antibacterials, antihypertensive agents, steroids, natural products as flavonoids, and others. Among these, 5-methyl-7-methoxyisoflavone which was incorrectly stated as methoxyvone in the Spectrum Collection library (see “Quality control” paragraph in validation section). Lovastatin, a drug that presented a mean value of ER = 1.23 ± 0.1 but with a demonstrated activity on OPC differentiation and myelin formation23,24 was added to this list to provide a convenient internal control for the screening. Moreover, statins are in an advanced phase of clinical evaluation in the therapy of multiple sclerosis alone or in combination with other drugs (www.clinicaltrials.gov).

Figure 1
Cell-based screen identifies molecules with myelin repair potential.
Table 1
Compound effect on OPC metabolic activity.

We then moved to test the direct effect of the 43 (42 + lovastatin) selected compounds on the proliferation and differentiation of mouse purified OPC using specific assays (screen step in Fig. 1b). By [3 H]thymidine incorporation experiments, we established that about half of the compounds significantly stimulated the proliferation of OPC in different experiments compared to vehicle alone, with an increase that ranged from 1.5 to 7 fold compared to untreated cultures (Table 2, Prolif raw data ± SE). In parallel experiments, real time RT-PCR was used to assess the effect of the compounds on the expression of transcripts for ceramide galactosyltransferase (CGT), the key enzyme in the biosynthesis of myelin cerebrosides, and the myelin basic protein (MBP), a major component of the myelin sheath. We focused our analysis on CGT and MBP since they are markers of early and late stage of oligodendrocyte differentiation, respectively. Several compounds induced both CGT and MBP transcripts compared to vehicle alone, ranging from a 1.5 to 23 fold increase (median 4.1), and from 2.4 to 60 fold (median 5.6), respectively (Table 2, CGT and MBP raw data ± SE).

Table 2
Raw data of the activity of the 43 selected compounds originating Figure 2.

To select the best hits among the compounds active on the different biological processes analyzed, we ranked the substances by assigning a score to each experimental setting (see score assignment in methods section for details). Note that we rewarded differentiation more than proliferation as remyelination frequently fails due to a declining efficiency in the later stages of OPC development25. The score contribution of each compound was derived by multiplying the assigned score to the raw data of the experimental setting it belongs to (Table 2). Finally, score contributions have been summed per-compound and all compounds have been compared according to their total and relative score (Fig. 2). Seven compounds appeared to emerge from the others. We therefore took them as potential hit compounds. All the selected hits share remarkable properties such as neuroprotective, anti-inflammatory and anti-oxidant activities.

Figure 2
Compound effect on OPC proliferation and differentiation.

Confirmation of hit compound myelin repair capability

To evaluate the remyelinating potential of the seven hits, we analyzed their activity in three biologically relevant systems of increasing complexity, sequentially discarding the compounds that were not able to achieve a significant effect in any one of the different in vitro/ex vivo tests (screen step in Fig. 1b).

Cultures of mixed glial cells containing astrocytes, oligodendrocyte lineage cells, and microglia were used to investigate compound effects on OPC developmental processes in the presence of mitogenic and differentiation factors secreted by the other glial cell populations. Cultures were incubated with each one of the seven compounds for 48 h, using T3 and T4 (30 and 40 ng/ml, respectively) as positive control (Table 3). Real time RT-PCR was used to analyze mRNA levels of CGT and MBP in three different experiments, showing that five compounds significantly up-regulated both transcripts (Table 3). The effect of the anti-inflammatory vulpinic acid and the anti-bacterial fenamisal did not reach statistical significance.

Table 3
Compound effect on OPC differentiation in cultures of mixed glial cells.

We next tested whether the five hit compounds could promote the differentiation of endogenous OPC in central nervous system (CNS) tissue. Cerebellar slices were generated from mice at postnatal day 7, a time that corresponds to the onset of myelination. Two hours after plating, slices were treated with compounds for 5 DIV and the expression MBP mRNA was evaluated by real time RT-PCR in 3–5 experiments. We focused on MBP transcripts since this myelin marker shows a wider range of expression than CGT during the developmental period analyzed. Progesterone (40 μM) was used as a positive control due to its myelin stimulating activity in this ex vivo model26. We found that only 5-methyl-7-methoxyisoflavone, edaravone, losartan and lovastatin significantly (p  0.05) increased MBP transcript levels relative to control slices, demonstrating their ability to promote OPC differentiation also in this experimental model (Fig. 3).

Figure 3
Compound induced MBP mRNA expression in cerebellar slices.

Subsequently, we evaluated the remyelinating potential of the four selected compounds in toxin-treated cerebellar slice cultures, an ex vivo model of myelin damage where limited remyelination occurs in basal culture conditions. Cerebellar slices from P10 mice were cultured for 7 DIV and then treated for 16 h with lysolecithin to induce demyelination, as previously described27. Immediately after toxin removal, slices were incubated with the compounds (20 μM) or DMSO (0.002% vehicle) for the indicated times, using progesterone (40 μM) as positive control28. By real time RT-PCR, we demonstrated that after 4 days of incubation edaravone, 5-methyl-7-methoxyisoflavone and lovastatin, but not losartan, significantly (p  0.05) stimulated MBP transcripts compared to control slices (Fig. 4a). By a method of image confocal analysis quantification developed in house, we also evaluated the ability of the three hit compounds to induce an increase in axonal remyelination. In particular, the co-localization of the myelin protein MBP and the axonal protein NFH was analyzed in lysolecithin demyelinated slices maintained in the absence or presence of the compounds (20 μM) for 7 days. Consistent with the MBP gene expression analysis, we found that edaravone, 5-methyl-7-methoxyisoflavone, and lovastatin significantly increased MBP+ myelin membranes relatively to NFH+ axons compared to control slices (Fig. 4 b and c).

Figure 4
Compound induced MBP expression and remyelination in lysolecithin-demyelinated cerebellar slices.

Validation of edaravone and 5-methyl-7-methoxyisoflavone biological activity and chemical structure

By applying established procedures of the modern drug discovery pipelines, we validated the biological activity and chemical structure of the antioxidant edaravone and the flavonoid 5-methyl-7-methoxyisoflavone. The hypolipidemic compound lovastatin was not included in this analysis since, as previously mentioned, it was used as an internal control of the screening for its documented remyelinating activity23.

Quality Control

The analytical characterization of purity of the two hit compounds was performed by Proton Nuclear Magnetic Resonance Spectroscopy (1H-NMR) and Liquid chromatography–mass spectrometry (LC/MS) techniques. On the basis of the results obtained, we identified with certainty that the active compound of the screening is 5-methyl-7-methoxyisoflavone, which is incorrectly identified in the library data set as the flavone methoxyvone (7 methoxy-5methyl 2 phenylcromen 4-one) as the two compounds differ only for the position of the phenyl group (isoflavones and flavones present the C ring in position 2′ and 3′ of the B ring, respectively). Overall we confirmed that the chemical compounds analyzed correspond to edaravone and 5-methyl-7-methoxyisoflavone and present a purity of 97% and 99%, respectively (Supplementary Fig. S1), demonstrating that the activity depends exclusively on compounds and not on impurities (less than 5%).

Retesting

To confirm the compound activity and avoid false positives, we purchased independent stocks of the two hit compounds from different suppliers and tested their remyelinating potential in organotypic demyelinated slices. The increased level of MBP transcripts in slices incubated with the new stocks of edaravone and 5-methyl-7-methoxyisoflavone (20 μM, 4 days) compared to control slices confirmed the remyelinating potential of the two hits, demonstrating that the results obtained with the two compounds are not false positives (Supplementary Fig. S2).

Dose-response curves

To assess the dose dependence of the assay’s readouts, we performed a dose-response analysis in cultures of mixed glial cells, as this is the most suitable test among those we had employed. Cells were incubated with edaravone and 5-methyl-7-methoxyisoflavone at five different concentrations (0.01–30 μM) for 48 h and expression of MBP mRNA was evaluated by real time RT-PCR. The results demonstrated that: i) the effect of both compounds was dose dependent and no activity was present at the dose of 0.01 μM, further confirming that there are no false positives; ii) the lowest and highest responses of 5-methyl-7-methoxyisoflavone were established at 0.1 and 1 μM concentration, respectively, with an evident decrease of the activity at 30 μM; whereas, edaravone showed a significant effect at the dose of 1 μM but it may not have reached its maximal efficacy at 30 μM (Fig. 5a and b). These results indicated that both compounds are active within a micromolar range but only edaravone did not show any toxic effect up to the concentration tested.

Figure 5
Dose-response curves of the two hit compounds in cultures of mixed glial cells.

Chemical analog selection and testing

To identify a relationship between the structure and activity of two hits, we selected and tested structurally related chemical analogs. We selected new chemical analogs from commercial chemical libraries by means of a public database (Zinc, PubChem and ChemSpider), finding four available analogs for edaravone (Figs 1a–d and 6a,b) and five for 5-methyl-7-methoxyisoflavone, of which only one was an isoflavone and four were flavones (Figs 2a, 3a–d and 6d,e,). The activity of the selected analogs (10 μM, 48 h) was evaluated in OPC purified cultures by MTT test, identifying two active analogs for edaravone (1b ER = 1.36 ± 0.029; and 1d ER = 1.44 ± 0.06) and none for 5-methyl-7-methoxyisoflavone. The biological activity of edaravone analogs was then confirmed in mixed glial cell cultures with dose-response curves at six different concentrations (0.001–30 μM; Supplementary Fig. S3) and in demyelinated cerebellar slices (Supplementary Fig. S4). By searching for chemical analogs in the Spectrum Collection library we found three analogs of edaravone (the antipyretic and analgesic drugs antipyrine, aminopyrine and ramifenazone; Fig. 6b) and 7 isoflavones (Fig. 6d) that did not overcome the MTT test in both primary and re-testing experiments. Since in the Spectrum Collection library many flavones and isoflavones are present, we selected only molecules with the same scaffold of the reference compound (Fig. 6d).

Figure 6
Activity of chemical analogs of edaravone and methoxyisoflavone and identification of edaravone chemical features.

Chemical functional feature identification

By using the chemical structure of the two hit compounds as reference and their commercial and library analogs, we were able to identify some functional groups indispensable for their activity. As shown in Fig. 6c, we identified four functional groups essential for the activity of edaravone: 2 H bond acceptors (position 1 and 3, in red) plus hydrophobic group (the methyl group in position 5, in green) and an aromatic ring (position 2, in orange). These four groups are all present in active analogs: i) 1b has a para-substitution on phenyl group in R, but the presence of a fluorine (a small and electronegative atom) did not prevent the activity; ii) 1d is different from edaravone only for the presence of an ethyl group, which is always a hydrophobic group, in R1. These substitutions do not change the features responsible for the activity. On the contrary, the following substitutions result in loss of activity: i) the presence of an isopropylamino or dimethylamino or a chlorine group in R2, as in ramifenazone, aminopyrine and 1c, respectively. Actually, the active compounds do not exhibit substituent in R2, suggesting that this position must remain empty; ii) the absence of H bond acceptor in position 1, as in antipyrine, ramifenazone and aminopyrine; iii) the replacement of the five-membered ring with one to six-membered, as in 1a, which increases the distance between the hydrophobic group in position 5 and the H bond acceptor in position 3.

The lack of active analogs of 5-methyl-7-methoxyisoflavone among the few commercially available did not allow for the identification of functional groups essential for its activity. However, we were able to define that: i) the presence of a methoxy group in position R2 may prevent the activity, as in 2a and in the library isoflavones biochanin A and formonetin (Fig. 6d); ii) the position of the ring in ’2 causes loss of activity; indeed, flavones 3a, 3b, 3c and 3d were inactive (Fig. 6e). Comparison between the chemical structure of 5-methyl-7-methoxyisoflavone and that of inactive isoflavones present in the library (Fig. 6d) showed that 5-methyl-7-methoxyisoflavone is the only compound that has a methyl group in R1, replacement with a methoxy group, as for 5–7 dimethoxyisoflavone (the two compounds differ only for this substituent), a hydroxyl group, as for 5–7 dihydroxyisoflavone, genistein, biochanin A, or a hydrogen, as for daidzein and formonetin, determines inactivity. These results suggest that the position in R1 must be occupied by a hydrophobic aliphatic substituent probably in combination with the absence of substitutions in R2. Finally, replacement in R with a hydroxyl group, as for genistein, daidzein, biochanin A or isopropoxy group, as for ipriflavone, are unfavorable for the activity.

Prediction of compound targets

We used Chemical Similarity Network Analysis Pulldown (CSNAP) as the computational method to look for a possible consensus target that could explain the effect produced by the hit compounds edaravone, 5-methyl-7-methoxyisoflavone, and lovastatin. We included in the analysis other top scoring compounds of our screening (chlormadinone acetate, fenamisal, losartan, vulpinic acid) and four putative remyelinating compounds from other phenotypic screenings (benztropine, clemastine, clobetasol, miconazole; 20,21,22). Biochanin A, a compound that did not exceed the screening selection, was used as a negative control. CSNAP highlighted Neuropeptide S Receptor 1 (NPSR1; Uniprot Q6W5P4) and isocitrate dehydrogenase (Uniprot O75874) as possible consensus targets, showing the highest cumulative Schwikowski-score (S-score) across our hit compounds with a meaningful pattern across the other top scoring and putative remyelinating compounds (Supplementary Dataset S1) as well. The NPS G-coupled receptor was selected as a target for in vitro validation since, besides having the highest score, it belongs to a class of receptors that are implicated in many diseases and are the targets of numerous therapeutic drugs29.

The investigation of the pharmacological action of edaravone, 5-methyl 7-metoxyisoflavone, lovastatin, benztropine, losartan and biochanin A at recombinant NPSR was performed in HEK293 cells stably expressing the mouse NPSR by a calcium mobilization assay30. The results of these experiments showed that NPS increased the intracellular calcium concentrations in a concentration dependent manner with pEC50 and Emax values of 8.75 (8.30–9.20) and 335 ± 32% over the basal values, respectively. The compounds did not stimulate calcium mobilization up to a concentration of 10 μM (Supplementary Fig. S5). These compounds were then assayed as NPSR antagonists vs the stimulatory action elicited by NPS. In this series of experiments the selective NPSR antagonist SHA 68 was used as a standard. At 100 nM, SHA 68 was able to elicit a rightward displacement of the concentration response curve to NPS with a pA2 value of 7.83 (7.65–8.01). At the concentration of 10 μM, compounds did not affect the concentration response curve to NPS (Supplementary Fig. S6). Therefore, the tested compounds do not interact with the NPSR, either as agonists or as antagonists.

Discussion

Multiple sclerosis pathophysiology is highly complex and difficult to reduce to in vitro or in vivo models where to investigate for new therapies. Nonetheless, high-throughput screenings for myelin repair are beginning to yield the first, consistent results20,21,22. However, because of the above complexity, false negatives and false positives are and will be inevitable. To limit them, comparisons of multiple screenings within and across studies are much needed. For the same reasons, and given the high numbers of compounds to be screened, it is important that these studies are structured as to balance inclusiveness and selectivity of the analysis.

CNS remyelination is the result of a specific sequence of events that are tightly regulated and recapitulate the process of oligodendrocyte generation that occurs during development25. Our staged series of assays, composed by multiple in vitro screens, allowed for the testing of compound effects on different aspects of OPC development and for the selection of the best hits by assessing their activity across the combined data set rather than their status within a single analysis. This multiple phenotypic screening highlighted a substantial number of compounds that are suitable for human use. These molecules belong to various classes of pharmaceutical drugs and represent useful information for further research in the field.

Three drugs, lovastatin, edaravone and 5-methyl-7-methoxyisoflavone, convincingly demonstrated activity in all the tiers of the screening process; their remyelinating potential arose from an increase in OPC proliferation and differentiation, which are essential components of the remyelination process. Based on the evidence that remyelination failure is more commonly associated with the impairment of the final stage of oligodendrocyte developmental program25, strategies to enhance remyelination are frequently designed to promote OPC differentiation into new myelinating oligodendrocytes31. However, several remyelinating molecules have been described to promote both proliferation and differentiation32,33,34 or proliferation alone35. These two processes are not mutually exclusive for the progression of oligodendrocyte lineage; indeed cAMP response element-binding protein (CREB), implicated in controlling myelin gene expression and oligodendrocyte generation, participates in the regulation of both OPC proliferation36 and differentiation37.

The biological plausibility of lovastatin as drugs with myelin repair potential confirms the validity of our approach and reinforces the rational for pursuing the clinical evaluation of lipophilic statins in the therapy of multiple sclerosis.

The hit compound edaravone is a free radical scavenger, approved as a neuroprotective agent for acute ischemic stroke in Japan, where it has been used for over 10 years. Surprisingly, it does not have market authorization in other countries38. It is a small bioactive molecule that easily crosses the blood brain barrier and has been indicated as a promising drug candidate for several neurodegenerative diseases39,40. Edaravone was also investigated in amyotrophic lateral sclerosis patients41 and recently received approval, again in Japan, for this disease. With respect to multiple sclerosis, encouraging studies in the model of experimental allergic encephalomyelitis already exist42. It is of interest that edaravone rescues oligodendrocyte differentiation after prolonged ischemic damage in mice43, reduces white matter injury in a mouse model of hypoxic-ischemia44, and has anti-inflammatory effects in activated microglia45. The identification of two active analogs of edaravone and the purity analysis validated its biological and structural activity, allowing us to conclude that the remyelinating activity of edaravone is not a false positive. The micromolar concentration range of activity and the lack of cytotoxicity shown by edaravone in our experimental models are consistent with previous studies performed on different cellular types, including all components of cerebrovascular units46,47. The weak activity of edaravone does not seem to be a major obstacle to its human use as the drug is registered not only for an acute disease (ischemic stroke), but also for a chronic condition (amyotrophic lateral sclerosis). Moreover, there are ongoing attempts to improve its poor oral bioavailability. If successful, these attempts may not only promote its use in chronic diseases but possibly also increase the commercial interest about its development as a licensed drug in different conditions48,49,50.

On these grounds, we suggest that it is urgent to initiate the clinical evaluation of edaravone in multiple sclerosis, considering also that based on our structure-activity relationship analysis, the molecule may be chemically modified to enhance its activity: four functional groups appear to be essential for the activity together with the presence of a five-membered ring that allows for the correct distance between the hydrophobic group in position 5 and the H bond acceptor in position 3.

The other hit of our screening, 5-methyl-7-methoxyisoflavone, is a semi-natural isoflavone. Isoflavones are the main group of flavonoids and are well characterized for their estrogenic and/or antioxidant properties. They are active ingredients of human diet since they are present in many daily-consumed foods, like vegetables and fruits. As plant-derived secondary metabolites, and in spite of some limitations about their final development as drugs51, they are relevant for drug discovery since they have a wide range of molecular targets. Their neuroprotective action may depend on estrogenic and other properties. 5-methyl-7-methoxyisoflavone has never been considered as a candidate for CNS repair. It shows some differences in chemical structure compared to traditional isoflavones: a methoxy group in R, found in some isoflavones, an unusual methyl group in R1 (Fig. 6d), but no hydroxyl groups, which are typical of flavonoids with mainly antioxidant activity. In addition, the lack of structure similarity to 17β-estradiol (the most potent estrogen in mammals) suggests that the activity of 5-methyl-7-methoxyisoflavone in our assays may not be related to the estrogenic properties typical of isoflavones such as genistein and daidzein, an interpretation that is supported - to some extent - by the inactivity of genistein and daidzein in our primary screening. Unfortunately, at variance with edaravone, the analogs of 5-methyl-7-methoxyisoflavone were not enough to allow the identification of a clear SAR. Driven synthetic modifications will be needed to improve pharmacological properties of this compound.

Overall, the possibility that the antioxidant action of the two hit compounds is not primarily responsible for their regenerative effect is supported by the fact that many other antioxidant compounds of the library did not show any activity in the screening process. Based on these considerations, the investigation of the molecular targets and pathways activated by edaravone and 5-methyl-7-methoxyisoflavone stands a good chance to uncover novel biological activities relevant to the regenerative process of remyelination. Apart from the elements of internal consistency of the data, and those in accord with published results from other groups, there are compounds that have demonstrated myelin repair potential in other experimental systems but whose efficacy did not emerge in our screening. In particular, the class of anticholinergic compounds with anti-muscarinic activity was not represented among our best hits in spite of the positive data from two different studies. While benztropine and clemastine had been identified as potential remyelinating compounds in screening assays performed at 1 μM final concentration20,21, they showed toxicity in our primary screening performed at 10 μM concentration. In MTT re-testing experiments, we confirmed the activity of both compounds at the lower concentration (1 μM, data not shown). This is consistent with the issue of low difference between active and toxic doses described for benztropine20.

The optimization of therapeutic algorithms, or the development of new molecules with better therapeutic indexes, would benefit from better knowledge about the targets that are relevant for the regenerative effects of these drugs. So far, and including this study, none of the screenings have unequivocally defined the therapeutic target(s) of the identified compounds. This is not uncommon (approximately, only the 20% of the small molecule drugs have an identified individual molecular target) and may reflect the need for a promiscuous mechanism of action for “neuroprotective” drugs to be effective. Interestingly, some of the recently identified drugs with myelin repair potential might share a common mechanisms of action; it has been recently shown that edaravone effects are mediated not only by the elimination of oxidative stress, but also by the increased production of brain-derived neurotrophic factor (BDNF) through phosphorylation of CREB52. Similarly, lipophilic statins may exert neurotrophic functions through the transcriptional activation of CREB and consequent BDNF increase53. Besides, the phosphorylation of CREB is also part of the cascade triggered during benztropine-induced oligodendrocyte differentiation20 and is involved in the regulation of OPC differentiation by GPR17, another potential target for the treatment of myelin damage54.

In conclusion, in spite of differences in the experimental platforms, and in the composition of the libraries of compounds, there are elements of consistency between the last screening studies on the remyelinating potential of small molecules. Moreover, all hits that emerged from the last four (20,21,22 and the present one) screenings are molecules that are registered for human use. This means that, while additional studies with different platforms are needed, there are already enough drugs to design experimental medicine trials to explore the clinical efficacy of the compounds (also in combination therapies) and increase our understanding of the biology of the disease. The results of this study and the clinical effects in two other neurological diseases strongly suggest that edaravone is one of the most interesting candidates. Future work will need to concentrate more on medicinal chemistry studies to optimize the compounds and increase the interest for industrial development. However, all these advancements also urge the support of political steps and regulatory solutions to facilitate the entire process of repurposing up to a full marketing authorization.

Materials and Methods

Experimental design

The goal of this study was to identify new therapeutic strategies for myelin repair that can be readily transferred to patients with demyelinating diseases of CNS. We carried out an unbiased in vitro and ex vivo screening in mammalian cells with the aim of repositioning drugs according to their remyelinating potential. The compound activity was first verified on the development of mouse OPC through multiple, sequential phenotypic screens and their myelin repair capability was then assessed in organotypic mouse brain slice cultures. The biological activity and chemical structure of hit compounds were validated, also allowing for the identification of their chemical functional features.

Animals

CD1 Swiss mice were purchased from Harlan Laboratories (San Pietro Al Natisone, Udine, Italy). The experimental procedures related to the use of CD1 Swiss mice were performed in accordance with the approved national law with the implementation of European Union Directive nr. 86/609/CEE - D.lgs 116/92 regarding the Protection of animals used for experimental research. The experiments were approved by the Service for Biotechnology and Animal Welfare of the Istituto Superiore di Sanità and by the Italian Ministry of Health (Authorization #271/SSA/2010 -18/03/2010).

Compound sources

The library tested was The Spectrum Collection (Micro- Source Discovery Inc., Groton, Conn.). An alphabetical list of the compounds is available at the Micro-Source Discovery website at www.msdiscovery.com/spectrum.html. The compounds are supplied as 10 mM solutions in dimethyl sulfoxide (DMSO). Edaravone and 5-methyl-7-methoxyisoflavone used in retesting experiments were purchased from Sigma-Aldrich and Ambinter (Orléans, France), respectively. Chemical analogs were purchased from Ambinter (2a and 3a), and from MolPort (Riga, Latvia; 1a–d and 3b,d).

Staged screening of the compounds

Primary screen

We screened 2,000 compounds of The Spectrum Collection library comprising primarily Food and Drug Administration (FDA)-approved compounds or natural products. Purified OPC obtained from neonatal mouse primary mixed glial cultures (for details see Supplementary methods) were incubated with or without drugs at 10 μM concentration in dimethyl sulphoxide (DMSO; 0.001% vehicle) for 48 h. PDGF(20 ng/ml, Preprotech, Rocky Hill, NJ) and T3 (30 ng/ml; Sigma-Aldrich, St Louis, MO) and T4 (40 ng/ml; Sigma-Aldrich) were used as positive controls. The 10 μM concentration was selected to obtain more primary hits and more chemical scaffold types while simultaneously removing compounds that may have problems of toxicity55. Compounds were first screened for their ability to reduce [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT)] using an automated microplate reader (Bio-rad, Hercules, CA, USA). Active compounds were then screened for OPC proliferation and differentiation assessing [3 H]thymidine incorporation and the expression of transcripts for the CGT and MBP markers of early and late stage of oligodendrocyte differentiation, respectively.

Score assignment

For the combined analysis of proliferation (Prolif), early differentiation (Diff_CGT) and late differentiation (Diff_MBP) data, we assigned a score to each of the experimental setting performed in purified OPC cultures: MTT score = 1, Prolif score = 10, Diff_CGT score = 3, Diff_MBP score = 5. Score numbers have been assigned with the aim of: (i) isolate MTT contribution [range 0–3]; (ii) keep Prolif scored contribution within the same range of Diff_ CGT and MBP [range 0–60]; (iii) reward both Diff_CGT and Diff_MBP more than proliferation since the final stage of oligodendrocyte developmental program is impaired during the remyelination process25; (iv) reward Diff_MBP more than Diff_CGT since MBP mRNA is subjected to a broader range of expression changes with respect to CGT mRNA during the in vitro developmental period analyzed [range 0–300] (Table 2).

Mixed glial cell cultures

Cells enzymatically dissociated from forebrains of newborn CD1 Swiss mice (1.2 × 105 cells/cm2) were grown in Dulbecco’s modified eagle medium containing 10% Foetal Bovine Serum (FBS), 2 mM glutamine and penicillin (50 μg/ml) and streptomycin (50 μg/ml); all media, sera and reagents were by GIBCO (Lifetechnologies, Grand Island, NY). After 9–10 DIV, when small groups of OPC started to proliferate on the top of the astroglia layer together with microglial cells, medium was replaced with a defined serum-free medium without thyroid hormones56. Cells were then incubated with the selected compounds (10 μM) or DMSO (0.001% vehicle) for 48 h. T3 (30 ng/ml; Sigma-Aldrich) and T4 (40 ng/ml; Sigma-Aldrich) were used as positive controls.

Organotypic cerebellar slice cultures

For myelination studies, organotypic cerebellar slices were prepared from mouse cerebellum as previously described for rat26. Briefly, slices were prepared from cerebellum of postnatal day 7 CD1 Swiss mice and grown in culture medium with compounds at 20 μM concentration or DMSO (0.002% vehicle) for 5 DIV (replacing fresh medium after 1 DIV and every 2–3 days). A compound dose two times higher than that used in monolayer cell cultures was used due to the thickness of the slices.

For remyelination studies, cerebellar slices were prepared from mouse cerebellum as previously described for rat27. Briefly, slices were prepared from cerebellum of postnatal day 10 CD1 Swiss mice and grown in culture medium for 7 DIV (replacing fresh medium after 1 DIV and every 2–3 days). Demyelination was then induced by addition of 0.5 mg/ml lysolecithin (Sigma-Aldrich) to the medium for 16 h. Two hours after toxin removal, cultures were incubated with the compounds (20 μM) or DMSO (0.002% vehicle) for 4 days (for gene expression experiments) or 7 days (for immunohistochemistry). Progesterone (40 μM, Sigma-Aldrich) was used as positive control.

Automated immunohistochemistry quantification of remyelination

We used a LSM 5 Pascal Laser Scanning Microscope to obtain stacks of photographs of MBP and NFH immunolabelling at 1 μm intervals in white matter areas at ×20 magnification in remyelinating slices with or without drugs. For each experimental point, 3 stacks were analyzed in three separate experiments in a blinded fashion by two different investigators. We developed an image analysis method to automatically quantify only MBP staining overlying axons. We implemented a statistical evaluator based on a sound and accepted methodology57 for pixel quantification which are both green and red above a defined intensity overlap, producing a mask of co-localization. The value of co-localization of each layer within a confocal stack was divided for the corresponding NFH staining to form the “Index of remyelination”. Our statistic combined three steps: (i) a split channel filter that isolates different color channels in a given image stack; (ii) a sub-color palette matching algorithm (e.g., this was used to sample the presence of the yellow color according to the palette); (iii) the assessment of biological expression co-localization58.

Analysis of purity

HPLC quantitative determination of the purity of edaravone and 5-methyl-7-methoxyisoflavone was performed by Varian Pro star (waters symmetry C18 25 × 4.6 5 μm, loop 10 μl, 254 nm). High resolution mass spectrum was obtained with the spectrometer Mariner Applied Biosystems ESI-TOF. Spectra NMR (DMSO-300K) were registered using a spectrometer Bruker AVANCE III 400 MHz Ultrashield Plus.

Chemical analog selection

Search in public Database (PubChem, Zinc and ChemSpider) of chemical analogs was performed with different approaches: bio-dimensional similarity using Tanimoto score calculated from 2D structure fingerprint and substructure search using main scaffolds.

Chemical Similarity Network Analysis

We used Chemical Similarity Network Analysis Pull-down (CSNAP; http://services.mbi.ucla.edu/CSNAP; University of California, Los Angeles) as a computational method to look for a possible consensus target of the hit compounds identified through the phenotype-based screening. The technique behind CSNAP utilizes chemical similarity networks for chemotype (consensus chemical pattern) recognition and drug target profiling. The result of this technique is presented via a Schwikowski-score that associates a compound to a list of targets (Supplementary Dataset S2) through the generation of ligand-target interaction fingerprint59. The objective of this analysis was the identification of robust interaction networks60 and the identification of a possible consensus target, i.e. a target shared by the majority of the candidate compounds that previously showed an effect when analyzed singularly in vitro and ex vivo.

Statistical analysis

For all the experiments in which statistical analysis was performed, a paired Student’s t-test was used. p  0.05 was considered significant.

Additional Information

How to cite this article: Eleuteri, C. et al. A staged screening of registered drugs highlights remyelinating drug candidates for clinical trials. Sci. Rep. 7, 45780; doi: 10.1038/srep45780 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Information:
Supplementary Dataset S1:
Supplementary Dataset S2:

Acknowledgments

This work has been supported by Fondazione Roma project (2009–2012). CENTERS is supported as a special project by the Italian Foundation for Multiple Sclerosis (FISM).The authors thank Dr. Matteo Floris for his support to database research and Dr. Adem Albayrak for final English language editing.

Footnotes

The authors declare no competing financial interests.

Author Contributions C.E. and C.V. performed all the laboratory experiments and contributed to data analysis and interpretation. S.O. was responsible for chemical analog selection and chemical functional feature identification. R.U. performed compound score assignment, image analysis and CSNAP. R.M. contributed to the performance of real time RT-PCR experiments. F.F. and G.C. performed calcium mobilization experiments. S.R. assisted in the conception of the studies, interpretation and analysis of the data. M.C.B. assisted in the conception and design of the project. G.R. assisted in the interpretation of the data and co-wrote the manuscript. M.S. conceived, designed the study and co-wrote the manuscript. C.A. conceived, coordinated the study and co-wrote the manuscript. All authors contributed to the critical review of the manuscript. All authors read and approved the final version of the manuscript.

References

  • Thompson A. J. A much-needed focus on progression in multiple sclerosis. Lancet Neurol. 14, 133–135 (2015). [PubMed]
  • Stankoff B. et al. . Repair strategies for multiple sclerosis: challenges, achievements and perspectives. Curr Opin Neurol. 29, 286–292 (2016). [PubMed]
  • Plenge R. M. Disciplined approach to drug discovery and early development. Sci Transl Med. 349, 1–5 (2016). [PubMed]
  • Bei F. et al. . Restoration of Visual Function by Enhancing Conduction in Regenerated Axons. Cell. 164, 219–232 (2016). [PMC free article] [PubMed]
  • Smith A. L. & Cohen J. A. Multiple sclerosis: Fingolimod failure in progressive MS INFORMS future trials. Nat Rev Neurol. 12, 253–254 (2016). [PubMed]
  • Salvetti M. et al. . Progressive MS: from pathophysiology to drug discovery. Mult Scler. 21, 1376–1384 (2015). [PubMed]
  • Dolgin E. Nonprofit disease groups earmark grants for drug repositioning. Nat Med. 17, 1027 (2011). [PubMed]
  • Nosengo N. New tricks for old drugs. Nature. 534, 314–316 (2016). [PubMed]
  • Reed J. C. et al. . The NIH’s role in accelerating translational sciences. Nat Biotech. 30, 16–19 (2012). [PubMed]
  • Ristori G. et al. . Use of Bacille Calmette-Guèrin (BCG) in multiple sclerosis. Neurology. 53, 1588–1589 (1999). [PubMed]
  • Ristori G. et al. . Riluzole in cerebellar ataxia: a randomized, double-blind, placebo-controlled pilot trial. Neurology. 74, 839–845 (2010). [PubMed]
  • Ristori G. et al. . Effects of Bacille Calmette Guèrin after the first demyelinating event in the CNS. Neurology. 82, 41–48 (2014). [PMC free article] [PubMed]
  • Barkhof F. et al. . Ibudilast in relapsing-remitting multiple sclerosis: a neuroprotectant? Neurology. 74, 1033–1040 (2010). [PubMed]
  • Arun T. et al. . Targeting ASIC in primary progressive multiple sclerosis: evidence of neuroprotection with amiloride. Brain. 136, 106–115 (2013). [PubMed]
  • Chataway J. et al. . Effect of high-dose simvastatin on brain atrophy and disability in secondary progressive multiple sclerosis (MS-STAT): a randomised, placebo-controlled, phase 2 trial. Lancet. 383, 2213–2221 (2014). [PubMed]
  • Romano S. et al. . Riluzole in patients with hereditary cerebellar ataxia: a randomised, double-blind, placebo-controlled trial. Lancet Neurol. 14, 985–991 (2015). [PubMed]
  • Giovannoni G., Baker D. & Schmierer K. The problem with repurposing: is there really an alternative to Big Pharma for developing new drugs for multiple sclerosis? Mult Scler Relat Disord. 4, 3–5 (2015). [PubMed]
  • Hurle M. R., Yang L., Xie Q., Rajpal D. K., Sanseau P. & Agarwal P. Computational drug repositioning: from data to therapeutics. Clin Pharmaco Ther. 93, 335–341 (2013). [PubMed]
  • Ricigliano V. A. et al. . Contribution of genome-wide association studies to scientific research: a pragmatic approach to evaluate their impact. PloS One. 8, e71198 (2013). [PMC free article] [PubMed]
  • Deshmukh V. A. et al. . A regenerative approach to the treatment of multiple sclerosis. Nature. 502, 327–332 (2013). [PMC free article] [PubMed]
  • Mei F. et al. . Micropillar arrays as a high throughput screening platform for therapeutics in multiple sclerosis. Nat Med. 20, 945–960 (2014). [PMC free article] [PubMed]
  • Najm F. J. et al. . Drug-based modulation of endogenous stem cells promotes functional remyelination in vivo. Nature. 522, 216–220 (2015). [PMC free article] [PubMed]
  • Paintlia A. S., Paintlia M. K., Singh A. K. & Singh I. Inhibition of rho family functions by lovastatin promotes myelin repair in ameliorating experimental autoimmune encephalomyelitis. Mol Pharmacol. 73, 1381–1393 (2008). [PMC free article] [PubMed]
  • Paintlia A. S., Paintlia M. K., Singh A. K., Orak J. K. & Singh I. Activation of PPAR-γ and PTEN cascade participates in lovastatin-mediated accelerated differentiation of oligodendrocyte progenitor cells. Glia. 58, 1669–1685 (2010). [PMC free article] [PubMed]
  • Franklin R. J. & Ffrench-Constant C. Remyelination in the CNS: from biology to therapy. Nat Rev Neurosci. 9, 839–855 (2008). [PubMed]
  • Ghoumari A. M. et al. . Progesterone and its metabolites increase myelin basic protein expression in organotypic slice cultures of rat cerebellum. J Neurochem. 86, 848–859 (2003). [PubMed]
  • Birgbauer E., Rao T. S. & Webb M. Lysolecithin induces demyelination in vitro in a cerebellar slice culture system. J Neurosci Res. 78, 157–166 (2004). [PubMed]
  • Hussain R. et al. . Progesterone and Nestorone facilitate axon remyelination: a role for progesterone receptors. Endocrinology. 152, 3820–3831 (2011). [PubMed]
  • Congreve M., Dias J. M. & Marshall F. H. Chapter One – Structure-Based Drug Design for G Protein-Coupled Receptors. Prog. Med. Chem. 53, 1–63 (2014). [PubMed]
  • Ruzza C. et al. . Further studies on the pharmacological profile of the neuropeptide S receptor antagonist SHA 68. Peptides. 31, 915–925 (2010). [PubMed]
  • Harlow D. E., Honce J. M. & Miravalle A. A. Remyelination Therapy in Multiple Sclerosis. Front Neurol. 6, 257 (2015). [PMC free article] [PubMed]
  • Deverman B. E. & Patterson P. H. Exogenous leukemia inhibitory factor stimulates oligodendrocyte progenitor cell proliferation and enhances hippocampal remyelination. J Neurosci. 32, 2100–2109 (2012). [PMC free article] [PubMed]
  • Zhang J. et al. . Fingolimod treatment promotes proliferation and differentiation of oligodendrocyte progenitor cells in mice with experimental autoimmune encephalomyelitis. Neurobiol Dis. 76, 57–66 (2015). [PubMed]
  • Matías-Guíu J., Oreja-Guevara C., Matias-Guiu J. A. & Gomez-Pinedo U. Vitamin D and remyelination in multiple sclerosis. Neurologia. S0213-4853, 30065–2 (2016).
  • Skripuletz T. et al. . Pivotal role of choline metabolites in remyelination. Brain. 138(Pt 2), 398–413 (2015). [PubMed]
  • Johnson J. R., Chu A. K. & Sato-Bigbee C. Possible role of CREB in the stimulation of oligodendrocyte precursor cell proliferation by neurotrophin-3. J Neurochem. 74,1409–17 (2000). [PubMed]
  • Syed Y. A. et al. . Inhibition of phosphodiesterase-4 promotes oligodendrocyte precursor cell differentiation and enhances CNS remyelination. EMBO Mol Med. 5, 1918–34 (2013). [PMC free article] [PubMed]
  • Wu S., Sena E., Egan K., Macleod M. & Mead G. Edaravone improves functional and structural outcomes in animal models of focal cerebral ischemia: a systematic review. Int J Stroke. 9, 101–6 (2014). [PubMed]
  • Yuan W. J. Neuroprotective effects of edaravone-administration on 6-OHDA-treated dopaminergic neurons. BMC Neurosci. 9, 75 (2008). [PMC free article] [PubMed]
  • Jiao S. S. Edaravone alleviates Alzheimer’s disease-type pathologies and cognitive deficits. Proc Natl Acad Sci USA. 112, 5225–30 (2015). [PubMed]
  • Abe K. et al. . Confirmatory double-blind, parallel-group, placebo-controlled study of efficacy and safety of edaravone (MCI-186) in amyotrophic lateral sclerosis patients. Amyotroph Lateral Scler Frontotemporal Degener. 15, 610–7 (2014). [PMC free article] [PubMed]
  • Moriya M. et al. . Edaravone, a free radical scavenger, ameliorates experimental autoimmune encephalomyelitis. Neurosci Lett. 440, 323–6 (2008). [PubMed]
  • Miyamoto N. et al. . Oxidative stress interferes with white matter renewal after prolonged cerebral hypoperfusion in mice. Stroke. 44, 3516–21 (2013). [PMC free article] [PubMed]
  • Shen Y., Liu X. B., Pleasure D. E. & Deng W. Axon-glia synapses are highly vulnerable to white matter injury in the developing brain. J Neurosci Res. 90, 105–121 (2012). [PMC free article] [PubMed]
  • Yuan H., Zha H., Rangarajan P., Ling E. A. & Wu C. Anti-inflammatory effects of Edaravone and Scutellarin in activated microglia in experimentally induced ischemia injury in rats and in BV-2 microglia. BMC Neurosci. 15, 125 (2014). [PMC free article] [PubMed]
  • Lee B. J., Egi Y., van Leyen K., Lo E. H. & Arai K. Edaravone, a free radical scavenger, protects components of the neurovascular unit against oxidative stress in vitro. Brain Res. 1307, 22–27 (2010). [PMC free article] [PubMed]
  • Yoshida H., Sasaki K., Namiki Y., Sato N. &Tada N. Edaravone, a novel radical scavenger, inhibits oxidative modification of low-density lipoprotein (LDL) and reverses oxidized LDL-mediated reduction in the expression of endothelial nitric oxide synthase. Atherosclerosis. 179, 97–102 (2005). [PubMed]
  • Sato T., Mizuno K. & Ishii F. A novel administration route of edaravone–II: mucosal absorption of edaravone from edaravone/hydroxypropyl-beta-cyclodextrin complex solution including L-cysteine and sodium hydrogen sulfite. Pharmacology. 85, 88–94 (2010). [PubMed]
  • Rong W. T. et al. . Hydroxypropyl-sulfobutyl-β-cyclodextrin improves the oral bioavailability of edaravone by modulating drug efflux pump of enterocytes. J Pharm Sci. 103, 730–742 (2014). [PubMed]
  • Parikh A., Kathawala K., Tan C. C., Garg S. & Zhou X. F. Development of a novel oral delivery system of edaravone for enhancing bioavailability. Int J Pharm. 515, 490–500 (2016). [PubMed]
  • Otake Y. & Walle T. Oxidation of the flavonoids galangin and kaempferide by human liver microsomes and CYP1A1, CYP1A2, and CYP2C9. Drug Metab Dispos. 30, 103–105 (2002). [PubMed]
  • Okuyama S. et al. . Edaravone enhances brain-derived neurotrophic factor production in the ischemic mouse brain. Pharmaceuticals (Basel). 8, 176–185 (2015). [PMC free article] [PubMed]
  • Roy A. et al. . HMG-CoA Reductase Inhibitors Bind to PPARα to Upregulate Neurotrophin Expression in the Brain and Improve Memory in Mice. Cell Metab. 22, 253–265 (2015). [PMC free article] [PubMed]
  • Daniele S. et al. . Does GRK-β arrestin machinery work as a “switch on” for GPR17-mediated activation of intracellular signaling pathways? Cell Signal. 26, 1310–1325 (2014). [PubMed]
  • Walters W. P. & Namchuk M. Designing screens: how to make your hits a hit. Nat Rev Drug Discov. 2, 259–266 (2013). [PubMed]
  • Agresti C., D’Urso D. & Levi G. Reversible inhibitory effects of interferon-gamma and tumour necrosis factor alpha on oligodendroglial lineage cell proliferation and differentiation in vitro. Eur J Neurosci. 8, 1106–1116 (1996). [PubMed]
  • Zhang H., Jarjour A. A., Boyd A. & Williams A. Central nervous system remyelination in culture: a tool for multiple sclerosis research. Exp Neurol. 230, 138–148 (2011). [PMC free article] [PubMed]
  • Lachmanovich E. et al. . Co-localization analysis of complex formation among membrane proteins by computerized fluorescence microscopy: application to immunofluorescence co-patching studies. J Microsc. 212, 122–131 (2003). [PubMed]
  • Schwikowski B., Uetz P. & Fields S. A network of protein-protein interactions in yeast. Nature Biotechnology. 18, 1257–1261 (2000). [PubMed]
  • Umeton R. et al. . Design of robust metabolic pathways. Proceedings of the 48th Design Automation Conference. 747–752 (2011).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group