Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

Computer-designed repurposing of chemical wastes into drugs

Abstract

As the chemical industry continues to produce considerable quantities of waste chemicals1,2, it is essential to devise ‘circular chemistry’3,4,5,6,7,8 schemes to productively back-convert at least a portion of these unwanted materials into useful products. Despite substantial progress in the degradation of some classes of harmful chemicals9, work on ‘closing the circle’—transforming waste substrates into valuable products—remains fragmented and focused on well known areas10,11,12,13,14,15. Comprehensive analyses of which valuable products are synthesizable from diverse chemical wastes are difficult because even small sets of waste substrates can, within few steps, generate millions of putative products, each synthesizable by multiple routes forming densely connected networks. Tracing all such syntheses and selecting those that also meet criteria of process and ‘green’ chemistries is, arguably, beyond the cognition of human chemists. Here we show how computers equipped with broad synthetic knowledge can help address this challenge. Using the forward-synthesis Allchemy platform16, we generate giant synthetic networks emanating from approximately 200 waste chemicals recycled on commercial scales, retrieve from these networks tens of thousands of routes leading to approximately 300 important drugs and agrochemicals, and algorithmically rank these syntheses according to the accepted metrics of sustainable chemistry17,18,19. Several of these routes we validate by experiment, including an industrially realistic demonstration on a ‘pharmacy on demand’ flow-chemistry platform20. Wide adoption of computerized waste-to-valuable algorithms can accelerate productive reuse of chemicals that would otherwise incur storage or disposal costs, or even pose environmental hazards.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: Examples of small molecules recycled from various types of industrial chemical waste.
Fig. 2: Generation and properties of ‘forward’ synthetic networks.
Fig. 3: Examples of highly ranked syntheses of more advanced drugs starting from waste substrates and few simple, auxiliary molecules used frequently in organic synthesis.
Fig. 4: Experimental validation of selected, computer-designed pathways in laboratory scale syntheses.
Fig. 5: Syntheses of COVID-19 intensive care unit medications or their intermediates performed on an automated, modular ODP platform.

Similar content being viewed by others

Data availability

All data in support of the findings of this study are available within the Article and its Supplementary Information. Syntheses of all drugs and selected most interesting agrochemicals identified in large-scale network searches are available for analysis and re-ranking at https://wasteresults.allchemy.net. User manuals are available in Supplementary Information section 2.

Code availability

The interactive Allchemy web application is freely available for academic users at https://waste.allchemy.net (given server capacity, to five concurrent academic users on a rolling basis and two-week slots). Reaction rules and the source code of Allchemy are proprietary.

