1
|
Kinose K, Shinoda K, Konishi T, Kawasaki H. Mutational analysis in Corynebacterium stationis MFS transporters for improving nucleotide bioproduction. Appl Microbiol Biotechnol 2024; 108:251. [PMID: 38436751 PMCID: PMC10912292 DOI: 10.1007/s00253-024-13080-y] [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: 10/04/2023] [Revised: 02/16/2024] [Accepted: 02/19/2024] [Indexed: 03/05/2024]
Abstract
Product secretion from an engineered cell can be advantageous for microbial cell factories. Extensive work on nucleotide manufacturing, one of the most successful microbial fermentation processes, has enabled Corynebacterium stationis to transport nucleotides outside the cell by random mutagenesis; however, the underlying mechanism has not been elucidated, hindering its applications in transporter engineering. Herein, we report the nucleotide-exporting major facilitator superfamily (MFS) transporter from the C. stationis genome and its hyperactive mutation at the G64 residue. Structural estimation and molecular dynamics simulations suggested that the activity of this transporter improved via two mechanisms: (1) enhancing interactions between transmembrane helices through the conserved "RxxQG" motif along with substrate binding and (2) trapping substrate-interacting residue for easier release from the cavity. Our results provide novel insights into how MFS transporters change their conformation from inward- to outward-facing states upon substrate binding to facilitate efflux and can contribute to the development of rational design approaches for efflux improvements in microbial cell factories. KEYPOINTS: • An MFS transporter from C. stationis genome and its mutation at residue G64 were assessed • It enhanced the transporter activity by strengthening transmembrane helix interactions and trapped substrate-interacting residues • Our results contribute to rational design approach development for efflux improvement.
Collapse
Affiliation(s)
- Keita Kinose
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
- Nagahama Institute for Biochemical Science, Oriental Yeast Co., Ltd., Nagahama, Shiga, Japan
| | - Keiko Shinoda
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
- Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo, Japan
- Research Organization of Information and Systems, The Institute of Statistical Mathematics, Tachikawa, Japan
| | - Tomoyuki Konishi
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan
| | - Hisashi Kawasaki
- Agro-Biotechnology Research Center, Graduate School of Agriculture and Life Sciences, The University of Tokyo, Tokyo, Japan.
- Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Tokyo, Japan.
| |
Collapse
|
2
|
Kondo HX, Nakamura H, Takano Y. Negative fragmentation approach for investigating the depolarization effect of neighboring residues on hydrogen bonds in π-helix. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140361] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
|
3
|
Quantum chemical studies on hydrogen bonds in helical secondary structures. Biophys Rev 2023; 14:1369-1378. [PMID: 36659988 PMCID: PMC9842822 DOI: 10.1007/s12551-022-01034-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 11/30/2022] [Indexed: 01/07/2023] Open
Abstract
We present a brief review of our recent computational studies of hydrogen bonds (H-bonds) in helical secondary structures of proteins, α-helix and 310-helix, using a Negative Fragmentation Approach with density functional theory. We found that the depolarized electronic structures of the carbonyl oxygen of the ith residue and the amide hydrogen of the (i + 4)th residue cause weaker H-bond in an α-helix than in an isolated H-bond. Our calculations showed that the H-bond energies in the 310-helix were also weaker than those of the isolated H-bonds. In the 310-helices, the adjacent N-H group at the (i + 1)th residue was closer to the C=O group of the H-bond pair than the adjacent C=O group in the 310-helices, whereas the adjacent C=O group at the (i + 1)th residue was close to the H-bond acceptor in α-helices. Therefore, the destabilization of the H-bond is attributed to the depolarization caused by the adjacent residue of the helical backbone connecting the H-bond donor and acceptor. The differences in the change in electron density revealed that such depolarizations were caused by the local electronic interactions in their neighborhood inside the helical structure and redistributed the electron density. We also present the improvements in the force field of classical molecular simulation, based on our findings. Supplementary Information The online version contains supplementary material available at 10.1007/s12551-022-01034-5.
Collapse
|
4
|
Depolarizing Effects in Hydrogen Bond Energy in 3 10-Helices Revealed by Quantum Chemical Analysis. Int J Mol Sci 2022; 23:ijms23169032. [PMID: 36012292 PMCID: PMC9409261 DOI: 10.3390/ijms23169032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2022] [Revised: 08/10/2022] [Accepted: 08/11/2022] [Indexed: 12/18/2022] Open
Abstract
Hydrogen-bond (H-bond) energies in 310-helices of short alanine peptides were systematically examined by precise DFT calculations with the negative fragmentation approach (NFA), a modified method based on the molecular tailoring approach. The contribution of each H-bond was evaluated in detail from the 310-helical conformation of total energies (whole helical model, WH3-10 model), and the results were compared with the property of H-bond in α-helix from our previous study. The H-bond energies of the WH3-10 model exhibited tendencies different from those exhibited by the α-helix in that they depended on the helical position of the relevant H-bond pair. H-bond pairs adjacent to the terminal H-bond pairs were observed to be strongly destabilized. The analysis of electronic structures indicated that structural characteristics cause the destabilization of the H-bond in 310-helices. We also found that the longer the helix length, the more stable the H-bond in the terminal pairs of the WH3-10 model, suggesting the action of H-bond cooperativity.
Collapse
|
5
|
Mahmood MI, Yamashita T. Influence of Lipid Bilayer on the GPCR Structure: Comparison of All-Atom Lipid Force Fields. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20210244] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Md. Iqbal Mahmood
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| |
Collapse
|
6
|
Nasrin SR, Ganser C, Nishikawa S, Kabir AMR, Sada K, Yamashita T, Ikeguchi M, Uchihashi T, Hess H, Kakugo A. Deformation of microtubules regulates translocation dynamics of kinesin. SCIENCE ADVANCES 2021; 7:eabf2211. [PMID: 34644102 PMCID: PMC10763888 DOI: 10.1126/sciadv.abf2211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2020] [Accepted: 08/20/2021] [Indexed: 06/13/2023]
Abstract
Microtubules, the most rigid components of the cytoskeleton, can be key transduction elements between external forces and the cellular environment. Mechanical forces induce microtubule deformation, which is presumed to be critical for the mechanoregulation of cellular events. However, concrete evidence is lacking. In this work, with high-speed atomic force microscopy, we unravel how microtubule deformation regulates the translocation of the microtubule-associated motor protein kinesin-1, responsible for intracellular transport. Our results show that the microtubule deformation by bending impedes the translocation dynamics of kinesins along them. Molecular dynamics simulation shows that the hindered translocation of kinesins can be attributed to an enhanced affinity of kinesins to the microtubule structural units in microtubules deformed by bending. This study advances our understanding of the role of cytoskeletal components in mechanotransduction.
Collapse
Affiliation(s)
| | - Christian Ganser
- Department of Creative Research, Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
| | - Seiji Nishikawa
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
| | | | - Kazuki Sada
- Faculty of Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, Tokyo 153-8904, Japan
| | - Mitsunori Ikeguchi
- Graduate School of Medical Life Science, Yokohama City University, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Takayuki Uchihashi
- Department of Creative Research, Exploratory Research Center on Life and Living Systems, National Institutes of Natural Sciences, Okazaki, Aichi 444-8787, Japan
- Department of Physics, Nagoya University, Chikusa-ku, Nagoya 464-8602, Japan
| | - Henry Hess
- Department of Biomedical Engineering, Columbia University, New York, NY 10027, USA
| | - Akira Kakugo
- Faculty of Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
- Graduate School of Chemical Sciences and Engineering, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
| |
Collapse
|
7
|
Okajima R, Hiraoka S, Yamashita T. Environmental Effects on Salt Bridge Stability in the Protein-Protein Interface: The Case of Hen Egg-White Lysozyme and Its Antibody, HyHEL-10. J Phys Chem B 2021; 125:1542-1549. [PMID: 33544613 DOI: 10.1021/acs.jpcb.0c09248] [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/29/2022]
Abstract
We studied the stability of two salt bridges between hen egg-white lysozyme (HEL) and its antibody, HyHEL-10, by using molecular dynamics simulations. It was observed that one salt bridge, D32H-K97Y, was stable, whereas the other, D99H-K97Y, was not. To understand this difference, we compared several reduced salt bridge models that incorporated the salt bridges and nearby residues. The results showed the importance of nearby residues, especially Y33H and W98H. Furthermore, to understand the effects of nearby salt bridges, we investigated two mutants, D32HA and D99HA. We found that the D32HA mutation considerably stabilized the D99H-K97Y salt bridge. The reduced model analysis indicated that this can be largely attributed to a conformational change of the main chain.
Collapse
Affiliation(s)
- Ryo Okajima
- Department of Basic Science, Graduate School of Arts and Sciences, the University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan.,Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Shuichi Hiraoka
- Department of Basic Science, Graduate School of Arts and Sciences, the University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| |
Collapse
|
8
|
TANIDA Y, SATO H, MANABE T, TERASHIMA C. Fast Conformation Search of Macrocyclic Peptides Using a Combination of Digital Annealer and REST2. JOURNAL OF COMPUTER CHEMISTRY-JAPAN 2021. [DOI: 10.2477/jccj.2021-0036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Yoshiaki TANIDA
- FUJITSU LIMITED, 1-5 Ohmiya-cho, Kawasaki-ku, Kawasaki, Kanagawa 212-0014
| | - Hiroyuki SATO
- FUJITSU LIMITED, 1-5 Ohmiya-cho, Kawasaki-ku, Kawasaki, Kanagawa 212-0014
| | - Toshio MANABE
- FUJITSU LIMITED, 1-5 Ohmiya-cho, Kawasaki-ku, Kawasaki, Kanagawa 212-0014
| | - Chieko TERASHIMA
- FUJITSU LIMITED, 1-5 Ohmiya-cho, Kawasaki-ku, Kawasaki, Kanagawa 212-0014
| |
Collapse
|
9
|
Qiu Y, Smith DGA, Stern CD, Feng M, Jang H, Wang LP. Driving torsion scans with wavefront propagation. J Chem Phys 2020; 152:244116. [PMID: 32610969 DOI: 10.1063/5.0009232] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
The parameterization of torsional/dihedral angle potential energy terms is a crucial part of developing molecular mechanics force fields. Quantum mechanical (QM) methods are often used to provide samples of the potential energy surface (PES) for fitting the empirical parameters in these force field terms. To ensure that the sampled molecular configurations are thermodynamically feasible, constrained QM geometry optimizations are typically carried out, which relax the orthogonal degrees of freedom while fixing the target torsion angle(s) on a grid of values. However, the quality of results and computational cost are affected by various factors on a non-trivial PES, such as dependence on the chosen scan direction and the lack of efficient approaches to integrate results started from multiple initial guesses. In this paper, we propose a systematic and versatile workflow called TorsionDrive to generate energy-minimized structures on a grid of torsion constraints by means of a recursive wavefront propagation algorithm, which resolves the deficiencies of conventional scanning approaches and generates higher quality QM data for force field development. The capabilities of our method are presented for multi-dimensional scans and multiple initial guess structures, and an integration with the MolSSI QCArchive distributed computing ecosystem is described. The method is implemented in an open-source software package that is compatible with many QM software packages and energy minimization codes.
Collapse
Affiliation(s)
- Yudong Qiu
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| | - Daniel G A Smith
- The Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Chaya D Stern
- Computational and Systems Biology Program, Sloan-Kettering Institute, New York, New York 10065, USA
| | - Mudong Feng
- Department of Chemistry and Biochemistry, UC San Diego, La Jolla, California 92093, USA
| | - Hyesu Jang
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| | - Lee-Ping Wang
- Department of Chemistry, UC Davis, Davis, California 95616, USA
| |
Collapse
|
10
|
Tanida Y, Matsuura A. Alchemical free energy calculations via metadynamics: Application to the theophylline-RNA aptamer complex. J Comput Chem 2020; 41:1804-1819. [PMID: 32449538 DOI: 10.1002/jcc.26221] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 04/03/2020] [Accepted: 04/26/2020] [Indexed: 11/08/2022]
Abstract
We propose a computational workflow for robust and accurate prediction of both binding poses and their affinities at early stage in designing drug candidates. Small, rigid ligands with few intramolecular degrees of freedom, for example, fragment-like molecules, have multiple binding poses, even at a single binding site, and their affinities are often close to each other. We explore various structures of ligand binding to a target through metadynamics using a small number of collective variables, followed by reweighting to obtain the atomic coordinates. After identifying each binding pose by cluster analysis, we perform alchemical free energy calculations on each structure to obtain the overall value. We applied this protocol in computing free energy of binding for the theophylline-RNA aptamer complex. Of the six (meta)stable structures found, the most favorable binding structure is consistent with the structure obtained by NMR. The overall free energy of binding reproduces the experimental values very well.
Collapse
|
11
|
Kamiya N, Kayanuma M, Fujitani H, Shinoda K. A New Lipid Force Field (FUJI). J Chem Theory Comput 2020; 16:3664-3676. [PMID: 32384238 DOI: 10.1021/acs.jctc.9b01195] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
To explore inhomogeneous and anisotropic systems such as lipid bilayers, the Lennard-Jones particle mesh Ewald (LJ-PME) method has been applied without a conventional isotropic dispersion correction. As the popular AMBER and CHARMM lipid force fields were developed using a cutoff scheme, their lipid bilayers unacceptably shrink when using the LJ-PME method. In this study, a new all-atom lipid force field (FUJI) was developed on the basis of the AMBER force-field scheme including the Lipid14 van der Waals parameters. Point charges were calculated using the restrained electrostatic potentials of many lipid conformers. Further, torsion energy profiles were calculated using high-level ab initio molecular orbitals (LCCSD(T)/Aug-cc-pVTZ//LMP2/Aug-cc-pVTZ), following which the molecular mechanical dihedral parameters were derived through a fast Fourier transform. By incorporation of these parameters into a new lipid force field without fitting experimental data, the desired lipid characteristics such as the area per lipid and lateral diffusion coefficients were obtained through GROMACS molecular dynamics simulations using the LJ-PME method and virtual hydrogen sites. The calculated area per lipid and lateral diffusion coefficients showed satisfactory agreement with experimental data. Furthermore, the electron-density profiles along the membrane normal were calculated for pure lipid bilayers, and the resulting membrane thicknesses agreed well with the experimental values. As the new lipid force field is compatible with FUJI for protein and small molecules, the new FUJI force field will offer accurate modeling for complex systems consisting of various membrane proteins and lipids.
Collapse
Affiliation(s)
- Nozomu Kamiya
- Fujitsu Limited Bio-IT R&D Office, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Megumi Kayanuma
- Research Center for Computational Design of Advanced Functional Materials, National Institute of Advanced Industrial Science and Technology, 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Keiko Shinoda
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| |
Collapse
|
12
|
Kondo HX, Kusaka A, Kitakawa CK, Onari J, Yamanaka S, Nakamura H, Takano Y. Hydrogen bond donors and acceptors are generally depolarized in α-helices as revealed by a molecular tailoring approach. J Comput Chem 2019; 40:2043-2052. [PMID: 31099907 PMCID: PMC6767508 DOI: 10.1002/jcc.25859] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2019] [Revised: 04/09/2019] [Accepted: 04/29/2019] [Indexed: 11/24/2022]
Abstract
Hydrogen-bond (H-bond) interaction energies in α-helices of short alanine peptides were systematically examined by precise density functional theory calculations, followed by a molecular tailoring approach. The contribution of each H-bond interaction in α-helices was estimated in detail from the entire conformation energies, and the results were compared with those in the minimal H-bond models, in which only H-bond donors and acceptors exist with the capping methyl groups. The former interaction energies were always significantly weaker than the latter energies, when the same geometries of the H-bond donors and acceptors were applied. The chemical origin of this phenomenon was investigated by analyzing the differences among the electronic structures of the local peptide backbones of the α-helices and those of the minimal H-bond models. Consequently, we found that the reduced H-bond energy originated from the depolarizations of both the H-bond donor and acceptor groups, due to the repulsive interactions with the neighboring polar peptide groups in the α-helix backbone. The classical force fields provide similar H-bond energies to those in the minimal H-bond models, which ignore the current depolarization effect, and thus they overestimate the actual H-bond energies in α-helices. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Hiroko X. Kondo
- School of Regional Innovation and Social Design EngineeringFaculty of Engineering, Kitami Institute of Technology, 165 Koen‐choKitamiHokkaido090‐8507Japan
| | - Ayumi Kusaka
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
| | - Colin K. Kitakawa
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Jinta Onari
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Shusuke Yamanaka
- Graduate School of ScienceOsaka University, 1‐1 MachikaneyamachoToyonakaOsaka560‐0043Japan
| | - Haruki Nakamura
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
| | - Yu Takano
- Institute for Protein ResearchOsaka University, 3‐2 YamadaokaSuitaOsaka565‐0871Japan
- Graduate School of Information SciencesHiroshima City University, 3‐4‐1 Ozuka‐Higashi Asa‐Minami‐KuHiroshima731‐3194Japan
| |
Collapse
|
13
|
Koyama T, Nakamoto M, Morishima K, Yamashita R, Yamashita T, Sasaki K, Kuruma Y, Mizuno N, Suzuki M, Okada Y, Ieda R, Uchino T, Tasumi S, Hosoya S, Uno S, Koyama J, Toyoda A, Kikuchi K, Sakamoto T. A SNP in a Steroidogenic Enzyme Is Associated with Phenotypic Sex in Seriola Fishes. Curr Biol 2019; 29:1901-1909.e8. [PMID: 31130458 DOI: 10.1016/j.cub.2019.04.069] [Citation(s) in RCA: 51] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2018] [Revised: 03/03/2019] [Accepted: 04/26/2019] [Indexed: 12/30/2022]
Abstract
Vertebrate sex development consists largely of two processes: "sex determination," the initial bifurcation of sexual identity, and "sex differentiation," which subsequently facilitates maleness or femaleness according to the sex determination signal. Steroid hormones promote multiple types of sexual dimorphism in eutherian mammals and avians [1-3], in which they are indispensable for proper sex differentiation. By contrast, in many poikilothermic vertebrates, steroid hormones have been proposed to be key players in sex determination as well as sex differentiation [4-8]. This hypothesis was introduced more than 50 years ago but has never been rigorously tested due to difficulties in discriminating the roles of steroids in sex determination and differentiation. We found that a missense SNP in the gene encoding the steroidogenic enzyme 17β-hydroxysteroid dehydrogenase 1 (Hsd17b1) was perfectly associated with ZZ/ZW sex determination in Seriola fishes. Biochemical analyses revealed that a glutamate residue present specifically in Z-type HSD17B1 attenuated interconversion between 17-keto and 17β-hydroxy steroids relative to the allelic product from the W chromosome, which harbors glycine at that position, by disrupting the hydrogen bond network between the steroid and the enzyme's catalytic residues. Hsd17b1 mRNA is constitutively expressed in undifferentiated and differentiating gonads of both genotypic sexes, whereas W-type mRNA is expressed only in genotypic females. Meanwhile, Cyp19a1 is predominantly expressed in differentiating ovary. We conclude that the combination of Hsd17b1 alleles determines sex by modulating endogenous estrogen levels in Seriola species. These findings strongly support the long-standing hypothesis on steroids in sex determination.
Collapse
Affiliation(s)
- Takashi Koyama
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Masatoshi Nakamoto
- Department of Marine Biosciences, Tokyo University of Marine Science and Technology, 4-5-7 Konan, Minato, Tokyo 108-8477, Japan
| | - Kagayaki Morishima
- Oita Marine Biological Technology Center, Nippon Suisan Kaisha, Ltd., 508-8 Ariakeura, Tsurumi, Saeki, Oita 876-1204, Japan
| | - Ryohei Yamashita
- Oita Marine Biological Technology Center, Nippon Suisan Kaisha, Ltd., 508-8 Ariakeura, Tsurumi, Saeki, Oita 876-1204, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro, Tokyo 153-8904, Japan
| | - Kohei Sasaki
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Komaba, Meguro, Tokyo 153-8904, Japan
| | - Yosuke Kuruma
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Naoki Mizuno
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Moe Suzuki
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Yoshiharu Okada
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Risa Ieda
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Tsubasa Uchino
- Department of Marine Biosciences, Tokyo University of Marine Science and Technology, 4-5-7 Konan, Minato, Tokyo 108-8477, Japan
| | - Satoshi Tasumi
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Sho Hosoya
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan
| | - Seiichi Uno
- Education and Research Center for Marine Resources and Environment, Faculty of Fisheries, Kagoshima University, 50-20 Shimoarata 4-Chome, Kagoshima 890-0056, Japan
| | - Jiro Koyama
- Education and Research Center for Marine Resources and Environment, Faculty of Fisheries, Kagoshima University, 50-20 Shimoarata 4-Chome, Kagoshima 890-0056, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, Center for Information Biology, National Institute of Genetics, 1111 Yata, Mishima, Shizuoka 411-8540, Japan
| | - Kiyoshi Kikuchi
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, 2971-4 Bentenjima, Maisaka, Hamamatsu, Shizuoka 431-0214, Japan.
| | - Takashi Sakamoto
- Department of Marine Biosciences, Tokyo University of Marine Science and Technology, 4-5-7 Konan, Minato, Tokyo 108-8477, Japan.
| |
Collapse
|
14
|
Yamashita T, Mizohata E, Nagatoishi S, Watanabe T, Nakakido M, Iwanari H, Mochizuki Y, Nakayama T, Kado Y, Yokota Y, Matsumura H, Kawamura T, Kodama T, Hamakubo T, Inoue T, Fujitani H, Tsumoto K. Affinity Improvement of a Cancer-Targeted Antibody through Alanine-Induced Adjustment of Antigen-Antibody Interface. Structure 2018; 27:519-527.e5. [PMID: 30595454 DOI: 10.1016/j.str.2018.11.002] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2017] [Revised: 08/13/2018] [Accepted: 11/01/2018] [Indexed: 12/19/2022]
Abstract
To investigate favorable single amino acid substitutions that improve antigen-antibody interactions, alanine (Ala) mutagenesis scanning of the interfacial residues of a cancer-targeted antibody, B5209B, was performed based on X-ray crystallography analysis. Two substitutions were shown to significantly enhance the binding affinity for the antigen, by up to 30-fold. One substitution improved the affinity by a gain of binding enthalpy, whereas the other substitution improved the affinity by a gain of binding entropy. Molecular dynamics simulations showed that the enthalpic improvement could be attributed to the stabilization of distant salt bridges located at the periphery of the antigen-antibody interface. The entropic improvement was due to the release of water molecules that were stably trapped in the antigen-antibody interface of the wild-type antibody. Importantly, these effects of the Ala substitutions were caused by subtle adjustments of the binding interface. These results will be helpful to design high-affinity antibodies with avoiding entropy-enthalpy compensation.
Collapse
Affiliation(s)
- Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Eiichi Mizohata
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Satoru Nagatoishi
- Department of Bioengineering, School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Takahiro Watanabe
- Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 4-6-1, Shirokanedai, Minato-ku, Tokyo 108-8639, Japan
| | - Makoto Nakakido
- Institute of Medical Science, The University of Tokyo, 4-6-1, Shirokanedai, Minato-ku, Tokyo 108-8639, Japan
| | - Hiroko Iwanari
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Yasuhiro Mochizuki
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Taisuke Nakayama
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yuji Kado
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Yuki Yokota
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Hiroyoshi Matsumura
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Takeshi Kawamura
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Tatsuhiko Kodama
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Takao Hamakubo
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Tsuyoshi Inoue
- Department of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamadaoka, Suita, Osaka 565-0871, Japan.
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1, Komaba, Meguro-ku, Tokyo 153-8904, Japan.
| | - Kouhei Tsumoto
- Department of Bioengineering, School of Engineering, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-8656, Japan; Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, 4-6-1, Shirokanedai, Minato-ku, Tokyo 108-8639, Japan; Institute of Medical Science, The University of Tokyo, 4-6-1, Shirokanedai, Minato-ku, Tokyo 108-8639, Japan.
| |
Collapse
|
15
|
Initiation of prolyl cis-trans isomerisation in the CDR-H3 loop of an antibody in response to antigen binding. Sci Rep 2017; 7:16964. [PMID: 29208911 PMCID: PMC5717248 DOI: 10.1038/s41598-017-16766-8] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Accepted: 11/16/2017] [Indexed: 01/12/2023] Open
Abstract
Proline cis-trans isomerisation is a regulatory mechanism used in a range of biological processes, and is related to various diseases such as Alzheimers disease and cancer. However, the details of the exact molecular mechanism by which it occurs are not known. Using X-ray crystallography, proline isomerisation has been shown to occur following formation of an antigen-antibody complex between the target epiregulin (EPR) and the antibody 9E5, at proline (Pro103), located in the third complementarity-determining region (CDR) of the heavy chain of 9E5. To obtain an accurate description of the pathway involved in cis-trans isomerisation in this system, we performed ten independent long molecular dynamics (MD) simulations starting at a stable transient bound structure obtained from many short binding MD simulations. As a result, we were able to describe the process by which cis-trans isomerisation is initiated, and suggest a catalysis mechanism for cis-trans isomerization in this antigen-antibody system. We found that Asp102, which is immediately adjacent to Pro103, rotates while changing its interacting partner residues in the light chain of 9E5, and at the same time EPR polar residues help to stabilise the intermediate states in the isomerisation process by interacting strongly with Asp102.
Collapse
|
16
|
Balaji GA, Nagendra HG, Balaji VN, Rao SN. Experimental conformational energy maps of proteins and peptides. Proteins 2017; 85:979-1001. [PMID: 28168743 DOI: 10.1002/prot.25266] [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: 10/05/2016] [Revised: 01/26/2017] [Accepted: 01/30/2017] [Indexed: 01/26/2023]
Abstract
We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X-ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (ϕ, ψ) plots in which global minima are predominantly observed either in the right-handed α-helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (ϕ,ψ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol-1 ) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X-ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979-1001. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Govardhan A Balaji
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - H G Nagendra
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - Vitukudi N Balaji
- Department of Biotechnology, Sir M Visvesvaraya Institute of Technology, Bangalore, 562157, India
| | - Shashidhar N Rao
- Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey, 08552
| |
Collapse
|
17
|
Jangra H, Haindl MH, Achrainer F, Hioe J, Gschwind RM, Zipse H. Conformational Preferences in Small Peptide Models: The Relevance ofcis/trans-Conformations. Chemistry 2016; 22:13328-35. [DOI: 10.1002/chem.201601828] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Indexed: 11/12/2022]
Affiliation(s)
- Harish Jangra
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| | - Michael H. Haindl
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Florian Achrainer
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| | - Johnny Hioe
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Ruth M. Gschwind
- Institut für Organische Chemie; Universität Regensburg; 93053 Regensburg Germany
| | - Hendrik Zipse
- Department of Chemistry; LMU München; Butenandstrasse 5-14 81377 München Germany
| |
Collapse
|
18
|
Sakano T, Mahamood MI, Yamashita T, Fujitani H. Molecular dynamics analysis to evaluate docking pose prediction. Biophys Physicobiol 2016; 13:181-194. [PMID: 27924273 PMCID: PMC5042163 DOI: 10.2142/biophysico.13.0_181] [Citation(s) in RCA: 35] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2015] [Accepted: 02/27/2016] [Indexed: 02/06/2023] Open
Abstract
The accurate prediction of a ligand-protein complex structure is important for computer-assisted drug development. Although many docking methods have been developed over the last three decades, the success of binding structure prediction remains greatly limited. The purpose of this study was to demonstrate the usefulness of molecular dynamics (MD) simulation in assessing a docking pose predicted using a docking program. If the predicted pose is not unstable in an aqueous environment, MD simulation equilibrates the system and removes the ligand from the predicted position. Here we investigated two proteins that are important potential therapeutic targets: β2 adrenergic receptor (β2AR) and PR-Set7. While β2AR is rigid and its ligands are very similar to the template ligand (carazolol), PR-Set7 is very flexible and its ligands vary greatly from the template ligand (histone H4 tail peptide). On an empirical basis, we usually expect that the docking prediction is accurate when the protein is rigid and its ligands are similar to the template ligand. The MD analyses in this study clearly suggested such a tendency. Furthermore, we discuss the possibility that the MD simulation can predict the binding pose of a ligand.
Collapse
Affiliation(s)
- Takako Sakano
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, Tokyo 153-8904, Japan
| | - Md Iqbal Mahamood
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, Tokyo 153-8904, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, Tokyo 153-8904, Japan
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, the University of Tokyo, Tokyo 153-8904, Japan
| |
Collapse
|
19
|
Kaminský J, Jensen F. Conformational Interconversions of Amino Acid Derivatives. J Chem Theory Comput 2016; 12:694-705. [PMID: 26691979 DOI: 10.1021/acs.jctc.5b00911] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.
Collapse
Affiliation(s)
- Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry, Flemingovo nám. 2, 166 10 Prague, Czech Republic
| | - Frank Jensen
- Department of Chemistry, Aarhus University , Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
20
|
Kado Y, Mizohata E, Nagatoishi S, Iijima M, Shinoda K, Miyafusa T, Nakayama T, Yoshizumi T, Sugiyama A, Kawamura T, Lee YH, Matsumura H, Doi H, Fujitani H, Kodama T, Shibasaki Y, Tsumoto K, Inoue T. Epiregulin Recognition Mechanisms by Anti-epiregulin Antibody 9E5: STRUCTURAL, FUNCTIONAL, AND MOLECULAR DYNAMICS SIMULATION ANALYSES. J Biol Chem 2015; 291:2319-30. [PMID: 26627827 DOI: 10.1074/jbc.m115.656009] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2015] [Indexed: 11/06/2022] Open
Abstract
Epiregulin (EPR) is a ligand of the epidermal growth factor (EGF) family that upon binding to its epidermal growth factor receptor (EGFR) stimulates proliferative signaling, especially in colon cancer cells. Here, we describe the three-dimensional structure of the EPR antibody (the 9E5(Fab) fragment) in the presence and absence of EPR. Among the six complementarity-determining regions (CDRs), CDR1-3 in the light chain and CDR2 in the heavy chain predominantly recognize EPR. In particular, CDR3 in the heavy chain dramatically moves with cis-trans isomerization of Pro(103). A molecular dynamics simulation and mutational analyses revealed that Arg(40) in EPR is a key residue for the specific binding of 9E5 IgG. From isothermal titration calorimetry analysis, the dissociation constant was determined to be 6.5 nm. Surface plasmon resonance analysis revealed that the dissociation rate of 9E5 IgG is extremely slow. The superimposed structure of 9E5(Fab)·EPR on the known complex structure of EGF·EGFR showed that the 9E5(Fab) paratope overlaps with Domains I and III on the EGFR, which reveals that the 9E5(Fab)·EPR complex could not bind to the EGFR. The 9E5 antibody will also be useful in medicine as a neutralizing antibody specific for colon cancer.
Collapse
Affiliation(s)
- Yuji Kado
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan, Interdisciplinary Program for Biomedical Sciences, Institute for Academic Initiatives, Osaka University, 1-1 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Eiichi Mizohata
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
| | - Satoru Nagatoishi
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of Tokyo, 4-6-1 Shirokanedai, Minato-ku Tokyo 108-8639, Japan, and
| | - Mariko Iijima
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Keiko Shinoda
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Takamitsu Miyafusa
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of Tokyo, 4-6-1 Shirokanedai, Minato-ku Tokyo 108-8639, Japan, and
| | - Taisuke Nakayama
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
| | - Takuma Yoshizumi
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
| | - Akira Sugiyama
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Takeshi Kawamura
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Young-Hun Lee
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Hiroyoshi Matsumura
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan
| | - Hirofumi Doi
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Tatsuhiko Kodama
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Yoshikazu Shibasaki
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, 4-6-1 Number 34 Komaba, Meguro-ku, Tokyo 153-8904, Japan
| | - Kouhei Tsumoto
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of Tokyo, 4-6-1 Shirokanedai, Minato-ku Tokyo 108-8639, Japan, and
| | - Tsuyoshi Inoue
- From the Division of Applied Chemistry, Graduate School of Engineering, Osaka University, 2-1 Yamada-Oka, Suita, Osaka 565-0871, Japan,
| |
Collapse
|
21
|
Cormanich RA, Bühl M, Rittner R. Understanding the conformational behaviour of Ac-Ala-NHMe in different media. A joint NMR and DFT study. Org Biomol Chem 2015. [PMID: 26219244 DOI: 10.1039/c5ob01296a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The conformational behaviour of Ac-Ala-NHMe was studied in the gas-phase and in solution by theoretical calculations (B3LYP-D3/aug-cc-pVDZ level) and experimental (1)H NMR. The conformational preferences of this compound were shown to result from a complex interplay between the strengths of possible intramolecular hydrogen bonds, steric interactions, hyperconjugation, entropy effects and the overall dipole moments. The Ac-Ala-N(Me)2 derivative was studied in addition, to design a system akin to Ac-Ala-NHMe, but with disrupted intramolecular hydrogen bonds involving the -NHMe group, mimicking the effect of polar protic solvents.
Collapse
Affiliation(s)
- Rodrigo A Cormanich
- Chemistry Institute, University of Campinas, Campinas, SP 13083-970, Brazil.
| | | | | |
Collapse
|
22
|
Robertson M, Tirado-Rives J, Jorgensen WL. Improved Peptide and Protein Torsional Energetics with the OPLSAA Force Field. J Chem Theory Comput 2015; 11:3499-509. [PMID: 26190950 PMCID: PMC4504185 DOI: 10.1021/acs.jctc.5b00356] [Citation(s) in RCA: 482] [Impact Index Per Article: 53.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2015] [Indexed: 11/28/2022]
Abstract
The development and validation of new peptide dihedral parameters are reported for the OPLS-AA force field. High accuracy quantum chemical methods were used to scan φ, ψ, χ1, and χ2 potential energy surfaces for blocked dipeptides. New Fourier coefficients for the dihedral angle terms of the OPLS-AA force field were fit to these surfaces, utilizing a Boltzmann-weighted error function and systematically examining the effects of weighting temperature. To prevent overfitting to the available data, a minimal number of new residue-specific and peptide-specific torsion terms were developed. Extensive experimental solution-phase and quantum chemical gas-phase benchmarks were used to assess the quality of the new parameters, named OPLS-AA/M, demonstrating significant improvement over previous OPLS-AA force fields. A Boltzmann weighting temperature of 2000 K was determined to be optimal for fitting the new Fourier coefficients for dihedral angle parameters. Conclusions are drawn from the results for best practices for developing new torsion parameters for protein force fields.
Collapse
Affiliation(s)
- Michael
J. Robertson
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - Julian Tirado-Rives
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| | - William L. Jorgensen
- Department of Chemistry, Yale University, New Haven, Connecticut 06520-8107, United States
| |
Collapse
|
23
|
Kurauchi R, Watanabe C, Fukuzawa K, Tanaka S. Novel type of virtual ligand screening on the basis of quantum-chemical calculations for protein–ligand complexes and extended clustering techniques. COMPUT THEOR CHEM 2015. [DOI: 10.1016/j.comptc.2015.02.016] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
24
|
Nakayama T, Mizohata E, Yamashita T, Nagatoishi S, Nakakido M, Iwanari H, Mochizuki Y, Kado Y, Yokota Y, Satoh R, Tsumoto K, Fujitani H, Kodama T, Hamakubo T, Inoue T. Structural features of interfacial tyrosine residue in ROBO1 fibronectin domain-antibody complex: Crystallographic, thermodynamic, and molecular dynamic analyses. Protein Sci 2015; 24:328-40. [PMID: 25492858 PMCID: PMC4353359 DOI: 10.1002/pro.2619] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2014] [Revised: 12/02/2014] [Accepted: 12/04/2014] [Indexed: 11/09/2022]
Abstract
ROBO1, fibronectin Type-III domain (Fn)-containing protein, is a novel immunotherapeutic target for hepatocellular carcinoma in humans. The crystal structure of the antigen-binding fragment (Fab) of B2212A, the monoclonal antibody against the third Fn domain (Fn3) of ROBO1, was determined in pursuit of antibody drug for hepatocellular carcinoma. This effort was conducted in the presence or absence of the antigen, with the chemical features being investigated by determining the affinity of the antibody using molecular dynamics (MD) and thermodynamics. The structural comparison of B2212A Fab between the complex and the free form revealed that the interfacial Tyr(L) 50 (superscripts L, H, and F stand for the residues in the light chain, heavy chain, and Fn3, respectively) played important roles in Fn3 recognition. That is, the aromatic ring of Tyr(L) 50 pivoted toward Phe(F) 68, forming a CH/π interaction and a new hydrogen bond with the carbonyl O atom of Phe(F) 68. MD simulations predicted that the Tyr(L) 50-Phe(F) 68 interaction almost entirely dominated Fab-Fn3 binding, and Ala-substitution of Tyr(L) 50 led to a reduced binding of the resultant complex. On the contrary, isothermal titration calorimetry experiments underscored that Ala-substitution of Tyr(L) 50 caused an increase of the binding enthalpy between B2212A and Fn3, but importantly, it induced an increase of the binding entropy, resulting in a suppression of loss in the Gibbs free energy in total. These results suggest that mutation analysis considering the binding entropy as well as the binding enthalpy will aid in the development of novel antibody drugs for hepatocellular carcinoma.
Collapse
Affiliation(s)
- Taisuke Nakayama
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| | - Eiichi Mizohata
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| | - Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Satoru Nagatoishi
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of TokyoTokyo, Japan
| | - Makoto Nakakido
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of TokyoTokyo, Japan
| | - Hiroko Iwanari
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Yasuhiro Mochizuki
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Yuji Kado
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| | - Yuki Yokota
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| | - Reiko Satoh
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| | - Kouhei Tsumoto
- Medical Proteomics Laboratory, The Institute of Medical Science, The University of TokyoTokyo, Japan
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Tatsuhiko Kodama
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Takao Hamakubo
- Department of Quantitative Biology and Medicine, Research Center for Advanced Science and Technology, The University of TokyoTokyo, Japan
| | - Tsuyoshi Inoue
- Structural Physical Chemistry, Division of Applied Chemistry, Graduate School of Engineering, Osaka UniversityOsaka, Japan
| |
Collapse
|
25
|
Yamashita T, Ueda A, Mitsui T, Tomonaga A, Matsumoto S, Kodama T, Fujitani H. The Feasibility of an Efficient Drug Design Method with High-Performance Computers. Chem Pharm Bull (Tokyo) 2015; 63:147-55. [DOI: 10.1248/cpb.c14-00596] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Affiliation(s)
- Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo
| | - Akihiko Ueda
- Bio-IT R&D Office, Next-Generation Healthcare Innovation Center, Fujitsu Limited
| | - Takashi Mitsui
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo
- Bio-IT R&D Office, Next-Generation Healthcare Innovation Center, Fujitsu Limited
| | - Atsushi Tomonaga
- Bio-IT R&D Office, Next-Generation Healthcare Innovation Center, Fujitsu Limited
| | - Shunji Matsumoto
- Bio-IT R&D Office, Next-Generation Healthcare Innovation Center, Fujitsu Limited
| | - Tatsuhiko Kodama
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo
| | - Hideaki Fujitani
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo
| |
Collapse
|
26
|
On accurate calculation of the potential of mean force between antigen and antibody: A case of the HyHEL-10-hen egg white lysozyme system. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.06.028] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
27
|
Nakamura M, Obata M, Morishita T, Oda T. An ab initio approach to free-energy reconstruction using logarithmic mean force dynamics. J Chem Phys 2014; 140:184110. [PMID: 24832256 DOI: 10.1063/1.4874654] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present an ab initio approach for evaluating a free energy profile along a reaction coordinate by combining logarithmic mean force dynamics (LogMFD) and first-principles molecular dynamics. The mean force, which is the derivative of the free energy with respect to the reaction coordinate, is estimated using density functional theory (DFT) in the present approach, which is expected to provide an accurate free energy profile along the reaction coordinate. We apply this new method, first-principles LogMFD (FP-LogMFD), to a glycine dipeptide molecule and reconstruct one- and two-dimensional free energy profiles in the framework of DFT. The resultant free energy profile is compared with that obtained by the thermodynamic integration method and by the previous LogMFD calculation using an empirical force-field, showing that FP-LogMFD is a promising method to calculate free energy without empirical force-fields.
Collapse
Affiliation(s)
- Makoto Nakamura
- Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan
| | - Masao Obata
- Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan
| | - Tetsuya Morishita
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba 305-8568, Japan
| | - Tatsuki Oda
- Graduate School of Natural Science and Technology, Kanazawa University, Kanazawa 920-1192, Japan
| |
Collapse
|
28
|
Tzanov AT, Cuendet MA, Tuckerman ME. How Accurately Do Current Force Fields Predict Experimental Peptide Conformations? An Adiabatic Free Energy Dynamics Study. J Phys Chem B 2014; 118:6539-52. [DOI: 10.1021/jp500193w] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Alexandar T. Tzanov
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Michel A. Cuendet
- Department of Chemistry, New York University, New York, New York 10003, United States
- Swiss Institute of Bioinformatics, UNIL Sorge, 1015 Lausanne, Switzerland
| | - Mark E. Tuckerman
- Department of Chemistry, New York University, New York, New York 10003, United States
- Courant Institute
of Mathematical Sciences, New York University, New York, New York 10003, United States
| |
Collapse
|
29
|
Kang YK, Park HS. Assessment of CCSD(T), MP2, DFT-D, CBS-QB3, and G4(MP2) methods for conformational study of alanine and proline dipeptides. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.03.067] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
30
|
Jiang F, Han W, Wu YD. The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development. Phys Chem Chem Phys 2013; 15:3413-28. [PMID: 23385383 DOI: 10.1039/c2cp43633g] [Citation(s) in RCA: 52] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.
Collapse
Affiliation(s)
- Fan Jiang
- Laboratory of Computational Chemistry and Drug Design, Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | | | | |
Collapse
|
31
|
Parchaňský V, Kapitán J, Kaminský J, Šebestík J, Bouř P. Ramachandran Plot for Alanine Dipeptide as Determined from Raman Optical Activity. J Phys Chem Lett 2013; 4:2763-2768. [PMID: 26706714 DOI: 10.1021/jz401366j] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Accessible values of the φ and ψ torsional angles determining peptide main chain conformation are traditionally displayed in the form of Ramachandran plots. The number of experimental methods making it possible to determine such conformational distribution is limited. In the present study, Raman optical activity (ROA) spectra of Ac-Ala-NHMe were measured and fit by theoretical curves. This revealed the most favored conformers and a large part of the potential energy surface (PES) of this model dipeptide. Such experimental PES compares well to quantum chemical computations, whereas molecular dynamics (MD) modeling reproduces it less faithfully. The surface shape is consistent with the temperature dependence of the spectra, as observed experimentally and predicted by MD. Despite errors associated with spectral modeling and the measurement, the results are likely to facilitate future applications of ROA spectroscopy.
Collapse
Affiliation(s)
- Václav Parchaňský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences , Flemingovo náměstí 2, 16610 Prague, Czech Republic
- Department of Analytical Chemistry, Institute of Chemical Technology , Technická 5, 16628 Prague, Czech Republic
| | - Josef Kapitán
- Department of Optics, Palacký University Olomouc , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jakub Kaminský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences , Flemingovo náměstí 2, 16610 Prague, Czech Republic
| | - Jaroslav Šebestík
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences , Flemingovo náměstí 2, 16610 Prague, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences , Flemingovo náměstí 2, 16610 Prague, Czech Republic
| |
Collapse
|
32
|
Okada O, Yamashita H, Takedomi K, Ono S, Sunada S, Kubodera H. Prediction of the binding affinity of compounds with diverse scaffolds by MP-CAFEE. Biophys Chem 2013; 180-181:119-26. [PMID: 23938954 DOI: 10.1016/j.bpc.2013.07.005] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2013] [Revised: 07/02/2013] [Accepted: 07/15/2013] [Indexed: 11/16/2022]
Abstract
Accurate methods to predict the binding affinities of compounds for target molecules are powerful tools in structure-based drug design (SBDD). A recently developed method called massively parallel computation of absolute binding free energy with a well-equilibrated system (MP-CAFEE) successfully predicted the binding affinities of compounds with relatively similar scaffolds. We investigate the applicability of MP-CAFEE for predicting the affinity of compounds having more diverse scaffolds for the target p38α, a mitogen-activated protein kinase. The calculated and experimental binding affinities correlate well, showing that MP-CAFEE can accurately rank the compounds with diverse scaffolds. We propose a method to determine the optimal number of sampling runs with respect to a predefined level of accuracy, which is established according to the stage in the SBDD process being considered. The optimal number of sampling runs for two key stages-lead identification and lead optimization-is estimated to be five and eight or more, respectively, in our model system using Cochrans sample size formula.
Collapse
Affiliation(s)
- Okimasa Okada
- Medicinal Chemistry Research Laboratories II, Mitsubishi Tanabe Pharma Corporation, Toda, Saitama 335-8505, Japan.
| | | | | | | | | | | |
Collapse
|
33
|
Hayashi T, Matsuura A, Sato H, Sakurai M. Full-Quantum chemical calculation of the absorption maximum of bacteriorhodopsin: a comprehensive analysis of the amino acid residues contributing to the opsin shift. Biophysics (Nagoya-shi) 2012; 8:115-25. [PMID: 27493528 PMCID: PMC4629650 DOI: 10.2142/biophysics.8.115] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2012] [Accepted: 06/14/2012] [Indexed: 12/01/2022] Open
Abstract
Herein, the absorption maximum of bacteriorhodopsin (bR) is calculated using our recently developed method in which the whole protein can be treated quantum mechanically at the level of INDO/S-CIS//ONIOM (B3LYP/6-31G(d,p): AMBER). The full quantum mechanical calculation is shown to reproduce the so-called opsin shift of bR with an error of less than 0.04 eV. We also apply the same calculation for 226 different bR mutants, each of which was constructed by replacing any one of the amino acid residues of the wild-type bR with Gly. This substitution makes it possible to elucidate the extent to which each amino acid contributes to the opsin shift and to estimate the inter-residue synergistic effect. It was found that one of the most important contributions to the opsin shift is the electron transfer from Tyr185 to the chromophore upon excitation. We also indicate that some aromatic (Trp86, Trp182) and polar (Ser141, Thr142) residues, located in the vicinity of the retinal polyene chain and the β-ionone ring, respectively, play an important role in compensating for the large blue-shift induced by both the counterion residues (Asp85, Asp212) and an internal water molecule (W402) located near the Schiff base linkage. In particular, the effect of Trp86 is comparable to that of Tyr185. In addition, Ser141 and Thr142 were found to contribute to an increase in the dipole moment of bR in the excited state. Finally, we provide a complete energy diagram for the opsin shift together with the contribution of the chromophore-protein steric interaction.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Center for Biological Resources and Informatics, Tokyo Institute of Technology, B-62 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | - Azuma Matsuura
- Fujitsu Laboratories, Ltd., 10-1 Morinosato-Wakamiya, Atsugi 243-0197, Japan
| | - Hiroyuki Sato
- Fujitsu Laboratories, Ltd., 10-1 Morinosato-Wakamiya, Atsugi 243-0197, Japan
| | - Minoru Sakurai
- Center for Biological Resources and Informatics, Tokyo Institute of Technology, B-62 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| |
Collapse
|
34
|
Free Energy Profile of APOBEC3G Protein Calculated by a Molecular Dynamics Simulation. BIOLOGY 2012; 1:245-59. [PMID: 24832225 PMCID: PMC4009775 DOI: 10.3390/biology1020245] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/17/2012] [Revised: 07/04/2012] [Accepted: 07/05/2012] [Indexed: 11/25/2022]
Abstract
The human APOBEC3G protein (A3G) is a single-stranded DNA deaminase that inhibits the replication of retrotransposons and retroviruses, including HIV-1. Atomic details of A3G’s catalytic mechanism have started to emerge, as the structure of its catalytic domain (A3Gctd) has been revealed by NMR and X-ray crystallography. The NMR and crystal structures are similar overall; however, differences are apparent for β2 strand (β2) and loops close to the catalytic site. To add some insight into these differences and to better characterize A3Gctd dynamics, we calculated its free energy profile by using the Generalized-Born surface area (GBSA) method accompanied with a molecular dynamics simulation. The GBSA method yielded an enthalpy term for A3Gctd’s free energy, and we developed a new method that takes into account the distribution of the protein’s dihedral angles to calculate its entropy term. The structure solved by NMR was found to have a lower energy than that of the crystal structure, suggesting that this conformation is dominant in solution. In addition, β2-loop-β2’ configuration was stable throughout a 20-ns molecular dynamics (MD) simulation. This finding suggests that in solution A3Gctd is not likely to adopt the continuous β2 strand configuration present in the APOBEC2 crystal structure. In the NMR structure, the solvent water accessibility of the catalytic Zn2+ was limited throughout the 20-ns MD simulation. This result explains previous observations in which A3G did not bind or catalyze single cytosine nucleotide, even when at excessive concentrations.
Collapse
|
35
|
Abstract
A new peptidomimetic is proposed, resulting from substitution of the C═O carbonyl group by a B-F bond at the amide linkage. The effects of such chemical alteration are theoretically investigated through comparative calculations on dimethyl-fluoro-aminoborane H(3)C-BF-NH-CH(3) and N-methylacetamide H(3)C-CO-NH-CH(3), the simplest model of a peptide linkage. While little difference is found regarding size, electronic structure, and plaque rigidity, substantial distinctions are, however, observed between the polarities and association energies of the two compounds, with a B-F···H-N hydrogen bond estimated to be about one-third as strong as the natural C═O···H-N one. The conformational maps of the corresponding dipeptide models exhibit similarities and distinctions, which partly account for helical oligomer properties. Although capable of a high level of organization, the chains made of fluoro-aminoborane units show overall less structuration and more plasticity than their peptidic counterparts. Contrasts with fluorine-containing peptidomimetic 2-fluoro-2-butene are further underlined.
Collapse
Affiliation(s)
- Simon Mathieu
- Laboratoire de Chimie et Physique Quantiques (CNRS, UMR-5626), IRSAMC, Université Paul-Sabatier, 31062 Toulouse Cedex, France
| | | |
Collapse
|
36
|
Buczek A, Wałęsa R, Broda MA. β-turn tendency in N-methylated peptides with dehydrophenylalanine residue: DFT study. Biopolymers 2012; 97:518-28. [DOI: 10.1002/bip.22034] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2010] [Revised: 02/02/2011] [Accepted: 02/04/2011] [Indexed: 01/10/2023]
|
37
|
|
38
|
Mitchell W, Matsumoto S. Large-scale integrated super-computing platform for next generation virtual drug discovery. Curr Opin Chem Biol 2011; 15:553-9. [PMID: 21723773 PMCID: PMC7108376 DOI: 10.1016/j.cbpa.2011.06.005] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2011] [Revised: 06/05/2011] [Accepted: 06/07/2011] [Indexed: 11/26/2022]
Abstract
Traditional drug discovery starts by experimentally screening chemical libraries to find hit compounds that bind to protein targets, modulating their activity. Subsequent rounds of iterative chemical derivitization and rescreening are conducted to enhance the potency, selectivity, and pharmacological properties of hit compounds. Although computational docking of ligands to targets has been used to augment the empirical discovery process, its historical effectiveness has been limited because of the poor correlation of ligand dock scores and experimentally determined binding constants. Recent progress in super-computing, coupled to theoretical insights, allows the calculation of the Gibbs free energy, and therefore accurate binding constants, for usually large ligand–receptor systems. This advance extends the potential of virtual drug discovery. A specific embodiment of the technology, integrating de novo, abstract fragment based drug design, sophisticated molecular simulation, and the ability to calculate thermodynamic binding constants with unprecedented accuracy, are discussed.
Collapse
Affiliation(s)
- Wayne Mitchell
- Experimental Therapeutics Centre, Agency for Science and Technology Research (A*STAR), 31 Biopolis Way, #03-01 Nanos, Biopolis, Singapore 138669, Singapore
| | | |
Collapse
|
39
|
Hioe J, Savasci G, Brand H, Zipse H. The stability of Cα peptide radicals: why glycyl radical enzymes? Chemistry 2011; 17:3781-9. [PMID: 21341321 DOI: 10.1002/chem.201002620] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2010] [Indexed: 11/06/2022]
Abstract
The conformational space of dipeptide models derived from glycine, alanine, phenylalanine, proline, tyrosine, and cysteine has been searched extensively and compared with the corresponding C(α) dipeptide radicals at the G3(MP2)-RAD level of theory. The results indicate that the (least-substituted) glycine dipeptide radical is the thermochemically most stable of these species. Analysis of the structural parameters indicates that this is due to repulsive interactions between the C(α) substituents and peptide units in the radical. A comparison of the conformational preferences of dipeptide radicals and their closed-shell parents also indicates that radical stability is a strongly conformation-dependent property.
Collapse
Affiliation(s)
- Johnny Hioe
- Department of Chemistry, LMU München, Butenandtstrasse 5-13, 81377 München, Germany
| | | | | | | |
Collapse
|
40
|
The DF-LCCSD(T0) correction of the φ/ψ force field dihedral parameters significantly influences the free energy profile of the alanine dipeptide. Chem Phys Lett 2011. [DOI: 10.1016/j.cplett.2011.01.030] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
41
|
Jiang F, Han W, Wu YD. Influence of side chain conformations on local conformational features of amino acids and implication for force field development. J Phys Chem B 2010; 114:5840-50. [PMID: 20392111 DOI: 10.1021/jp909088e] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Statistical analysis of coil regions in protein structures has been used to obtain the local backbone phi, psi preferences of amino acids, which agree well with the NMR experiments of unfolded peptides and proteins. We analyzed the conformational features of amino acid residues in a restricted coil library of 4220 high-resolution protein crystal structures. In addition to Gly, Ala, and Pro, the phi, psi distribution (Ramachandran plot) of each amino acid is analyzed with respect to three side chain conformers: g+ (chi(1) approximately -60 degrees), g- (chi(1) approximately 60 degrees), and t (chi(1) approximately 180 degrees). The statistical study indicates that the effect of side chain conformations on phi, psi distributions is even greater than the effect of amino acid types. On the basis of the chi(1), phi, psi conformational preferences, the amino acids in addition to Gly, Pro, and Ala can be divided into five types: (1) ordinary amino acids, (2) Ser, (3) Asp and Asn, (4) Val and Ile, and (5) Thr, each with distinguished chi(1) rotamers. The alpha-helix, beta-sheet, and type-I beta-turn preferences of the different rotamers of various amino acid types can be captured by their intrinsic phi, psi preferences from our coil library. Molecular dynamics simulations of dipeptide Ac-X-NHMe and tetrapeptide Ac-A-X-A-NHMe models give nearly the same side chain rotamer distributions. However, for many amino acids, both OPLS-AA/L and AMBER-FF03 force fields give very different chi(1) rotamer distributions from the coil library. This may partially explain why dipeptide models sometimes cannot reproduce those of protein structures well. The current coil library analysis may be valuable in improving the force field for protein simulations.
Collapse
Affiliation(s)
- Fan Jiang
- Laboratory of Chemical Genomics, Shenzhen Graduate School of Peking University, Shenzhen 518055, China
| | | | | |
Collapse
|
42
|
Kang YK, Byun BJ. Assessment of density functionals with long-range and/or empirical dispersion corrections for conformational energy calculations of peptides. J Comput Chem 2010; 31:2915-23. [DOI: 10.1002/jcc.21587] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
43
|
Docking study and binding free energy calculation of poly (ADP-ribose) polymerase inhibitors. J Mol Model 2010; 17:383-9. [DOI: 10.1007/s00894-010-0728-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2010] [Accepted: 04/21/2010] [Indexed: 12/15/2022]
|
44
|
Matsuura A, Hayashi T, Sato H, Takahashi A, Sakurai M. Theoretical study on the absorption maxima of real GFPs. Chem Phys Lett 2010. [DOI: 10.1016/j.cplett.2009.11.074] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|