1
|
Maruyama Y, Yoshida N. RISMiCal: A software package to perform fast RISM/3D-RISM calculations. J Comput Chem 2024; 45:1470-1482. [PMID: 38472097 DOI: 10.1002/jcc.27340] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2023] [Revised: 02/23/2024] [Accepted: 02/28/2024] [Indexed: 03/14/2024]
Abstract
Solvent plays an essential role in a variety of chemical, physical, and biological processes that occur in the solution phase. The reference interaction site model (RISM) and its three-dimensional extension (3D-RISM) serve as powerful computational tools for modeling solvation effects in chemical reactions, biological functions, and structure formations. We present the RISM integrated calculator (RISMiCal) program package, which is based on RISM and 3D-RISM theories with fast GPU code. RISMiCal has been developed as an integrated RISM/3D-RISM program that has interfaces with external programs such as Gaussian16, GAMESS, and Tinker. Fast 3D-RISM programs for single- and multi-GPU codes written in CUDA would enhance the availability of these hybrid methods because they require the performance of many computationally expensive 3D-RISM calculations. We expect that our package can be widely applied for chemical and biological processes in solvent. The RISMiCal package is available at https://rismical-dev.github.io.
Collapse
Affiliation(s)
- Yutaka Maruyama
- Data Science Center for Creative Design and Manufacturing, The Institute of Statistical Mathematics, Tachikawa, Tokyo, Japan
- Department of Physics, School of Science and Technology, Meiji University, Kawasaki-shi, Kanagawa, Japan
| | - Norio Yoshida
- Graduate School of Informatics, Nagoya University, Chikusa, Nagoya, Japan
| |
Collapse
|
2
|
Fedotova MV, Chuev GN. The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands. Int J Mol Sci 2024; 25:3697. [PMID: 38612508 PMCID: PMC11011341 DOI: 10.3390/ijms25073697] [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: 03/09/2024] [Revised: 03/22/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.
Collapse
Affiliation(s)
- Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., 142290 Pushchino, Russia
| |
Collapse
|
3
|
Sk MF, Samanta S, Poddar S, Kar P. Microsecond dynamics of H10N7 influenza neuraminidase reveals the plasticity of loop regions and drug resistance due to the R292K mutation. J Comput Chem 2024; 45:247-263. [PMID: 37787086 DOI: 10.1002/jcc.27234] [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: 03/22/2023] [Revised: 08/12/2023] [Accepted: 09/14/2023] [Indexed: 10/04/2023]
Abstract
At the beginning of the last century, multiple pandemics caused by influenza (flu) viruses severely impacted public health. Despite the development of vaccinations and antiviral medications to prevent and control impending flu outbreaks, unforeseen novel strains and continuously evolving old strains continue to represent a serious threat to human life. Therefore, the recently identified H10N7, for which not much data is available for rational structure-based drug design, needs to be further explored. Here, we investigated the structural dynamics of neuraminidase N7 upon binding of inhibitors, and the drug resistance mechanisms against the oseltamivir (OTV) and laninamivir (LNV) antivirals due to the crucial R292K mutation on the N7 using the computational microscope, molecular dynamics (MD) simulations. In this study, each system underwent long 2 × 1 μs MD simulations to answer the conformational changes and drug resistance mechanisms. These long time-scale dynamics simulations and free energy landscapes demonstrated that the mutant systems showed a high degree of conformational variation compared to their wildtype (WT) counterparts, and the LNV-bound mutant exhibited an extended 150-loop conformation. Further, the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) calculation and MM/GBSA free energy decomposition were used to characterize the binding of OTV and LNV with WT, and R292K mutated N7, revealing the R292K mutation as drug-resistant, facilitated by a decline in binding interaction and a reduction in the dehydration penalty. Due to the broader binding pocket cavity of the smaller K292 mutant residue relative to the wildtype, the drug carboxylate to K292 hydrogen bonding was lost, and the area surrounding the K292 residue was more accessible to water molecules. This implies that drug resistance could be reduced by strengthening the hydrogen bond contacts between N7 inhibitors and altered N7, creating inhibitors that can form a hydrogen bond to the mutant K292, or preserving the closed cavity conformations.
Collapse
Affiliation(s)
- Md Fulbabu Sk
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
- Theoretical and Computational Biophysics Group, NIH Resource for Macromolecular Modeling and Visualization, Beckman Institute for Advanced Science and Technology, University of Illinois Urbana-Champaign, Urbana, Illinois, USA
| | - Sunanda Samanta
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - Sayan Poddar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Indore, India
| |
Collapse
|
4
|
Nguyen DV, Hengphasatporn K, Danova A, Suroengrit A, Boonyasuppayakorn S, Fujiki R, Shigeta Y, Rungrotmongkol T, Chavasiri W. Structure-yeast α-glucosidase inhibitory activity relationship of 9-O-berberrubine carboxylates. Sci Rep 2023; 13:18865. [PMID: 37914757 PMCID: PMC10620162 DOI: 10.1038/s41598-023-45116-0] [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: 06/24/2023] [Accepted: 10/16/2023] [Indexed: 11/03/2023] Open
Abstract
Thirty-five 9-O-berberrubine carboxylate derivatives were synthesized and evaluated for yeast α-glucosidase inhibitory activity. All compounds demonstrated better inhibitory activities than the parent compounds berberine (BBR) and berberrubine (BBRB), and a positive control, acarbose. The structure-activity correlation study indicated that most of the substituents on the benzoate moiety such as methoxy, hydroxy, methylenedioxy, benzyloxy, halogen, trifluoromethyl, nitro and alkyl can contribute to the activities except multi-methoxy, fluoro and cyano. In addition, replacing benzoate with naphthoate, cinnamate, piperate or diphenylacetate also led to an increase in inhibitory activities except with phenyl acetate. 9, 26, 27, 28 and 33 exhibited the most potent α-glucosidase inhibitory activities with the IC50 values in the range of 1.61-2.67 μM. Kinetic study revealed that 9, 26, 28 and 33 interacted with the enzyme via competitive mode. These four compounds were also proved to be not cytotoxic at their IC50 values. The competitive inhibition mechanism of these four compounds against yeast α-glucosidase was investigated using molecular docking and molecular dynamics simulations. The binding free energy calculations suggest that 26 exhibited the strongest binding affinity, and its binding stability is supported by hydrophobic interactions with D68, F157, F158 and F177. Therefore, 9, 26, 28 and 33 would be promising candidates for further studies of antidiabetic activity.
Collapse
Affiliation(s)
- Duy Vu Nguyen
- Department of Chemistry, Faculty of Science, Center of Excellence in Natural Products Chemistry, Chulalongkorn University, Pathumwan, Bangkok, 10330, Thailand
| | - Kowit Hengphasatporn
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan
| | - Ade Danova
- Department of Chemistry, Faculty of Science, Center of Excellence in Natural Products Chemistry, Chulalongkorn University, Pathumwan, Bangkok, 10330, Thailand
- Organic Chemistry Division, Department of Chemistry, Faculty of Mathematics and Natural Sciences, Institut Teknologi Bandung, Bandung, West Java, 40132, Indonesia
| | - Aphinya Suroengrit
- Department of Microbiology, Faculty of Medicine, Center of Excellence in Applied Medical Virology, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Siwaporn Boonyasuppayakorn
- Department of Microbiology, Faculty of Medicine, Center of Excellence in Applied Medical Virology, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Ryo Fujiki
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki, 305-8577, Japan
| | - Thanyada Rungrotmongkol
- Bioinformatics and Computational Biology Program, Graduated School, Chulalongkorn University, Bangkok, 10330, Thailand
- Department of Biochemistry, Faculty of Science, Center of Excellence in Biocatalyst and Sustainable Biotechnology, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Warinthorn Chavasiri
- Department of Chemistry, Faculty of Science, Center of Excellence in Natural Products Chemistry, Chulalongkorn University, Pathumwan, Bangkok, 10330, Thailand.
| |
Collapse
|
5
|
Hirata F. Structural Fluctuation, Relaxation, and Folding of Protein: An Approach Based on the Combined Generalized Langevin and RISM/3D-RISM Theories. Molecules 2023; 28:7351. [PMID: 37959769 PMCID: PMC10647392 DOI: 10.3390/molecules28217351] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2023] [Revised: 09/25/2023] [Accepted: 09/28/2023] [Indexed: 11/15/2023] Open
Abstract
In 2012, Kim and Hirata derived two generalized Langevin equations (GLEs) for a biomolecule in water, one for the structural fluctuation of the biomolecule and the other for the density fluctuation of water, by projecting all the mechanical variables in phase space onto the two dynamic variables: the structural fluctuation defined by the displacement of atoms from their equilibrium positions, and the solvent density fluctuation. The equation has an expression similar to the classical Langevin equation (CLE) for a harmonic oscillator, possessing terms corresponding to the restoring force proportional to the structural fluctuation, as well as the frictional and random forces. However, there is a distinct difference between the two expressions that touches on the essential physics of the structural fluctuation, that is, the force constant, or Hessian, in the restoring force. In the CLE, this is given by the second derivative of the potential energy among atoms in a protein. So, the quadratic nature or the harmonicity is only valid at the minimum of the potential surface. On the contrary, the linearity of the restoring force in the GLE originates from the projection of the water's degrees of freedom onto the protein's degrees of freedom. Taking this into consideration, Kim and Hirata proposed an ansatz for the Hessian matrix. The ansatz is used to equate the Hessian matrix with the second derivative of the free-energy surface or the potential of the mean force of a protein in water, defined by the sum of the potential energy among atoms in a protein and the solvation free energy. Since the free energy can be calculated from the molecular mechanics and the RISM/3D-RISM theory, one can perform an analysis similar to the normal mode analysis (NMA) just by diagonalizing the Hessian matrix of the free energy. This method is referred to as the Generalized Langevin Mode Analysis (GLMA). This theory may be realized to explore a variety of biophysical processes, including protein folding, spectroscopy, and chemical reactions. The present article is devoted to reviewing the development of this theory, and to providing perspective in exploring life phenomena.
Collapse
Affiliation(s)
- Fumio Hirata
- Institute for Molecular Science, National Institute of Natural Sciences, Okazaki 444-8585, Japan
| |
Collapse
|
6
|
Hirata F. Generalized Langevin Mode Analysis (GLMA) for Local Density Fluctuation of Water in an Inhomogeneous Field of a Biomolecule. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023]
|
7
|
Osaki K, Ekimoto T, Yamane T, Ikeguchi M. 3D-RISM-AI: A Machine Learning Approach to Predict Protein-Ligand Binding Affinity Using 3D-RISM. J Phys Chem B 2022; 126:6148-6158. [PMID: 35969673 PMCID: PMC9421647 DOI: 10.1021/acs.jpcb.2c03384] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2022] [Revised: 07/27/2022] [Indexed: 11/30/2022]
Abstract
Hydration free energy (HFE) is a key factor in improving protein-ligand binding free energy (BFE) prediction accuracy. The HFE itself can be calculated using the three-dimensional reference interaction model (3D-RISM); however, the BFE predictions solely evaluated using 3D-RISM are not correlated to the experimental BFE for abundant protein-ligand pairs. In this study, to predict the BFE for multiple sets of protein-ligand pairs, we propose a machine learning approach incorporating the HFEs obtained using 3D-RISM, termed 3D-RISM-AI. In the learning process, structural metrics, intra-/intermolecular energies, and HFEs obtained via 3D-RISM of ∼4000 complexes in the PDBbind database (ver. 2018) were used. The BFEs predicted using 3D-RISM-AI were well correlated to the experimental data (Pearson's correlation coefficient of 0.80 and root-mean-square error of 1.91 kcal/mol). As important factors for the prediction, the difference in the solvent accessible surface area between the bound and unbound structures and the hydration properties of the ligands were detected during the learning process.
Collapse
Affiliation(s)
- Kazu Osaki
- Graduate
School of Medical Life Science, Yokohama
City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Toru Ekimoto
- Graduate
School of Medical Life Science, Yokohama
City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Tsutomu Yamane
- Center
for Computational Science, RIKEN, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| | - Mitsunori Ikeguchi
- Graduate
School of Medical Life Science, Yokohama
City University, 1-7-29 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
- Center
for Computational Science, RIKEN, 1-7-22 Suehiro-cho, Tsurumi-ku, Yokohama 230-0045, Japan
| |
Collapse
|
8
|
Sugita M, Onishi I, Irisa M, Yoshida N, Hirata F. Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory. Molecules 2021; 26:E271. [PMID: 33430461 PMCID: PMC7826681 DOI: 10.3390/molecules26020271] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2020] [Revised: 12/23/2020] [Accepted: 12/28/2020] [Indexed: 11/18/2022] Open
Abstract
There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.
Collapse
Affiliation(s)
- Masatake Sugita
- Department of Computer Science, School of Computing, Tokyo Institute of Technology, W8-76, 2-12-1, Ookayama Meguro-ku, Tokyo 152-8550, Japan;
| | - Itaru Onishi
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka 820-8502, Japan; (I.O.); (M.I.)
| | - Masayuki Irisa
- Department of Bioscience and Bioinformatics, Kyushu Institute of Technology, Iizuka, Fukuoka 820-8502, Japan; (I.O.); (M.I.)
| | - Norio Yoshida
- Department of Chemistry, Kyushu University, Fukuoka, Fukuoka 812-8581, Japan;
| | - Fumio Hirata
- Theoretical and Computational Molecular Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
9
|
Verma K, Mahalapbutr P, Auepattanapong A, Khaikate O, Kuhakarn C, Takahashi K, Rungrotmongkol T. Molecular dynamics simulations of sulfone derivatives in complex with DNA topoisomerase IIα ATPase domain. J Biomol Struct Dyn 2020; 40:1692-1701. [PMID: 33089727 DOI: 10.1080/07391102.2020.1831961] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Human topoisomerase II alpha (TopoIIα) is a crucial enzyme involved in maintaining genomic integrity during the process of DNA replication and mitotic division. It is a vital therapeutic target for designing novel anticancer agents in targeted cancer therapy. Sulfones, members of organosulfur compounds, have been reported to possess various biological activities such as antimicrobial, anti-inflammatory, anti-HIV, anticancer, and antimalarial properties. In the present study, a series of sulfones was selected to evaluate their inhibitory activity against TopoIIα using computational approaches. Molecular docking results revealed that several sulfone analogs bind efficiently to the ATPase domain of TopoIIα. Among them, sulfones 18a, 60a, *4 b, *8 b, *3c, and 8c exhibit higher binding affinity than the known TopoII inhibitor, salvicine. Molecular dynamics simulations and free energy calculations based on MM/PB(GB)SA method demonstrated that sulfone *8 b strongly interacts with amino acid residues in the ATP-binding pocket (E87, N91, D94, I125, I141, F142, S149, G161, and A167), driven mainly by an electrostatic attraction and a strong H-bond formation at G161 residue. Altogether, the obtained results predicted that sulfones could have a high potential to be a lead molecule for targeting TopoIIα.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Kanika Verma
- Biocatalyst and Environmental Biotechnology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Panupong Mahalapbutr
- Department of Biochemistry, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
| | - Atima Auepattanapong
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok, Thailand
| | - Onnicha Khaikate
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok, Thailand
| | - Chutima Kuhakarn
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok, Thailand
| | - Kaito Takahashi
- Institute of Atomic and Sciences, Academia Sinica, Taipei, Taiwan
| | - Thanyada Rungrotmongkol
- Biocatalyst and Environmental Biotechnology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand.,Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
| |
Collapse
|
10
|
Nutho B, Pengthaisong S, Tankrathok A, Lee VS, Ketudat Cairns JR, Rungrotmongkol T, Hannongbua S. Structural Basis of Specific Glucoimidazole and Mannoimidazole Binding by Os3BGlu7. Biomolecules 2020; 10:biom10060907. [PMID: 32549280 PMCID: PMC7356692 DOI: 10.3390/biom10060907] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Revised: 06/09/2020] [Accepted: 06/11/2020] [Indexed: 01/07/2023] Open
Abstract
β-Glucosidases and β-mannosidases hydrolyze substrates that differ only in the epimer of the nonreducing terminal sugar moiety, but most such enzymes show a strong preference for one activity or the other. Rice Os3BGlu7 and Os7BGlu26 β-glycosidases show a less strong preference, but Os3BGlu7 and Os7BGlu26 prefer glucosides and mannosides, respectively. Previous studies of crystal structures with glucoimidazole (GIm) and mannoimidazole (MIm) complexes and metadynamic simulations suggested that Os7BGlu26 hydrolyzes mannosides via the B2,5 transition state (TS) conformation preferred for mannosides and glucosides via their preferred 4H3/4E TS conformation. However, MIm is weakly bound by both enzymes. In the present study, we found that MIm was not bound in the active site of crystallized Os3BGlu7, but GIm was tightly bound in the -1 subsite in a 4H3/4E conformation via hydrogen bonds with the surrounding residues. One-microsecond molecular dynamics simulations showed that GIm was stably bound in the Os3BGlu7 active site and the glycone-binding site with little distortion. In contrast, MIm initialized in the B2,5 conformation rapidly relaxed to a E3/4H3 conformation and moved out into a position in the entrance of the active site, where it bound more stably despite making fewer interactions. The lack of MIm binding in the glycone site in protein crystals and simulations implies that the energy required to distort MIm to the B2,5 conformation for optimal active site residue interactions is sufficient to offset the energy of those interactions in Os3BGlu7. This balance between distortion and binding energy may also provide a rationale for glucosidase versus mannosidase specificity in plant β-glycosidases.
Collapse
Affiliation(s)
- Bodee Nutho
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand;
| | - Salila Pengthaisong
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand; (S.P.); (A.T.)
- Center for Biomolecular Structure, Function and Application, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
| | - Anupong Tankrathok
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand; (S.P.); (A.T.)
- Center for Biomolecular Structure, Function and Application, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
| | - Vannajan Sanghiran Lee
- Department of Chemistry, Faculty of Science, University of Malaya, Kuala Lumpur 50603, Malaysia;
| | - James R. Ketudat Cairns
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand; (S.P.); (A.T.)
- Center for Biomolecular Structure, Function and Application, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand
- Correspondence: (J.R.K.C.); (T.R.); (S.H.); Tel.: +66-4422-4304 (J.R.K.C.); +66-2218-5426 (T.R.); +66-2218-7602 (S.H.)
| | - Thanyada Rungrotmongkol
- Biocatalyst and Environmental Biotechnology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok 10330, Thailand
- Correspondence: (J.R.K.C.); (T.R.); (S.H.); Tel.: +66-4422-4304 (J.R.K.C.); +66-2218-5426 (T.R.); +66-2218-7602 (S.H.)
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand;
- Correspondence: (J.R.K.C.); (T.R.); (S.H.); Tel.: +66-4422-4304 (J.R.K.C.); +66-2218-5426 (T.R.); +66-2218-7602 (S.H.)
| |
Collapse
|
11
|
Kato Y, Takahashi K, Ito F, Suzuki S, Fukui K, Mimaki M, Suzuki K. Novel oseltamivir-resistant mutations distant from the active site of influenza B neuraminidase. J Biomol Struct Dyn 2020; 39:3491-3500. [PMID: 32406800 DOI: 10.1080/07391102.2020.1765872] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We performed a neuraminidase sequence analysis of thirty-two pediatric patients with influenza B who visited Teikyo University Hospital from January 2016 to March 2017, and found oseltamivir-resistant samples belonging to the Yamagata and Victoria lineages. Comparison with the neuraminidase sequence of oseltamivir-susceptible B/Brisbane/60/2008 revealed 5 common amino acid substitutions in many of these samples. According to the binding free energy calculation, the N340D and E358K substitutions reduced the affinity of oseltamivir to neuraminidase. Unexpectedly, these substitutions were located distant from the oseltamivir-binding site in neuraminidase. According to the molecular dynamics simulations, the N340D substitution rearranged complicated hydrogen bond networks in an extensive surface region of neuraminidase. The E358K substitution extensively altered the electrostatic potential map of the overall neuraminidase structure. Through these novel mechanisms, the N340D and E358K substitutions indirectly influenced the affinity reduction. These results may be useful for designing drugs for the treatment of oseltamivir-resistant virus infections.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Yusuke Kato
- Institute for Enzyme Research, Tokushima University, Tokushima, Japan.,Asia International Institute of Infectious Disease Control, Teikyo University, Tokyo, Japan
| | - Kazuhiro Takahashi
- Asia International Institute of Infectious Disease Control, Teikyo University, Tokyo, Japan.,Department of Pediatrics, Teikyo University School of Medicine, Tokyo, Japan
| | - Fuyu Ito
- Asia International Institute of Infectious Disease Control, Teikyo University, Tokyo, Japan
| | - Shoichi Suzuki
- Asia International Institute of Infectious Disease Control, Teikyo University, Tokyo, Japan
| | - Kiyoshi Fukui
- Institute for Enzyme Research, Tokushima University, Tokushima, Japan
| | - Masakazu Mimaki
- Department of Pediatrics, Teikyo University School of Medicine, Tokyo, Japan
| | - Kazuo Suzuki
- Asia International Institute of Infectious Disease Control, Teikyo University, Tokyo, Japan
| |
Collapse
|
12
|
Sanachai K, Mahalapbutr P, Choowongkomon K, Poo-arporn RP, Wolschann P, Rungrotmongkol T. Insights into the Binding Recognition and Susceptibility of Tofacitinib toward Janus Kinases. ACS OMEGA 2020; 5:369-377. [PMID: 31956784 PMCID: PMC6964278 DOI: 10.1021/acsomega.9b02800] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/29/2019] [Accepted: 12/04/2019] [Indexed: 05/06/2023]
Abstract
Janus kinases (JAKs) are enzymes involved in signaling pathways that affect hematopoiesis and immune cell functions. JAK1, JAK2, and JAK3 play different roles in numerous diseases of the immune system and have also been considered as potential targets for cancer therapy. In the present study, the susceptibility of the oral JAK inhibitor tofacitinib against these three JAKs was elucidated using the 500-ns molecular dynamics (MD) simulations and free energy calculations based on MM-PB(GB)SA, QM/MM-GBSA (PM3 and SCC-DFTB), and SIE methods. The obtained results revealed that tofacitinib could interact with all JAKs at the ATP-binding site via electrostatic attraction, hydrogen bond formation, and in particular van der Waals interaction. The conserved glutamate and leucine residues (E957 and L959 of JAK1, E930 and L932 of JAK2, and E903 and L905 of JAK3) located in the hinge region stabilized tofacitinib binding through strongly formed hydrogen bonds. Complexation with the incoming tofacitinib led to a closed conformation of the ATP-binding site and a decreased protein fluctuation at the glycine loop of the JAK protein. The binding affinities of tofacitinib/JAKs were ranked in the order of JAK3 > JAK2 ∼ JAK1, which are in line with the reported experimental data.
Collapse
Affiliation(s)
- Kamonpan Sanachai
- Structural
and Computational Biology Research Unit, Department of
Biochemistry, Faculty of Science and Program in Bioinformatics and Computational
Biology, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
| | - Panupong Mahalapbutr
- Structural
and Computational Biology Research Unit, Department of
Biochemistry, Faculty of Science and Program in Bioinformatics and Computational
Biology, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
| | - Kiattawee Choowongkomon
- Department
of Biochemistry, Faculty of Science, Kasetsart
University, Bangkok 10900, Thailand
| | - Rungtiva P. Poo-arporn
- Biological
Engineering Program, Faculty of Engineering, King Mongkut’s University of Technology Thonburi, Bangkok 10140, Thailand
| | - Peter Wolschann
- Department of Pharmaceutical Chemistry,
Faculty of Life Sciences and Institute of Theoretical
Chemistry, University of Vienna, Vienna 1090, Austria
| | - Thanyada Rungrotmongkol
- Structural
and Computational Biology Research Unit, Department of
Biochemistry, Faculty of Science and Program in Bioinformatics and Computational
Biology, Faculty of Science, Chulalongkorn
University, Bangkok 10330, Thailand
- E-mail: , . Tel: +66 2 2185426. Fax: +66 22185418
| |
Collapse
|
13
|
Source of oseltamivir resistance due to single E119D and double E119D/H274Y mutations in pdm09H1N1 influenza neuraminidase. J Comput Aided Mol Des 2019; 34:27-37. [PMID: 31773463 DOI: 10.1007/s10822-019-00251-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2019] [Accepted: 11/09/2019] [Indexed: 12/24/2022]
Abstract
Influenza epidemics are responsible for an average of 3-5 millions of severe cases and up to 500,000 deaths around the world. One of flu pandemic types is influenza A(H1N1)pdm09 virus (pdm09H1N1). Oseltamivir is the antiviral drug used to treat influenza targeting at neuraminidase (NA) located on the viral surface. Influenza virus undergoes high mutation rates and leads to drug resistance, and thus the development of more efficient drugs is required. In the present study, all-atom molecular dynamics simulations were applied to understand the oseltamivir resistance caused by the single E119D and double E119D/H274Y mutations on NA. The obtained results in terms of binding free energy and intermolecular interactions in the ligand-protein interface showed that the oseltamivir could not be well accommodated in the binding pocket of both NA mutants and the 150-loop moves out from oseltamivir as an "open" state. A greater number of water molecules accessible to the binding pocket could disrupt the oseltamivir binding with NA target as seen be high mobility of oseltamivir at the active site. Additionally, our finding could guide to the design and development of novel NA inhibitor drugs.
Collapse
|
14
|
First principle determination of dielectric constant of polar liquids by the extended RISM theory with the bridge diagrams concerning dipole-dipole interactions. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111567] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
15
|
Cao S, Konovalov KA, Unarta IC, Huang X. Recent Developments in Integral Equation Theory for Solvation to Treat Density Inhomogeneity at Solute–Solvent Interface. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201900049] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Siqin Cao
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
| | - Kirill A. Konovalov
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
| | - Ilona Christy Unarta
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
- Bioengineering Graduate Programthe Hong Kong University of Science and TechnologyHong Kong of Chinese National EngineeringResearch Center for Tissue Restoration and Reconstructionthe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
| | - Xuhui Huang
- Department of Chemistrythe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- Center of System Biology and Human HealthState Key Laboratory of Molecular Neuroscience, Hong Kong Branch Clear Water Bay Kowloon Hong Kong
- Bioengineering Graduate Programthe Hong Kong University of Science and TechnologyHong Kong of Chinese National EngineeringResearch Center for Tissue Restoration and Reconstructionthe Hong Kong University of Science and Technology Clear Water Bay Kowloon Hong Kong
- HKUST‐Shenzhen Research Institute Hi‐Tech Park, Nanshan Shenzhen 518057 China
| |
Collapse
|
16
|
Kongkaew S, Rungrotmongkol T, Punwong C, Noguchi H, Takeuchi F, Kungwan N, Wolschann P, Hannongbua S. Interactions of HLA-DR and Topoisomerase I Epitope Modulated Genetic Risk for Systemic Sclerosis. Sci Rep 2019; 9:745. [PMID: 30679605 PMCID: PMC6345791 DOI: 10.1038/s41598-018-37038-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2018] [Accepted: 11/29/2018] [Indexed: 12/11/2022] Open
Abstract
The association of systemic sclerosis with anti-Topoisomerase 1 antibody (ATASSc) with specific alleles of human leukocyte antigen (HLA)-DR has been observed among various ethnics. The anti-Topoisomerase 1 antibody is a common autoantibody in SSc with diffuse cutaneous scleroderma, which is one of the clinical subtypes of SSc. On the other hand, an immunodominant peptide of topoisomerase 1 (Top1) self-protein (residues 349-368) was reported to have strong association with ATASSc. In this study, molecular dynamics simulation was performed on the complexes of Top1 peptide with various HLA-DR subtypes divided into ATASSc-associated alleles (HLA-DRB1*08:02, HLA-DRB1*11:01 and HLA-DRB1*11:04), suspected allele (HLA-DRB5*01:02), and non-associated allele (HLA-DRB1*01:01). The unique interaction for each system was compared to the others in terms of dynamical behaviors, binding free energies and solvation effects. Our results showed that three HLA-DR/Top1 complexes of ATASSc association mostly exhibited high protein stability and increased binding efficiency without solvent interruption, in contrast to non-association. The suspected case (HLA-DRB5*01:02) binds Top1 as strongly as the ATASSc association case, which implied a highly possible risk for ATASSc development. This finding might support ATASSc development mechanism leading to a guideline for the treatment and avoidance of pathogens like Top1 self-peptide risk for ATASSc.
Collapse
Affiliation(s)
- Sirilak Kongkaew
- Program in Biotechnology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.,The Center of Excellence in Computational Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Thanyada Rungrotmongkol
- Biocatalyst and Environmental Biotechnology Research unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand. .,Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Chutintorn Punwong
- Department of Physics, Faculty of Science, Prince of Songkla University, Hat Yai, Songkhla, 90110, Thailand
| | - Hiroshi Noguchi
- School of Pharmacy, Nihon Pharmaceutical University, Saitama, 361-0806, Japan.,School of Pharmaceutical Sciences, University of Shizuoka, Shizuoka, 422-8526, Japan
| | - Fujio Takeuchi
- School of Pharmaceutical Sciences, University of Shizuoka, Shizuoka, 422-8526, Japan.,Faculty of Health and Nutrition, Tokyo Seiei University, Tokyo, 124-8530, Japan
| | - Nawee Kungwan
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand.,Center of Excellence in Materials Science and Technology, Chiang Mai University, Chiang Mai, 50200, Thailand
| | - Peter Wolschann
- The Center of Excellence in Computational Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.,Department of Pharmaceutical Chemistry, University of Vienna, Vienna, 1090, Austria.,Institute of Theoretical Chemistry, University of Vienna, Vienna, 1090, Austria
| | - Supot Hannongbua
- The Center of Excellence in Computational Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| |
Collapse
|
17
|
Glycan binding and specificity of viral influenza neuraminidases by classical molecular dynamics and replica exchange molecular dynamics simulations. J Biomol Struct Dyn 2018; 37:3354-3365. [DOI: 10.1080/07391102.2018.1514326] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|
18
|
Panman W, Nutho B, Chamni S, Dokmaisrijan S, Kungwan N, Rungrotmongkol T. Computational screening of fatty acid synthase inhibitors against thioesterase domain. J Biomol Struct Dyn 2018; 36:4114-4125. [PMID: 29161996 DOI: 10.1080/07391102.2017.1408496] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2017] [Accepted: 11/10/2017] [Indexed: 12/17/2022]
Abstract
Thioesterase (TE) domain of fatty acid synthase (FAS) is an attractive therapeutic target for design and development of anticancer drugs. In this present work, we search for the potential FAS inhibitors of TE domain from the ZINC database based on similarity search using three natural compounds as templates, including flavonoids, terpenoids, and phenylpropanoids. Molecular docking was used to predict the interaction energy of each screened ligand compared to the reference compound, which is methyl γ-linolenylfluorophosphonate (MGLFP). Based on this computational technique, rosmarinic acid and its eight analogs were observed as a new series of potential FAS inhibitors, which showed a stronger binding affinity than MGLFP. Afterward, nine docked complexes were studied by molecular dynamics simulations for investigating protein-ligand interactions and binding free energies using MM-PB(GB)SA, MM-3DRISM-KH, and QM/MM-GBSA methods. The binding free energy calculation indicated that the ZINC85948835 (R34) displayed the strongest binding efficiency against the TE domain of FAS. There are eight residues (S2308, I2250, E2251, Y2347, Y2351, F2370, L2427, and E2431) mainly contributed for the R34 binding. Moreover, R34 could directly form hydrogen bonds with S2308, which is one of the catalytic triad of TE domain. Therefore, our finding suggested that R34 could be a potential candidate as a novel FAS-TE inhibitor for further drug design.
Collapse
Affiliation(s)
- Wanwisa Panman
- a Multidisciplinary Program of Petrochemistry and Polymer Science, Faculty of Science , Chulalongkorn University , Bangkok 10330 , Thailand
| | - Bodee Nutho
- b Program in Biotechnology, Faculty of Science , Chulalongkorn University , Bangkok 10330 , Thailand
| | - Supakarn Chamni
- c Faculty of Pharmaceutical Sciences, Department of Pharmacognosy and Pharmaceutical Botany , Chulalongkorn University , Bangkok 10330 , Thailand
| | - Supaporn Dokmaisrijan
- d Division of Chemistry, School of Science , Walailak University , Nakon Si Thammarat 80161 , Thailand
| | - Nawee Kungwan
- e Faculty of Science, Department of Chemistry , Chiang Mai University , 239 Huay Kaew Road, Muang District, Chiang Mai 50200 , Thailand
- f Research Center on Chemistry for Development of Health Promoting Products from Northern Resources , Chiang Mai University , Chiang Mai , 50200 , Thailand
| | - Thanyada Rungrotmongkol
- g Faculty of Science, Structural and Computational Biology Research Group, Department of Biochemistry , Chulalongkorn University , Bangkok 10330 , Thailand
- h Faculty of Science, Program in Bioinformatics and Computational Biology , Chulalongkorn University , Bangkok 10330 , Thailand
| |
Collapse
|
19
|
Matter H, Güssregen S. Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands. Bioorg Med Chem Lett 2018; 28:2343-2352. [PMID: 29880400 DOI: 10.1016/j.bmcl.2018.05.061] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2018] [Revised: 05/28/2018] [Accepted: 05/30/2018] [Indexed: 11/18/2022]
Abstract
Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.
Collapse
Affiliation(s)
- Hans Matter
- Sanofi-Aventis Deutschland GmbH, Integrated Drug Discovery (IDD), Synthetic Molecular Design, Building G838, Industriepark Höchst, 65926 Frankfurt am Main, Germany.
| | - Stefan Güssregen
- Sanofi-Aventis Deutschland GmbH, Integrated Drug Discovery (IDD), Synthetic Molecular Design, Building G838, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| |
Collapse
|
20
|
Hirata F, Sugita M, Yoshida M, Akasaka K. Perspective: Structural fluctuation of protein and Anfinsen's thermodynamic hypothesis. J Chem Phys 2018; 148:020901. [PMID: 29331129 DOI: 10.1063/1.5013104] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
The thermodynamics hypothesis, casually referred to as "Anfinsen's dogma," is described theoretically in terms of a concept of the structural fluctuation of protein or the first moment (average structure) and the second moment (variance and covariance) of the structural distribution. The new theoretical concept views the unfolding and refolding processes of protein as a shift of the structural distribution induced by a thermodynamic perturbation, with the variance-covariance matrix varying. Based on the theoretical concept, a method to characterize the mechanism of folding (or unfolding) is proposed. The transition state, if any, between two stable states is interpreted as a gap in the distribution, which is created due to an extensive reorganization of hydrogen bonds among back-bone atoms of protein and with water molecules in the course of conformational change. Further perspective to applying the theory to the computer-aided drug design, and to the material science, is briefly discussed.
Collapse
Affiliation(s)
- Fumio Hirata
- Toyota Physical and Chemical Research Institute, Nagakute, Aichi 480-1192, Japan
| | - Masatake Sugita
- College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga 525-8577, Japan
| | - Masasuke Yoshida
- Faculty of Life Sciences, Kyoto Sangyo University, Kyoto, Kyoto 603-8555, Japan
| | - Kazuyuki Akasaka
- Kyoto Prefectural University of Medicine, Kyoto, Kyoto 602-8566, Japan
| |
Collapse
|
21
|
Guan S, Xu Y, Qiao Y, Kuai Z, Qian M, Jiang X, Wang S, Zhang H, Kong W, Shan Y. A novel small molecule displays two different binding modes during inhibiting H1N1 influenza A virus neuraminidases. J Struct Biol 2018; 202:142-149. [DOI: 10.1016/j.jsb.2017.12.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Revised: 12/08/2017] [Accepted: 12/27/2017] [Indexed: 02/06/2023]
|
22
|
Yoshida N, Higashi M, Motoki H, Hirota S. Theoretical analysis of the domain-swapped dimerization of cytochrome c: An MD and 3D-RISM approach. J Chem Phys 2018; 148:025102. [DOI: 10.1063/1.5009785] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
Affiliation(s)
- Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishiku, Fukuoka 819-0395, Japan
| | - Masahiro Higashi
- Department of Chemistry, Biology, and Marine Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Hideyoshi Motoki
- Department of Chemistry, Biology, and Marine Science, University of the Ryukyus, 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | - Shun Hirota
- Graduate School of Materials Science, Nara Institute of Science and Technology, 8916-5 Takayama, Ikoma, Nara 630-0192, Japan
| |
Collapse
|
23
|
Zheng G, Yang F, Fu T, Tu G, Chen Y, Yao X, Xue W, Zhu F. Computational characterization of the selective inhibition of human norepinephrine and serotonin transporters by an escitalopram scaffold. Phys Chem Chem Phys 2018; 20:29513-29527. [DOI: 10.1039/c8cp06232c] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
Selective inhibition of human norepinephrine and serotonin transporters has been studied by computational approaches. 4 warm spots in hNET and 4 in hSERT were found to exert a pronounced effect on inhibition by the studied ligands.
Collapse
Affiliation(s)
- Guoxun Zheng
- College of Pharmaceutical Sciences
- Zhejiang University
- Hangzhou 310058
- China
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
| | - Fengyuan Yang
- College of Pharmaceutical Sciences
- Zhejiang University
- Hangzhou 310058
- China
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
| | - Tingting Fu
- College of Pharmaceutical Sciences
- Zhejiang University
- Hangzhou 310058
- China
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
| | - Gao Tu
- College of Pharmaceutical Sciences
- Zhejiang University
- Hangzhou 310058
- China
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
| | - Yuzong Chen
- Bioinformatics and Drug Design Group
- Department of Pharmacy
- National University of Singapore
- Singapore 117543
- Singapore
| | - Xiaojun Yao
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry
- Lanzhou University
- Lanzhou 730000
- China
| | - Weiwei Xue
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
- Chongqing University
- Chongqing 401331
- China
| | - Feng Zhu
- College of Pharmaceutical Sciences
- Zhejiang University
- Hangzhou 310058
- China
- School of Pharmaceutical Sciences and Collaborative Innovation Center for Brain Science
| |
Collapse
|
24
|
Mahalapbutr P, Chusuth P, Kungwan N, Chavasiri W, Wolschann P, Rungrotmongkol T. Molecular recognition of naphthoquinone-containing compounds against human DNA topoisomerase IIα ATPase domain: A molecular modeling study. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.10.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
25
|
Hasegawa T, Sugita M, Kikuchi T, Hirata F. A Systematic Analysis of the Binding Affinity between the Pim-1 Kinase and Its Inhibitors Based on the MM/3D-RISM/KH Method. J Chem Inf Model 2017; 57:2789-2798. [PMID: 29019402 DOI: 10.1021/acs.jcim.7b00158] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
A systematic study of the binding affinities of 16 lead compounds targeting the Pim-1 kinase based on the 3D-RISM/KH theory and MD simulations is reported. The results show a correlation coefficient R = 0.69 between the theoretical and experimental values of the binding free energy. This demonstrates that the method is applicable to the problem of compound screening and lead optimization, for which relative values of the free energy among the compounds have significance. We elucidate the contribution of the ligand fragments to the binding free energy. Our results indicate that the interactions between the residues and the triazolo[4,3-b]pyridazine scaffold as well as the phenyl ring of the ligand molecule make significant contributions to stabilization of the complex. Using the 3D-RISM/KH theory, we further analyze the probability distribution of a ligand fragment around the protein-ligand complex in which the substituent around the phenyl ring is removed from the ligand. The results demonstrate that the 3D-RISM/KH theory is capable of predicting the position of substitution on a ligand that has a higher affinity to a target protein.
Collapse
Affiliation(s)
- Takeshi Hasegawa
- Department of Bioinformatics, College of Life Science, Ritsumeikan University , 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Masatake Sugita
- Department of Bioinformatics, College of Life Science, Ritsumeikan University , 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Takeshi Kikuchi
- Department of Bioinformatics, College of Life Science, Ritsumeikan University , 1-1-1 Noji-higashi, Kusatsu, Shiga 525-8577, Japan
| | - Fumio Hirata
- Toyota Physical and Chemical Research Institute , 41-1 Yokomichi, Nagakute, Aichi 480-1192, Japan
| |
Collapse
|
26
|
Yoshida N. Role of Solvation in Drug Design as Revealed by the Statistical Mechanics Integral Equation Theory of Liquids. J Chem Inf Model 2017; 57:2646-2656. [DOI: 10.1021/acs.jcim.7b00389] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Norio Yoshida
- Department of Chemistry,
Graduate School of Science, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395 Japan
| |
Collapse
|
27
|
Wongpituk P, Nutho B, Panman W, Kungwan N, Wolschann P, Rungrotmongkol T, Nunthaboot N. Structural dynamics and binding free energy of neral-cyclodextrins inclusion complexes: molecular dynamics simulation. MOLECULAR SIMULATION 2017. [DOI: 10.1080/08927022.2017.1356458] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Peerapong Wongpituk
- Faculty of Science, Department of Chemistry and Center of Excellence for Innovation in Chemistry, Mahasarakham University, Mahasarakham, Thailand
| | - Bodee Nutho
- Faculty of Science, Program in Biotechnology, Chulalongkorn University, Bangkok, Thailand
| | - Wanwisa Panman
- Faculty of Science, Multidisciplinary Program of Petrochemistry and Polymer Science, Chulalongkorn University, Bangkok, Thailand
| | - Nawee Kungwan
- Faculty of Science, Department of Chemistry, Chiang Mai University, Chiang Mai, Thailand
| | - Peter Wolschann
- Faculty of Science, Structural and Computational Biology Research Group, Department of Biochemistry, Chulalongkorn University, Bangkok, Thailand
- Department of Pharmaceutical Technology and Biopharmaceutics, University of Vienna, Vienna, Austria
- Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Thanyada Rungrotmongkol
- Faculty of Science, Structural and Computational Biology Research Group, Department of Biochemistry, Chulalongkorn University, Bangkok, Thailand
- Faculty of Science, Ph.D. Program in Bioinformatics and Computational Biology, Chulalongkorn University, Bangkok, Thailand
- Faculty of Science, Molecular Sensory Science Center, Chulalongkorn University, Bangkok, Thailand
| | - Nadtanet Nunthaboot
- Faculty of Science, Department of Chemistry and Center of Excellence for Innovation in Chemistry, Mahasarakham University, Mahasarakham, Thailand
| |
Collapse
|
28
|
Meeprasert A, Hannongbua S, Kungwan N, Rungrotmongkol T. Effect of D168V mutation in NS3/4A HCV protease on susceptibilities of faldaprevir and danoprevir. MOLECULAR BIOSYSTEMS 2017; 12:3666-3673. [PMID: 27731877 DOI: 10.1039/c6mb00610h] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Hepatitis C virus (HCV) is a serious cause of liver inflammation, cirrhosis and the development of hepatocellular carcinoma. Its NS3/4A serine protease functions to cleave a specific peptide bond, which is an important step in HCV replication. Thus the NS3/4A protease has become one of the main drug-targets in the design and development of anti-HCV agents. Unfortunately, high mutation rates in HCV have been reported due to the lack of RNA proofreading activity resulting in drug resistance. Herein, all-atom molecular dynamics simulations were employed to understand and illustrate the effects of the NS3/4A D168V mutation on faldaprevir (FDV) and danoprevir (DNV) binding efficiency. The D168V mutation was shown to interrupt the hydrogen bonding network of Q80R155D168R123 embedded in the extended S2 and partial S4 subsites of the NS3 protein and as a result the R123 side chain was displaced and moved out from the binding pocket. By means of MM/PBSA and MM/GBSA binding free energy calculations, the FDV and DNV binding affinities were shown to be significantly reduced by ∼10-15 kcal mol-1 and ∼4-9 kcal mol-1 relative to the wild-type complexes, respectively, which somewhat agrees with the experimental resistance folds.
Collapse
Affiliation(s)
- Arthitaya Meeprasert
- Structural and Computational Biology Research Group, Department of Biochemistry, Faculty of Science Chulalongkorn University, 254, Phayathai Road, Pathumwan, Bangkok 10330, Thailand
| | - Supot Hannongbua
- Computational Chemistry Unit Cell, Department of Chemistry, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Pathumwan, Bangkok 10330, Thailand
| | - Nawee Kungwan
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
| | - Thanyada Rungrotmongkol
- Structural and Computational Biology Research Group, Department of Biochemistry, Faculty of Science Chulalongkorn University, 254, Phayathai Road, Pathumwan, Bangkok 10330, Thailand and PhD Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand.
| |
Collapse
|
29
|
Hirata F. A molecular theory of the structural dynamics of protein induced by a perturbation. J Chem Phys 2016; 145:234106. [DOI: 10.1063/1.4971799] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023] Open
Affiliation(s)
- Fumio Hirata
- Toyota Physical and Chemical Research Institute (Toyota Riken), Nagakute, Aichi 480-1192, Japan
| |
Collapse
|
30
|
Phanich J, Rungrotmongkol T, Kungwan N, Hannongbua S. Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study. J Comput Aided Mol Des 2016; 30:917-926. [PMID: 27714494 DOI: 10.1007/s10822-016-9981-5] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2016] [Accepted: 09/27/2016] [Indexed: 12/12/2022]
Abstract
The H7N9 avian influenza virus is a novel re-assortment from at least four different strains of virus. Neuraminidase, which is a glycoprotein on the surface membrane, has been the target for drug treatment. However, some H7N9 strains that have been isolated from patient after drug treatment have a R292K mutation in neuraminidase. This substitution was found to facilitate drug resistance using protein- and virus- assays, in particular it gave a high resistance to the most commonly used drug, oseltamivir. The aim of this research is to understand the source of oseltamivir resistance using MD simulations and the MM/PB(GB)SA binding free energy approaches. Both methods can predict the reduced susceptibility of oseltamivir in good agreement to the IC 50 binding energy, although MM/GBSA underestimates this prediction compared to the MM/PBSA calculation. Electrostatic interaction is the main contribution for oseltamivir binding in terms of both interaction and solvation. We found that the source of the drug resistance is a decrease in the binding interaction combined with the reduction of the dehydration penalty. The smaller K292 mutated residue has a larger binding pocket cavity compared to the wild-type resulting in the loss of drug carboxylate-K292 hydrogen bonding and an increased accessibility for water molecules around the K292 mutated residue. In addition, oseltamivir does not bind well to the R292K mutant complex as shown by the high degree of fluctuation in ligand RMSD during the simulation and the change in angular distribution of bulky side chain groups.
Collapse
Affiliation(s)
- Jiraphorn Phanich
- Computational Chemistry Unit Cell, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Thanyada Rungrotmongkol
- Structural and Computational Biology Research Group, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
- Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Nawee Kungwan
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand
| | - Supot Hannongbua
- Computational Chemistry Unit Cell, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| |
Collapse
|
31
|
Sugita M, Hirata F. Predicting the binding free energy of the inclusion process of 2-hydroxypropyl-β-cyclodextrin and small molecules by means of the MM/3D-RISM method. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:384002. [PMID: 27452185 DOI: 10.1088/0953-8984/28/38/384002] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
A protocol to calculate the binding free energy of a host-guest system is proposed based on the MM/3D-RISM method, taking cyclodextrin derivatives and their ligands as model systems. The protocol involves the procedure to identify the most probable binding mode (MPBM) of receptors and ligands by means of the umbrella sampling method. The binding free energies calculated by the MM/3D-RISM method for the complexes of the seven ligands with the MPBM of the cyclodextrin, and with the fluctuated structures around it, are in agreement with the corresponding experimental data in a semi-quantitative manner. It suggests that the protocol proposed here is promising for predicting the binding affinity of a small ligand to a relatively rigid receptor such as cyclodextrin.
Collapse
Affiliation(s)
- Masatake Sugita
- Department of Bioinformatics, College of Life Sciences, Ritsumeikan University, Kusatsu, Shiga, Japan
| | | |
Collapse
|
32
|
Tanimoto S, Higashi M, Yoshida N, Nakano H. The ion dependence of carbohydrate binding of CBM36: an MD and 3D-RISM study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:344005. [PMID: 27366974 DOI: 10.1088/0953-8984/28/34/344005] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
The molecular recognition process of the carbohydrate-binding module family 36 (CBM36) was examined theoretically. The mechanism of xylan binding by CBM36 and the role of Ca(2+) were investigated by the combined use of molecular dynamics simulations and the 3D reference interaction site model method. The CBM36 showed affinity for xylan after Ca(2+) binding, but not after Mg(2+) binding. Free-energy component analysis of the xylan-binding process revealed that the major factor for xylan-binding affinity is the electrostatic interaction between the Ca(2+) and the hydroxyl oxygens of xylan. The van der Waals interaction between the hydrophobic side chain of CBM36 and the glucopyranose ring of xylan also contributes to the stabilization of the xylan-binding state. Dehydration on the formation of the complex has the opposite effect on these interactions. The affinity of CBM36 for xylan results from a balance of the interactions between the binding ion and solvents, hydrophilic residues around xylan, and the hydroxyl oxygens of xylan. When CBM binds Ca(2+), these interactions are well balanced; in contrast, when CBM binds Mg(2+), the dehydration penalty is excessively large.
Collapse
Affiliation(s)
- Shoichi Tanimoto
- Department of Chemistry, Graduate School of Science, Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan
| | | | | | | |
Collapse
|
33
|
Bodnarchuk MS. Water, water, everywhere… It's time to stop and think. Drug Discov Today 2016; 21:1139-46. [DOI: 10.1016/j.drudis.2016.05.009] [Citation(s) in RCA: 62] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Revised: 04/15/2016] [Accepted: 05/13/2016] [Indexed: 12/11/2022]
|
34
|
Schaduangrat N, Phanich J, Rungrotmongkol T, Lerdsamran H, Puthavathana P, Ubol S. The significance of naturally occurring neuraminidase quasispecies of H5N1 avian influenza virus on resistance to oseltamivir: a point of concern. J Gen Virol 2016; 97:1311-1323. [PMID: 26935590 DOI: 10.1099/jgv.0.000444] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Viral adaptability and survival arise due to the presence of quasispecies populations that are able to escape the immune response or produce drug-resistant variants. However, the presence of H5N1 virus with natural mutations acquired without any drug selection pressure poses a great threat. Cloacal samples collected from the 2004-2005 epidemics in Thailand from Asian open-billed storks revealed one major and several minor quasispecies populations with mutations on the oseltamivir (OTV)-binding site of the neuraminidase gene (NA) without prior exposure to a drug. Therefore, this study investigated the binding between the NA-containing novel mutations and OTV drug using molecular dynamic simulations and plaque inhibition assay. The results revealed that the mutant populations, S236F mutant, S236F/C278Y mutant, A250V/V266A/P271H/G285S mutant and C278Y mutant, had a lower binding affinity with OTV as compared with the WT virus due to rearrangement of amino acid residues and increased flexibility in the 150-loop. This result was further emphasized through the IC50 values obtained for the major population and WT virus, 104.74 nM and 18.30 nM, respectively. Taken together, these data suggest that H5N1 viruses isolated from wild birds have already acquired OTV-resistant point mutations without any exposure to a drug.
Collapse
Affiliation(s)
- Nalini Schaduangrat
- Department of Microbiology, Faculty of Science, Mahidol University, 272 Rama 6 Road, Ratchatewi, Bangkok 10400, Thailand
| | - Jiraphorn Phanich
- Computational Chemistry Unit Cell, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Thanyada Rungrotmongkol
- Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand.,Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Hatairat Lerdsamran
- Department of Microbiology, Faculty of Medicine, Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - Pilaipan Puthavathana
- Department of Microbiology, Faculty of Medicine, Siriraj Hospital, Mahidol University, Bangkok, Thailand
| | - Sukathida Ubol
- Department of Microbiology, Faculty of Science, Mahidol University, 272 Rama 6 Road, Ratchatewi, Bangkok 10400, Thailand
| |
Collapse
|