1
|
Wakabayashi T, Oide M, Nakasako M. CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase. Sci Rep 2024; 14:11165. [PMID: 38750092 PMCID: PMC11096400 DOI: 10.1038/s41598-024-61793-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 05/09/2024] [Indexed: 05/18/2024] Open
Abstract
Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.
Collapse
Grants
- JPMJPR22E2 Japan Science and Technology Agency
- jp13480214 Japan Society for the Promotion of Science
- jp19204042 Japan Society for the Promotion of Science
- jp22244054 Japan Society for the Promotion of Science
- jp21H01050 Japan Society for the Promotion of Science
- jp26800227 Japan Society for the Promotion of Science
- 18J11653 Japan Society for the Promotion of Science
- jp15076210 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp20050030 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp22018027 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp23120525, jp25120725 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp15H01647 Ministry of Education, Culture, Sports, Science and Technology of Japan
- jp17H05891 Ministry of Education, Culture, Sports, Science and Technology of Japan
Collapse
Affiliation(s)
- Taiki Wakabayashi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan
| | - Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan
- PRESTO, Japan Science and Technology Agency, Chiyoda-Ku, Tokyo, 102-0076, Japan
- Protein Research Institute, Osaka University, Yamadaoka, Suita, Osaka, 565-0871, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-Ku, Yokohama, Kanagawa, 223-8522, Japan.
- RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-Cho, Sayo-Gun, Hyogo, 679-5148, Japan.
| |
Collapse
|
2
|
Zhang L, Dai W, Rong S, Schwaneberg U, Xu G, Ni Y. Engineering diaryl alcohol dehydrogenase KpADH reveals importance of retaining hydration shell in organic solvent tolerance. Protein Sci 2024; 33:e4933. [PMID: 38501647 PMCID: PMC10949390 DOI: 10.1002/pro.4933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 01/27/2024] [Accepted: 02/01/2024] [Indexed: 03/20/2024]
Abstract
Alcohol dehydrogenases (ADHs) are synthetically important biocatalysts for the asymmetric synthesis of chiral alcohols. The catalytic performance of ADHs in the presence of organic solvents is often important since most prochiral ketones are highly hydrophobic. Here, the organic solvent tolerance of KpADH from Kluyveromyces polyspora was semi-rationally evolved. Using tolerant variants obtained, meticulous experiments and computational studies were conducted to explore properties including stability, activity and kinetics in the presence of various organic solvents. Compared with WT, variant V231D exhibited 1.9-fold improvement in ethanol tolerance, while S237G showed a 6-fold increase in catalytic efficiency, a higherT 50 15 $$ {\mathrm{T}}_{50}^{15} $$ , as well as 15% higher tolerance in 7.5% (v/v) ethanol. Based on 3 × 100 ns MD simulations, the increased tolerance of V231D and S237G against ethanol may be ascribed to their enhanced ability in retaining water molecules and repelling ethanol molecules. Moreover, 6.3-fold decreased KM value of V231D toward hydrophilic ketone substrate confirmed its capability of retaining hydration shell. Our results suggest that retaining hydration shell surrounding KpADH is critical for its tolerance to organic solvents, as well as catalytic performance. This study provides useful guidance for engineering organic solvent tolerance of KpADH and other ADHs.
Collapse
Affiliation(s)
- Lu Zhang
- Key Laboratory of Industrial Biotechnology, Ministry of EducationSchool of Biotechnology, Jiangnan UniversityWuxiChina
| | - Wei Dai
- Key Laboratory of Industrial Biotechnology, Ministry of EducationSchool of Biotechnology, Jiangnan UniversityWuxiChina
| | - Shuo Rong
- Key Laboratory of Industrial Biotechnology, Ministry of EducationSchool of Biotechnology, Jiangnan UniversityWuxiChina
| | | | - Guochao Xu
- Key Laboratory of Industrial Biotechnology, Ministry of EducationSchool of Biotechnology, Jiangnan UniversityWuxiChina
| | - Ye Ni
- Key Laboratory of Industrial Biotechnology, Ministry of EducationSchool of Biotechnology, Jiangnan UniversityWuxiChina
| |
Collapse
|
3
|
Wakabayashi T, Oide M, Kato T, Nakasako M. Coenzyme-binding pathway on glutamate dehydrogenase suggested from multiple-binding sites visualized by cryo-electron microscopy. FEBS J 2023; 290:5514-5535. [PMID: 37682540 DOI: 10.1111/febs.16951] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2022] [Revised: 08/10/2023] [Accepted: 09/05/2023] [Indexed: 09/09/2023]
Abstract
The structure of hexameric glutamate dehydrogenase (GDH) in the presence of the coenzyme nicotinamide adenine dinucleotide phosphate (NADP) was visualized using cryogenic transmission electron microscopy to investigate the ligand-binding pathways to the active site of the enzyme. Each subunit of GDH comprises one hexamer-forming core domain and one nucleotide-binding domain (NAD domain), which spontaneously opens and closes the active-site cleft situated between the two domains. In the presence of NADP, the potential map of GDH hexamer, assuming D3 symmetry, was determined at a resolution of 2.4 Å, but the NAD domain was blurred due to the conformational variety. After focused classification with respect to the NAD domain, the potential maps interpreted as NADP molecules appeared at five different sites in the active-site cleft. The subunits associated with NADP molecules were close to one of the four metastable conformations in the unliganded state. Three of the five binding sites suggested a pathway of NADP molecules to approach the active-site cleft for initiating the enzymatic reaction. The other two binding modes may rarely appear in the presence of glutamate, as demonstrated by the reaction kinetics. Based on the visualized structures and the results from the enzymatic kinetics, we discussed the binding modes of NADP to GDH in the absence and presence of glutamate.
Collapse
Grants
- JPMJPR22E2 Japan Science and Technology Agency
- 18J11653 Japan Society for the Promotion of Science
- jp13480214 Japan Society for the Promotion of Science
- jp19204042 Japan Society for the Promotion of Science
- jp21H01050 Japan Society for the Promotion of Science
- jp22244054 Japan Society for the Promotion of Science
- jp26800227 Japan Society for the Promotion of Science
- jp15076210 Ministry of Education, Culture, Sports, Science and Technology
- jp15H01647 Ministry of Education, Culture, Sports, Science and Technology
- jp17H05891 Ministry of Education, Culture, Sports, Science and Technology
- jp20050030 Ministry of Education, Culture, Sports, Science and Technology
- jp22018027 Ministry of Education, Culture, Sports, Science and Technology
- jp23120525 Ministry of Education, Culture, Sports, Science and Technology
- jp25120725 Ministry of Education, Culture, Sports, Science and Technology
- 0436 Japan Agency for Medical Research and Development
Collapse
Affiliation(s)
- Taiki Wakabayashi
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan
- RIKEN SPring-8 Center, Sayo-gun, Hyogo, Japan
- RIKEN Cluster for Pioneering Research, Wako, Japan
| | - Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan
- RIKEN SPring-8 Center, Sayo-gun, Hyogo, Japan
- RIKEN Cluster for Pioneering Research, Wako, Japan
- PRESTO, Japan Science and Technology Agency, Tokyo, Japan
| | - Takayuki Kato
- Protein Research Institute, Osaka University, Suita, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan
- RIKEN SPring-8 Center, Sayo-gun, Hyogo, Japan
| |
Collapse
|
4
|
Sato K, Oide M, Nakasako M. Prediction of hydrophilic and hydrophobic hydration structure of protein by neural network optimized using experimental data. Sci Rep 2023; 13:2183. [PMID: 36750742 PMCID: PMC9905073 DOI: 10.1038/s41598-023-29442-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2022] [Accepted: 02/06/2023] [Indexed: 02/09/2023] Open
Abstract
The hydration structures of proteins, which are necessary for their folding, stability, and functions, were visualized using X-ray and neutron crystallography and transmission electron microscopy. However, complete visualization of hydration structures over the entire protein surface remains difficult. To compensate for this incompleteness, we developed a three-dimensional convolutional neural network to predict the probability distribution of hydration water molecules on the hydrophilic and hydrophobic surfaces, and in the cavities of proteins. The neural network was optimized using the distribution patterns of protein atoms around the hydration water molecules identified in the high-resolution X-ray crystal structures. We examined the feasibility of the neural network using water sites in the protein crystal structures that were not included in the datasets. The predicted distribution covered most of the experimentally identified hydration sites, with local maxima appearing in their vicinity. This computational approach will help to highlight the relevance of hydration structures to the biological functions and dynamics of proteins.
Collapse
Affiliation(s)
- Kochi Sato
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, Kanagawa, 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan
| | - Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, Kanagawa, 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan.,PRESTO, Japan Science and Technology Agency, Chiyoda-ku, Tokyo, 102-0076, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, Kanagawa, 223-8522, Japan. .,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo, 679-5148, Japan.
| |
Collapse
|
5
|
Sultan A, Ali R, Ishrat R, Ali S. Anti-HIV and anti-HCV small molecule protease inhibitors in-silico repurposing against SARS-CoV-2 M pro for the treatment of COVID-19. J Biomol Struct Dyn 2022; 40:12848-12862. [PMID: 34569411 DOI: 10.1080/07391102.2021.1979097] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
The COVID-19 pandemic is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). It is a global health emergency warranting development and implementation of targeted treatment. The enzyme main protease (Mpro; also known as 3C-like protease) is emerging as an attractive drug target. This enzyme plays an indispensable role in processing the translated polyproteins of viral RNA. Inhibiting the activity of Mpro would wedge viral replication. To facilitate the discovery of targeted therapy for COVID-19, we carried out the structure-assisted repurposing of existing protease inhibiting small molecules to target SARS-CoV-2 Mpro. Based on the structure of SARS-CoV-2 Mpro, here we report the small drug molecule namely saquinavir as its potent inhibitor. Findings support the premise that this promising antiviral protease inhibiting small drug molecule can be validated and implemented for the treatment and clinical management of COVID-19 pandemic disease.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Armiya Sultan
- Department of Biosciences, Jamia Millia Islamia (A Central University), New Delhi, India
| | - Rafat Ali
- Computational Laboratory, Center for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia (A Central University), New Delhi, India
| | - Romana Ishrat
- Computational Laboratory, Center for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia (A Central University), New Delhi, India
| | - Sher Ali
- Department of Life Sciences, Sharda University, Greater Noida, UP, India
| |
Collapse
|
6
|
Chang HY, Wu KY, Chen WC, Weng JT, Chen CY, Raj A, Hamaguchi HO, Chuang WT, Wang X, Wang CL. Water-Induced Self-Assembly of Amphiphilic Discotic Molecules for Adaptive Artificial Water Channels. ACS NANO 2021; 15:14885-14890. [PMID: 34410689 DOI: 10.1021/acsnano.1c04994] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Inspired by the induced-fit mechanism in nature, we developed the process of water-induced self-assembly (WISA) to make water an active substrate that regulates the self-assembly and function of amphiphilic discotic molecules (ADMs). The ADM is an isotropic liquid that self-assembles only when in contact with water. Characterization results indicate that water fits into the hydrophilic core of the ADMs and induces the formation of a hexagonal columnar phase (Colh), where each column contains a hydrated artificial water channel (AWC). The hydrated AWCs are adaptive rather than static; the dynamic incorporation/removal of water results in the reversible assembly/disassembly of the adaptive AWCs (aAWCs). Furthermore, its dynamic characteristics can enable water to act as an orientation-directional guest molecule that controls the growth direction of the aAWCs. Well-aligned aAWC arrays that showed the ability of water transport were obtained via a "directional WISA" method. In WISA, water thus governs the supramolecular chemistry and function of synthetic molecules as it does with natural materials. By making water an active component in adaptive chemistry and enabling host molecules to dynamically interact with water, this adaptive aquatic material may motivate the development of synthetic molecules further toward biomaterials.
Collapse
Affiliation(s)
- Hsi-Yen Chang
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Kuan-Yi Wu
- National Synchrotron Radiation Research Center, 101 Hsin-Ann Road, Hsinchu, 30076, Taiwan
| | - Wei-Chun Chen
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Jing-Ting Weng
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Chin-Yi Chen
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Ankit Raj
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Hiro-O Hamaguchi
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| | - Wei-Tsung Chuang
- National Synchrotron Radiation Research Center, 101 Hsin-Ann Road, Hsinchu, 30076, Taiwan
| | - Xiaosong Wang
- Department of Chemistry, Waterloo Institute for Nanotechnology, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Chien-Lung Wang
- Department of Applied Chemistry, National Yang Ming Chiao Tung University, 1001 Ta Hsueh Road, Hsinchu, 30010, Taiwan
| |
Collapse
|
7
|
Valadares VS, Martins LC, Roman EA, Valente AP, Cino EA, Moraes AH. Conformational dynamics of Tetracenomycin aromatase/cyclase regulate polyketide binding and enzyme aggregation propensity. Biochim Biophys Acta Gen Subj 2021; 1865:129949. [PMID: 34139289 DOI: 10.1016/j.bbagen.2021.129949] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2021] [Revised: 06/01/2021] [Accepted: 06/10/2021] [Indexed: 11/25/2022]
Abstract
BACKGROUND The N-terminal domain of Tetracenomycin aromatase/cyclase (TcmN), an enzyme derived from Streptomyces glaucescens, is involved in polyketide cyclization, aromatization, and folding. Polyketides are a diverse class of secondary metabolites produced by certain groups of bacteria, fungi, and plants with various pharmaceutical applications. Examples include antibiotics, such as tetracycline, and anticancer drugs, such as doxorubicin. Because TcmN is a promising enzyme for in vitro production of polyketides, it is important to identify conditions that enhance its thermal resistance and optimize its function. METHODS TcmN unfolding, stability, and dynamics were evaluated by fluorescence spectroscopy, circular dichroism, nuclear magnetic resonance 15N relaxation experiments, and microsecond molecular dynamics (MD) simulations. RESULTS TcmN thermal resistance was enhanced at low protein and high salt concentrations, was pH-dependent, and denaturation was irreversible. Conformational dynamics on the μs-ms timescale were detected for residues in the substrate-binding cavity, and two predominant conformers representing opened and closed cavity states were observed in the MD simulations. CONCLUSION Based on the results, a mechanism was proposed in which the thermodynamics and kinetics of the TcmN conformational equilibrium modulate enzyme function by favoring ligand binding and avoiding aggregation. GENERAL SIGNIFICANCE Understanding the principles underlying TcmN stability and dynamics may help in designing mutants with optimal properties for biotechnological applications.
Collapse
Affiliation(s)
- Veronica S Valadares
- Department of Biochemistry and Immunology, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
| | - Luan C Martins
- Graduate Program in Bioinformatics, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
| | - Ernesto A Roman
- Facultad de Farmacia y Bioquimica, Universidad de Buenos Aires, Buenos Aires, Argentina
| | - Ana Paula Valente
- National Center of Nuclear Magnetic Resonance, Leopoldo de Meis Institute of Medical Biochemistry, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Elio A Cino
- Department of Biochemistry and Immunology, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil.
| | - Adolfo H Moraes
- Departamento de Química, Instituto de Ciências Exatas, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil; Department of NMR-based Structural Biology, Max-Planck Institute for Biophysical Chemistry, Göttingen, Germany.
| |
Collapse
|
8
|
Mapping the Intramolecular Communications among Different Glutamate Dehydrogenase States Using Molecular Dynamics. Biomolecules 2021; 11:biom11060798. [PMID: 34072154 PMCID: PMC8228935 DOI: 10.3390/biom11060798] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2021] [Revised: 05/18/2021] [Accepted: 05/25/2021] [Indexed: 12/23/2022] Open
Abstract
Glutamate dehydrogenase (GDH) is a ubiquitous enzyme that catalyzes the reversible oxidative deamination of glutamate to α-ketoglutarate. It acts as an important branch-point enzyme between carbon and nitrogen metabolisms. Due to the multifaceted roles of GDH in cancer, hyperinsulinism/hyperammonemia, and central nervous system development and pathologies, tight control of its activity is necessitated. To date, several GDH structures have been solved in its closed form; however, intrinsic structural information in its open and apo forms are still deficient. Moreover, the allosteric communications and conformational changes taking place in the three different GDH states are not well studied. To mitigate these drawbacks, we applied unbiased molecular dynamic simulations (MD) and network analysis to three different GDH states i.e., apo, active, and inactive forms, for investigating their modulatory mechanisms. In this paper, based on MD and network analysis, crucial residues important for signal transduction, conformational changes, and maps of information flow among the different GDH states were elucidated. Moreover, with the recent findings of allosteric modulators, an allosteric wiring illustration of GDH intramolecular signal transductions would be of paramount importance to obtain the process of this enzyme regulation. The structural insights gained from this study will pave way for large-scale screening of GDH regulators and could support researchers in the design and development of new and potent GDH ligands.
Collapse
|
9
|
Ren W, Dokainish HM, Shinobu A, Oshima H, Sugita Y. Unraveling the Coupling between Conformational Changes and Ligand Binding in Ribose Binding Protein Using Multiscale Molecular Dynamics and Free-Energy Calculations. J Phys Chem B 2021; 125:2898-2909. [PMID: 33728914 PMCID: PMC10954230 DOI: 10.1021/acs.jpcb.0c11600] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Conformational changes of proteins upon ligand binding are usually explained in terms of several mechanisms including the induced fit, conformational selection, or their mixtures. Due to the slow time scales, conventional molecular dynamics (cMD) simulations based on the atomistic models cannot easily simulate the open-to-closed conformational transition in proteins. In our previous study, we have developed an enhanced sampling scheme (generalized replica exchange with solute tempering selected surface charged residues: gREST_SSCR) for multidomain proteins and applied it to ligand-mediated conformational changes in the G134R mutant of ribose-binding protein (RBPG134R) in solution. The free-energy landscape (FEL) of RBPG134R in the presence of a ribose at the binding site included the open and closed states and two intermediates, open-like and closed-like forms. Only the open and open-like forms existed in the FEL without a ribose. In the current study, the coupling between the conformational changes and ligand binding is further investigated using coarse-grained MD, multiple atomistic cMD, and free-energy calculations. The ribose is easily dissociated from the binding site of wild-type RBP and RBPG134R in the cMD simulations starting from the open and open-like forms. In contrast, it is stable at the binding site in the simulations from the closed and closed-like forms. The free-energy calculations provide the binding affinities of different structures, supporting the results of cMD simulations. Importantly, cMD simulations from the closed-like structures reveal transitions toward the closed one in the presence of a bound ribose. On the basis of the computational results, we propose a molecular mechanism in which conformational selection and induced fit happen in the first and second halves of the open-to-closed transition in RBP, respectively.
Collapse
Affiliation(s)
- Weitong Ren
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Hisham M. Dokainish
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
| | - Ai Shinobu
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Hiraku Oshima
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| | - Yuji Sugita
- Theoretical
Molecular Science Laboratory, RIKEN Cluster
for Pioneering Research, 2-1 Hirosawa, Wako, Saitama 351-0198, Japan
- Computational
Biophysics Research Team, RIKEN Center for
Computational Science, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
- Laboratory
for Biomolecular Function Simulation, RIKEN
Center for Biosystems Dynamics Research, Integrated Innovation Building 7F, 6-7-1 minatojima-minamimachi,
Chuo-ku, Kobe, Hyogo 650-0047, Japan
| |
Collapse
|
10
|
Cui H, Zhang L, Eltoukhy L, Jiang Q, Korkunç SK, Jaeger KE, Schwaneberg U, Davari MD. Enzyme Hydration Determines Resistance in Organic Cosolvents. ACS Catal 2020. [DOI: 10.1021/acscatal.0c03233] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Haiyang Cui
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| | - Lingling Zhang
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| | - Lobna Eltoukhy
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| | - Qianjia Jiang
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| | - Seval Kübra Korkunç
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| | - Karl-Erich Jaeger
- Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Wilhelm Johnen Strasse, Jülich 52426, Germany
- Institute of Bio-and Geosciences IBG 1: Biotechnology, Forschungszentrum Jülich GmbH, Wilhelm Johnen Strasse, Jülich 52426, Germany
| | - Ulrich Schwaneberg
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
- DWI-Leibniz Institute for Interactive Materials, Forckenbeckstrasse 50, Aachen 52074, Germany
| | - Mehdi D. Davari
- Institute of Biotechnology, RWTH Aachen University, Worringerweg 3, Aachen 52074, Germany
| |
Collapse
|
11
|
Oroguchi T, Oide M, Wakabayashi T, Nakasako M. Assessment of Force Field Accuracy Using Cryogenic Electron Microscopy Data of Hyper-thermostable Glutamate Dehydrogenase. J Phys Chem B 2020; 124:8479-8494. [PMID: 32841031 DOI: 10.1021/acs.jpcb.0c04464] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Molecular dynamics (MD) simulations in biophysically relevant time scales of microseconds is a powerful tool for studying biomolecular processes, but results often display force field dependency. Therefore, assessment of force field accuracy using experimental data of biomolecules in solution is essential for simulation studies. Here, we propose the use of structural models obtained via cryo-electron microscopy (cryoEM), which provides biomolecular structures in vitreous ice mimicking the environment in solution. The accuracy of the AMBER (ff99SB-ILDN-NMR, ff14SB, ff15ipq, and ff15FB) and CHARMM (CHARMM22 and CHARMM36m) force fields was assessed by comparing their MD trajectories with the cryoEM data of thermostable hexameric glutamate dehydrogenase (GDH), which included a cryoEM map at a resolution of approximately 3 Å and structure models of subunits reflecting metastable conformations in domain motion occurring in GDH. In the assessment, we validated the force fields with respect to the reproducibility and stability of secondary structures and intersubunit interactions in the cryoEM data. Furthermore, we evaluated the force fields regarding the reproducibility of the energy landscape in the domain motion expected from the cryoEM data. As a result, among the six force fields, ff15FB and ff99SB-ILDN-NMR displayed good agreement with the experiment. The present study demonstrated the advantages of the high-resolution cryoEM map and suggested the optimal force field to reproduce experimentally observed protein structures.
Collapse
Affiliation(s)
- Tomotaka Oroguchi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-ku, Yokohama, Kanagawa 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan
| | - Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-ku, Yokohama, Kanagawa 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan
| | - Taiki Wakabayashi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-ku, Yokohama, Kanagawa 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoko-ku, Yokohama, Kanagawa 223-8522, Japan.,RIKEN SPring-8 Center, 1-1-1 Kouto, Sayo-cho, Sayo-gun, Hyogo 679-5148, Japan
| |
Collapse
|
12
|
Setny P. Conserved internal hydration motifs in protein kinases. Proteins 2020; 88:1578-1591. [DOI: 10.1002/prot.25977] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2019] [Revised: 05/18/2020] [Accepted: 07/05/2020] [Indexed: 01/14/2023]
Affiliation(s)
- Piotr Setny
- Centre of New Technologies University of Warsaw Warsaw Poland
| |
Collapse
|
13
|
Cui H, Stadtmüller THJ, Jiang Q, Jaeger K, Schwaneberg U, Davari MD. How to Engineer Organic Solvent Resistant Enzymes: Insights from Combined Molecular Dynamics and Directed Evolution Study. ChemCatChem 2020. [DOI: 10.1002/cctc.202000422] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Haiyang Cui
- Lehrstuhl für Biotechnologie RWTH Aachen University Worringerweg 3 52074 Aachen Germany
| | - Tom H. J. Stadtmüller
- Lehrstuhl für Biotechnologie RWTH Aachen University Worringerweg 3 52074 Aachen Germany
| | - Qianjia Jiang
- Lehrstuhl für Biotechnologie RWTH Aachen University Worringerweg 3 52074 Aachen Germany
| | - Karl‐Erich Jaeger
- Institute of Molecular Enzyme Technology Heinrich Heine University Düsseldorf and Research Center Jülich Wilhelm Johnen Strasse 52426 Jülich Germany
| | - Ulrich Schwaneberg
- Lehrstuhl für Biotechnologie RWTH Aachen University Worringerweg 3 52074 Aachen Germany
- DWI-Leibniz Institute for Interactive Materials Forckenbeckstraße 50 52074 Aachen Germany
| | - Mehdi D. Davari
- Lehrstuhl für Biotechnologie RWTH Aachen University Worringerweg 3 52074 Aachen Germany
| |
Collapse
|
14
|
Hydrogen molecules can modulate enzymatic activity and structural properties of pepsin in vitro. Colloids Surf B Biointerfaces 2020; 189:110856. [DOI: 10.1016/j.colsurfb.2020.110856] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Revised: 02/03/2020] [Accepted: 02/07/2020] [Indexed: 01/19/2023]
|
15
|
Oide M, Kato T, Oroguchi T, Nakasako M. Energy landscape of domain motion in glutamate dehydrogenase deduced from cryo-electron microscopy. FEBS J 2020; 287:3472-3493. [PMID: 31976609 DOI: 10.1111/febs.15224] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Revised: 11/09/2019] [Accepted: 01/21/2020] [Indexed: 11/28/2022]
Abstract
Analysis of the conformational changes of protein is important to elucidate the mechanisms of protein motions correlating with their function. Here, we studied the spontaneous domain motion of unliganded glutamate dehydrogenase from Thermococcus profundus using cryo-electron microscopy and proposed a novel method to construct free-energy landscape of protein conformations. Each subunit of the homo-hexameric enzyme comprises nucleotide-binding domain (NAD domain) and hexamer-forming core domain. A large active-site cleft is situated between the two domains and varies from open to close according to the motion of a NAD domain. A three-dimensional map reconstructed from all cryo-electron microscopy images displayed disordered volumes of NAD domains, suggesting that NAD domains in the collected images adopted various conformations in domain motion. Focused classifications on NAD domain of subunits provided several maps of possible conformations in domain motion. To deduce what kinds of conformations appeared in EM images, we developed a novel analysis method that describe the EM maps as a linear combination of representative conformations appearing in a 200-ns molecular dynamics simulation as reference. The analysis enabled us to estimate the appearance frequencies of the representative conformations, which illustrated a free-energy landscape in domain motion. In the open/close domain motion, two free-energy basins hindered the direct transformation from open to closed state. Structure models constructed for representative EM maps in classifications demonstrated the correlation between the energy landscape and conformations in domain motion. Based on the results, the domain motion in glutamate dehydrogenase and the analysis method to visualize conformational changes and free-energy landscape were discussed. DATABASE: The EM maps of the four conformations were deposited to Electron Microscopy Data Bank (EMDB) as accession codes EMD-9845 (open), EMD-9846 (half-open1), EMD-9847 (half-open2), and EMD-9848 (closed), respectively. In addition, the structural models built for the four conformations were deposited to the Protein Data Bank (PDB) as accession codes 6JN9 (open), 6JNA (half-open1), 6JNC (half-open2), and 6JND (closed), respectively.
Collapse
Affiliation(s)
- Mao Oide
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| | - Takayuki Kato
- Graduate School of Frontier Biosciences, Osaka University, Suita, Japan
| | - Tomotaka Oroguchi
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, Yokohama, Japan.,RIKEN SPring-8 Center, Sayo-gun, Japan
| |
Collapse
|
16
|
He P, Sarkar S, Gallicchio E, Kurtzman T, Wickstrom L. Role of Displacing Confined Solvent in the Conformational Equilibrium of β-Cyclodextrin. J Phys Chem B 2019; 123:8378-8386. [PMID: 31509409 DOI: 10.1021/acs.jpcb.9b07028] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
This study investigates the role of hydration and its relationship to the conformational equilibrium of the host molecule β-cyclodextrin. Molecular dynamics simulations indicate that the unbound β-cyclodextrin exhibits two state behavior in explicit solvent due to the opening and closing of its cavity. In implicit solvent, these transitions are not observed, and there is one dominant conformation of β-cyclodextrin with an open cavity. Based on these observations, we investigate the hypothesis that the expulsion of thermodynamically unfavorable water molecules into the bulk plays an important role in controlling the accessibility of the closed macrostate at room temperature. We compare the results of the molecular mechanics analytical generalized Born plus nonpolar solvation approach to those obtained through grid inhomogeneous solvation theory analysis with explicit solvation to elucidate the thermodynamic forces at play. The work illustrates the use of continuum solvent models to tease out solvation effects related to the inhomogeneity and the molecular nature of water and demonstrates the key role of the thermodynamics of enclosed hydration in driving the conformational equilibrium of molecules in solution.
Collapse
Affiliation(s)
- Peng He
- Center for Biophysics & Computational Biology/ICMS, Department of Chemistry , Temple University , Philadelphia , Pennsylvania 19122 , United States
| | - Sheila Sarkar
- Department of Science , Borough of Manhattan Community College, The City University of New York , New York , New York 10007 , United States
| | - Emilio Gallicchio
- Department of Chemistry , Brooklyn College, The City University of New York , Brooklyn , New York 11210 , United States.,Ph.D. Programs in Chemistry & Biochemistry , The Graduate Center of the City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Tom Kurtzman
- Department of Chemistry , Lehman College, The City University of New York , Bronx , New York 10468 , United States.,Ph.D. Programs in Chemistry & Biochemistry , The Graduate Center of the City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Lauren Wickstrom
- Department of Science , Borough of Manhattan Community College, The City University of New York , New York , New York 10007 , United States
| |
Collapse
|
17
|
Pal RK, Gadhiya S, Ramsey S, Cordone P, Wickstrom L, Harding WW, Kurtzman T, Gallicchio E. Inclusion of enclosed hydration effects in the binding free energy estimation of dopamine D3 receptor complexes. PLoS One 2019; 14:e0222902. [PMID: 31568493 PMCID: PMC6768453 DOI: 10.1371/journal.pone.0222902] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Accepted: 08/30/2019] [Indexed: 01/04/2023] Open
Abstract
Confined hydration and conformational flexibility are some of the challenges encountered for the rational design of selective antagonists of G-protein coupled receptors. We present a set of C3-substituted (-)-stepholidine derivatives as potent binders of the dopamine D3 receptor. The compounds are characterized biochemically, as well as by computer modeling using a novel molecular dynamics-based alchemical binding free energy approach which incorporates the effect of the displacement of enclosed water molecules from the binding site. The free energy of displacement of specific hydration sites is obtained using the Hydration Site Analysis method with explicit solvation. This work underscores the critical role of confined hydration and conformational reorganization in the molecular recognition mechanism of dopamine receptors and illustrates the potential of binding free energy models to represent these key phenomena.
Collapse
Affiliation(s)
- Rajat Kumar Pal
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, NY 11210, United States of America
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
| | - Satishkumar Gadhiya
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Steven Ramsey
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Lehman College, 250 Bedford Park Blvd. West, Bronx, NY 10468, United States of America
| | - Pierpaolo Cordone
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Lauren Wickstrom
- Department of Science, Borough of Manhattan Community College, 199 Chambers Street, New York, NY 10007, United States of America
| | - Wayne W. Harding
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Hunter College, 695 Park Avenue, NY 10065, United States of America
| | - Tom Kurtzman
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- Department of Chemistry, Lehman College, 250 Bedford Park Blvd. West, Bronx, NY 10468, United States of America
| | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College, 2900 Bedford Avenue, Brooklyn, NY 11210, United States of America
- PhD Program in Biochemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- PhD Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, United States of America
- * E-mail:
| |
Collapse
|
18
|
Fukuda A, Oroguchi T, Nakasako M. Dipole-dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme. Biochim Biophys Acta Gen Subj 2019; 1864:129406. [PMID: 31377191 DOI: 10.1016/j.bbagen.2019.07.015] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/28/2018] [Revised: 07/24/2019] [Accepted: 07/29/2019] [Indexed: 10/26/2022]
Abstract
The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320-370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole-dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.
Collapse
Affiliation(s)
- Asahi Fukuda
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokihama 223-8522, Japan
| | - Tomotaka Oroguchi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokihama 223-8522, Japan
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokihama 223-8522, Japan.
| |
Collapse
|
19
|
Cametti M, Valkonen A, Rissanen K. Entrapment of a linear water pentamer into a uranyl-salophen dimer in the solid state. Supramol Chem 2019. [DOI: 10.1080/10610278.2019.1632458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Massimo Cametti
- Department of Chemistry, Materials and Chemical Engineering"Giulio Natta", Politecnico di Milano, Milano, Italy
| | - Arto Valkonen
- Department of Chemistry, University of Jyväskylä, Jyväskylä, Finland
| | - Kari Rissanen
- Department of Chemistry, University of Jyväskylä, Jyväskylä, Finland
| |
Collapse
|
20
|
Dahanayake JN, Shahryari E, Roberts KM, Heikes ME, Kasireddy C, Mitchell-Koch KR. Protein Solvent Shell Structure Provides Rapid Analysis of Hydration Dynamics. J Chem Inf Model 2019; 59:2407-2422. [PMID: 30865440 DOI: 10.1021/acs.jcim.9b00009] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The solvation layer surrounding a protein is clearly an intrinsic part of protein structure-dynamics-function, and our understanding of how the hydration dynamics influences protein function is emerging. We have recently reported simulations indicating a correlation between regional hydration dynamics and the structure of the solvation layer around different regions of the enzyme Candida antarctica lipase B, wherein the radial distribution function (RDF) was used to calculate the pairwise entropy, providing a link between dynamics (diffusion) and thermodynamics (excess entropy) known as Rosenfeld scaling. Regions with higher RDF values/peaks in the hydration layer (the first peak, within 6 Å of the protein surface) have faster diffusion in the hydration layer. The finding thus hinted at a handle for rapid evaluation of hydration dynamics at different regions on the protein surface in molecular dynamics simulations. Such an approach may move the analysis of hydration dynamics from a specialized venture to routine analysis, enabling an informatics approach to evaluate the role of hydration dynamics in biomolecular function. This paper first confirms that the correlation between regional diffusive dynamics and hydration layer structure (via water center of mass around protein side-chain atom RDF) is observed as a general relationship across a set of proteins. Second, it seeks to devise an approach for rapid analysis of hydration dynamics, determining the minimum amount of information and computational effort required to get a reliable value of hydration dynamics from structural data in MD simulations based on the protein-water RDF. A linear regression model using the integral of the hydration layer in the water-protein RDF was found to provide statistically equivalent apparent diffusion coefficients at the 95% confidence level for a set of 92 regions within five different proteins. In summary, RDF analysis of 10 ns of data after simulation convergence is sufficient to accurately map regions of fast and slow hydration dynamics around a protein surface. Additionally, it is anticipated that a quick look at protein-water RDFs, comparing peak heights, will be useful to provide a qualitative ranking of regions of faster and slower hydration dynamics at the protein surface for rapid analysis when investigating the role of solvent dynamics in protein function.
Collapse
Affiliation(s)
- Jayangika N Dahanayake
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Elaheh Shahryari
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Kirsten M Roberts
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Micah E Heikes
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Chandana Kasireddy
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| | - Katie R Mitchell-Koch
- Department of Chemistry , Wichita State University , 1845 Fairmount Street , Wichita , Kansas 67260-0051 , United States
| |
Collapse
|
21
|
Kitazawa S, Aoshima Y, Wakamoto T, Kitahara R. Water-Protein Interactions Coupled with Protein Conformational Transition. Biophys J 2018; 115:981-987. [PMID: 30146267 PMCID: PMC6139601 DOI: 10.1016/j.bpj.2018.08.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2018] [Revised: 07/16/2018] [Accepted: 08/03/2018] [Indexed: 10/28/2022] Open
Abstract
Conformational fluctuations of proteins are crucially important for their functions. However, changes in the location and dynamics of hydrated water in many proteins accompanied by the conformational transition have not been fully understood. Here, we used phase-modulated clean chemical exchange NMR approach to investigate pressure-induced changes in water-to-amide proton exchange occurring at sub-second time scale. With the transition of ubiquitin from its native conformation (N1) to an alternative conformation (N2) at 250 MPa, proton exchange rates of residues 32-35, 40-41, and 71, which are located at the C-terminal side of the protein, were significantly increased. These observations can be explained by the destabilization of the hydrogen bonds in the backbone and partial exposure of those amide groups to solvent in N2. We conclude that phase-modulated clean chemical exchange NMR approach coupled with pressure perturbation will be a useful tool for investigations of more open and hydrated protein structures.
Collapse
Affiliation(s)
| | - Yu Aoshima
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | - Takuro Wakamoto
- Graduate School of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | | |
Collapse
|
22
|
Capdevila DA, Edmonds KA, Campanello GC, Wu H, Gonzalez-Gutierrez G, Giedroc DP. Functional Role of Solvent Entropy and Conformational Entropy of Metal Binding in a Dynamically Driven Allosteric System. J Am Chem Soc 2018; 140:9108-9119. [PMID: 29953213 PMCID: PMC6425489 DOI: 10.1021/jacs.8b02129] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Allostery is a regulatory phenomenon whereby ligand binding to one site influences the binding of the same or a different ligand to another site on a macromolecule. The physical origins of allosteric regulation remain under intense investigation. In general terms, ligand-induced structural changes, perturbations of residue-specific dynamics, and surrounding solvent molecules all potentially contribute to the global energetics of allostery. While the role of solvent is generally well understood in regulatory events associated with major protein structural rearrangements, the degree to which protein dynamics impact solvent degrees of freedom is unclear, particularly in cases of dynamically driven allostery. With the aid of new crystal structures, extensive calorimetric and residue-specific dynamics studies over a range of time scales and temperatures, we dissect for the first time the relative degree to which changes in solvent entropy and residue-specific dynamics impact dynamically driven, allosteric inhibition of DNA binding by Zn in the zinc efflux repressor, CzrA (chromosomal zinc-regulated repressor). We show that non-native residue-specific dynamics in allosterically impaired CzrA mutants are accompanied by significant perturbations in solvent entropy that cannot be predicted from crystal structures. We conclude that functional dynamics are not necessarily restricted to protein residues but involve surface water molecules that may be responding to ligand (Zn)-mediated perturbations in protein internal motions that define the conformational ensemble, rather than major structural rearrangements.
Collapse
Affiliation(s)
- Daiana A. Capdevila
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102 United States
| | - Katherine A. Edmonds
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102 United States
| | - Gregory C. Campanello
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102 United States
| | - Hongwei Wu
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102 United States
| | - Giovanni Gonzalez-Gutierrez
- Department of Molecular and Cellular Biochemistry, Indiana University, Bloomington, Indiana 47405 United States
| | - David P. Giedroc
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405-7102 United States
- Department of Molecular and Cellular Biochemistry, Indiana University, Bloomington, Indiana 47405 United States
| |
Collapse
|
23
|
Dahanayake JN, Mitchell-Koch KR. How Does Solvation Layer Mobility Affect Protein Structural Dynamics? Front Mol Biosci 2018; 5:65. [PMID: 30057902 PMCID: PMC6053501 DOI: 10.3389/fmolb.2018.00065] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2018] [Accepted: 06/20/2018] [Indexed: 11/18/2022] Open
Abstract
Solvation is critical for protein structural dynamics. Spectroscopic studies have indicated relationships between protein and solvent dynamics, and rates of gas binding to heme proteins in aqueous solution were previously observed to depend inversely on solution viscosity. In this work, the solvent-compatible enzyme Candida antarctica lipase B, which functions in aqueous and organic solvents, was modeled using molecular dynamics simulations. Data was obtained for the enzyme in acetonitrile, cyclohexane, n-butanol, and tert-butanol, in addition to water. Protein dynamics and solvation shell dynamics are characterized regionally: for each α-helix, β-sheet, and loop or connector region. Correlations are seen between solvent mobility and protein flexibility. So, does local viscosity explain the relationship between protein structural dynamics and solvation layer dynamics? Halle and Davidovic presented a cogent analysis of data describing the global hydrodynamics of a protein (tumbling in solution) that fits a model in which the protein's interfacial viscosity is higher than that of bulk water's, due to retarded water dynamics in the hydration layer (measured in NMR τ2 reorientation times). Numerous experiments have shown coupling between protein and solvation layer dynamics in site-specific measurements. Our data provides spatially-resolved characterization of solvent shell dynamics, showing correlations between regional solvation layer dynamics and protein dynamics in both aqueous and organic solvents. Correlations between protein flexibility and inverse solvent viscosity (1/η) are considered across several protein regions and for a rather disparate collection of solvents. It is seen that the correlation is consistently higher when local solvent shell dynamics are considered, rather than bulk viscosity. Protein flexibility is seen to correlate best with either the local interfacial viscosity or the ratio of the mobility of an organic solvent in a regional solvation layer relative to hydration dynamics around the same region. Results provide insight into the function of aqueous proteins, while also suggesting a framework for interpreting and predicting enzyme structural dynamics in non-aqueous solvents, based on the mobility of solvents within the solvation layer. We suggest that Kramers' theory may be used in future work to model protein conformational transitions in different solvents by incorporating local viscosity effects.
Collapse
|
24
|
Abstract
Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.
Collapse
Affiliation(s)
- R. Bryn Fenwick
- The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla CA 92037
| | - David Oyen
- The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla CA 92037
| | - H. Jane Dyson
- The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla CA 92037
| | - Peter E. Wright
- The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla CA 92037
| |
Collapse
|
25
|
Oroguchi T, Nakasako M. Influences of lone-pair electrons on directionality of hydrogen bonds formed by hydrophilic amino acid side chains in molecular dynamics simulation. Sci Rep 2017; 7:15859. [PMID: 29158598 PMCID: PMC5696464 DOI: 10.1038/s41598-017-16203-w] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Accepted: 11/08/2017] [Indexed: 11/08/2022] Open
Abstract
The influence of lone-pair electrons on the directionality of hydrogen bonds that are formed by oxygen and nitrogen atoms in the side chains of nine hydrophilic was investigated using molecular dynamics simulations. The simulations were conducted using two types of force fields; one incorporated lone-pair electrons placed at off-atom sites and the other did not. The density distributions of the hydration water molecules around the oxygen and nitrogen atoms were calculated from the simulation trajectories, and were compared with the empirical hydration distribution functions, which were constructed from a large number of hydration water molecules found in the crystal structures of proteins. Only simulations using the force field explicitly incorporating lone-pair electrons reproduced the directionality of hydrogen bonds that is observed in the empirical distribution functions for the deprotonated oxygen and nitrogen atoms in the sp 2-hybridization. The amino acids that include such atoms are functionally important glutamate, aspartate, and histidine. Therefore, a set of force field that incorporates lone-pair electrons as off-atom charge sites would be effective for considering hydrogen bond formation by these amino acids in molecular dynamics simulation studies.
Collapse
Affiliation(s)
- Tomotaka Oroguchi
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan.
- RIKEN SPring-8 Center, 1-1-1 Kohto, Sayo, Sayo-gun, Hyogo, 679-5148, Japan.
| | - Masayoshi Nakasako
- Department of Physics, Faculty of Science and Technology, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan
- RIKEN SPring-8 Center, 1-1-1 Kohto, Sayo, Sayo-gun, Hyogo, 679-5148, Japan
| |
Collapse
|
26
|
Abstract
Szent-Győrgi called water the "matrix of life" and claimed that there was no life without it. This statement is true, as far as we know, on our planet, but it is not clear whether it must hold throughout the cosmos. To evaluate that question requires a close consideration of the many varied and subtle roles that water plays in living cells-a consideration that must be free of both an assumed essentialism that gives water an almost mystical life-giving agency and a traditional tendency to see it as a merely passive solvent. Water is a participant in the "life of the cell," and here I describe some of the features of that active agency. Water's value for molecular biology comes from both the structural and dynamic characteristics of its status as a complex, structured liquid as well as its nature as a polar, protic, and amphoteric reagent. Any discussion of water as life's matrix must, however, begin with an acknowledgment that our understanding of it as both a liquid and a solvent is still incomplete.
Collapse
|
27
|
Samoila F, Sirghi L. Disjoining Pressure in Partial Wetting on the Nanoscale. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:5188-5196. [PMID: 28485609 DOI: 10.1021/acs.langmuir.7b01156] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Partial wetting on the nanoscale may result in the formation of sessile liquid nanodroplets on flat substrates. In this case, the molecular forces generate a strong interaction between nanodroplet interfaces. This interaction is expressed in the mean-field approximation by the disjoining pressure and determines an important deviation from the spherical cap shape of the nanodroplets. This deviation is observed on the atomic force microscopy images of sessile nanodroplets of oleic acid on glass. The disjoining pressure was manipulated by hydroxylation of the glass surface. This surface modification generated a strong negative disjoining pressure due to structural forces arising from the orientation of oleic acid molecules with their polar heads toward the substrate. As a result, the shape of oleic acid nanodroplets showed large deviations from the spherical cap shape, with the liquid-vapor interface tilting angle with respect to the plane substrate having a maximum (herein considered to be the contact angle) a certain distance from the substrate, followed by its decrease to zero at the droplet edge. The integration of the augmented Young-Laplace equation, where the dependence of the negative structural disjoining pressure on the interface separation distance was assumed to be an exponential decay, yielded height profiles of droplets in good agreement with the experiment.
Collapse
Affiliation(s)
- Florentina Samoila
- Iasi Plasma Advanced Research Center (IPARC), Faculty of Physics, "Alexandru Ioan Cuza", University of Iasi , Iasi-700506, Romania
| | - Lucel Sirghi
- Iasi Plasma Advanced Research Center (IPARC), Faculty of Physics, "Alexandru Ioan Cuza", University of Iasi , Iasi-700506, Romania
| |
Collapse
|
28
|
Custer GS, Das P, Matysiak S. Interplay between Conformational Heterogeneity and Hydration in the Folding Landscape of a Designed Three-Helix Bundle. J Phys Chem B 2017; 121:2731-2738. [PMID: 28282142 DOI: 10.1021/acs.jpcb.6b12286] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Water is known to play a critical role in protein folding and stability. Here we develop and employ a coarse-grained (CG) model to directly explore the role of water in shaping the conformational landscape explored during protein folding. For this purpose, we simulate a designed sequence with binary patterning of neutral and hydrophobic residues, which is capable of folding to a three-helix bundle in explicit water. We find two folded states of this sequence, with rotation of the helices occurring to trade between hydrophobic packing and water expulsion from the core. This work provides insight into the role of water and hydrophobicity in generating competing folded states for a protein.
Collapse
Affiliation(s)
- Gregory S Custer
- Fischell Department of Bioengineering, University of Maryland , College Park, Maryland 20742, United States
| | - Payel Das
- IBM Thomas J. Watson Research Center , Yorktown Heights, New York 10598, United States
| | - Silvina Matysiak
- Fischell Department of Bioengineering, University of Maryland , College Park, Maryland 20742, United States
| |
Collapse
|
29
|
High anisotropy and frustration: the keys to regulating protein function efficiently in crowded environments. Curr Opin Struct Biol 2017; 42:50-58. [DOI: 10.1016/j.sbi.2016.10.014] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 09/16/2016] [Accepted: 10/19/2016] [Indexed: 11/17/2022]
|
30
|
Sohn WY, Habka S, Gloaguen E, Mons M. Unifying the microscopic picture of His-containing turns: from gas phase model peptides to crystallized proteins. Phys Chem Chem Phys 2017. [DOI: 10.1039/c7cp03058d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The presence in crystallized proteins of a local anchoring between the side chain of a His residue, located in the central position of a γ- or β-turn, and its local main chain environment, is assessed by the comparison of protein structures with relevant isolated model peptides.
Collapse
Affiliation(s)
- Woon Yong Sohn
- LIDYL
- CEA
- CNRS
- Université Paris-Saclay
- 91191 Gif-sur-Yvette Cedex
| | - Sana Habka
- LIDYL
- CEA
- CNRS
- Université Paris-Saclay
- 91191 Gif-sur-Yvette Cedex
| | - Eric Gloaguen
- LIDYL
- CEA
- CNRS
- Université Paris-Saclay
- 91191 Gif-sur-Yvette Cedex
| | - Michel Mons
- LIDYL
- CEA
- CNRS
- Université Paris-Saclay
- 91191 Gif-sur-Yvette Cedex
| |
Collapse
|
31
|
Modestova Y, Ugarova NN. Color-shifting mutations in the C-domain of L. mingrelica firefly luciferase provide new information about the domain alternation mechanism. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2016; 1864:1818-1826. [DOI: 10.1016/j.bbapap.2016.09.007] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/01/2016] [Revised: 09/09/2016] [Accepted: 09/14/2016] [Indexed: 02/07/2023]
|