References

  1. Yavrom, D. An Overview of Hazardous Waste Generation (EPA, accessed 28 April 2021); https://rcrapublic.epa.gov/rcra-public-web/action/posts/2

  2. Production-related Waste Managed by Chemical (EPA, accessed 1 July 2021); https://www.epa.gov/trinationalanalysis/waste-managed-chemical-and-industry

  3. Stahel, W. R. The circular economy. Nature 531, 435–438 (2016).

    Article  CAS  ADS  Google Scholar 

  4. Ellen MacArthur Foundation, World Economic Forum & McKinsey & Company. The New Plastics Economy: Rethinking the Future of Plastics (Ellen MacArthur Foundation, 2016).

  5. Winans, K., Kendall, A. & Deng, H. The history and current applications of the circular economy concept. Renew. Sust. Ener. Rev. 68, 825–833 (2017).

    Article  Google Scholar 

  6. Keijer, T., Bakker, V. & Slootweg, J. C. Circular chemistry to enable a circular economy. Nat. Chem. 11, 190–195 (2019).

    Article  CAS  Google Scholar 

  7. Kümmerer, K., Clark, J. H. & Zuin, V. G. Rethinking chemistry for a circular economy. Science 367, 369–370 (2020).

    Article  ADS  Google Scholar 

  8. Kümmerer, K. Sustainable chemistry: a future guiding principle. Angew. Chem. Int. Ed. 56, 16420–16421 (2017).

    Article  Google Scholar 

  9. Tullo, A. H. Plastic has a problem; is chemistry the solution? Chem. Eng. News 97, 29–34 (2019).

    Google Scholar 

  10. Zeng, H. & Li, C.-J. Conversion of lignin into high value chemical products. In Green Chemistry and Chemical Engineering (eds Han, B. & Wu, T.) 385–403 (Springer, 2018).

  11. Sun, Z., Balint, F., de Santi, A., Saravanakumar, E. & Barta, K. Bright side of lignin depolymerization: toward new platform chemicals. Chem. Rev. 118, 614–678 (2018).

    Article  CAS  Google Scholar 

  12. Park, C. & Lee, J. Recent achievements in CO2-assisted and CO2-catalyzed biomass conversion reactions. Green Chem. 22, 2628–2642 (2020).

    Article  CAS  Google Scholar 

  13. Antonetti, C., Licursi, D., Fulignati, S., Valentini, G. & Raspolli Galletti, A. M. New frontiers in the catalytic synthesis of levulinic acid: from sugars to raw and waste biomass as starting feedstock. Catalysts 6, 196 (2016).

    Article  Google Scholar 

  14. Dabral, S. & Schaub, T. The use of carbon dioxide (CO2) as a building block in organic synthesis from an industrial perspective. Adv. Synth. Catal. 361, 223–246 (2018).

    Article  Google Scholar 

  15. Zhang, F. et al. Polyethylene upcycling to long-chain alkylaromatics by tandem hydrogenolysis/aromatization. Science 370, 437–441 (2020).

    Article  CAS  ADS  Google Scholar 

  16. Wołos, A. et al. Synthetic connectivity, emergence, and autocatalysis in the network of prebiotic chemistry. Science 369, eaaw1955 (2020).

    Article  Google Scholar 

  17. Adams, J. P. et al. Development of GSK’s reagent guides – embedding sustainability into reagent selection. Green Chem. 15, 1542 (2013).

    Article  CAS  Google Scholar 

  18. Henderson, R. K., Hill, A. P., Redman, A. M. & Sneddon, H. F. Development of GSK’s acid and base selection guides. Green Chem. 17, 945–949 (2015).

    Article  CAS  Google Scholar 

  19. Henderson, R. K. et al. Expanding GSK’s solvent selection guide – embedding sustainability into solvent selection starting at medicinal chemistry. Green Chem. 13, 854 (2011).

    Article  CAS  Google Scholar 

  20. Rogers, L. et al. Continuous production of five active pharmaceutical ingredients in flexible plug-and-play modules: a demonstration campaign. Org. Process Res. Dev. 24, 2183–2196 (2020).

    Article  CAS  Google Scholar 

  21. Brown, D. G. & Boström, J. Analysis of past and present synthetic methodologies on medicinal chemistry: where have all the new reactions gone? J. Med. Chem. 59, 4443–4458 (2015).

    Article  Google Scholar 

  22. Roughley, S. D. & Jordan, A. M. The medicinal chemist’s toolbox: an analysis of reactions used in the pursuit of drug candidates. J. Med. Chem. 54, 3451–3479 (2011).

    Article  CAS  Google Scholar 

  23. Molga, K., Gajewska, E. P., Szymkuć, S. & Grzybowski, B. A. The logic of translating chemical knowledge into machine-processable forms: a modern playground for physical-organic chemistry. React. Chem. Eng. 4, 1506–1521 (2019).

    Article  CAS  Google Scholar 

  24. Segler, M. H. S., Preuss, M. & Waller, M. P. Planning chemical syntheses with deep neural networks and symbolic AI. Nature 555, 604–610 (2018).

    Article  CAS  ADS  Google Scholar 

  25. Szymkuć, S. et al. Computer-assisted synthetic planning: the end of the beginning. Angew. Chem. Int. Ed. 55, 5904–5937 (2016).

    Article  Google Scholar 

  26. Klucznik, T. et al. Efficient syntheses of diverse, medicinally relevant targets planned by computer and executed in the laboratory. Chem 4, 522–532 (2018).

    Article  CAS  Google Scholar 

  27. Gajewska, E. P. et al. Algorithmic discovery of tactical combinations for advanced organic syntheses. Chem 6, 280–293 (2020).

    Article  CAS  Google Scholar 

  28. Mikulak-Klucznik, B. et al. Computational planning of the synthesis of complex natural products. Nature 588, 83–88 (2020).

    Article  CAS  ADS  Google Scholar 

  29. Molga, K., Dittwald, P. & Grzybowski, B. A. Computational design of syntheses leading to compound libraries or isotopically labelled targets. Chem. Sci. 10, 9219–9232 (2019).

    Article  CAS  Google Scholar 

  30. Electronic Code of Federal Regulations (eCFR, accessed 1 July 2021); https://www.ecfr.gov/cgi-bin/textidx?SID=2b4d2d375e73ebc5c93d8b2fe632cb6f&mc=true&node=pt40.28.355&rgn=div

  31. Candidate List of Substances of Very High Concern for Authorisation (ECHA, accessed 1 September 2021); https://echa.europa.eu/candidate-list-table

  32. Trost, B. M. Atom economy—a challenge for organic synthesis. Angew. Chem. Int. Ed. Eng. 34, 259–281 (1995).

    Article  CAS  Google Scholar 

  33. Benson, S. W. & Buss, J. H. Additivity rules for the estimation of molecular properties. Thermodynamic properties. J. Chem. Phys. 29, 546–572 (1958).

    Article  CAS  ADS  Google Scholar 

  34. Bishop, K. J. M., Klajn, R. & Grzybowski, B. A. The core and most useful molecules in organic chemistry. Angew. Chem. Int. Ed. 45, 5348–5354 (2006).

    Article  CAS  Google Scholar 

  35. Fialkowski, M., Bishop, K. J. M., Chubukov, V. A., Campbell, C. J. & Grzybowski, B. A. Architecture and evolution of organic chemistry. Angew. Chem. Int. Ed. 44, 7263–7269 (2005).

    Article  CAS  Google Scholar 

  36. Wishart, D. S. et al. DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucl. Acids Res. 36, D901–D906 (2007).

    Article  Google Scholar 

  37. Rogers, D. J. & Tanimoto, T. T. A computer program for classifying plants. Science 132, 1115–1118 (1960).

    Article  CAS  ADS  Google Scholar 

  38. Skoraczyński, G. et al. Predicting the outcomes of organic reactions via machine learning: are current descriptors sufficient? Sci. Rep. 7, 3582 (2017).

    Article  ADS  Google Scholar 

  39. Li, J., Albrecht, J., Borovika, A. & Eastgate, M. D. Evolving green chemistry metrics into predictive tools for decision making and benchmarking analytics. ACS Sustain. Chem. Eng. 6, 1121–1132 (2017).

    Article  Google Scholar 

  40. Borovika, A. et al. The PMI Predictor app to enable green-by-design chemical synthesis. Nat. Sustain. 2, 1034–1040 (2019).

    Article  Google Scholar 

  41. Kane, S. P. The Top 300 of 2021 (ClinCalc, accessed 1 July 2021); https://clincalc.com/DrugStats/Top300Drugs.aspx

  42. Resilient Drug Supply Project (Center for Infectious Disease Research and Policy, the University of Minnesota, accessed 1 July 2021); https://www.cidrap.umn.edu/sites/default/files/public/downloads/cidrap-rds-drug-shortages.pdf

  43. Rogers, L., et. al.  Continuous production of five active pharmaceutical ingredients in flexible plug-and-play modules: A demonstration campaign, Org. Proc. Res. Dev.24, 2183-2196 (2020).

    CAS  Google Scholar 

  44. Tengzhou Wutong Spice Co. Ltd. Reaction kettle device applicable to producing veratraldehyde and derivatives thereof. Chinese patent 203170325U (2013).

  45. Guilin Teachers Technical College. Preparation method for aryl acetic acid derivative. Chinese patent 102070433A (2013).

  46. Paterson, J. Poddutoori, P. & Romakh, V. Mechanism for production of biobased products from plant lignin. W.O. patent 2013/173316A1 (2013).

  47. Dunn, R. O. Separation of chloronitrobenzene isomers by crystallization and fractionation. US patent 3311666A (1967).

  48. Liaoning Shuntong Chemical Co. Ltd. A kind of preparation method of parachloroanilinum hydrochloride. Chinese patent 110467533A (2019).

  49. Choi, W. J., Byun, J. W., Ahn, J. H., Ha, Y. W. & Seo, J.-H. Process of biologically producing a p-hydroxybenzoic acid. US patent 9206449B2 (2015).

  50. Sheldon, R. A. The E factor 25 years on: the rise of green chemistry and sustainability. Green Chem. 19, 18–43 (2017).

    Article  CAS  Google Scholar 

  51. Turner, N. J. & O’Reilly, E. Biocatalytic retrosynthesis. Nat. Chem. Biol. 9, 285–288 (2013).

    Article  CAS  Google Scholar 

  52. Sheldon, R. A. & Woodley, J. M. Role of biocatalysis in sustainable chemistry. Chem. Rev. 118, 801–838 (2018).

    Article  CAS  Google Scholar 

  53. Oakes, F. T. & Leonard, N. J. Broadened scope of translocative rearrangements. Substituted 1,2,3-triazolo[1,5-a]-1,3,5-triazines. J. Org. Chem. 50, 4986–4989 (1985).

    Article  CAS  Google Scholar 

  54. Chatterjee, M., Ishizaka, T. & Kawanami, H. Reductive amination of furfural to furfurylamine using aqueous ammonia solution and molecular hydrogen: an environmentally friendly approach. Green Chem. 18, 487–496 (2016).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Development of the medicinal chemistry modules within the Allchemy platform (by A.W., R.R., S.S., M.M. and B.A.G.) has been supported by internal funds of Allchemy, Inc. D.K. and R.O. gratefully acknowledge funding from the National Science Centre, Poland (award 2016/23/B/ST5/03307) which supported the laboratory-scale syntheses described in this paper. Analysis of pathways and writing of the paper by B.A.G. was supported by the Institute for Basic Science, Korea (project code IBS-R020-D1). Development of the manufacturing processes for COVID-19 medications was supported by DARPA & CARES Act grant (HR0011-16-2-0029). The views, opinions and/or findings expressed are those of the author(s) and should not be interpreted as representing the official views or policies of the Department of Defense or the US Government.

