1
|
Kench T, Rahardjo A, Terrones GG, Bellamkonda A, Maher TE, Storch M, Kulik HJ, Vilar R. A Semi-Automated, High-Throughput Approach for the Synthesis and Identification of Highly Photo-Cytotoxic Iridium Complexes. Angew Chem Int Ed Engl 2024; 63:e202401808. [PMID: 38404222 DOI: 10.1002/anie.202401808] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Revised: 02/20/2024] [Accepted: 02/21/2024] [Indexed: 02/27/2024]
Abstract
The discovery of new compounds with pharmacological properties is usually a lengthy, laborious and expensive process. Thus, there is increasing interest in developing workflows that allow for the rapid synthesis and evaluation of libraries of compounds with the aim of identifying leads for further drug development. Herein, we apply combinatorial synthesis to build a library of 90 iridium(III) complexes (81 of which are new) over two synthesise-and-test cycles, with the aim of identifying potential agents for photodynamic therapy. We demonstrate the power of this approach by identifying highly active complexes that are well-tolerated in the dark but display very low nM phototoxicity against cancer cells. To build a detailed structure-activity relationship for this class of compounds we have used density functional theory (DFT) calculations to determine some key electronic parameters and study correlations with the experimental data. Finally, we present an optimised semi-automated synthesise-and-test protocol to obtain multiplex data within 72 hours.
Collapse
Affiliation(s)
- Timothy Kench
- Department of Chemistry, Imperial College London, White City Campus, W12 0BZ, London, UK
| | - Arielle Rahardjo
- Department of Chemistry, Imperial College London, White City Campus, W12 0BZ, London, UK
| | - Gianmarco G Terrones
- Department of Chemical Engineering, Massachusetts Institute of Technology, 02139, Cambridge, MA, USA
| | | | - Thomas E Maher
- Department of Chemistry, Imperial College London, White City Campus, W12 0BZ, London, UK
- Institute of Chemical Biology, Imperial College London, White City Campus, W12 0BZ, London, UK
| | - Marko Storch
- Department of Infectious Disease, Imperial College London, South Kensington Campus, SW7 2AZ, London, UK
- London Biofoundry, Imperial College Translation and Innovation Hub, W12 0BZ, London, UK
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, 02139, Cambridge, MA, USA
- Department of Chemistry, Massachusetts Institute of Technology, 02139, Cambridge, MA, USA
| | - Ramon Vilar
- Department of Chemistry, Imperial College London, White City Campus, W12 0BZ, London, UK
- Institute of Chemical Biology, Imperial College London, White City Campus, W12 0BZ, London, UK
| |
Collapse
|
2
|
Vennelakanti V, Kilic IB, Terrones GG, Duan C, Kulik HJ. Machine Learning Prediction of the Experimental Transition Temperature of Fe(II) Spin-Crossover Complexes. J Phys Chem A 2024; 128:204-216. [PMID: 38148525 DOI: 10.1021/acs.jpca.3c07104] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
Spin-crossover (SCO) complexes are materials that exhibit changes in the spin state in response to external stimuli, with potential applications in molecular electronics. It is challenging to know a priori how to design ligands to achieve the delicate balance of entropic and enthalpic contributions needed to tailor a transition temperature close to room temperature. We leverage the SCO complexes from the previously curated SCO-95 data set [Vennelakanti et al. J. Chem. Phys. 159, 024120 (2023)] to train three machine learning (ML) models for transition temperature (T1/2) prediction using graph-based revised autocorrelations as features. We perform feature selection using random forest-ranked recursive feature addition (RF-RFA) to identify the features essential to model transferability. Of the ML models considered, the full feature set RF and recursive feature addition RF models perform best, achieving moderate correlation to experimental T1/2 values. We then compare ML T1/2 predictions to those from three previously identified best-performing density functional approximations (DFAs) which accurately predict SCO behavior across SCO-95, finding that the ML models predict T1/2 more accurately than the best-performing DFAs. In addition, we study ML model predictions for a set of 18 SCO complexes for which only estimated T1/2 values are available. Upon excluding outliers from this set, the RF-RFA RF model shows a strong correlation to estimated T1/2 values with a Pearson's r of 0.82. In contrast, DFA-predicted T1/2 values have large errors and show no correlation to estimated T1/2 values over the same set of complexes. Overall, our study demonstrates slightly superior performance of ML models in comparison with some of the best-performing DFAs, and we expect ML models to improve further as larger data sets of SCO complexes are curated and become available for model training.
Collapse
Affiliation(s)
- Vyshnavi Vennelakanti
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Irem B Kilic
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Gianmarco G Terrones
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
3
|
Duan C, Nandy A, Terrones GG, Kastner DW, Kulik HJ. Active Learning Exploration of Transition-Metal Complexes to Discover Method-Insensitive and Synthetically Accessible Chromophores. JACS Au 2023; 3:391-401. [PMID: 36873700 PMCID: PMC9976347 DOI: 10.1021/jacsau.2c00547] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Revised: 11/15/2022] [Accepted: 11/16/2022] [Indexed: 06/18/2023]
Abstract
Transition-metal chromophores with earth-abundant transition metals are an important design target for their applications in lighting and nontoxic bioimaging, but their design is challenged by the scarcity of complexes that simultaneously have well-defined ground states and optimal target absorption energies in the visible region. Machine learning (ML) accelerated discovery could overcome such challenges by enabling the screening of a larger space but is limited by the fidelity of the data used in ML model training, which is typically from a single approximate density functional. To address this limitation, we search for consensus in predictions among 23 density functional approximations across multiple rungs of "Jacob's ladder". To accelerate the discovery of complexes with absorption energies in the visible region while minimizing the effect of low-lying excited states, we use two-dimensional (2D)efficient global optimization to sample candidate low-spin chromophores from multimillion complex spaces. Despite the scarcity (i.e., ∼0.01%) of potential chromophores in this large chemical space, we identify candidates with high likelihood (i.e., >10%) of computational validation as the ML models improve during active learning, representing a 1000-fold acceleration in discovery. Absorption spectra of promising chromophores from time-dependent density functional theory verify that 2/3 of candidates have the desired excited-state properties. The observation that constituent ligands from our leads have demonstrated interesting optical properties in the literature exemplifies the effectiveness of our construction of a realistic design space and active learning approach.
Collapse
Affiliation(s)
- Chenru Duan
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Chemistry, Massachusetts Institute of
Technology, Cambridge, Massachusetts 02139, United States
| | - Aditya Nandy
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Chemistry, Massachusetts Institute of
Technology, Cambridge, Massachusetts 02139, United States
| | - Gianmarco G. Terrones
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
| | - David W. Kastner
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Biological Engineering, Massachusetts
Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department
of Chemical Engineering, Massachusetts Institute
of Technology, Cambridge, Massachusetts 02139, United States
- Department
of Chemistry, Massachusetts Institute of
Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
4
|
Terrones GG, Duan C, Nandy A, Kulik HJ. Low-cost machine learning prediction of excited state properties of iridium-centered phosphors. Chem Sci 2023; 14:1419-1433. [PMID: 36794185 PMCID: PMC9906783 DOI: 10.1039/d2sc06150c] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Accepted: 01/05/2023] [Indexed: 01/07/2023] Open
Abstract
Prediction of the excited state properties of photoactive iridium complexes challenges ab initio methods such as time-dependent density functional theory (TDDFT) both from the perspective of accuracy and of computational cost, complicating high-throughput virtual screening (HTVS). We instead leverage low-cost machine learning (ML) models and experimental data for 1380 iridium complexes to perform these prediction tasks. We find the best-performing and most transferable models to be those trained on electronic structure features from low-cost density functional tight binding calculations. Using artificial neural network (ANN) models, we predict the mean emission energy of phosphorescence, the excited state lifetime, and the emission spectral integral for iridium complexes with accuracy competitive with or superseding that of TDDFT. We conduct feature importance analysis to determine that high cyclometalating ligand ionization potential correlates to high mean emission energy, while high ancillary ligand ionization potential correlates to low lifetime and low spectral integral. As a demonstration of how our ML models can be used for HTVS and the acceleration of chemical discovery, we curate a set of novel hypothetical iridium complexes and use uncertainty-controlled predictions to identify promising ligands for the design of new phosphors while retaining confidence in the quality of the ANN predictions.
Collapse
Affiliation(s)
- Gianmarco G. Terrones
- Department of Chemical Engineering, Massachusetts Institute of TechnologyCambridgeMA 02139USA
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of TechnologyCambridgeMA 02139USA,Department of Chemistry, Massachusetts Institute of TechnologyCambridgeMA 02139USA
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of TechnologyCambridgeMA 02139USA,Department of Chemistry, Massachusetts Institute of TechnologyCambridgeMA 02139USA
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of TechnologyCambridgeMA 02139USA,Department of Chemistry, Massachusetts Institute of TechnologyCambridgeMA 02139USA
| |
Collapse
|