1
|
Albayati SH, Nezhad NG, Taki AG, Rahman RNZRA. Efficient and easible biocatalysts: Strategies for enzyme improvement. A review. Int J Biol Macromol 2024; 276:133978. [PMID: 39038570 DOI: 10.1016/j.ijbiomac.2024.133978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Revised: 06/19/2024] [Accepted: 07/16/2024] [Indexed: 07/24/2024]
Abstract
Owing to the environmental friendliness and vast advantages that enzymes offer in the biotechnology and industry fields, biocatalysts are a prolific investigation field. However, the low catalytic activity, stability, and specific selectivity of the enzyme limit the range of the reaction enzymes involved in. A comprehensive understanding of the protein structure and dynamics in terms of molecular details enables us to tackle these limitations effectively and enhance the catalytic activity by enzyme engineering or modifying the supports and solvents. Along with different strategies including computational, enzyme engineering based on DNA recombination, enzyme immobilization, additives, chemical modification, and physicochemical modification approaches can be promising for the wide spread of industrial enzyme usage. This is attributed to the successful application of biocatalysts in industrial and synthetic processes requires a system that exhibits stability, activity, and reusability in a continuous flow process, thereby reducing the production cost. The main goal of this review is to display relevant approaches for improving enzyme characteristics to overcome their industrial application.
Collapse
Affiliation(s)
- Samah Hashim Albayati
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia; Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
| | - Nima Ghahremani Nezhad
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia; Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia
| | - Anmar Ghanim Taki
- Department of Radiology Techniques, Health and Medical Techniques College, Alnoor University, Mosul, Iraq
| | - Raja Noor Zaliha Raja Abd Rahman
- Enzyme and Microbial Technology Research Centre, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia; Department of Microbiology, Faculty of Biotechnology and Biomolecular Sciences, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia; Institute Bioscience, Universiti Putra Malaysia, Serdang 43400, Selangor, Malaysia.
| |
Collapse
|
2
|
Eschenbach P, Neugebauer J. Subsystem density-functional theory: A reliable tool for spin-density based properties. J Chem Phys 2022; 157:130902. [PMID: 36209003 DOI: 10.1063/5.0103091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Subsystem density-functional theory compiles a set of features that allow for efficiently calculating properties of very large open-shell radical systems such as organic radical crystals, proteins, or deoxyribonucleic acid stacks. It is computationally less costly than correlated ab initio wave function approaches and can pragmatically avoid the overdelocalization problem of Kohn-Sham density-functional theory without employing hard constraints on the electron-density. Additionally, subsystem density-functional theory calculations commonly start from isolated fragment electron densities, pragmatically preserving a priori specified subsystem spin-patterns throughout the calculation. Methods based on subsystem density-functional theory have seen a rapid development over the past years and have become important tools for describing open-shell properties. In this Perspective, we address open questions and possible developments toward challenging future applications in connection with subsystem density-functional theory for spin-dependent properties.
Collapse
Affiliation(s)
- Patrick Eschenbach
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Simulation, Westfälische Wilhelms-Universität Münster, Corrensstraße 36, 48149 Münster, Germany
| |
Collapse
|
3
|
Zhang Y, Zhu X, Zhang H, Yan J, Xu P, Wu P, Wu S, Bai C. Mechanism Study of Proteins under Membrane Environment. MEMBRANES 2022; 12:membranes12070694. [PMID: 35877897 PMCID: PMC9322369 DOI: 10.3390/membranes12070694] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/12/2022] [Revised: 07/01/2022] [Accepted: 07/05/2022] [Indexed: 11/24/2022]
Abstract
Membrane proteins play crucial roles in various physiological processes, including molecule transport across membranes, cell communication, and signal transduction. Approximately 60% of known drug targets are membrane proteins. There is a significant need to deeply understand the working mechanism of membrane proteins in detail, which is a challenging work due to the lack of available membrane structures and their large spatial scale. Membrane proteins carry out vital physiological functions through conformational changes. In the current study, we utilized a coarse-grained (CG) model to investigate three representative membrane protein systems: the TMEM16A channel, the family C GPCRs mGlu2 receptor, and the P4-ATPase phospholipid transporter. We constructed the reaction pathway of conformational changes between the two-end structures. Energy profiles and energy barriers were calculated. These data could provide reasonable explanations for TMEM16A activation, the mGlu2 receptor activation process, and P4-ATPase phospholipid transport. Although they all belong to the members of membrane proteins, they behave differently in terms of energy. Our work investigated the working mechanism of membrane proteins and could give novel insights into other membrane protein systems of interest.
Collapse
Affiliation(s)
- Yue Zhang
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, China
- Warshel Institute for Computational Biology, Shenzhen 518172, China
| | - Xiaohong Zhu
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, China
- Warshel Institute for Computational Biology, Shenzhen 518172, China
| | - Honghui Zhang
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- Warshel Institute for Computational Biology, Shenzhen 518172, China
| | - Junfang Yan
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- Warshel Institute for Computational Biology, Shenzhen 518172, China
| | - Peiyi Xu
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- Warshel Institute for Computational Biology, Shenzhen 518172, China
| | - Peng Wu
- School of Biomedical Engineering, Health Science Center, Shenzhen University, Shenzhen 518055, China;
| | - Song Wu
- South China Hospital, Health Science Center, Shenzhen University, Shenzhen 518116, China
- Correspondence: (S.W.); (C.B.)
| | - Chen Bai
- School of Life and Health Sciences, School of Medicine, The Chinese University of Hong Kong, Shenzhen 518172, China; (Y.Z.); (X.Z.); (H.Z.); (J.Y.); (P.X.)
- Warshel Institute for Computational Biology, Shenzhen 518172, China
- Chenzhu Biotechnology Co., Ltd., Hangzhou 310005, China
- Correspondence: (S.W.); (C.B.)
| |
Collapse
|
4
|
Shi D, An K, Zhang H, Xu P, Bai C. Application of Coarse-Grained (CG) Models to Explore Conformational Pathway of Large-Scale Protein Machines. ENTROPY 2022; 24:e24050620. [PMID: 35626506 PMCID: PMC9140642 DOI: 10.3390/e24050620] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Revised: 04/16/2022] [Accepted: 04/27/2022] [Indexed: 12/14/2022]
Abstract
Protein machines are clusters of protein assemblies that function in order to control the transfer of matter and energy in cells. For a specific protein machine, its working mechanisms are not only determined by the static crystal structures, but also related to the conformational transition dynamics and the corresponding energy profiles. With the rapid development of crystallographic techniques, the spatial scale of resolved structures is reaching up to thousands of residues, and the concomitant conformational changes become more and more complicated, posing a great challenge for computational biology research. Previously, a coarse-grained (CG) model aiming at conformational free energy evaluation was developed and showed excellent ability to reproduce the energy profiles by accurate electrostatic interaction calculations. In this study, we extended the application of the CG model to a series of large-scale protein machine systems. The spike protein trimer of SARS-CoV-2, ATP citrate lyase (ACLY) tetramer, and P4-ATPases systems were carefully studied and discussed as examples. It is indicated that the CG model is effective to depict the energy profiles of the conformational pathway between two endpoint structures, especially for large-scale systems. Both the energy change and energy barrier between endpoint structures provide reasonable mechanism explanations for the associated biological processes, including the opening of receptor binding domain (RBD) of spike protein, the phospholipid transportation of P4-ATPase, and the loop translocation of ACLY. Taken together, the CG model provides a suitable alternative in mechanistic studies related to conformational change in large-scale protein machines.
Collapse
Affiliation(s)
- Danfeng Shi
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen 518172, China; (D.S.); (K.A.); (H.Z.); (P.X.)
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, China
| | - Ke An
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen 518172, China; (D.S.); (K.A.); (H.Z.); (P.X.)
- School of Chemistry and Materials Science, University of Science and Technology of China, Hefei 230026, China
| | - Honghui Zhang
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen 518172, China; (D.S.); (K.A.); (H.Z.); (P.X.)
| | - Peiyi Xu
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen 518172, China; (D.S.); (K.A.); (H.Z.); (P.X.)
| | - Chen Bai
- Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong (Shenzhen), Shenzhen 518172, China; (D.S.); (K.A.); (H.Z.); (P.X.)
- Correspondence:
| |
Collapse
|
5
|
Planas-Iglesias J, Marques SM, Pinto GP, Musil M, Stourac J, Damborsky J, Bednar D. Computational design of enzymes for biotechnological applications. Biotechnol Adv 2021; 47:107696. [PMID: 33513434 DOI: 10.1016/j.biotechadv.2021.107696] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 01/12/2021] [Accepted: 01/13/2021] [Indexed: 12/14/2022]
Abstract
Enzymes are the natural catalysts that execute biochemical reactions upholding life. Their natural effectiveness has been fine-tuned as a result of millions of years of natural evolution. Such catalytic effectiveness has prompted the use of biocatalysts from multiple sources on different applications, including the industrial production of goods (food and beverages, detergents, textile, and pharmaceutics), environmental protection, and biomedical applications. Natural enzymes often need to be improved by protein engineering to optimize their function in non-native environments. Recent technological advances have greatly facilitated this process by providing the experimental approaches of directed evolution or by enabling computer-assisted applications. Directed evolution mimics the natural selection process in a highly accelerated fashion at the expense of arduous laboratory work and economic resources. Theoretical methods provide predictions and represent an attractive complement to such experiments by waiving their inherent costs. Computational techniques can be used to engineer enzymatic reactivity, substrate specificity and ligand binding, access pathways and ligand transport, and global properties like protein stability, solubility, and flexibility. Theoretical approaches can also identify hotspots on the protein sequence for mutagenesis and predict suitable alternatives for selected positions with expected outcomes. This review covers the latest advances in computational methods for enzyme engineering and presents many successful case studies.
Collapse
Affiliation(s)
- Joan Planas-Iglesias
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Sérgio M Marques
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Gaspar P Pinto
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Milos Musil
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic; IT4Innovations Centre of Excellence, Faculty of Information Technology, Brno University of Technology, 61266 Brno, Czech Republic
| | - Jan Stourac
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic.
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5/A13, 625 00 Brno, Czech Republic.
| |
Collapse
|
6
|
Multistate Reactive Molecular Dynamics Simulations of Proton Diffusion in Water Clusters and in the Bulk. J Phys Chem B 2019; 123:9846-9861. [PMID: 31647873 DOI: 10.1021/acs.jpcb.9b03258] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
The molecular mechanics with proton transfer (MMPT) force field is combined with multistate adiabatic reactive molecular dynamics (MS-ARMD) to describe proton transport in the condensed phase. Parametrization for small protonated water clusters based on electronic structure calculations at the MP2/6-311+G(2d,2p) level of theory and refinement by comparing with infrared spectra for a protonated water tetramer yields a force field which faithfully describes the minimum energy structures of small protonated water clusters. In protonated water clusters up to (H2O)100H+, the proton hopping rate is around 100 hops/ns. This rate converges for 21 ≤ n ≤ 31, and no further speedup in bulk water is found. This indicates that bulklike behavior requires the solvation of a Zundel motif by ∼25 water molecules, which corresponds to the second solvation sphere. For smaller cluster sizes, the number of available states (i.e., the number of proton acceptors) is too small and slows down proton-transfer rates. The cluster simulations confirm that the excess proton is typically located on the surface. The free-energy surface as a function of the weights of the two lowest states and a configurational parameter suggests that the "special pair" plays a central role in rapid proton transport. The barriers between this minimum-energy structure and the Zundel and Eigen minima are sufficiently low (∼1 kcal/mol, consistent with recent experiments and commensurate with a hopping rate of ∼100/ns or 1 every 10 ps), leading to a highly dynamic environment. These findings are also consistent with recent experiments which find that Zundel-type hydration geometries are prevalent in bulk water.
Collapse
|
7
|
Bruno A, Costantino G, Sartori L, Radi M. The In Silico Drug Discovery Toolbox: Applications in Lead Discovery and Optimization. Curr Med Chem 2019; 26:3838-3873. [PMID: 29110597 DOI: 10.2174/0929867324666171107101035] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2017] [Revised: 09/27/2017] [Accepted: 09/28/2017] [Indexed: 01/04/2023]
Abstract
BACKGROUND Discovery and development of a new drug is a long lasting and expensive journey that takes around 20 years from starting idea to approval and marketing of new medication. Despite R&D expenditures have been constantly increasing in the last few years, the number of new drugs introduced into market has been steadily declining. This is mainly due to preclinical and clinical safety issues, which still represent about 40% of drug discontinuation. To cope with this issue, a number of in silico techniques are currently being used for an early stage evaluation/prediction of potential safety issues, allowing to increase the drug-discovery success rate and reduce costs associated with the development of a new drug. METHODS In the present review, we will analyse the early steps of the drug-discovery pipeline, describing the sequence of steps from disease selection to lead optimization and focusing on the most common in silico tools used to assess attrition risks and build a mitigation plan. RESULTS A comprehensive list of widely used in silico tools, databases, and public initiatives that can be effectively implemented and used in the drug discovery pipeline has been provided. A few examples of how these tools can be problem-solving and how they may increase the success rate of a drug discovery and development program have been also provided. Finally, selected examples where the application of in silico tools had effectively contributed to the development of marketed drugs or clinical candidates will be given. CONCLUSION The in silico toolbox finds great application in every step of early drug discovery: (i) target identification and validation; (ii) hit identification; (iii) hit-to-lead; and (iv) lead optimization. Each of these steps has been described in details, providing a useful overview on the role played by in silico tools in the decision-making process to speed-up the discovery of new drugs.
Collapse
Affiliation(s)
- Agostino Bruno
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Gabriele Costantino
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| | - Luca Sartori
- Experimental Therapeutics Unit, IFOM - The FIRC Institute for Molecular Oncology Foundation, Via Adamello 16 - 20139 Milano, Italy
| | - Marco Radi
- Dipartimento di Scienze degli Alimenti e del Farmaco, Universita degli Studi di Parma, Viale delle Scienze, 27/A, 43124 Parma, Italy
| |
Collapse
|
8
|
Kokkonen P, Bednar D, Pinto G, Prokop Z, Damborsky J. Engineering enzyme access tunnels. Biotechnol Adv 2019; 37:107386. [PMID: 31026496 DOI: 10.1016/j.biotechadv.2019.04.008] [Citation(s) in RCA: 107] [Impact Index Per Article: 21.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2018] [Revised: 04/16/2019] [Accepted: 04/18/2019] [Indexed: 12/14/2022]
Abstract
Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.
Collapse
Affiliation(s)
- Piia Kokkonen
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic
| | - David Bednar
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Gaspar Pinto
- International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Zbynek Prokop
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic
| | - Jiri Damborsky
- Loschmidt Laboratories, Department of Experimental Biology and RECETOX, Faculty of Science, Masaryk University, Kamenice 5, 625 00 Brno, Czech Republic; International Clinical Research Center, St. Anne's University Hospital Brno, Pekarska 53, 656 91 Brno, Czech Republic.
| |
Collapse
|
9
|
Yin F, Tang C, Tang Y, Gui Y, Zhao Z. Reactive Molecular Dynamics Study of Effects of Small-Molecule Organic Acids on PMIA Thermal Decomposition. J Phys Chem B 2018; 122:10384-10392. [PMID: 30362742 DOI: 10.1021/acs.jpcb.8b09343] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Reactive molecular dynamics was used to investigate the atomic-level mechanism of formic acid-accelerated deterioration of meta-aramid (PMIA) fibers. The simulation results showed that formic acid promoted PMIA decomposition. The activation energy of a composite system (PF) consisting of formic acid and PMIA was 106.94 kJ/mol at 2000-3000 K, which is 11.95% lower than that of pure PMIA. The main small-molecule products of the PF system were H/C/O-containing molecules (H2O, CO, and CO2), hydrocarbon molecules (e.g., CH4, •C2H, C2H4, and C3H4), N-containing molecules (N2, NH3, and HCN), H2, and various free radicals. Formic acid can promote the production of small molecules such as CO, CO2, and H2O. The N-H bonds, C-N bonds and the amide C═O double bond of PMIA were vulnerable to CO, H ions, and free radicals produced by formic acid decomposition, and this decreased the PMIA stability. Temperature is an important factor in the thermal decomposition of PMIA and can accelerate reactions in the PF system. The initial reaction rate of PMIA at 3000 K was 8.1 times that at 2000 K, and the intermediate reaction rate was 6.2 times that at 2200 K; temperature also affects the types of pyrolysis products, for example, hydrocarbons are high-temperature products.
Collapse
Affiliation(s)
- Fei Yin
- College of Engineering and Technology , Southwest University , Chongqing 400715 , China
| | - Chao Tang
- College of Engineering and Technology , Southwest University , Chongqing 400715 , China
| | - Yujing Tang
- College of Engineering and Technology , Southwest University , Chongqing 400715 , China
| | - Yingang Gui
- College of Engineering and Technology , Southwest University , Chongqing 400715 , China
| | - Zhongyong Zhao
- College of Engineering and Technology , Southwest University , Chongqing 400715 , China
| |
Collapse
|
10
|
Affiliation(s)
- Vito Genna
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Elisa Donati
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163, Genoa, Italy
| |
Collapse
|
11
|
Yin F, Tang C, Wang Q, Liu X, Tang Y. Molecular Dynamics Simulations on the Thermal Decomposition of Meta-Aramid Fibers. Polymers (Basel) 2018; 10:E691. [PMID: 30960616 PMCID: PMC6403965 DOI: 10.3390/polym10070691] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2018] [Revised: 06/17/2018] [Accepted: 06/19/2018] [Indexed: 11/17/2022] Open
Abstract
The thermal decomposition mechanism of a meta-aramid fiber was simulated at the atomic level using the ReaxFF reactive force field. The simulation results indicated that the main initial decomposition positions of the meta-aramid fiber elements were Caromatic ring⁻N and C=O, which could be used as targets for the modification of meta-aramid fibers. The meta-aramid fiber elements first decomposed into C6⁻C13 and then into smaller segments and micromolecular gases. The temperature was shown to be the key factor affecting the thermal decomposition of the meta-aramid fibers. More complex compositions and stable gases were produced at high temperatures than at lower temperatures. HCN was a decomposition product at high temperature, suggesting that its presence could be used for detecting thermal faults in meta-aramid fibers. Generation path tracing of the thermal decomposition products NH₃ and H₂O was also performed. NH₃ was produced when the NH₂ group captured an H atom adjacent to the system. H₂O was formed after a carbonyl group captured an H atom, became a hydroxyl group, with subsequent intramolecular dehydration or intermolecular hydrogen abstraction.
Collapse
Affiliation(s)
- Fei Yin
- College of Engineering and Technology, Southwest University, Chongqing 400715, China.
| | - Chao Tang
- College of Engineering and Technology, Southwest University, Chongqing 400715, China.
| | - Qian Wang
- Electric Power Research Institute of State Grid Chongqing Electric Power Company, Chongqing 401123, China.
| | - Xiong Liu
- Electric Power Research Institute of State Grid Chongqing Electric Power Company, Chongqing 401123, China.
| | - Yujing Tang
- College of Engineering and Technology, Southwest University, Chongqing 400715, China.
| |
Collapse
|
12
|
Manjari SR, Banavali NK. Structural Articulation of Biochemical Reactions Using Restrained Geometries and Topology Switching. J Chem Inf Model 2018; 58:453-463. [PMID: 29357231 DOI: 10.1021/acs.jcim.7b00699] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A strategy named "restrained geometries and topology switching" (RGATS) is presented to obtain detailed trajectories for complex biochemical reactions using molecular mechanics (MM) methods. It enables prediction of realistic dynamical pathways for chemical reactions, especially for accurately characterizing the structural adjustments of highly complex environments to any proximal biochemical reaction. It can be used to generate reactive conformations, model stepwise or concerted reactions in complex environments, and probe the influence of changes in the environment. Its ability to take reactively nonoptimal conformations and generate favorable starting conformations for a biochemical reaction is illustrated for a proton transfer between two model compounds. Its ability to study concerted reactions in explicit solvent is illustrated using proton transfers between an ammonium ion and two conserved histidines in an ammonia transporter channel embedded in a lipid membrane. Its ability to characterize the changes induced by subtle differences in the active site environment is illustrated using nucleotide addition by a DNA polymerase in the presence of two versus three Mg2+ ions. RGATS can be employed within any MM program and requires no additional software implementation. This allows the full assortment of computational methods implemented in all available MM programs to be used to tackle virtually any question about biochemical reactions that is answerable without using a quantum mechanical (QM) model. It can also be applied to generate reasonable starting structures for more detailed and expensive QM or QM/MM methods. In particular, this strategy enables rapid prediction of reactant, intermediary, or product state structures in any macromolecular context, with the only requirement being that the structure in any one of these states is either known or can be accurately modeled.
Collapse
Affiliation(s)
- Swati R Manjari
- Laboratory of Computational and Structural Biology, Division of Genetics, NYS Department of Health , CMS 2008, Biggs Laboratory, Wadsworth Center, Empire State Plaza, Albany, New York 12201-0509, United States
| | - Nilesh K Banavali
- Laboratory of Computational and Structural Biology, Division of Genetics, NYS Department of Health , CMS 2008, Biggs Laboratory, Wadsworth Center, Empire State Plaza, Albany, New York 12201-0509, United States.,Department of Biomedical Sciences, School of Public Health, State University of New York at Albany , Albany, New York 12222, United States
| |
Collapse
|
13
|
Gillet N, Lévy B, Moliner V, Demachy I, de la Lande A. Theoretical estimation of redox potential of biological quinone cofactors. J Comput Chem 2017; 38:1612-1621. [PMID: 28470751 DOI: 10.1002/jcc.24802] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2016] [Revised: 03/09/2017] [Accepted: 03/10/2017] [Indexed: 11/10/2022]
Abstract
Redox potentials are essential to understand biological cofactor reactivity and to predict their behavior in biological media. Experimental determination of redox potential in biological system is often difficult due to complexity of biological media but computational approaches can be used to estimate them. Nevertheless, the quality of the computational methodology remains a key issue to validate the results. Instead of looking to the best absolute results, we present here the calibration of theoretical redox potential for quinone derivatives in water coupling QM + MM or QM/MM scheme. Our approach allows using low computational cost theoretical level, ideal for long simulations in biological systems, and determination of the uncertainties linked to the calculations. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Natacha Gillet
- Laboratoire de Chimie-Physique, Université Paris Sud, CNRS, UMR 8000. 15, rue Jean Perrin, 91405 Orsay, CEDEX, France.,Departament de Química Física i Analítica, Universitat Jaume I, Castellón, 12071, Spain
| | - Bernard Lévy
- Laboratoire de Chimie-Physique, Université Paris Sud, CNRS, UMR 8000. 15, rue Jean Perrin, 91405 Orsay, CEDEX, France
| | - Vicent Moliner
- Departament de Química Física i Analítica, Universitat Jaume I, Castellón, 12071, Spain.,Department of Chemistry, University of Bath, Bath BA2 7AY, United Kingdom
| | - Isabelle Demachy
- Laboratoire de Chimie-Physique, Université Paris Sud, CNRS, UMR 8000. 15, rue Jean Perrin, 91405 Orsay, CEDEX, France
| | - Aurélien de la Lande
- Laboratoire de Chimie-Physique, Université Paris Sud, CNRS, UMR 8000. 15, rue Jean Perrin, 91405 Orsay, CEDEX, France
| |
Collapse
|
14
|
Poberžnik M, Purg M, Repič M, Mavri J, Vianello R. Empirical Valence Bond Simulations of the Hydride-Transfer Step in the Monoamine Oxidase A Catalyzed Metabolism of Noradrenaline. J Phys Chem B 2016; 120:11419-11427. [PMID: 27734680 DOI: 10.1021/acs.jpcb.6b09011] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines, such as dopamine, serotonin, and noradrenaline (NA), which is why they have been extensively implicated in the etiology and course of various neurodegenerative disorders and, accordingly, used as primary pharmacological targets to treat these debilitating cognitive diseases. The precise chemical mechanism through which MAOs regulate the amine concentration, which is vital for the development of novel inhibitors, is still not unambiguously determined in the literature. In this work, we present atomistic empirical valence bond simulations of the rate-limiting step of the MAO-A-catalyzed NA (norepinephrine) degradation, involving hydride transfer from the substrate α-methylene group to the flavin moiety of the flavin adenine dinucleotide prosthetic group, employing the full dimensionality and thermal fluctuations of the hydrated enzyme, with extensive configurational sampling. We show that MAO-A lowers the free energy of activation by 14.3 kcal mol-1 relative to that of the same reaction in aqueous solution, whereas the calculated activation free energy of ΔG‡ = 20.3 ± 1.6 kcal mol-1 is found to be in reasonable agreement with the correlated experimental value of 16.5 kcal mol-1. The results presented here strongly support the fact that both MAO-A and MAO-B isoforms function by the same hydride-transfer mechanism. We also considered a few point mutations of the "aromatic cage" tyrosine residue (Tyr444Phe, Tyr444Leu, Tyr444Trp, Tyr444His, and Tyr444Glu), and the calculated changes in the reaction barriers are in agreement with the experimental values, thus providing further support to the proposed mechanism.
Collapse
Affiliation(s)
- Matic Poberžnik
- Department of Physical and Organic Chemistry, Jožef Stefan Institute , Jamova cesta 39, SI-1000 Ljubljana, Slovenia
| | - Miha Purg
- Department of Cell and Molecular Biology, Uppsala Biomedical Centre , Husargatan 3, S-75124 Uppsala, Sweden
| | - Matej Repič
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry , Hajdrihova ulica 19, SI-1000 Ljubljana, Slovenia
| | - Janez Mavri
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry , Hajdrihova ulica 19, SI-1000 Ljubljana, Slovenia
| | - Robert Vianello
- Computational Organic Chemistry and Biochemistry Group, Ruđer Bošković Institute , Bijenička cesta 54, HR-10000 Zagreb, Croatia
| |
Collapse
|
15
|
Litman YE, Videla PE, Rodriguez J, Laria D. Positional Isotope Exchange in HX·(H2O)n (X = F, I) Clusters at Low Temperatures. J Phys Chem A 2016; 120:7213-24. [DOI: 10.1021/acs.jpca.6b06681] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yair E. Litman
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Pablo E. Videla
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
| | - Javier Rodriguez
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
- ECyT, UNSAM, Martín
de Irigoyen 3100, 1650 San Martín, Pcia. de Buenos Aires, Argentina
| | - Daniel Laria
- Departamento
de Química Inorgánica, Analítica y Química-Física
and INQUIMAE-CONICET, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires Ciudad Universitaria, Pabellón II, 1428 Buenos Aires, Argentina
- Departamento
de Física de la Materia Condensada, Comisión Nacional de Energía Atómica, Avenida Libertador 8250, 1429 Buenos Aires, Argentina
| |
Collapse
|
16
|
Sousa SF, Ribeiro AJM, Neves RPP, Brás NF, Cerqueira NMFSA, Fernandes PA, Ramos MJ. Application of quantum mechanics/molecular mechanics methods in the study of enzymatic reaction mechanisms. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2016. [DOI: 10.1002/wcms.1281] [Citation(s) in RCA: 100] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Sérgio Filipe Sousa
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - António J. M. Ribeiro
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - Rui P. P. Neves
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - Natércia F. Brás
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - Nuno M. F. S. A. Cerqueira
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - Pedro A. Fernandes
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| | - Maria João Ramos
- Departamento de Química e Bioquímica, Faculdade de Ciências UCIBIO, REQUIMTE, Universidade do Porto Porto Portugal
| |
Collapse
|
17
|
Repič M, Vianello R, Purg M, Duarte F, Bauer P, Kamerlin SCL, Mavri J. Empirical valence bond simulations of the hydride transfer step in the monoamine oxidase B catalyzed metabolism of dopamine. Proteins 2014; 82:3347-55. [PMID: 25220264 DOI: 10.1002/prot.24690] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2014] [Accepted: 09/10/2014] [Indexed: 01/05/2023]
Abstract
Monoamine oxidases (MAOs) A and B are flavoenzymes responsible for the metabolism of biogenic amines such as dopamine, serotonin and noradrenaline. In this work, we present a comprehensive study of the rate-limiting step of dopamine degradation by MAO B, which consists in the hydride transfer from the methylene group of the substrate to the flavin moiety of the FAD prosthetic group. This article builds on our previous quantum chemical study of the same reaction using a cluster model (Vianello et al., Eur J Org Chem 2012; 7057), but now considering the full dimensionality of the hydrated enzyme with extensive configurational sampling. We show that MAO B is specifically tuned to catalyze the hydride transfer step from the substrate to the flavin moiety of the FAD prosthetic group and that it lowers the activation barrier by 12.3 kcal mol⁻¹ compared to the same reaction in aqueous solution, a rate enhancement of more than nine orders of magnitude. Taking into account the deprotonation of the substrate prior to the hydride transfer reaction, the activation barrier in the enzyme is calculated to be 16.1 kcal mol⁻¹, in excellent agreement with the experimental value of 16.5 kcal mol⁻¹. Additionally, we demonstrate that the protonation state of the active site residue Lys296 does not have an influence on the hydride transfer reaction.
Collapse
Affiliation(s)
- Matej Repič
- Laboratory for Biocomputing and Bioinformatics, National Institute of Chemistry, Hajdrihova 19, SI-1000, Ljubljana, Slovenia
| | | | | | | | | | | | | |
Collapse
|
18
|
Mlýnský V, Banáš P, Šponer J, van der Kamp MW, Mulholland AJ, Otyepka M. Comparison of ab Initio, DFT, and Semiempirical QM/MM Approaches for Description of Catalytic Mechanism of Hairpin Ribozyme. J Chem Theory Comput 2014; 10:1608-22. [PMID: 26580373 DOI: 10.1021/ct401015e] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
We have analyzed the capability of state-of-the-art multiscale computational approaches to provide atomic-resolution electronic structure insights into possible catalytic scenarios of the hairpin ribozyme by evaluating potential and free energy surfaces of the reactions by various hybrid QM/MM methods. The hairpin ribozyme is a unique catalytic RNA that achieves rate acceleration similar to other small self-cleaving ribozymes but without direct metal ion participation. Guanine 8 (G8) and adenine 38 (A38) have been identified as the catalytically essential nucleobases. However, their exact catalytic roles are still being investigated. In line with the available experimental data, we considered two reaction scenarios involving protonated A38H(+) as a general acid which is further assisted by either canonical G8 or deprotonated G8(-) forms. We used the spin-component scaled Møller-Plesset (SCS-MP2) method at the complete basis set limit as the reference method. The semiempirical AM1/d-PhoT and SCC-DFTBPR methods provided acceptable activation barriers with respect to the SCS-MP2 data but predicted significantly different reaction pathways. DFT functionals (BLYP and MPW1K) yielded the same reaction pathway as the SCS-MP2 method. The activation barriers were slightly underestimated by the GGA BLYP functional, although with accuracy comparable to the semiempirical methods. The SCS-MP2 method and hybrid MPW1K functional gave activation barriers that were closest to those derived from experimentally measured rate constants.
Collapse
Affiliation(s)
- Vojtěch Mlýnský
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics , Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
- CEITEC-Central European Institute of Technology, Masaryk University , Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Marc W van der Kamp
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, U.K
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Cantock's Close, Bristol BS8 1TS, U.K
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| |
Collapse
|
19
|
Chakraborty S. A quantitative measure of electrostatic perturbation in holo and apo enzymes induced by structural changes. PLoS One 2013; 8:e59352. [PMID: 23516628 PMCID: PMC3597595 DOI: 10.1371/journal.pone.0059352] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2012] [Accepted: 02/13/2013] [Indexed: 11/19/2022] Open
Abstract
Biological pathways are subject to subtle manipulations that achieve a wide range of functional variation in differing physiological niches. In many instances, changes in the structure of an enzyme on ligand binding germinate electrostatic perturbations that form the basis of its changed catalytic or transcriptional efficiency. Computational methods that seek to gain insights into the electrostatic changes in enzymes require expertise to setup and computing prowess. In the current work, we present a fast, easy and reliable methodology to compute electrostatic perturbations induced by ligand binding (MEPP). The theoretical foundation of MEPP is the conserved electrostatic potential difference (EPD) in cognate pairs of active site residues in proteins with the same functionality. Previously, this invariance has been used to unravel promiscuous serine protease and metallo-β-lactamase scaffolds in alkaline phosphatases. Given that a similarity in EPD is significant, we expect differences in the EPD to be significant too. MEPP identifies residues or domains that undergo significant electrostatic perturbations, and also enumerates residue pairs that undergo significant polarity change. The gain in a certain polarity of a residue with respect to neighboring residues, or the reversal of polarity between two residues might indicate a change in the preferred ligand. The methodology of MEPP has been demonstrated on several enzymes that employ varying mechanisms to perform their roles. For example, we have attributed the change in polarity in residue pairs to be responsible for the loss of metal ion binding in fructose 1,6-bisphosphatases, and corroborated the pre-organized state of the active site of the enzyme with respect to functionally relevant changes in electric fields in ketosteroid isomerases.
Collapse
Affiliation(s)
- Sandeep Chakraborty
- Department of Biological Sciences, Tata Institute of Fundamental Research, Mumbai, India.
| |
Collapse
|
20
|
Labas A, Szabo E, Mones L, Fuxreiter M. Optimization of reorganization energy drives evolution of the designed Kemp eliminase KE07. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2013; 1834:908-17. [PMID: 23380188 DOI: 10.1016/j.bbapap.2013.01.005] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/16/2012] [Revised: 01/06/2013] [Accepted: 01/08/2013] [Indexed: 10/27/2022]
Abstract
Understanding enzymatic evolution is essential to engineer enzymes with improved activities or to generate enzymes with tailor-made activities. The computationally designed Kemp eliminase KE07 carries out an unnatural reaction by converting of 5-nitrobenzisoxazole to cyanophenol, but its catalytic efficiency is significantly lower than those of natural enzymes. Three series of designed Kemp eliminases (KE07, KE70, KE59) were shown to be evolvable with considerable improvement in catalytic efficiency. Here we use the KE07 enzyme as a model system to reveal those forces, which govern enzymatic evolution and elucidate the key factors for improving activity. We applied the Empirical Valence Bond (EVB) method to construct the free energy pathway of the reaction in the original KE07 design and the evolved R7 1/3H variant. We analyzed catalytic effect of residues and demonstrated that not all mutations in evolution are favorable for activity. In contrast to the small decrease in the activation barrier, in vitro evolution significantly reduced the reorganization energy. We developed an algorithm to evaluate group contributions to the reorganization energy and used this approach to screen for KE07 variants with potential for improvement. We aimed to identify those mutations that facilitate enzymatic evolution, but might not directly increase catalytic efficiency. Computational results in accord with experimental data show that all mutations, which appear during in vitro evolution were either neutral or favorable for the reorganization energy. These results underscore that distant mutations can also play role in optimizing efficiency via their contribution to the reorganization energy. Exploiting this principle could be a promising strategy for computer-aided enzyme design. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.
Collapse
Affiliation(s)
- A Labas
- DE OEC-Momentum Laboratory of Protein dynamics, University of Debrecen, Debrecen, Hungary
| | | | | | | |
Collapse
|
21
|
Mones L, Kulhánek P, Simon I, Laio A, Fuxreiter M. The energy gap as a universal reaction coordinate for the simulation of chemical reactions. J Phys Chem B 2009; 113:7867-73. [PMID: 19432459 DOI: 10.1021/jp9000576] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The selection of a proper reaction coordinate is a major bottleneck in simulations of chemical reactions in complex systems. Increasing the number of variables that are used to bias the reaction largely affects the convergence and leads to an unbearable increase in computational price. This problem can be overcome by employing a complex reaction coordinate that depends on many geometrical variables of the system, such as the energy gap (EGAP) in the empirical valence bond (EVB) method. EGAP depends on all of the coordinates of the system, and its robustness has been demonstrated for a variety of enzymatic reactions. In this work, we demonstrate that EGAP, derived from a classical representation, can be used as a reaction coordinate in systems described with any quantum chemistry Hamiltonian. Benefits of using EGAP as a reaction coordinate as compared to a traditional geometrical variable are illustrated in the case of a symmetric nucleophilic substitution reaction in water solution. EGAP is shown to provide a significantly more efficient sampling and allows a better localization of the transition state as compared to a geometrical reaction coordinate.
Collapse
Affiliation(s)
- Letif Mones
- Institute of Enzymology, Biological Research Center, Hungarian Academy of Sciences, P.O. Box 7, H-1518 Budapest, Hungary
| | | | | | | | | |
Collapse
|
22
|
Gao J, Wong KY, Major DT, Cembran A, Song L, Lin YL, Fan Y, Ma S. Kinetic Isotope Effects from Hybrid Classical and Quantum Path Integral Computations. QUANTUM TUNNELLING IN ENZYME-CATALYSED REACTIONS 2009. [DOI: 10.1039/9781847559975-00105] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Jiali Gao
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Kin-Yiu Wong
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Dan T. Major
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Alessandro Cembran
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Lingchun Song
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Yen-lin Lin
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Yao Fan
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| | - Shuhua Ma
- Department of Chemistry, Digital Technology Center and Supercomputer Institute, University of Minnesota 207 Pleasant Street S.E. Minneapolis MN 55455-0431 USA
| |
Collapse
|
23
|
Abstract
This review discusses methods for the incorporation of quantum mechanical effects into enzyme kinetics simulations in which the enzyme is an explicit part of the model. We emphasize three aspects: (a) use of quantum mechanical electronic structure methods such as molecular orbital theory and density functional theory, usually in conjunction with molecular mechanics; (b) treating vibrational motions quantum mechanically, either in an instantaneous harmonic approximation, or by path integrals, or by a three-dimensional wave function coupled to classical nuclear motion; (c) incorporation of multidimensional tunneling approximations into reaction rate calculations.
Collapse
Affiliation(s)
- Jiali Gao
- Department of Chemistry and Supercomputer Institute, University of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431, USA.
| | | |
Collapse
|
24
|
Robson B. The design of biologically active polypeptides. CRC CRITICAL REVIEWS IN BIOCHEMISTRY 1983; 14:273-96. [PMID: 6357633 DOI: 10.3109/10409238309102796] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
On the one hand considerable advances are being made in the computer-aided calculation of polypeptide structure and behaviour, while on the other even very simple reasoning has lead to some apparent successes in the design and synthesis of biologically active peptides. While this looks extremely promising, there is a case for arguing that progress made in both these areas is to some degree illusory. Such criticisms as can be made, however, turn out to be highly constructive, suggesting a profitable and still workable approach to the problem of rational polypeptide design.
Collapse
|