Machine Learning and Molecular Modelling: A Synergistic Approach to Rapid Reactivity Prediction
Lead Research Organisation:
UNIVERSITY OF BATH
Department Name: Chemistry
Abstract
The computational design of new chemical reactions is regarded as one of the "Holy Grails" of computational organic chemistry and biochemistry. Accurate and fast computational approaches to predicting chemical reactivity would provide cost-effective alternatives to time-consuming experimental approaches, and in some cases animal testing, in drug design, toxicology and chemical synthesis. Of great importance are mechanism-based prediction models because they are much more likely to reach general acceptance compared to computational "black-box" models which offer no insight into how and why predictions are made. Providing such insight is especially important for models to gain regulatory acceptance in toxicology and drug design. However, no current computational approach (molecular modelling or machine learning (ML)) to reactivity prediction offers the combination of fast, accurate predictions with clear mechanistic insight; one or more of these desirable characteristics must be sacrificed in pursuit of the others. This project will develop a novel, synergistic molecular modelling and ML approach to rapid, high-accuracy and mechanism-based reactivity prediction for use in toxicology, drug design and chemical synthesis and thus help realise the "Holy Grail".
We will train and validate ML models on large datasets (~10,000 compounds) that can correct energy barriers obtained from rapid molecular modelling techniques to those derived from prohibitively slow, high-accuracy methods. Our synergistic approach to reaction modelling will thus be to derive mechanistic insight from these rapid molecular modelling techniques and use our ML models to obtain fast and accurate reaction barriers. Models for C-N bond-forming reactions will be developed for use in covalent drug design (targeting lysine), toxicology (predicting mutagenicity and respiratory sensitisation) and pharmaceutical drug synthesis planning. To demonstrate the broad utility of our synergistic approach, we will use it to rationalise experimental reactivity data of biologically and synthetically relevant systems for which the use of current modelling approaches would be prohibitively slow. Rather than requiring a supercomputer, predictions will be possible even on a laptop which will represent a paradigm shift in reaction modelling.
We will train and validate ML models on large datasets (~10,000 compounds) that can correct energy barriers obtained from rapid molecular modelling techniques to those derived from prohibitively slow, high-accuracy methods. Our synergistic approach to reaction modelling will thus be to derive mechanistic insight from these rapid molecular modelling techniques and use our ML models to obtain fast and accurate reaction barriers. Models for C-N bond-forming reactions will be developed for use in covalent drug design (targeting lysine), toxicology (predicting mutagenicity and respiratory sensitisation) and pharmaceutical drug synthesis planning. To demonstrate the broad utility of our synergistic approach, we will use it to rationalise experimental reactivity data of biologically and synthetically relevant systems for which the use of current modelling approaches would be prohibitively slow. Rather than requiring a supercomputer, predictions will be possible even on a laptop which will represent a paradigm shift in reaction modelling.
Organisations
Publications
Beaglehole I
(2025)
Machine Learning Transition State Geometries and Applications in Reaction Property Prediction
in WIREs Computational Molecular Science
Beaglehole I
(2024)
Machine learning for transition state geometry prediction
Beaglehole I
(2024)
Machine learning for transition state geometry prediction
Espley S
(2024)
Distortion/interaction analysis via machine learning
in Digital Discovery
Espley S
(2023)
Machine learning reaction barriers in low data regimes: a horizontal and diagonal transfer learning approach
in Digital Discovery
Espley S
(2025)
Correction: Distortion/interaction analysis via machine learning
in Digital Discovery
Farrar EHE
(2022)
Noncovalent Interactions in the Oxazaborolidine-Catalyzed Enantioselective Mukaiyama Aldol.
in The Journal of organic chemistry
| Description | This award focused on the development of reaction modelling methodology that is fast, accurate, and provides mechanistic insight from transition state (TS) geometries. Such computational tools enable the rational design of new reactions thus helping to reduce the need for trial-and-error experimentation in chemical synthesis and drug discovery. We have reported a hybrid machine learning (ML)-modelling approach which meets all of the above criteria for reaction barrier prediction and thus serves as a high-throughput method for rapid reactivity design. In this work, we trained ML models that corrected rapid but approximate semi-empirical quantum mechanical (SQM) reaction barriers to DFT-quality barriers for reactions of importance to the pharmaceutical industry. This approach enabled barrier predictions to be made in seconds with errors below the chemical accuracy threshold of 1 kcal mol-1. Rapid mechanistic insight was available from the SQM TSs which were found to be good approximations of the DFT geometries. We also showed that the amount of expensive data needed to develop these reactivity prediction models can be reduced from thousands to just a few tens of data points which was particularly important when working with large and complex systems of high synthetic importance. |
| Exploitation Route | We have published our results open access to reach the broadest possible audience. It is hoped that our methods will be used widely within the academic and industrial chemical research communities to aid in reactivity design and thus reduce trial-and-error experimentation. It is likely that our low data approaches will be of particular importance because they offer the cheapest route to the routine use of computational methods in the design of experimental work. |
| Sectors | Chemicals Digital/Communication/Information Technologies (including Software) Pharmaceuticals and Medical Biotechnology |
| Description | The award focused on developing computational tools to reduce the need for trial-and-error experimentation in chemical synthesis and drug discovery. The work has begun to influence the UK pharmaceutical industry; this is best evidenced by the recruitment of a team member into a job in this industry as a result of the skills and knowledge they developed during the award. |
| First Year Of Impact | 2024 |
| Sector | Pharmaceuticals and Medical Biotechnology |
| Impact Types | Economic |
| Title | Dataset for "Distortion/Interaction Analysis via Machine Learning" |
| Description | Machine learning (ML) has previously been applied to predict reaction barriers for a variety of different chemical reactions. This is seen as the end point for this type of study however, post-reaction barrier analysis/energy decomposition approaches can provide insight into chemical reactivity. One such approach that has previously been used to provide information on chemical reactivity, for cycloaddition reactions in particular, is distortion/interaction-activation strain analysis (DIAS). We demonstrate that ML can be coupled with cheap and rapid semi-empirical quantum mechanical methods (SQM) to predict distortion and interaction energies at a fraction of the computational cost associated with running density functional theory (DFT) calculations. This dataset includes all the structural data in the form of Gaussian16 (Revision A.03 and C.01) output files for the four datasets used in this work and, the literature dataset reactions. |
| Type Of Material | Database/Collection of data |
| Year Produced | 2025 |
| Provided To Others? | Yes |
| Impact | It supports our paper that is mentioned in the title |
| URL | https://researchdata.bath.ac.uk/id/eprint/1480 |
| Title | Dataset for "Distortion/Interaction Analysis via Machine Learning" |
| Description | Machine learning (ML) has previously been applied to predict reaction barriers for a variety of different chemical reactions. This is seen as the end point for this type of study however, post-reaction barrier analysis/energy decomposition approaches can provide insight into chemical reactivity. One such approach that has previously been used to provide information on chemical reactivity, for cycloaddition reactions in particular, is distortion/interaction-activation strain analysis (DIAS). We demonstrate that ML can be coupled with cheap and rapid semi-empirical quantum mechanical methods (SQM) to predict distortion and interaction energies at a fraction of the computational cost associated with running density functional theory (DFT) calculations. This dataset includes all the structural data in the form of Gaussian16 (Revision A.03 and C.01) output files for the four datasets used in this work and, the literature dataset reactions. |
| Type Of Material | Database/Collection of data |
| Year Produced | 2024 |
| Provided To Others? | Yes |
| Impact | It supports our paper that is mentioned in the title |
| URL | https://researchdata.bath.ac.uk/id/eprint/1398 |
| Title | Dataset for "Machine learning reaction barriers in low data regimes: a horizontal and diagonal transfer learning approach" |
| Description | Machine learning (ML) has previously been used to predict density functional theory (DFT) free energy reaction barriers on a variety of different reactions from semi-empirical quantum mechanical (SQM) inputs. These models can require expensive dataset curation and can struggle with generalisability outside of the datasets immediate chemical space. One such approach that can drastically lower the number of required training points is transfer learning (TL). We demonstrate that various TL techniques can be used to provide highly accurate results with a fraction of the training points required for standard ML, thus lowering the overall computational cost of barrier predictions. This dataset includes all the structural data in the form of Gaussian16 (Revision A.03 and C.01) output files for the Diels-Alder and [3+2] cycloaddition reactions used for this ML/TL analysis. This data archive also includes exemplar code for performing some standard ML from the manuscript. |
| Type Of Material | Database/Collection of data |
| Year Produced | 2023 |
| Provided To Others? | Yes |
| Impact | It supports our paper that is mentioned in the title |
| URL | https://researchdata.bath.ac.uk/id/eprint/1229 |
| Title | Dataset for "Reformulating Reactivity Design for Data-Efficient Machine Learning" |
| Description | This dataset contains the Gaussian 16 output files for the dataset of aza-Michael addition reactions used in the publication "Fast Identification of Reactions with Desired Barriers by Reformulating Machine Learning Activation Energies". The structures of the methylamine nucleophile, the 1000 Michael acceptor electrophiles and their 1000 transition states were all optimised at the wB97X-D/def2-TZVP level of theory with the IEFPCM(water) implicit solvent model. Before optimisation all Michael acceptors and transition states were conformationally searched using the MMFF force field in Schrödinger's MacroModel software and the lowest energy conformer was selected for DFT calculation. This dataset also contains the Gaussian 16 output files for the SVWN/def2-SVP single-point energy calculations on the dihydrogen activation catalyst and transition state structures. |
| Type Of Material | Database/Collection of data |
| Year Produced | 2023 |
| Provided To Others? | Yes |
| Impact | It supports our paper that is mentioned in the title |
| URL | https://researchdata.bath.ac.uk/id/eprint/1240 |
| Title | Dataset for "The Catalytic Enantioselective [1,2]-Wittig Rearrangement Cascade of Allylic Ethers" |
| Description | This data set includes output files from the quantum chemical calculations run with Gaussian16 (Revision C.01) that support our computational mechanistic study of the enantioselective [1,2]-Wittig rearrangement of allylic ethers. It also contains three sets of in situ reaction monitoring data (collected by University of St Andrews contributors) and a Python script that fits the rate constants of a first-order kinetics model to the experimental data. |
| Type Of Material | Database/Collection of data |
| Year Produced | 2025 |
| Provided To Others? | Yes |
| Impact | It supports our paper that is mentioned in the title |
| URL | https://researchdata.bath.ac.uk/id/eprint/1337 |
| Description | Interdisciplinary Computer Aided Synthesis (ICAS) Symposium |
| Form Of Engagement Activity | Participation in an activity, workshop or similar |
| Part Of Official Scheme? | No |
| Geographic Reach | National |
| Primary Audience | Industry/Business |
| Results and Impact | We organised a two-day symposium for interdisciplinary computer-aided synthesis (ICAS) in September 2023 in Bath, hosting 23 synthetic and computational chemists from academic and industrial institutions across the UK. The aim of this symposium was to bring together chemists from several different disciplines and backgrounds to facilitate discussion, generate ideas for academic-industrial collaboration, and to understand different perspectives on computer-aided synthesis. During this time, the work associated with the grant was presented and discussed, as well as many other related pieces of work from other academic and industrial institutions. The event was an overall success and received overwhelmingly positive feedback from participants. The research insights obtained from the event regarding this research project have allowed us to understand more about the impact of our work, as well as the best and most efficient way to achieve it. It is also very beneficial to understand where others institutions may be working towards common goals and to educate others on what we are trying to achieve and how that could potentially be useful for them. Finally, the event has also allowed us to expand our scientific network, affording further opportunities to develop our research and collaroate with others in the future. |
| Year(s) Of Engagement Activity | 2023 |
