1
|
Kim SM, Kang SH, Jeon BW, Kim YH. Tunnel engineering of gas-converting enzymes for inhibitor retardation and substrate acceleration. BIORESOURCE TECHNOLOGY 2024; 394:130248. [PMID: 38158090 DOI: 10.1016/j.biortech.2023.130248] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Revised: 12/20/2023] [Accepted: 12/21/2023] [Indexed: 01/03/2024]
Abstract
Carbon monoxide dehydrogenase (CODH), formate dehydrogenase (FDH), hydrogenase (H2ase), and nitrogenase (N2ase) are crucial enzymatic catalysts that facilitate the conversion of industrially significant gases such as CO, CO2, H2, and N2. The tunnels in the gas-converting enzymes serve as conduits for these low molecular weight gases to access deeply buried catalytic sites. The identification of the substrate tunnels is imperative for comprehending the substrate selectivity mechanism underlying these gas-converting enzymes. This knowledge also holds substantial value for industrial applications, particularly in addressing the challenges associated with separation and utilization of byproduct gases. In this comprehensive review, we delve into the emerging field of tunnel engineering, presenting a range of approaches and analyses. Additionally, we propose methodologies for the systematic design of enzymes, with the ultimate goal of advancing protein engineering strategies.
Collapse
Affiliation(s)
- Suk Min Kim
- School of Energy and Chemical Engineering, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea.
| | - Sung Heuck Kang
- School of Energy and Chemical Engineering, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Byoung Wook Jeon
- School of Energy and Chemical Engineering, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Yong Hwan Kim
- School of Energy and Chemical Engineering, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea; Graduate School of Carbon Neutrality, Ulsan National Institute of Science and Technology (UNIST), 50 UNIST-gil, Ulsan 44919, Republic of Korea.
| |
Collapse
|
2
|
Kurt F, Filiz E, Yildiz K, Akbudak MA. Genome-Wide Identification, Characterization and Expression Profiling of Potato ( Solanum tuberosum) Frataxin ( FH) Gene. Genes (Basel) 2023; 14:468. [PMID: 36833395 PMCID: PMC9957314 DOI: 10.3390/genes14020468] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2023] [Revised: 02/06/2023] [Accepted: 02/10/2023] [Indexed: 02/15/2023] Open
Abstract
Frataxin (FH) plays a crucial role in the biogenesis of mitochondria and the regulation of iron in the cells of various organisms. However, there has been very little research on FH in plants. In this study, the potato FH gene (StFH) was identified and characterized using a genome-wide approach, and its sequence was compared to those of FH genes from Arabidopsis, rice, and maize. The FH genes were found to have a lineage-specific distribution and were more conserved in monocots than in dicots. While multiple copies of FH genes have been reported in some species, including plants, only one isoform of FH was found in potato. The expression of StFH in leaves and roots was analyzed under two different abiotic stress conditions, and the results showed that StFH was upregulated more in leaves and that its expression levels increased with the severity of the stress. This is the first study to examine the expression of an FH gene under abiotic stress conditions.
Collapse
Affiliation(s)
- Firat Kurt
- Department of Plant Production and Technologies, Faculty of Applied Sciences, Mus Alparslan University, 49250 Mus, Turkey
| | - Ertugrul Filiz
- Department of Crop and Animal Production, Cilimli Vocational School, Duzce University, Cilimli, 81750 Duzce, Turkey
| | - Kubra Yildiz
- Department of Agricultural Biotechnology, Akdeniz University, 07058 Antalya, Turkey
| | - M. Aydın Akbudak
- Department of Agricultural Biotechnology, Akdeniz University, 07058 Antalya, Turkey
| |
Collapse
|
3
|
Yadahalli S, Jayanthi LP, Gosavi S. A Method for Assessing the Robustness of Protein Structures by Randomizing Packing Interactions. Front Mol Biosci 2022; 9:849272. [PMID: 35832734 PMCID: PMC9271847 DOI: 10.3389/fmolb.2022.849272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Accepted: 04/27/2022] [Indexed: 12/02/2022] Open
Abstract
Many single-domain proteins are not only stable and water-soluble, but they also populate few to no intermediates during folding. This reduces interactions between partially folded proteins, misfolding, and aggregation, and makes the proteins tractable in biotechnological applications. Natural proteins fold thus, not necessarily only because their structures are well-suited for folding, but because their sequences optimize packing and fit their structures well. In contrast, folding experiments on the de novo designed Top7 suggest that it populates several intermediates. Additionally, in de novo protein design, where sequences are designed for natural and new non-natural structures, tens of sequences still need to be tested before success is achieved. Both these issues may be caused by the specific scaffolds used in design, i.e., some protein scaffolds may be more tolerant to packing perturbations and varied sequences. Here, we report a computational method for assessing the response of protein structures to packing perturbations. We then benchmark this method using designed proteins and find that it can identify scaffolds whose folding gets disrupted upon perturbing packing, leading to the population of intermediates. The method can also isolate regions of both natural and designed scaffolds that are sensitive to such perturbations and identify contacts which when present can rescue folding. Overall, this method can be used to identify protein scaffolds that are more amenable to whole protein design as well as to identify protein regions which are sensitive to perturbations and where further mutations should be avoided during protein engineering.
Collapse
|
4
|
Cammarota RC, Liu W, Bacsa J, Davies HML, Sigman MS. Mechanistically Guided Workflow for Relating Complex Reactive Site Topologies to Catalyst Performance in C–H Functionalization Reactions. J Am Chem Soc 2022; 144:1881-1898. [DOI: 10.1021/jacs.1c12198] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Ryan C. Cammarota
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| | - Wenbin Liu
- Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| | - John Bacsa
- Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| | - Huw M. L. Davies
- Department of Chemistry, Emory University, 1515 Dickey Drive, Atlanta, Georgia 30322, United States
| | - Matthew S. Sigman
- Department of Chemistry, University of Utah, 315 South 1400 East, Salt Lake City, Utah 84112, United States
| |
Collapse
|
5
|
Zhu C, Zhang AM, Li Y, Li HX, Qian Y, Fu Y, Wu X, Li Y. A biomimetic metal–organic framework with cuboid inner cavities for enantioselective separation. Inorg Chem Front 2022. [DOI: 10.1039/d2qi00152g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A biomimetic metal–organic framework with cuboid inner cavities and multiple recognition sites was constructed from a phenylalanine-derived ligand. It can enantioselectively separate various racemic alcohols, diols and epoxides with ee up to 99.5%.
Collapse
Affiliation(s)
- Chengfeng Zhu
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - A-Mei Zhang
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Ying Li
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Han-Xue Li
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Yijian Qian
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Yanming Fu
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Xiang Wu
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| | - Yougui Li
- Anhui Province Key Laboratory of Advanced Catalytic Materials and Reaction Engineering, School of Chemistry and Chemical Engineering, Hefei University of Technology, Hefei 230009, P. R. China
| |
Collapse
|
6
|
Şterbuleac D. Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators. RSC Med Chem 2021; 12:1503-1518. [PMID: 34671734 PMCID: PMC8459385 DOI: 10.1039/d1md00140j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2021] [Accepted: 07/21/2021] [Indexed: 01/10/2023] Open
Abstract
Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.
Collapse
Affiliation(s)
- Daniel Şterbuleac
- Department of Health and Human Development, "Ştefan cel Mare" University of Suceava Str. Universităţii 13, 720229, E Building Suceava Romania
- Department of Forestry and Environmental Protection, "Ştefan cel Mare" University of Suceava Str. Universităţii 13, 720229, E Building Suceava Romania
- Integrated Center for Research, Development and Innovation in Advanced Materials, Nanotechnologies and Distributed Systems for Fabrication and Control (MANSiD), "Ştefan cel Mare" University of Suceava Str. Universităţii 13 720229 Suceava Romania
| |
Collapse
|
7
|
Hollander M, Rasp D, Aziz M, Helms V. ProPores2: Web Service and Stand-Alone Tool for Identifying, Manipulating, and Visualizing Pores in Protein Structures. J Chem Inf Model 2021; 61:1555-1559. [PMID: 33844545 DOI: 10.1021/acs.jcim.1c00154] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Surface pockets, cavities, and tunnels in the 3D structures of proteins play integral functional roles such as enabling enzymatic catalysis, ligand binding, or transport of ions or small molecules across biomembranes. ProPores2 facilitates understanding and analysis of these processes by identifying pores and lining residues, determining their axes, and opening closed connections via side-chain rotation. The fast stand-alone tool introduces a novel mode for pore identification, improved axis determination, and additional features such as parallel batch processing and a graphical user interface. The new web service features an integrated and customizable protein viewer with an option to analyze and view more than one structure at once. This feature facilitates side-by-side comparisons of pores in different conformations of the same protein or of identified pores before and after opening gates within the same protein. ProPores2 is freely and publicly available at https://service.bioinformatik.uni-saarland.de/propores.
Collapse
Affiliation(s)
- Markus Hollander
- Center for Bioinformatics, Saarland Informatics Campus, Saarland University, D-66041 Saarbrücken, Germany
| | - David Rasp
- Center for Bioinformatics, Saarland Informatics Campus, Saarland University, D-66041 Saarbrücken, Germany
| | - Moomal Aziz
- Center for Bioinformatics, Saarland Informatics Campus, Saarland University, D-66041 Saarbrücken, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland Informatics Campus, Saarland University, D-66041 Saarbrücken, Germany
| |
Collapse
|
8
|
Meng S, An R, Li Z, Schwaneberg U, Ji Y, Davari MD, Wang F, Wang M, Qin M, Nie K, Liu L. Tunnel engineering for modulating the substrate preference in cytochrome P450 BsβHI. BIORESOUR BIOPROCESS 2021; 8:26. [PMID: 38650198 PMCID: PMC10992877 DOI: 10.1186/s40643-021-00379-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2020] [Accepted: 03/25/2021] [Indexed: 01/07/2023] Open
Abstract
An active site is normally located inside enzymes, hence substrates should go through a tunnel to access the active site. Tunnel engineering is a powerful strategy for refining the catalytic properties of enzymes. Here, P450BsβHI (Q85H/V170I) derived from hydroxylase P450Bsβ from Bacillus subtilis was chosen as the study model, which is reported as a potential decarboxylase. However, this enzyme showed low decarboxylase activity towards long-chain fatty acids. Here, a tunnel engineering campaign was performed for modulating the substrate preference and improving the decarboxylation activity of P450BsβHI. The finally obtained BsβHI-F79A variant had a 15.2-fold improved conversion for palmitic acid; BsβHI-F173V variant had a 3.9-fold improved conversion for pentadecanoic acid. The study demonstrates how the substrate preference can be modulated by tunnel engineering strategy.
Collapse
Affiliation(s)
- Shuaiqi Meng
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Ruipeng An
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Zhongyu Li
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Ulrich Schwaneberg
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
- DWI-Leibniz Institute for Interactive Materials, Forckenbeckstraße 50, 52074, Aachen, Germany
| | - Yu Ji
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Mehdi D Davari
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, 52074, Aachen, Germany
| | - Fang Wang
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Meng Wang
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Meng Qin
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Kaili Nie
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China
| | - Luo Liu
- Beijing Bioprocess Key Laboratory, Beijing University of Chemical Technology, Beijing, 100029, People's Republic of China.
| |
Collapse
|
9
|
Chwastyk M, Panek EA, Malinowski J, Jaskólski M, Cieplak M. Properties of Cavities in Biological Structures-A Survey of the Protein Data Bank. Front Mol Biosci 2020; 7:591381. [PMID: 33240933 PMCID: PMC7677499 DOI: 10.3389/fmolb.2020.591381] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2020] [Accepted: 10/06/2020] [Indexed: 11/13/2022] Open
Abstract
We performed a PDB-wide survey of proteins to assess their cavity content, using the SPACEBALL algorithm to calculate the cavity volumes. In addition, we determined the hydropathy character of the cavities. We demonstrate that the cavities of most proteins are hydrophilic, but smaller proteins tend to have cavities with hydrophobic walls. We propose criteria for distinguishing between cavities and pockets, and single out proteins with the largest cavities.
Collapse
Affiliation(s)
- Mateusz Chwastyk
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| | - Ewa A Panek
- Department of Biometry, Warsaw University of Life Sciences, Warsaw, Poland
| | - Jan Malinowski
- Faculty of Physics, University of Warsaw, Warsaw, Poland
| | - Mariusz Jaskólski
- Center for Biocrystallographic Research, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznań, Poland.,Department of Crystallography, Faculty of Chemistry, Adam Mickiewicz University, Poznań, Poland
| | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Warsaw, Poland
| |
Collapse
|
10
|
Xiao F, Song X, Tian P, Gan M, Verkhivker GM, Hu G. Comparative Dynamics and Functional Mechanisms of the CYP17A1 Tunnels Regulated by Ligand Binding. J Chem Inf Model 2020; 60:3632-3647. [PMID: 32530640 DOI: 10.1021/acs.jcim.0c00447] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
As an important member of cytochrome P450 (CYP) enzymes, CYP17A1 is a dual-function monooxygenase with a critical role in the synthesis of many human steroid hormones, making it an attractive therapeutic target. The emerging structural information about CYP17A1 and the growing number of inhibitors for these enzymes call for a systematic strategy to delineate and classify mechanisms of ligand transport through tunnels that control catalytic activity. In this work, we applied an integrated computational strategy to different CYP17A1 systems with a panel of ligands to systematically study at the atomic level the mechanism of ligand-binding and tunneling dynamics. Atomistic simulations and binding free energy computations identify the dynamics of dominant tunnels and characterize energetic properties of critical residues responsible for ligand binding. The common transporting pathways including S, 3, and 2c tunnels were identified in CYP17A1 binding systems, while the 2c tunnel is a newly formed pathway upon ligand binding. We employed and integrated several computational approaches including the analysis of functional motions and sequence conservation, atomistic modeling of dynamic residue interaction networks, and perturbation response scanning analysis to dissect ligand tunneling mechanisms. The results revealed the hinge-binding and sliding motions as main functional modes of the tunnel dynamic, and a group of mediating residues as key regulators of tunnel conformational dynamics and allosteric communications. We have also examined and quantified the mutational effects on the tunnel composition, conformational dynamics, and long-range allosteric behavior. The results of this investigation are fully consistent with the experimental data, providing novel rationale to the experiments and offering valuable insights into the relationships between the structure and function of the channel networks and a robust atomistic model of activation mechanisms and allosteric interactions in CYP enzymes.
Collapse
Affiliation(s)
- Fei Xiao
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Xingyu Song
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Peiyi Tian
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Mi Gan
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| | - Gennady M Verkhivker
- Department of Computational and Data Sciences, Chapman University, One University Drive, Orange, California 92866, United States.,Department of Biomedical and Pharmaceutical Sciences, Chapman University Pharmacy School, 9401 Jeronimo Rd, Irvine, California 92618, United States
| | - Guang Hu
- Center for Systems Biology, Department of Bioinformatics, School of Biology and Basic Medical Sciences, Soochow University, Suzhou 215123, China
| |
Collapse
|
11
|
Mitusińska K, Raczyńska A, Bzówka M, Bagrowska W, Góra A. Applications of water molecules for analysis of macromolecule properties. Comput Struct Biotechnol J 2020; 18:355-365. [PMID: 32123557 PMCID: PMC7036622 DOI: 10.1016/j.csbj.2020.02.001] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2019] [Revised: 01/26/2020] [Accepted: 02/01/2020] [Indexed: 01/12/2023] Open
Abstract
Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.
Collapse
Affiliation(s)
| | | | | | | | - Artur Góra
- Tunneling Group, Biotechnology Centre, Silesian University of Technology, Krzywoustego 8, Gliwice, Poland
| |
Collapse
|
12
|
Manak M. Voronoi-based detection of pockets in proteins defined by large and small probes. J Comput Chem 2019; 40:1758-1771. [PMID: 30932214 DOI: 10.1002/jcc.25828] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2018] [Revised: 03/04/2019] [Accepted: 03/05/2019] [Indexed: 01/06/2023]
Abstract
The function of enzymatic proteins is given by their ability to bind specific small molecules into their active sites. These sites can often be found in pockets on a hypothetical boundary between the protein and its environment. Detection, analysis, and visualization of pockets find its use in protein engineering and drug discovery. Many definitions of pockets and algorithms for their computation have been proposed. Kawabata and Go defined them as the regions of empty space into which a small spherical probe can enter but a large probe cannot and developed programs that can compute their approximate shape. In this article, this definition was slightly modified in order to capture the existence of large internal holes, and a Voronoi-based method for the computation of the exact shape of these modified regions is introduced. The method first puts a finite number of large probes on the protein exterior surface and then, considering both large probes and atomic balls as obstacles for the small probe, the method computes the exact shape of the regions for the small probe. This is all achieved with Voronoi diagrams, which help with the safe navigation of spherical probes among spherical obstacles. Detected regions are internally represented as graphs of vertices and edges describing possible movements of the center of the small probe on Voronoi edges. The surface bounding each region is obtained from this representation and used for visualization, volume estimation, and comparison with other approaches. © 2019 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Martin Manak
- New Technologies for the Information Society (NTIS), Faculty of Applied Sciences, University of West Bohemia, Univerzitni 8, 30614, Pilsen, Czech Republic
| |
Collapse
|
13
|
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
|
14
|
Four Major Channels Detected in the Cytochrome P450 3A4: A Step toward Understanding Its Multispecificity. Int J Mol Sci 2019; 20:ijms20040987. [PMID: 30823507 PMCID: PMC6412807 DOI: 10.3390/ijms20040987] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Revised: 02/12/2019] [Accepted: 02/20/2019] [Indexed: 12/27/2022] Open
Abstract
We computed the network of channels of the 3A4 isoform of the cytochrome P450 (CYP) on the basis of 16 crystal structures extracted from the Protein Data Bank (PDB). The calculations were performed with version 2 of the CCCPP software that we developed for this research project. We identified the minimal cost paths (MCPs) output by CCCPP as probable ways to access to the buried active site. The algorithm of calculation of the MCPs is presented in this paper, with its original method of visualization of the channels. We found that these MCPs constitute four major channels in CYP3A4. Among the many channels proposed by Cojocaru et al. in 2007, we found that only four of them open in 3A4. We provide a refined description of these channels together with associated quantitative data.
Collapse
|
15
|
Weiß RG, Chudoba R, Setny P, Dzubiella J. Affinity, kinetics, and pathways of anisotropic ligands binding to hydrophobic model pockets. J Chem Phys 2018; 149:094902. [DOI: 10.1063/1.5025118] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Affiliation(s)
- R. Gregor Weiß
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Laboratory of Physical Chemistry, ETH Zürich, CH-8093 Zürich, Switzerland
| | - Richard Chudoba
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Research Group Simulations of Energy Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Straße 3, D-79104 Freiburg, Germany
| | - Piotr Setny
- Centre of New Technologies, University of Warsaw, Stefana Banacha 2c, 00-927 Warsaw, Poland
| | - Joachim Dzubiella
- Institut für Physik, Humboldt-Universität zu Berlin, Newtonstraße 15, D-12489 Berlin, Germany
- Research Group Simulations of Energy Materials, Helmholtz-Zentrum Berlin, Hahn-Meitner-Platz 1, D-14109 Berlin, Germany
- Physikalisches Institut, Albert-Ludwigs-Universität Freiburg, Hermann-Herder Straße 3, D-79104 Freiburg, Germany
| |
Collapse
|
16
|
Ligand Access Channels in Cytochrome P450 Enzymes: A Review. Int J Mol Sci 2018; 19:ijms19061617. [PMID: 29848998 PMCID: PMC6032366 DOI: 10.3390/ijms19061617] [Citation(s) in RCA: 50] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2018] [Revised: 05/27/2018] [Accepted: 05/28/2018] [Indexed: 12/22/2022] Open
Abstract
Quantitative structure-activity relationships may bring invaluable information on structural elements of both enzymes and substrates that, together, govern substrate specificity. Buried active sites in cytochrome P450 enzymes are connected to the solvent by a network of channels exiting at the distal surface of the protein. This review presents different in silico tools that were developed to uncover such channels in P450 crystal structures. It also lists some of the experimental evidence that actually suggest that these predicted channels might indeed play a critical role in modulating P450 functions. Amino acid residues at the entrance of the channels may participate to a first global ligand recognition of ligands by P450 enzymes before they reach the buried active site. Moreover, different P450 enzymes show different networks of predicted channels. The plasticity of P450 structures is also important to take into account when looking at how channels might play their role.
Collapse
|
17
|
Simões T, Lopes D, Dias S, Fernandes F, Pereira J, Jorge J, Bajaj C, Gomes A. Geometric Detection Algorithms for Cavities on Protein Surfaces in Molecular Graphics: A Survey. COMPUTER GRAPHICS FORUM : JOURNAL OF THE EUROPEAN ASSOCIATION FOR COMPUTER GRAPHICS 2017; 36:643-683. [PMID: 29520122 PMCID: PMC5839519 DOI: 10.1111/cgf.13158] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
Detecting and analyzing protein cavities provides significant information about active sites for biological processes (e.g., protein-protein or protein-ligand binding) in molecular graphics and modeling. Using the three-dimensional structure of a given protein (i.e., atom types and their locations in 3D) as retrieved from a PDB (Protein Data Bank) file, it is now computationally viable to determine a description of these cavities. Such cavities correspond to pockets, clefts, invaginations, voids, tunnels, channels, and grooves on the surface of a given protein. In this work, we survey the literature on protein cavity computation and classify algorithmic approaches into three categories: evolution-based, energy-based, and geometry-based. Our survey focuses on geometric algorithms, whose taxonomy is extended to include not only sphere-, grid-, and tessellation-based methods, but also surface-based, hybrid geometric, consensus, and time-varying methods. Finally, we detail those techniques that have been customized for GPU (Graphics Processing Unit) computing.
Collapse
Affiliation(s)
- Tiago Simões
- Instituto de Telecomunicações, Portugal
- Universidade da Beira Interior, Portugal
| | | | - Sérgio Dias
- Instituto de Telecomunicações, Portugal
- Universidade da Beira Interior, Portugal
| | | | - João Pereira
- INESC-ID Lisboa, Portugal
- Instituto Superior Técnico, Universidade de Lisboa, Portugal
| | - Joaquim Jorge
- INESC-ID Lisboa, Portugal
- Instituto Superior Técnico, Universidade de Lisboa, Portugal
| | | | - Abel Gomes
- Instituto de Telecomunicações, Portugal
- Universidade da Beira Interior, Portugal
| |
Collapse
|
18
|
Cerisier N, Regad L, Triki D, Petitjean M, Flatters D, Camproux AC. Statistical Profiling of One Promiscuous Protein Binding Site: Illustrated by Urokinase Catalytic Domain. Mol Inform 2017; 36. [PMID: 28696518 DOI: 10.1002/minf.201700040] [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: 02/27/2017] [Accepted: 06/26/2017] [Indexed: 12/21/2022]
Abstract
While recent literature focuses on drug promiscuity, the characterization of promiscuous binding sites (ability to bind several ligands) remains to be explored. Here, we present a proteochemometric modeling approach to analyze diverse ligands and corresponding multiple binding sub-pockets associated with one promiscuous binding site to characterize protein-ligand recognition. We analyze both geometrical and physicochemical profile correspondences. This approach was applied to examine the well-studied druggable urokinase catalytic domain inhibitor binding site, which results in a large number of complex structures bound to various ligands. This approach emphasizes the importance of jointly characterizing pocket and ligand spaces to explore the impact of ligand diversity on sub-pocket properties and to establish their main profile correspondences. This work supports an interest in mining available 3D holo structures associated with a promiscuous binding site to explore its main protein-ligand recognition tendency.
Collapse
Affiliation(s)
- Natacha Cerisier
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| | - Leslie Regad
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| | - Dhoha Triki
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| | - Michel Petitjean
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| | - Delphine Flatters
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| | - Anne-Claude Camproux
- INSERM, UMRS-973, MTi,35, rue Hélène Brion, 75205, PARIS CEDEX 13.,University Paris Diderot, Sorbonne Paris Cité, UMRS-973, MTi
| |
Collapse
|
19
|
Benkaidali L, André F, Moroy G, Tangour B, Maurel F, Petitjean M. The Cytochrome P450 3A4 has three Major Conformations: New Clues to Drug Recognition by this Promiscuous Enzyme. Mol Inform 2017; 36. [PMID: 28685969 DOI: 10.1002/minf.201700044] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2017] [Accepted: 06/28/2017] [Indexed: 01/02/2023]
Abstract
We computed the channels of the 3A4 isoform of the cytochrome P450 3A4 (CYP) on the basis of 24 crystal structures extracted from the Protein Data Bank (PDB). We identified three major conformations (denoted C, O1 and O2) using an enhanced version of the CCCPP software that we developed for the present work, while only two conformations (C and O2 ) are considered in the literature. We established the flowchart of definition of these three conformations in function of the structural and physicochemical parameters of the ligand. The channels are characterized with qualitative and quantitative parameters, and not only with their surrounding secondary structures as it is usually done in the literature.
Collapse
Affiliation(s)
- Lydia Benkaidali
- ITODYS, CNRS UMR 7086, Université Paris Diderot, France.,Université de Carthage, Faculté des Sciences, Département de Chimie, Bizerte, Tunisie
| | - François André
- CEA/I2BC, CNRS UMR 9198, Université Paris-Saclay, France
| | - Gautier Moroy
- MTi, INSERM UMR-S 973, Université Paris Diderot, France
| | - Bahoueddine Tangour
- Unité de Recherche de Modélisation en Sciences Fondamentales et Didactique, BP244, Université de Tunis El Manar, 2092, Tunis, Tunisie
| | | | - Michel Petitjean
- MTi, INSERM UMR-S 973, Université Paris Diderot, France.,Epôle de génoinformatique, CNRS UMR 7592, Institut Jacques Monod, Paris, France
| |
Collapse
|
20
|
Cerisier N, Regad L, Triki D, Camproux AC, Petitjean M. Cavity Versus Ligand Shape Descriptors: Application to Urokinase Binding Pockets. J Comput Biol 2017; 24:1134-1137. [PMID: 28570103 PMCID: PMC5684670 DOI: 10.1089/cmb.2017.0061] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
We analyzed 78 binding pockets of the human urokinase plasminogen activator (uPA) catalytic domain extracted from a data set of crystallized uPA-ligand complexes. These binding pockets were computed with an original geometric method that does NOT involve any arbitrary parameter, such as cutoff distances, angles, and so on. We measured the deviation from convexity of each pocket shape with the pocket convexity index (PCI). We defined a new pocket descriptor called distributional sphericity coefficient (DISC), which indicates to which extent the protein atoms of a given pocket lie on the surface of a sphere. The DISC values were computed with the freeware PCI. The pocket descriptors and their high correspondences with ligand descriptors are crucial for polypharmacology prediction. We found that the protein heavy atoms lining the urokinases binding pockets are either located on the surface of their convex hull or lie close to this surface. We also found that the radii of the urokinases binding pockets and the radii of their ligands are highly correlated (r = 0.9).
Collapse
Affiliation(s)
- Natacha Cerisier
- 1 MTi, INSERM UMR-S 973, Université Paris Diderot , Paris, France
| | - Leslie Regad
- 1 MTi, INSERM UMR-S 973, Université Paris Diderot , Paris, France
| | - Dhoha Triki
- 1 MTi, INSERM UMR-S 973, Université Paris Diderot , Paris, France
| | | | - Michel Petitjean
- 1 MTi, INSERM UMR-S 973, Université Paris Diderot , Paris, France .,2 Epôle de Génoinformatique, Institut Jacques Monod, CNRS, UMR7592, Université Paris Diderot , Paris, France
| |
Collapse
|
21
|
Caumes G, Borrel A, Abi Hussein H, Camproux AC, Regad L. Investigating the Importance of the Pocket-estimation Method in Pocket-based Approaches: An Illustration Using Pocket-ligand Classification. Mol Inform 2017; 36. [PMID: 28452177 DOI: 10.1002/minf.201700025] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2017] [Accepted: 04/06/2017] [Indexed: 11/12/2022]
Abstract
Small molecules interact with their protein target on surface cavities known as binding pockets. Pocket-based approaches are very useful in all of the phases of drug design. Their first step is estimating the binding pocket based on protein structure. The available pocket-estimation methods produce different pockets for the same target. The aim of this work is to investigate the effects of different pocket-estimation methods on the results of pocket-based approaches. We focused on the effect of three pocket-estimation methods on a pocket-ligand (PL) classification. This pocket-based approach is useful for understanding the correspondence between the pocket and ligand spaces and to develop pharmacological profiling models. We found pocket-estimation methods yield different binding pockets in terms of boundaries and properties. These differences are responsible for the variation in the PL classification results that can have an impact on the detected correspondence between pocket and ligand profiles. Thus, we highlighted the importance of the pocket-estimation method choice in pocket-based approaches.
Collapse
Affiliation(s)
- Géraldine Caumes
- Molécules thérapeutiques In silico (MTi), INSERM UMR-S973, University Paris Diderot, 35 rue Hélène Brion, 75013, Paris Cedex, France.,IMPMC, UMR 7590, Equipe de Géobiologie, Université Pierre et Marie Curie, 4 place Jussieu, 75252, Paris Cedex, France
| | - Alexandre Borrel
- Molécules thérapeutiques In silico (MTi), INSERM UMR-S973, University Paris Diderot, 35 rue Hélène Brion, 75013, Paris Cedex, France.,Division of Pharmaceutical Chemistry and Technology, Faculty of Pharmacy, University of Helsinki, P.O. Box 56, FI-00014, Helsinki, Finland
| | - Hiba Abi Hussein
- Molécules thérapeutiques In silico (MTi), INSERM UMR-S973, University Paris Diderot, 35 rue Hélène Brion, 75013, Paris Cedex, France
| | - Anne-Claude Camproux
- Molécules thérapeutiques In silico (MTi), INSERM UMR-S973, University Paris Diderot, 35 rue Hélène Brion, 75013, Paris Cedex, France
| | - Leslie Regad
- Molécules thérapeutiques In silico (MTi), INSERM UMR-S973, University Paris Diderot, 35 rue Hélène Brion, 75013, Paris Cedex, France
| |
Collapse
|
22
|
Petitjean M. A Fast Algorithm to Compute Conical Pockets in Proteins. Application to the Structural Characterization of γ-Carbonic Anhydrases. Mol Inform 2017; 36. [PMID: 28402608 DOI: 10.1002/minf.201600155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Accepted: 03/30/2017] [Indexed: 11/12/2022]
Abstract
Some major proteins families, such as carbonic anhydrases (CAs), have a conical cavity at the active site. No algorithm was available to compute conical cavities, so we needed to design one. The fast algorithm we designed let us show on a set of 717 CAs extracted from the PDB database that γ-CAs are characterized by active site cavity cone angles significantly larger than those of α-CAs and β-CAs: the generatrix-axis angles are greater than 60° for the γ-CAs while they are smaller than 50° for the other CAs. Free binaries of the CONICA software implementing the algorithm are available through a software repository at http://petitjeanmichel.free.fr/itoweb.petitjean.freeware.html.
Collapse
Affiliation(s)
- Michel Petitjean
- Université Paris Diderot, MTi (INSERM UMR-S 973), 35 rue Hélène Brion, F-75013, Paris, France.,Université Paris Diderot, E-pôle de génoinformatique (IJM, CNRS UMR 7592), 15 rue Hélène Brion, F-75013 Paris, France
| |
Collapse
|
23
|
Miteva MA, Villoutreix BO. Computational Biology and Chemistry in MTi: Emphasis on the Prediction of Some ADMET Properties. Mol Inform 2017; 36. [DOI: 10.1002/minf.201700008] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2017] [Accepted: 02/03/2017] [Indexed: 12/21/2022]
Affiliation(s)
- Maria A. Miteva
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico , Inserm UMR−S 973; 35 rue Helene Brion 75013 Paris France
- INSERM, U973; F-75205 Paris France
| | - Bruno O. Villoutreix
- Université Paris Diderot, Sorbonne Paris Cité, Molécules Thérapeutiques In Silico , Inserm UMR−S 973; 35 rue Helene Brion 75013 Paris France
- INSERM, U973; F-75205 Paris France
| |
Collapse
|
24
|
Global vision of druggability issues: applications and perspectives. Drug Discov Today 2016; 22:404-415. [PMID: 27939283 DOI: 10.1016/j.drudis.2016.11.021] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2016] [Revised: 10/10/2016] [Accepted: 11/25/2016] [Indexed: 02/04/2023]
Abstract
During the preliminary stage of a drug discovery project, the lack of druggability information and poor target selection are the main causes of frequent failures. Elaborating on accurate computational druggability prediction methods is a requirement for prioritizing target selection, designing new drugs and avoiding side effects. In this review, we describe a survey of recently reported druggability prediction methods mainly based on networks, statistical pocket druggability predictions and virtual screening. An application for a frequent mutation of p53 tumor suppressor is presented, illustrating the complementarity of druggability prediction approaches, the remaining challenges and potential new drug development perspectives.
Collapse
|
25
|
Analytical algorithms for ligand cone angles calculations. Application to triphenylphosphine palladium complexes. CR CHIM 2015. [DOI: 10.1016/j.crci.2015.04.004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
26
|
Access channels to the buried active site control substrate specificity in CYP1A P450 enzymes. Biochim Biophys Acta Gen Subj 2015; 1850:696-707. [DOI: 10.1016/j.bbagen.2014.12.015] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2014] [Revised: 12/03/2014] [Accepted: 12/11/2014] [Indexed: 12/22/2022]
|
27
|
Juers DH, Ruffin J. MAP_CHANNELS: a computation tool to aid in the visualization and characterization of solvent channels in macromolecular crystals. J Appl Crystallogr 2014; 47:2105-2108. [PMID: 25484846 PMCID: PMC4248570 DOI: 10.1107/s160057671402281x] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2014] [Accepted: 10/17/2014] [Indexed: 11/10/2022] Open
Abstract
A computation tool is described that facilitates visualization and characterization of solvent channels or pores within macromolecular crystals. A scalar field mapping the shortest distance to protein surfaces is calculated on a grid covering the unit cell and is written as a map file. The map provides a multiscale representation of the solvent channels, which when viewed in standard macromolecular crystallographic software packages gives an intuitive sense of the solvent channel architecture. The map is analysed to yield descriptors of the topology and the morphology of the solvent channels, including bottleneck radii, tortuosity, width variation and anisotropy.
Collapse
Affiliation(s)
- Douglas H. Juers
- Department of Physics, Whitman College, Walla Walla, WA 99362, USA
| | - Jon Ruffin
- Department of Physics, Whitman College, Walla Walla, WA 99362, USA
| |
Collapse
|
28
|
Principal Component Analysis reveals correlation of cavities evolution and functional motions in proteins. J Mol Graph Model 2014; 55:13-24. [PMID: 25424655 DOI: 10.1016/j.jmgm.2014.10.011] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2014] [Revised: 10/16/2014] [Accepted: 10/18/2014] [Indexed: 11/24/2022]
Abstract
Protein conformation has been recognized as the key feature determining biological function, as it determines the position of the essential groups specifically interacting with substrates. Hence, the shape of the cavities or grooves at the protein surface appears to drive those functions. However, only a few studies describe the geometrical evolution of protein cavities during molecular dynamics simulations (MD), usually with a crude representation. To unveil the dynamics of cavity geometry evolution, we developed an approach combining cavity detection and Principal Component Analysis (PCA). This approach was applied to four systems subjected to MD (lysozyme, sperm whale myoglobin, Dengue envelope protein and EF-CaM complex). PCA on cavities allows us to perform efficient analysis and classification of the geometry diversity explored by a cavity. Additionally, it reveals correlations between the evolutions of the cavities and structures, and can even suggest how to modify the protein conformation to induce a given cavity geometry. It also helps to perform fast and consensual clustering of conformations according to cavity geometry. Finally, using this approach, we show that both carbon monoxide (CO) location and transfer among the different xenon sites of myoglobin are correlated with few cavity evolution modes of high amplitude. This correlation illustrates the link between ligand diffusion and the dynamic network of internal cavities.
Collapse
|