Author information

Authors and Affiliations

Authors

Contributions

A.W., R.R., S.S., M.M. and B.A.G. designed and developed the Allchemy platform and performed the analyses and calculations described in the paper. D.K. performed the syntheses described in Fig. 4, with supervision from R.O. B.T.H. and J.S. developed the process for cisatracurium intermediate; G.B., J.M.M. and B.T.H. developed the process for propofol; and J.A.M.L., D.T.M. and L.R. developed the process for midazolam; all these are described in Fig. 5. B.A.G. conceived and supervised the research and wrote the paper with help from other authors.

Corresponding author

Correspondence to Bartosz A. Grzybowski.

Ethics declarations

Competing interests

A.W., R.R., S.S., M.M. and B.A.G. are consultants and/or stakeholders of Allchemy, Inc. Allchemy software is the property of Allchemy, Inc., USA. On Demand Pharmaceuticals (ODP) is planning to file for US Food and Drug Administration (FDA) approvals on the flow processes described in this study. All queries about access options to Allchemy, including academic collaborations, should be sent to saraszymkuc@allchemy.net. On Demand Pharmaceuticals inquiries can be sent to lrogers@ondemandpharma.com.

Peer review

Peer review information

Nature thanks Fabrice Gallou, Frank Roschanger and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Additional information

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

