1
|
Bukowy T, Brown ML, Popelier PLA. Toward Gaussian Process Regression Modeling of a Urea Force Field. J Phys Chem A 2024; 128:8551-8560. [PMID: 39303098 DOI: 10.1021/acs.jpca.4c04117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/22/2024]
Abstract
FFLUX is a next-generation, machine-learnt force field built on three cornerstones: quantum chemical topology, Gaussian process regression, and (high-rank) multipolar electrostatics. It is capable of performing molecular dynamics with near-quantum accuracy at a lower computational cost than standard ab initio molecular dynamics. Previous work with FFLUX was concerned with water and formamide. In this study, we go one step further and challenge FFLUX to model urea, a larger and more flexible system. In result, we have trained urea models at the B3LYP/aug-cc-pVTZ level of theory, with a mean absolute error of 0.4 kJ mol-1 and a maximum prediction error below 7.0 kJ mol-1. To test their performance in molecular dynamics simulations, two sets of FFLUX geometry optimizations were carried out: 5 dimers corresponding to energy minima and 75 random dimers. The 5 dimers were recovered with a root-mean-square deviation below 0.1 Å with respect to their ab initio references. Out of the 75 random dimers, 68% converged to the qualitatively same dimer as those obtained at the ab initio level. Furthermore, we have ranked the 5 FFLUX-optimized dimers in the order of their relative FFLUX single-point energies and compared them with the ab initio method. The energy ranking fully agreed but for one crossover between two successive minima. Finally, we have demonstrated the importance of geometry-dependent (i.e., flexible) multipole moments, showing that the lack of multipole moment flexibility can lead to average errors in the total intermolecular electrostatic energy of more than 2 orders of magnitude.
Collapse
Affiliation(s)
- Tomasz Bukowy
- Department of Chemistry, University of Manchester, Manchester M13 9PL, Great Britain
| | - Matthew L Brown
- Department of Chemistry, University of Manchester, Manchester M13 9PL, Great Britain
| | - Paul L A Popelier
- Department of Chemistry, University of Manchester, Manchester M13 9PL, Great Britain
| |
Collapse
|
2
|
Brown M, Isamura BK, Skelton JM, Popelier PLA. Incorporating Noncovalent Interactions in Transfer Learning Gaussian Process Regression Models for Molecular Simulations. J Chem Theory Comput 2024; 20:5994-6008. [PMID: 38981081 PMCID: PMC11270819 DOI: 10.1021/acs.jctc.4c00402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2024] [Revised: 05/30/2024] [Accepted: 06/18/2024] [Indexed: 07/11/2024]
Abstract
FFLUX is a quantum chemical topology-based multipolar force field that uses Gaussian process regression machine learning models to predict atomic energies and multipole moments on the fly for fast and accurate molecular dynamics simulations. These models have previously been trained on monomers, meaning that many-body effects, for example, intermolecular charge transfer, are missed in simulations. Moreover, dispersion and repulsion have been modeled using Lennard-Jones potentials, necessitating careful parametrization. In this work, we take an important step toward addressing these shortcomings and show that models trained on clusters, in this case, a dimer, can be used in FFLUX simulations by preparing and benchmarking a formamide dimer model. To mitigate the computational costs associated with training higher-dimensional models, we rely on the transfer of hyperparameters from a smaller source model to a larger target model, enabling an order of magnitude faster training than with a direct learning approach. The dimer model allows for simulations that account for two-body effects, including intermolecular polarization and charge penetration, and that do not require nonbonded potentials. We show that addressing these limitations allows for simulations that are closer to quantum mechanics than previously possible with the monomeric models.
Collapse
Affiliation(s)
- Matthew
L. Brown
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| | - Bienfait K. Isamura
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| | - Jonathan M. Skelton
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| | - Paul L. A. Popelier
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United
Kingdom
| |
Collapse
|
3
|
Wang L, Wu Y, Hu J, Yin D, Wei W, Wen J, Chen X, Gao C, Zhou Y, Liu J, Hu G, Li X, Wu J, Zhou Z, Liu L, Song W. Unlocking the function promiscuity of old yellow enzyme to catalyze asymmetric Morita-Baylis-Hillman reaction. Nat Commun 2024; 15:5737. [PMID: 38982157 PMCID: PMC11233575 DOI: 10.1038/s41467-024-50141-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/29/2024] [Accepted: 07/02/2024] [Indexed: 07/11/2024] Open
Abstract
Exploring the promiscuity of native enzymes presents a promising strategy for expanding their synthetic applications, particularly for catalyzing challenging reactions in non-native contexts. In this study, we explore the promiscuous potential of old yellow enzymes (OYEs) to facilitate the Morita-Baylis-Hillman reaction (MBH reaction), leveraging substrate similarities between MBH reaction and reduction reaction. Using mass spectrometry and spectroscopic techniques, we confirm promiscuity of GkOYE in both MBH and reduction reactions. By blocking H- and H+ transfer pathways, we engineer GkOYE.8, which loses its reduction ability but enhances its MBH activity. The structural basis of MBH reaction catalyzed by GkOYE.8 is obtained through mutation studies and kinetic simulations. Furthermore, enantiocomplementary mutants GkOYE.11 and GkOYE.13 are obtained by directed evolution, exhibiting the ability to accept various aromatic aldehydes and alkenes as substrates. This study demonstrates the potential of leveraging substrate similarities to unlock enzyme functionalities, enabling the catalysis of new-to-nature reactions.
Collapse
Affiliation(s)
- Lei Wang
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Yaoyun Wu
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Jun Hu
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Dejing Yin
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
| | - Wanqing Wei
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Jian Wen
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Xiulai Chen
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Cong Gao
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Yiwen Zhou
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Jia Liu
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Guipeng Hu
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Xiaomin Li
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Jing Wu
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Zhi Zhou
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China
| | - Liming Liu
- School of Biotechnology, Jiangnan University, Wuxi, 214122, China
- Key Laboratory of Industrial Biotechnology of Ministry of Education, Jiangnan University, Wuxi, 214122, China
| | - Wei Song
- School of Life Sciences and Health Engineering, Jiangnan University, Wuxi, 214122, China.
| |
Collapse
|
4
|
Wu Y, Cui Y, Song W, Wei W, He Z, Tao J, Yin D, Chen X, Gao C, Liu J, Liu L, Wu J. Reprogramming the Transition States to Enhance C-N Cleavage Efficiency of Rhodococcus opacusl-Amino Acid Oxidase. JACS AU 2024; 4:557-569. [PMID: 38425913 PMCID: PMC10900486 DOI: 10.1021/jacsau.3c00672] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/30/2023] [Revised: 12/21/2023] [Accepted: 12/26/2023] [Indexed: 03/02/2024]
Abstract
l-Amino acid oxidase (LAAO) is an important biocatalyst used for synthesizing α-keto acids. LAAO from Rhodococcus opacus (RoLAAO) has a broad substrate spectrum; however, its low total turnover number limits its industrial use. To overcome this, we aimed to employ crystal structure-guided density functional theory calculations and molecular dynamic simulations to investigate the catalytic mechanism. Two key steps were identified: S → [TS1] in step 1 and Int1 → [TS2] in step 2. We reprogrammed the transition states [TS1] and [TS2] to reduce the identified energy barrier and obtain a RoLAAO variant capable of catalyzing 19 kinds of l-amino acids to the corresponding high-value α-keto acids with a high total turnover number, yield (≥95.1 g/L), conversion rate (≥95%), and space-time yields ≥142.7 g/L/d in 12-24 h, in a 5 L reactor. Our results indicated the promising potential of the developed RoLAAO variant for use in the industrial production of α-keto acids while providing a potential catalytic-mechanism-guided protein design strategy to achieve the desired physical and catalytic properties of enzymes.
Collapse
Affiliation(s)
- Yaoyun Wu
- School
of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
- School
of Biotechnology, Jiangnan University, Wuxi 214122, China
| | - Yaozhong Cui
- Jiangsu
Xishan Senior High School, Wuxi 214174, China
| | - Wei Song
- School
of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China
| | - Wanqing Wei
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
| | - Zhizhen He
- School
of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China
| | - Jinyang Tao
- School
of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China
| | - Dejing Yin
- School
of Biotechnology, Jiangnan University, Wuxi 214122, China
| | - Xiulai Chen
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
| | - Cong Gao
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
| | - Jia Liu
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
| | - Liming Liu
- State
Key Laboratory of Food Science and Resources, Jiangnan University, Wuxi 214122, China
| | - Jing Wu
- School
of Life Sciences and Health Engineering, Jiangnan University, Wuxi 214122, China
| |
Collapse
|
5
|
Brown M, Skelton JM, Popelier PLA. Application of the FFLUX Force Field to Molecular Crystals: A Study of Formamide. J Chem Theory Comput 2023; 19:7946-7959. [PMID: 37847867 PMCID: PMC10653110 DOI: 10.1021/acs.jctc.3c00578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Indexed: 10/19/2023]
Abstract
In this work, we present the first application of the quantum chemical topology force field FFLUX to the solid state. FFLUX utilizes Gaussian process regression machine learning models trained on data from the interacting quantum atom partitioning scheme to predict atomic energies and flexible multipole moments that change with geometry. Here, the ambient (α) and high-pressure (β) polymorphs of formamide are used as test systems and optimized using FFLUX. Optimizing the structures with increasing multipolar ranks indicates that the lattice parameters of the α phase differ by less than 5% to the experimental structure when multipole moments up to the quadrupole are used. These differences are found to be in line with the dispersion-corrected density functional theory. Lattice dynamics calculations are also found to be possible using FFLUX, yielding harmonic phonon spectra comparable to dispersion-corrected DFT while enabling larger supercells to be considered than is typically possible with first-principles calculations. These promising results indicate that FFLUX can be used to accurately determine properties of molecular solids that are difficult to access using DFT, including the structural dynamics, free energies, and properties at finite temperature.
Collapse
Affiliation(s)
- Matthew
L. Brown
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, Britain
| | - Jonathan M. Skelton
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, Britain
| | - Paul L. A. Popelier
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, Britain
| |
Collapse
|
6
|
Burn M, Popelier PLA. Gaussian Process Regression Models for Predicting Atomic Energies and Multipole Moments. J Chem Theory Comput 2023; 19:1370-1380. [PMID: 36757024 PMCID: PMC9979601 DOI: 10.1021/acs.jctc.2c00731] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
Developing a force field is a difficult task because its design is typically pulled in opposite directions by speed and accuracy. FFLUX breaks this trend by utilizing Gaussian process regression (GPR) to predict, at ab initio accuracy, atomic energies and multipole moments as obtained from the quantum theory of atoms in molecules (QTAIM). This work demonstrates that the in-house FFLUX training pipeline can generate successful GPR models for six representative molecules: peptide-capped glycine and alanine, glucose, paracetamol, aspirin, and ibuprofen. The molecules were sufficiently distorted to represent configurations from an AMBER-GAFF2 molecular dynamics run. All internal degrees of freedom were covered corresponding to 93 dimensions in the case of the largest molecule ibuprofen (33 atoms). Benefiting from active learning, the GPR models contain only about 2000 training points and return largely sub-kcal mol-1 prediction errors for the validation sets. A proof of concept has been reached for transferring the model produced through active learning on one atomic property to that of the remaining atomic properties. The prediction of electrostatic interaction can be assessed at the intermolecular level, and the vast majority of interactions have a root-mean-square error of less than 0.1 kJ mol-1 with a maximum value of ∼1 kJ mol-1 for a glycine and paracetamol dimer.
Collapse
|
7
|
Brown M, Skelton JM, Popelier PLA. Construction of a Gaussian Process Regression Model of Formamide for Use in Molecular Simulations. J Phys Chem A 2023; 127:1702-1714. [PMID: 36756842 PMCID: PMC9969515 DOI: 10.1021/acs.jpca.2c06566] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/10/2023]
Abstract
FFLUX, a novel force field based on quantum chemical topology, can perform molecular dynamics simulations with flexible multipole moments that change with geometry. This is enabled by Gaussian process regression machine learning models, which accurately predict atomic energies and multipole moments up to the hexadecapole. We have constructed a model of the formamide monomer at the B3LYP/aug-cc-pVTZ level of theory capable of sub-kJ mol-1 accuracy, with the maximum prediction error for the molecule being 0.8 kJ mol-1. This model was used in FFLUX simulations along with Lennard-Jones parameters to successfully optimize the geometry of formamide dimers with errors smaller than 0.1 Å compared to those obtained with D3-corrected B3LYP/aug-cc-pVTZ. Comparisons were also made to a force field constructed with static multipole moments and Lennard-Jones parameters. FFLUX recovers the expected energy ranking of dimers compared to the literature, and changes in C═O and C-N bond lengths associated with hydrogen bonding were found to be consistent with density functional theory.
Collapse
|
8
|
Symons BCB, Popelier PLA. Application of Quantum Chemical Topology Force Field FFLUX to Condensed Matter Simulations: Liquid Water. J Chem Theory Comput 2022; 18:5577-5588. [PMID: 35939826 PMCID: PMC9476653 DOI: 10.1021/acs.jctc.2c00311] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
We present here the first application of the quantum
chemical topology
force field FFLUX to condensed matter simulations. FFLUX offers many-body
potential energy surfaces learnt exclusively from ab initio data using Gaussian process regression. FFLUX also includes high-rank,
polarizable multipole moments (up to quadrupole moments in this work)
that are learnt from the same ab initio calculations
as the potential energy surfaces. Many-body effects (where a body
is an atom) and polarization are captured by the machine learning
models. The choice to use machine learning in this way allows the
force field’s representation of reality to be improved (e.g., by including higher order many-body effects) with
little detriment to the computational scaling of the code. In this
manner, FFLUX is inherently future-proof. The “plug and play”
nature of the machine learning models also ensures that FFLUX can
be applied to any system of interest, not just liquid water. In this
work we study liquid water across a range of temperatures and compare
the predicted bulk properties to experiment as well as other state-of-the-art
force fields AMOEBA(+CF), HIPPO, MB-Pol and SIBFA21. We find that
FFLUX finds a place amongst these.
Collapse
Affiliation(s)
- Benjamin C B Symons
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Paul L A Popelier
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|