Assessing the validity of DLPNO‐CCSD (T) in the calculation of activation and reaction energies of ubiquitous enzymatic reactions

2020 ◽  
Vol 41 (29) ◽  
pp. 2459-2468
Author(s):  
Pedro Paiva ◽  
Maria J. Ramos ◽  
Pedro A. Fernandes
Keyword(s):  
Enzymatic Reactions ◽  
Reaction Energies

scholarly journals Towards a barrier height benchmark set for biologically relevant systems

2016 ◽  
Author(s):  
Jimmy C Kromann ◽  
Anders S Christensen ◽  
Qiang Cui ◽  
Jan H. Jensen
Keyword(s):  
Structural Model ◽  
Geometry Optimization ◽  
System Size ◽  
Model Systems ◽  
Maximum Deviation ◽  
Enzymatic Reactions ◽  
Data Set ◽  
Biologically Relevant ◽  
Enzyme Model ◽  
Reaction Energies

We have collected computed barrier heights and reaction energies (and associated model structures) for five enzymes from studies published by Himo and co-workers. Using this data, obtained at the B3LYP/6- 311+G(2d,2p)[LANL2DZ]//B3LYP/6-31G(d,p) level of theory, we then benchmark PM6, PM7, PM7-TS, and DFTB3 and discuss the influence of system size, bulk solvation, and geometry re-optimization on the error. The mean absolute differences (MADs) observed for these five enzyme model systems are similar to those observed for PM6 and PM7 for smaller systems (10-15 kcal/mol), while DFTB results in a MAD that is significantly lower (6 kcal/mol). The MADs for PMx and DFTB3 are each dominated by large errors for a single system and if the system is disregarded the MADs fall to 4-5 kcal/mol. Overall, results for the condensed phase are neither more or less accurate relative to B3LYP than those in the gas phase. With the exception of PM7-TS, the MAD for small and large structural models are very similar, with a maximum deviation of 3 kcal/mol for PM6. Geometry optimization with PM6 shows that for one system this method predicts a different mechanism compared to B3LYP/6-31G(d,p). For the remaining systems geometry optimization of the large structural model increases the MAD relative to single points, by 2.5 and 1.8 kcal/mol for barriers and reaction energies. For the small structural model the corresponding MADs decrease by 0.4 and 1.2 kcal/mol, respectively. However, despite these small changes, significant changes in the structures are observed for some systems, such as proton transfer and hydrogen bonding rearrangements. The paper represents the first step in the process of creating a benchmark set of barriers computed for systems that are relatively large and representative of enzymatic reactions, a considerable challenge for any one research group but possible through a concerted effort by the community. We end by outlining steps needed to expand and improve the data set and how other researchers can contribute to the process.

scholarly journals Towards a barrier height benchmark set for biologically relevant systems

2016 ◽  
Author(s):  
Jimmy C Kromann ◽  
Anders S Christensen ◽  
Qiang Cui ◽  
Jan H. Jensen
Keyword(s):  
Structural Model ◽  
Geometry Optimization ◽  
System Size ◽  
Model Systems ◽  
Maximum Deviation ◽  
Enzymatic Reactions ◽  
Data Set ◽  
Biologically Relevant ◽  
Enzyme Model ◽  
Reaction Energies