Extended data figures and tables

Extended Data Fig. 1 Statistics of the numbers of synthetic pathways.

a, Distribution of the numbers of synthetic paths for 69 drugs obtained from 189 waste molecules (as in Extended Data Fig. 4) within 7 synthetic generations. b, Distribution of number of synthetic paths for 127 drugs obtained from 189 waste molecules and 1,000 auxiliary, popular reagents (as in Fig. 3) within 8 synthetic generations. c, Histogram plotting the average number of pathways for drug targets of given molecular weights (mass ranges without drugs are omitted for clarity). Note the logarithmic vertical scale.

Extended Data Fig. 2 Arcs tracing 52 syntheses of the allyl alcohol agrochemical.

For large number of syntheses, the arc representation becomes messy (for example, some arcs extend beyond the window view)—nevertheless, this viewing modality enables rapid assessment of key intermediates (that is, hub nodes shared between different arcs).

Extended Data Fig. 3 Ranking of pathways according to process criteria.

Allchemy identified 42 pathways, originating from 11 waste substrates and leading to acetaminophen. a, Screenshot shows the settings of the Cost function ranking these syntheses only for the overall length of synthesis (that is, assigning high cost for the execution of each step). b, The shortest route, but it entails Friedel–Crafts acylation using toxic and environmentally problematic DCM solvent (in red font) for which no ‘greener’ replacement is suggested. c, Another scoring scheme, this time assigning more penalty for problematic solvents and reagents (the checkboxes in the bottom assign even higher penalties, +200, to assure that any pathways violating these conditions will be at the very bottom of the rankings; here extremely toxic solvents are penalized). d, The top-ranking synthesis is one step longer but does not involve Friedel–Crafts acylation in DCM. Similar rankings are available for all syntheses generated by the web app (https://waste.allchemy.net) and for results of large-scale calculations from the text available at https://wasteresults.allchemy.net. For the use of these software programmes, please see the user manual in Supplementary Information section 2 and watch Supplementary Video 1.

Extended Data Fig. 4 High-ranking syntheses of nine drugs and four agrochemicals entailing three or more steps and starting solely from recycled ‘waste’ molecules.

The waste substrates are shown in red in the inner circle. Small stars indicate geographical locations (Europe, Asia, North America; see Fig. 1 and Supplementary Information section 1) at which companies producing and/or recycling these substrates are located. Larger stars next to some of the drugs and agrochemicals indicate that they can be synthesized from waste substrates available at the same geographical location. For instance, the drug hexoprenaline can be made from waste recycled solely in Europe (both of its substrates, muconic acid and 3,4-dihydroxyphenylglycol, are denoted by small light-blue stars) and is thus marked with a large light-blue star. The agrochemical triethanolamine is denoted by two large stars (light blue, Europe; orange, North America) because it can be made from the same-continent wastes either in Europe or in North America (one of the waste substrates, ethylene glycol, is recycled on all three continents and, accordingly, is denoted by three small stars; the other substrate, formaldehyde, is recycled at industrial scales only in Europe and North America and thus has two small stars). Synthetic pathways to different targets are differentiated by colours. Within each pathway, reaction arrows for steps already reported in the literature are coloured whereas those without literature precedent are in black. Hazardous substances30,31 are marked by yellow ellipses (for example, formaldehyde). Details of all syntheses shown in this figure as well as other routes of each target are available at https://wasteresults.allchemy.net. Note that all pathways are top-scoring with exception of that for dapsone for which the software also found a two-step route starting from chlorobenzene waste, sulfuric acid and ammonia; however, this pathway is already known and patented and is not shown.

Extended Data Fig. 5 Additional details of highly ranked syntheses of more advanced drugs starting from waste substrates and few simple, auxiliary molecules used frequently in organic synthesis.

The figure accompanies and extends Fig. 3. All colour-coding schemes are the same. Details of all syntheses shown in this figure as well as other routes of each target are available at https://wasteresults.allchemy.net. HNPhth, phthalimide.

Extended Data Fig. 6 Two syntheses of Cysview top-ranked by different Cost functions.

Route traced by blue arrows was scored for overall length—it is concise but entails several harmful reagents (in red font) and intermediates (allyl alcohol in a yellow ellipse). Path traced by violet arrows is less convergent but uses greener reaction conditions (for details, see text).

Supplementary information

Supplementary Information

This file contains Supplementary Figures S1–S67, Supplementary Tables S1–S11 and Supplementary References.

Supplementary Video 1

The video illustrates how to perform calculations and analyse results in the web application available at https://waste.allchemy.net. This web application is freely available to academic users as specified under Code availability. To register a new account, send an e-mail to admin@allchemy.net from your academic address.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wołos, A., Koszelewski, D., Roszak, R. et al. Computer-designed repurposing of chemical wastes into drugs. Nature 604, 668–676 (2022). https://doi.org/10.1038/s41586-022-04503-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-022-04503-9

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing