1
|
Jestilä JS, Wu N, Priante F, Foster AS. Accelerated Lignocellulosic Molecule Adsorption Structure Determination. J Chem Theory Comput 2024; 20:2297-2312. [PMID: 38408381 PMCID: PMC10939001 DOI: 10.1021/acs.jctc.3c01292] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 02/07/2024] [Accepted: 02/08/2024] [Indexed: 02/28/2024]
Abstract
Here, we present a study combining Bayesian optimization structural inference with the machine learning interatomic potential Neural Equivariant Interatomic Potential (NequIP) to accelerate and enable the study of the adsorption of the conformationally flexible lignocellulosic molecules β-d-xylose and 1,4-β-d-xylotetraose on a copper surface. The number of structure evaluations needed to map out the relevant potential energy surfaces are reduced by Bayesian optimization, while NequIP minimizes the time spent on each evaluation, ultimately resulting in cost-efficient and reliable sampling of large systems and configurational spaces. Although the applicability of Bayesian optimization for the conformational analysis of the more flexible xylotetraose molecule is restricted by the sample complexity bottleneck, the latter can be effectively bypassed with external conformer search tools, such as the Conformer-Rotamer Ensemble Sampling Tool, facilitating the subsequent lower-dimensional global minimum adsorption structure determination. Finally, we demonstrate the applicability of the described approach to find adsorption structures practically equivalent to the density functional theory counterparts at a fraction of the computational cost.
Collapse
Affiliation(s)
- Joakim S. Jestilä
- Department
of Applied Physics, Aalto University, 00076 Aalto, Espoo, Finland
| | - Nian Wu
- Department
of Applied Physics, Aalto University, 00076 Aalto, Espoo, Finland
| | - Fabio Priante
- Department
of Applied Physics, Aalto University, 00076 Aalto, Espoo, Finland
| | - Adam S. Foster
- Department
of Applied Physics, Aalto University, 00076 Aalto, Espoo, Finland
- Nano
Life Science Institute (WPI-NanoLSI), Kanazawa
University, 920-1192 Kanazawa, Japan
| |
Collapse
|
2
|
Ma L, Li G, Liu Y, Li Z, Miao Y, Wan Q, Liu D, Zhang R. Investigating the effect of substrate binding on the catalytic activity of xylanase. Appl Microbiol Biotechnol 2023; 107:6873-6886. [PMID: 37715802 DOI: 10.1007/s00253-023-12774-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2023] [Revised: 08/28/2023] [Accepted: 09/03/2023] [Indexed: 09/18/2023]
Abstract
XynAF1 from Aspergillus fumigatus Z5 is an efficient thermophilic xylanase belonging to glycoside hydrolase family 10 (GH10). The non-catalytic amino acids N179 and R246 in its catalytic center formed one and three intermolecular H-bonds with the substrate in the aglycone region, respectively. Here we purified XynAF1-N179S and XynAF1-R246K, and obtained the protein-product complex structures by X-ray diffraction. The snapshots indicated that mutations at N179 and R246 had decreased the substrate-binding ability in the aglycone region. XynAF1-N179S, XynAF1-R246K, and XynAF1-N179S-R246K lost one, three, and four H-bonds with the substrate in comparison with the wild-type XynAF1, respectively, but this had little influence on the protein structure. As expected, N179S, R246K, and N179S-R246K led to a gradual decrease of substrate affinity of XynAF1. Interestingly, the enzyme assay showed that N179S increased catalytic efficiency, while both R246K and N179S-R246K had decreased catalytic efficiency. KEY POINTS: • The non-catalytic amino acids of XynAF1 could form H-bonds with the substrate. • The protein-product complex structures were obtained by X-ray diffraction. • The enzyme-substrate-binding capacity could affect enzyme catalytic efficiency.
Collapse
Affiliation(s)
- Lei Ma
- College of Life Sciences and Engineering, Henan University of Urban Construction, Pingdingshan, 467000, Henan, People's Republic of China
| | - Guangqi Li
- Key Laboratory of Microbial Resources Collection and Preservation, Ministry of Agriculture, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, 100081, People's Republic of China
| | - Yunpeng Liu
- Key Laboratory of Microbial Resources Collection and Preservation, Ministry of Agriculture, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, 100081, People's Republic of China
| | - Zhihong Li
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, 210095, People's Republic of China
| | - Youzhi Miao
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, 210095, People's Republic of China
| | - Qun Wan
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, 210095, People's Republic of China
| | - Dongyang Liu
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, 210095, People's Republic of China
| | - Ruifu Zhang
- Key Laboratory of Microbial Resources Collection and Preservation, Ministry of Agriculture, Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing, 100081, People's Republic of China.
- Jiangsu Provincial Key Lab for Organic Solid Waste Utilization, National Engineering Research Center for Organic-based Fertilizers, Jiangsu Collaborative Innovation Center for Solid Organic Waste Resource Utilization, Nanjing Agricultural University, Nanjing, 210095, People's Republic of China.
| |
Collapse
|
3
|
Ardèvol A, Rovira C. Reaction Mechanisms in Carbohydrate-Active Enzymes: Glycoside Hydrolases and Glycosyltransferases. Insights from ab Initio Quantum Mechanics/Molecular Mechanics Dynamic Simulations. J Am Chem Soc 2015; 137:7528-47. [DOI: 10.1021/jacs.5b01156] [Citation(s) in RCA: 151] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Albert Ardèvol
- Departament
de Química Orgànica and Institut de Química Teòrica
i Computacional (IQTCUB), Universitat de Barcelona, Martí
i Franquès 1, 08028 Barcelona, Spain
| | - Carme Rovira
- Departament
de Química Orgànica and Institut de Química Teòrica
i Computacional (IQTCUB), Universitat de Barcelona, Martí
i Franquès 1, 08028 Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Passeig Lluís Companys, 23, 08018 Barcelona, Spain
| |
Collapse
|
4
|
Iglesias-Fernández J, Raich L, Ardèvol A, Rovira C. The complete conformational free energy landscape of β-xylose reveals a two-fold catalytic itinerary for β-xylanases. Chem Sci 2015; 6:1167-1177. [PMID: 29560204 PMCID: PMC5811086 DOI: 10.1039/c4sc02240h] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2014] [Accepted: 10/27/2014] [Indexed: 01/28/2023] Open
Abstract
Unraveling the conformational catalytic itinerary of glycoside hydrolases (GHs) is a growing topic of interest in glycobiology, with major impact in the design of GH inhibitors. β-xylanases are responsible for the hydrolysis of glycosidic bonds in β-xylans, a group of hemicelluloses of high biotechnological interest that are found in plant cell walls. The precise conformations followed by the substrate during catalysis in β-xylanases have not been unambiguously resolved, with three different pathways being proposed from structural analyses. In this work, we compute the conformational free energy landscape (FEL) of β-xylose to predict the most likely catalytic itineraries followed by β-xylanases. The calculations are performed by means of ab initio metadynamics, using the Cremer-Pople puckering coordinates as collective variables. The computed FEL supports only two of the previously proposed itineraries, 2SO → [2,5B]ǂ → 5S1 and 1S3 → [4H3]ǂ → 4C1, which clearly appear in low energy regions of the FEL. Consistently, 2SO and 1S3 are conformations preactivated for catalysis in terms of free energy/anomeric charge and bond distances. The results however exclude the OE → [OS2]ǂ → B2,5 itinerary that has been recently proposed for a family 11 xylanase. Classical and ab initio QM/MM molecular dynamics simulations reveal that, in this case, the observed OE conformation has been enforced by enzyme mutation. These results add a word of caution on using modified enzymes to inform on catalytic conformational itineraries of glycoside hydrolases.
Collapse
Affiliation(s)
- Javier Iglesias-Fernández
- Departament de Química Orgànica and Institut de Química Teòrica i Computacional (IQTCUB) , Universitat de Barcelona , Martí i Franquès 1 , 08028 Barcelona , Spain .
| | - Lluís Raich
- Departament de Química Orgànica and Institut de Química Teòrica i Computacional (IQTCUB) , Universitat de Barcelona , Martí i Franquès 1 , 08028 Barcelona , Spain .
| | - Albert Ardèvol
- Department of Chemistry and Applied Biosciences , ETH Zürich , USI Campus , 6900 Lugano , Switzerland
| | - Carme Rovira
- Departament de Química Orgànica and Institut de Química Teòrica i Computacional (IQTCUB) , Universitat de Barcelona , Martí i Franquès 1 , 08028 Barcelona , Spain .
- Institució Catalana de Recerca i Estudis Avançats (ICREA) , Passeig Lluís Companys , 23 , 08018 Barcelona , Spain
| |
Collapse
|
5
|
QM/MM study of catalytic mechanism of Xylanase Cex from Cellulomonas fimi. J Mol Graph Model 2012; 37:67-76. [DOI: 10.1016/j.jmgm.2012.04.005] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2012] [Revised: 04/01/2012] [Accepted: 04/17/2012] [Indexed: 12/13/2022]
|
6
|
Amorim L, Marcelo F, Rousseau C, Nieto L, Jiménez‐Barbero J, Marrot J, Rauter AP, Sollogoub M, Bols M, Blériot Y. Direct Experimental Evidence for the High Chemical Reactivity of α‐ and β‐Xylopyranosides Adopting a
2,5
B
Conformation in Glycosyl Transfer. Chemistry 2011; 17:7345-56. [DOI: 10.1002/chem.201003251] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2010] [Indexed: 12/12/2022]
Affiliation(s)
- Luis Amorim
- UPMC Univ Paris 06, Institut Parisien de Chimie Moléculaire (UMR 7201), FR 2769, C181, 4 place Jussieu, 75005 Paris (France)
| | - Filipa Marcelo
- UPMC Univ Paris 06, Institut Parisien de Chimie Moléculaire (UMR 7201), FR 2769, C181, 4 place Jussieu, 75005 Paris (France)
- Carbohydrate Chemistry Group, CQB‐FCUL, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa, 1749‐016 Lisbon (Portugal)
| | - Cyril Rousseau
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Kbh Ø (Denmark)
- Present address: Université d'Artois, IUT de Béthune, UCCS Artois, UMR 8181, 1230 rue de l'Université, BP 819, 62408 Béthune cedex (France)
| | - Lidia Nieto
- Centro de Investigaciones Biológicas, CSIC, 28040 Madrid (Spain)
| | | | - Jérôme Marrot
- Institut Lavoisier, Université de Versailles‐Saint‐Quentin, UMR 8180, 78035 Versailles (France)
| | - Amélia P. Rauter
- Carbohydrate Chemistry Group, CQB‐FCUL, Departamento de Química e Bioquímica, Faculdade de Ciências da Universidade de Lisboa, 1749‐016 Lisbon (Portugal)
| | - Matthieu Sollogoub
- UPMC Univ Paris 06, Institut Parisien de Chimie Moléculaire (UMR 7201), FR 2769, C181, 4 place Jussieu, 75005 Paris (France)
| | - Mikael Bols
- Department of Chemistry, University of Copenhagen, Universitetsparken 5, 2100 Kbh Ø (Denmark)
| | - Yves Blériot
- UPMC Univ Paris 06, Institut Parisien de Chimie Moléculaire (UMR 7201), FR 2769, C181, 4 place Jussieu, 75005 Paris (France)
- Present address: Université de Poitiers, UMR 6514, Laboratoire “Synthèse et Réactivité des Substances Naturelles”, 40 avenue du Recteur Pineau, 86022 Poitiers Cedex (France)
| |
Collapse
|
7
|
Soliman MES, Ruggiero GD, Pernía JJR, Greig IR, Williams IH. Computational mutagenesis reveals the role of active-site tyrosine in stabilising a boat conformation for the substrate: QM/MM molecular dynamics studies of wild-type and mutant xylanases. Org Biomol Chem 2009; 7:460-8. [DOI: 10.1039/b814695k] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
8
|
Ionescu AR, Whitfield DM, Zgierski MZ, Nukada T. Investigations into the role of oxacarbenium ions in glycosylation reactions by ab initio molecular dynamics. Carbohydr Res 2006; 341:2912-20. [PMID: 17069777 DOI: 10.1016/j.carres.2006.09.027] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2006] [Revised: 09/14/2006] [Accepted: 09/21/2006] [Indexed: 10/24/2022]
Abstract
We present a constrained ab initio molecular dynamics method that allows the modeling of the conformational interconversions of glycopyranosyl oxacarbenium ions. The model was successfully tested by estimating the barriers to ring inversion for two 4-substituted tetrahydropyranosyl oxacarbenium ions. The model was further extended to predict the pathways that connect the (4)H(3) half-chair conformation of 2,3,4,6-tetra-O-methyl-d-glucopyranosyl cation to its inverted (5)S(1) conformation and the (4)H(3) half-chair conformation of 2,3,4,6-tetra-O-methyl-d-mannopyranosyl cation to its inverted (3)E conformation. The modeled interconversion pathways reconcile a large body of experimental work on the acid-catalyzed hydrolysis of glycosides and the mechanisms of a number of glucosidases and mannosidases.
Collapse
Affiliation(s)
- Andrei R Ionescu
- Steacie Institute for Molecular Sciences, NRC Canada, 100 Sussex Drive, Ottawa, Ont., Canada
| | | | | | | |
Collapse
|
9
|
Biarnés X, Nieto J, Planas A, Rovira C. Substrate Distortion in the Michaelis Complex of Bacillus 1,3–1,4-β-Glucanase. J Biol Chem 2006; 281:1432-41. [PMID: 16260784 DOI: 10.1074/jbc.m507643200] [Citation(s) in RCA: 77] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
The structure and dynamics of the enzyme-substrate complex of Bacillus 1,3-1,4-beta-glucanase, one of the most active glycoside hydrolases, is investigated by means of Car-Parrinello molecular dynamics simulations (CPMD) combined with force field molecular dynamics (QM/MM CPMD). It is found that the substrate sugar ring located at the -1 subsite adopts a distorted 1S3 skew-boat conformation upon binding to the enzyme. With respect to the undistorted 4C1 chair conformation, the 1S3 skew-boat conformation is characterized by: (a) an increase of charge at the anomeric carbon (C1), (b) an increase of the distance between C1 and the leaving group, and (c) a decrease of the intraring O5-C1 distance. Therefore, our results clearly show that the distorted conformation resembles both structurally and electronically the transition state of the reaction in which the substrate acquires oxocarbenium ion character, and the glycosidic bond is partially broken. Together with analysis of the substrate conformational dynamics, it is concluded that the main determinants of substrate distortion have a structural origin. To fit into the binding pocket, it is necessary that the aglycon leaving group is oriented toward the beta region, and the skew-boat conformation naturally fulfills this premise. Only when the aglycon is removed from the calculation the substrate recovers the all-chair conformation, in agreement with the recent determination of the enzyme product structure. The QM/MM protocol developed here is able to predict the conformational distortion of substrate binding in glycoside hydrolases because it accounts for polarization and charge reorganization at the -1 sugar ring. It thus provides a powerful tool to model E.S complexes for which experimental information is not yet available.
Collapse
Affiliation(s)
- Xevi Biarnés
- Centre especial de Recerca en Química Teòrica, Parc Científic de Barcelona, Josep Samitier 1-5, 08028 Barcelona, Spain
| | | | | | | |
Collapse
|