We have collected computed barrier heights and reaction energies (and associated model structures) for five enzymes from studies published by Himo and co-workers. Using this data, obtained at the B3LYP/6- 311+G(2d,2p)[LANL2DZ]//B3LYP/6-31G(d,p) level of theory, we then benchmark PM6, PM7, PM7-TS, and DFTB3 and discuss the influence of system size, bulk solvation, and geometry re-optimization on the error. The mean absolute differences (MADs) observed for these five enzyme model systems are similar to those observed for PM6 and PM7 for smaller systems (10-15 kcal/mol), while DFTB results in a MAD that is significantly lower (6 kcal/mol). The MADs for PMx and DFTB3 are each dominated by large errors for a single system and if the system is disregarded the MADs fall to 4-5 kcal/mol. Overall, results for the condensed phase are neither more or less accurate relative to B3LYP than those in the gas phase. With the exception of PM7-TS, the MAD for small and large structural models are very similar, with a maximum deviation of 3 kcal/mol for PM6. Geometry optimization with PM6 shows that for one system this method predicts a different mechanism compared to B3LYP/6-31G(d,p). For the remaining systems geometry optimization of the large structural model increases the MAD relative to single points, by 2.5 and 1.8 kcal/mol for barriers and reaction energies. For the small structural model the corresponding MADs decrease by 0.4 and 1.2 kcal/mol, respectively. However, despite these small changes, significant changes in the structures are observed for some systems, such as proton transfer and hydrogen bonding rearrangements. The paper represents the first step in the process of creating a benchmark set of barriers computed for systems that are relatively large and representative of enzymatic reactions, a considerable challenge for any one research group but possible through a concerted effort by the community. We end by outlining steps needed to expand and improve the data set and how other researchers can contribute to the process.

scholarly journals Towards a barrier height benchmark set for biologically relevant systems

PeerJ ◽  
2016 ◽  
Vol 4 ◽  
pp. e1994 ◽  
Author(s):  
Jimmy C. Kromann ◽  
Anders S. Christensen ◽  
Qiang Cui ◽  
Jan H. Jensen
Keyword(s):  
Structural Model ◽  
Geometry Optimization ◽  
System Size ◽  
Model Systems ◽  
Maximum Deviation ◽  
Enzymatic Reactions ◽  
Data Set ◽  
Biologically Relevant ◽  
Enzyme Model ◽  
Reaction Energies

We have collected computed barrier heights and reaction energies (and associated model structures) for five enzymes from studies published by Himo and co-workers. Using this data, obtained at the B3LYP/6- 311+G(2d,2p)[LANL2DZ]//B3LYP/6-31G(d,p) level of theory, we then benchmark PM6, PM7, PM7-TS, and DFTB3 and discuss the influence of system size, bulk solvation, and geometry re-optimization on the error. The mean absolute differences (MADs) observed for these five enzyme model systems are similar to those observed for PM6 and PM7 for smaller systems (10–15 kcal/mol), while DFTB results in a MAD that is significantly lower (6 kcal/mol). The MADs for PMx and DFTB3 are each dominated by large errors for a single system and if the system is disregarded the MADs fall to 4–5 kcal/mol. Overall, results for the condensed phase are neither more or less accurate relative to B3LYP than those in the gas phase. With the exception of PM7-TS, the MAD for small and large structural models are very similar, with a maximum deviation of 3 kcal/mol for PM6. Geometry optimization with PM6 shows that for one system this method predicts a different mechanism compared to B3LYP/6-31G(d,p). For the remaining systems, geometry optimization of the large structural model increases the MAD relative to single points, by 2.5 and 1.8 kcal/mol for barriers and reaction energies. For the small structural model, the corresponding MADs decrease by 0.4 and 1.2 kcal/mol, respectively. However, despite these small changes, significant changes in the structures are observed for some systems, such as proton transfer and hydrogen bonding rearrangements. The paper represents the first step in the process of creating a benchmark set of barriers computed for systems that are relatively large and representative of enzymatic reactions, a considerable challenge for any one research group but possible through a concerted effort by the community. We end by outlining steps needed to expand and improve the data set and how other researchers can contribute to the process.

scholarly journals Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification

2021 ◽  
Author(s):  
Felix Brandt ◽  
Christoph Jacob
Keyword(s):  
Uncertainty Quantification ◽  
Computational Chemistry ◽  
Point Charge ◽  
Enzymatic Reactions ◽  
Variation Analysis ◽  
Computationally Efficient ◽  
Sources Of Uncertainty ◽  
Charge Variation ◽  
Reaction Energies ◽  
Selection Of

While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step towards quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region, and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol \textit{O}-methyltransferase, and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.

scholarly journals Towards a barrier height benchmark set for biologically relevant systems

2016 ◽  
Author(s):  
Jimmy C Kromann ◽  
Anders S Christensen ◽  
Jan H. Jensen
Keyword(s):  
Structural Model ◽  
Geometry Optimization ◽  
System Size ◽  
Model Systems ◽  
Maximum Deviation ◽  
Enzymatic Reactions ◽  
Data Set ◽  
Biologically Relevant ◽  
Enzyme Model ◽  
Reaction Energies

We have collected computed barrier heights and reaction energies (and associated model structures) for five enzymes from studies published by Himo and co-workers. Using this data, obtained at the B3LYP/6- 311+G(2d,2p)[LANL2DZ]//B3LYP/6-31G(d,p) level of theory, we then benchmark PM6, PM7, PM7-TS, and DFTB3 and discuss the influence of system size, bulk solvation, and geometry re-optimization on the error. The mean absolute differences (MADs) observed for these five enzyme model systems are similar to those observed for PM6 and PM7 for smaller systems (10-15 kcal/mol), while DFTB results in a MAD that is significantly lower (6 kcal/mol). The MADs for PMx and DFTB3 are each dominated by large errors for a single system and if the system is disregarded the MADs fall to 4-5 kcal/mol. Overall, results for the condensed phase are neither more or less accurate relative to B3LYP than those in the gas phase. With the exception of PM7-TS, the MAD for small and large structural models are very similar, with a maximum deviation of 3 kcal/mol for PM6. Geometry optimization with PM6 shows that for one system this method predicts a different mechanism compared to B3LYP/6-31G(d,p). For the remaining systems geometry optimization of the large structural model increases the MAD relative to single points, by 2.5 and 1.8 kcal/mol for barriers and reaction energies. For the small structural model the corresponding MADs decrease by 0.4 and 1.2 kcal/mol, respectively. However, despite these small changes, significant changes in the structures are observed for some systems, such as proton transfer and hydrogen bonding rearrangements. The paper represents the first step in the process of creating a benchmark set of barriers computed for systems that are relatively large and representative of enzymatic reactions, a considerable challenge for any one research group but possible through a concerted effort by the community. We end by outlining steps needed to expand and improve the data set and how other researchers can contribute to the process.

scholarly journals Systematic QM Region Construction in QM/MM Calculations Based on Uncertainty Quantification

2022 ◽  
Author(s):  
Felix Brandt ◽  
Christoph Jacob
Keyword(s):  
Uncertainty Quantification ◽  
Computational Chemistry ◽  
Point Charge ◽  
Enzymatic Reactions ◽  
Variation Analysis ◽  
Computationally Efficient ◽  
Sources Of Uncertainty ◽  
Charge Variation ◽  
Reaction Energies ◽  
Selection Of

While QM/MM studies of enzymatic reactions are widely used in computational chemistry, the results of such studies are subject to numerous sources of uncertainty, and the effect of different choices by the simulation scientist that are required when setting up QM/MM calculations is often unclear. In particular, the selection of the QM region is crucial for obtaining accurate and reliable results. Simply including amino acids by their distance to the active site is mostly not sufficient as necessary residues are missing or unimportant residues are included without evidence. Here, we take a first step towards quantifying uncertainties in QM/MM calculations by assessing the sensitivity of QM/MM reaction energies with respect to variations of the MM point charges. We show that such a point charge variation analysis (PCVA) can be employed to judge the accuracy of QM/MM reaction energies obtained with a selected QM region, and devise a protocol to systematically construct QM regions that minimize this uncertainty. We apply such a PCVA to the example of catechol \textit{O}-methyltransferase, and demonstrate that it provides a simple and reliable approach for the construction of the QM region. Our PCVA-based scheme is computationally efficient and requires only calculations for a system with a minimal QM region. Our work highlights the promise of applying methods of uncertainty quantification in computational chemistry.

Development of 96 Multiple Injection-GC-MS Technique and Its Application in Protein Engineering of Natural and Non-Natural Enzymatic Reactions

2019 ◽  
Author(s):  
Anja Knorrscheidt ◽  
Pascal Püllmann ◽  
Eugen Schell ◽  
Dominik Homann ◽  
Erik Freier ◽  
...  
Keyword(s):  
Microtiter Plate ◽  
Axial Ligand ◽  
Cell System ◽  
Enzymatic Reactions ◽  
The Novel ◽  
Internal Standards ◽  
E Coli ◽  
Enzyme Screening ◽  
Plate Analysis ◽  
Microtiter Plate Analysis

Directed evolution requires the screening of enzyme libraries in biological matrices. Available assays are mostly substrate or enzyme specific. Chromatographic techniques like LC and GC overcome this limitation, but require long analysis times. The herein developed multiple injections in a single experimental run (MISER) using GC coupled to MS allows the injection of samples every 33 s resulting in 96-well microtiter plate analysis within 50 min. This technique is implementable in any GC-MS system with autosampling. Since the GC-MS is far less prone to ion suppression than LCMS, no chromatographic separation is required. This allows the utilisation of an internal standards and the detection of main and side-product. To prove the feasibility of the system in enzyme screening, two libraries were assessed: i) YfeX library in an E. coli whole cell system for the carbene-transfer reaction on indole revealing the novel axial ligand tryptophan, ii) a library of 616 chimeras of fungal unspecific peroxygenase (UPO) in S. cerevisiae supernatant for hydroxylation of tetralin resulting in novel constructs. The data quality and representation are automatically assessed by a new R-script.

Modular Simulation to Determine the Optimal Operating Policy of a Batch Reactor for the Enzymatic Fructose Reduction to Mannitol with the in situ Continuous Enzymatic Regeneration of the NAD Cofactor

2017 ◽  
Vol 68 (9) ◽  
pp. 2196-2203 ◽  
Author(s):  
Mara Crisan ◽  
Gheorghe Maria
Keyword(s):  
Batch Reactor ◽  
Production Capacity ◽  
Optimization Procedure ◽  
Enzymatic Reactions ◽  
Mannitol Dehydrogenase ◽  
Optimal Operating ◽  
Operating Policy ◽  
Enzymatic Reduction ◽  
Optimal Operating Policy

Novel coupled enzymatic systems reported important applications in the industrial bio-catalysis. Multi-enzymatic reactions can successfully replace complex chemical syntheses, using milder reaction conditions, and generating less waste. For such systems acting simultaneously, the model-based engineering calculations (design, reactor operation optimization) are difficult tasks, because they must account for interacting reactions, differences in enzymes optimal activity domains and deactivation kinetics. The determination of the optimal operating mode (enzyme ratios, enzyme feeding policy, temperature, pH) often turns into a difficult multi-objective optimization problem with multiple constraints to be solved for every particular system. The paper focuses on applying a modular screening procedure that can identify the optimal operating policy of an enzymatic reactor, which minimizes the enzyme consumption, given the process kinetic model, and an imposed production capacity. Following an optimization procedure, the process effectiveness is evaluated in a systematic approach, by including simple batch reactor (BR), batch with intermittent addition of the key-enzyme following certain optimal policies (BRP). Exemplification is made for the case of the enzymatic reduction of D-fructose to mannitol by using suspended MDH (mannitol dehydrogenase) and NADH (Nicotinamide adenine dinucleotide) cofactor, with the in-situ continuous regeneration of the cofactor by the expense of formate degradation in the presence of suspended FDH (Formate dehydrogenase).

Synthesis of Rare Pentoses Using Microbial and Enzymatic Reactions

2016 ◽  
Vol 20 (14) ◽  
pp. 1456-1464 ◽  
Author(s):  
Zijie Li ◽  
Xiao-Dong Gao ◽  
Li Cai
Keyword(s):  
Enzymatic Reactions
Sign in / Sign up

Export Citation Format

Share Document