1
|
Yasuda N, Ali S, Aman A, Krusong K, Herfindo N, Chavasiri W, Choowongkomon K, Wolschann P, Mahalapbutr P, Rungrotmongkol T, Hannongbua S. In vitro and in silico studies of the inclusion complexation of 8-bromobaicalein with β-cyclodextrins. J Mol Graph Model 2024; 132:108840. [PMID: 39128382 DOI: 10.1016/j.jmgm.2024.108840] [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: 05/18/2024] [Revised: 07/29/2024] [Accepted: 07/29/2024] [Indexed: 08/13/2024]
Abstract
Baicalein, a flavone derived from Scutellaria baicalensis Georgi, exhibits potent anti-inflammatory, antiviral, and anticancer properties. Its derivative, known as 8-bromobaicalein (BB), has been found to have strong cytotoxic effect on MCF-7 human breast cancer cells. However, its limited solubility in water has hindered its potential for wider applications. To address this issue, we investigated the use of cyclodextrins specifically βCD, 2,6-di-O-methyl-β-cyclodextrin (DMβCD), and hydroxypropyl-β-cyclodextrin (HPβCD) to improve the solubility of BB through inclusion complexation. During 250 ns molecular dynamics simulations, it was found that BB can form inclusion complexes with all βCDs. These complexes exhibit two distinct orientations: chromone group insertion (C-form) and phenyl group insertion (P-form). The formation of these complexes is primarily driven by van der Waals interactions. DMβCD has the highest number of atom contacts with BB and the lowest solvent accessibility in the hydrophobic cavity. These results coincide with the highest binding affinity from the MM/GBSA-based free energy calculation method. Experimental phase solubility diagrams revealed a 1:1 stoichiometric ratio (AL type) between BB and βCDs, in which BB/DMβCD showed the highest stability. The formation of inclusion complexes was confirmed by differential scanning calorimetry and scanning electron microscope methods. Additionally, the BB/DMβCD inclusion complex demonstrated significantly higher anticancer activity against MCF-7 human breast cancer cells compared to BB alone. These findings underscore the potential of DMβCD for formulating BB in pharmaceutical and medical applications.
Collapse
Affiliation(s)
- Noriyuki Yasuda
- Program in Biotechnology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Saba Ali
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Aamir Aman
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Kuakarun Krusong
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Noval Herfindo
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Warinthorn Chavasiri
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Kiattawee Choowongkomon
- Department of Biochemistry, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand
| | - Peter Wolschann
- Institute of Theoretical Chemistry, University of Vienna, 1090, Vienna, Austria
| | - Panupong Mahalapbutr
- Department of Biochemistry, Center for Translational Medicine, Faculty of Medicine, Khon Kaen University, Khon Kaen, 40002, Thailand.
| | - Thanyada Rungrotmongkol
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand; Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| |
Collapse
|
2
|
Ali S, Aman A, Hengphasatporn K, Oopkaew L, Todee B, Fujiki R, Harada R, Shigeta Y, Krusong K, Choowongkomon K, Chavasiri W, Wolschann P, Mahalapbutr P, Rungrotmongkol T. Evaluating solubility, stability, and inclusion complexation of oxyresveratrol with various β-cyclodextrin derivatives using advanced computational techniques and experimental validation. Comput Biol Chem 2024; 112:108111. [PMID: 38879954 DOI: 10.1016/j.compbiolchem.2024.108111] [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: 01/24/2024] [Revised: 05/10/2024] [Accepted: 05/23/2024] [Indexed: 06/18/2024]
Abstract
Oxyresveratrol (OXY), a natural stilbenoid in mulberry fruits, is known for its diverse pharmacological properties. However, its clinical use is hindered by low water solubility and limited bioavailability. In the present study, the inclusion complexes of OXY with β-cyclodextrin (βCD) and its three analogs, dimethyl-β-cyclodextrin (DMβCD), hydroxypropyl-β-cyclodextrin (HPβCD) and sulfobutylether-β-cyclodextrin (SBEβCD), were investigated using in silico and in vitro studies. Molecular docking revealed two binding orientations of OXY, namely, 4',6'-dihydroxyphenyl (A-form) and 5,7-benzenediol ring (B-form). Molecular Dynamics simulations suggested the formation of inclusion complexes with βCDs through two distinct orientations, with OXY/SBEβCD exhibiting maximum atom contacts and the lowest solvent-exposed area in the hydrophobic cavity. These results corresponded well with the highest binding affinity observed in OXY/SBEβCD when assessed using the MM/GBSA method. Beyond traditional simulation methods, Ligand-binding Parallel Cascade Selection Molecular Dynamics method was employed to investigate how the drug enters and accommodates within the hydrophobic cavity. The in silico results aligned with stability constants: SBEβCD (2060 M-1), HPβCD (1860 M-1), DMβCD (1700 M-1), and βCD (1420 M-1). All complexes exhibited a 1:1 binding mode (AL type), with SBEβCD enhancing OXY solubility (25-fold). SEM micrographs, DSC thermograms, FT-IR and 1H NMR spectra confirm the inclusion complex formation, revealing novel surface morphologies, distinctive thermal behaviors, and new peaks. Notably, the inhibitory impact on the proliferation of breast cancer cell lines, MCF-7, exhibited by inclusion complexes particularly OXY/DMβCD, OXY/HPβCD, and OXY/SBEβCD were markedly superior compared to that of OXY alone.
Collapse
Affiliation(s)
- Saba Ali
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Aamir Aman
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok 10330, Thailand
| | - Kowit Hengphasatporn
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| | - Lipika Oopkaew
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Bunyaporn Todee
- Department of Chemistry, Faculty of Science, Mahidol University, Rama VI Rd, Bangkok 10400, Thailand
| | - Ryo Fujiki
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| | - Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan
| | - Kuakarun Krusong
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Kiattawee Choowongkomon
- Department of Biochemistry, Faculty of Science, Kasetsart University, Bangkok 10900, Thailand
| | - Warinthorn Chavasiri
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok 10330, Thailand
| | - Peter Wolschann
- Institute of Theoretical Chemistry, University of Vienna, Vienna 1090, Austria
| | - Panupong Mahalapbutr
- Department of Biochemistry, Center for Translational Medicine, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand.
| | - Thanyada Rungrotmongkol
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand; Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok 10330, Thailand.
| |
Collapse
|
3
|
Nagaraj K, Karuppiah C, Wadaan MA, Maity P, Kaliyaperumal R, Vaishnavi E, Rajaraman D, Abhijith SM, Ramaraj SK, Mathivanan I. Synthesis, characterization, molecular modeling, binding energies of β-cyclodextrin-inclusion complexes of quercetin: Modification of photo physical behavior upon β-CD complexation. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 313:124091. [PMID: 38447439 DOI: 10.1016/j.saa.2024.124091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2023] [Revised: 02/22/2024] [Accepted: 02/26/2024] [Indexed: 03/08/2024]
Abstract
We prepared a naturally occurring flavanoid namely quercetin from tea leaves and analyzed by Absorption, Emission, FT-IR, 1H, 13C nmr spectra and ESI-MS analysis. The inclusion behavior of quercetin in cyclodextrins like α-, β-, γ-, per-6-ABCD and mono-6-ABCD cavities were supported such as UV-vis., Emission, FT-IR and ICD spectra and energy minimization studies. From the absorption and emission results, the type of complexes formed were found to depend on stoichiometry of Host:Guest. FT-IR data of CD complexes of quercetin supported inclusion complex formation of the substrate with α-, β- and γ-CDs. The inclusion of host-guest complexation of quercetin with α-, β-, γ-CDs, per-6-ABCD and mono-6-ABCDs provides very valuable information about the CD:quercetin complexes, the study also shows that β-CD complexation improves water solubility, chemical stability and bioavailability of quercetin. Besides, phase solubility studies also supported the formation of 1:1 drug-CD soluble complexes. All these spectral results provide insight into the binding behavior of substrate into CD cavity in the order per-6-ABCD > Mono-6-ABCD > γ-CD > β-CD > α-CD. The proposed model also finds strong support from the fact with excess CD this exciton coupling disappears indicates the formation of only 1:1 complex.
Collapse
Affiliation(s)
- Karuppiah Nagaraj
- School of Pharmacy, National Forensic Sciences University, 6M56+XP8, Police Bhavan Rd, Sector 9, Gandhinagar, Gujarat 382007, India.
| | - Chelladurai Karuppiah
- Battery Research Center for Green Energy, Ming Chi University of Technology, New Taipei City 24301, Taiwan; PG & Research Department of Chemistry, Thiagarajar College, Madurai, Tamil Nadu, India
| | - Mohammad Ahmad Wadaan
- Department of Zoology, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
| | - Prasenjit Maity
- School of Engineering and Technology, National Forensic Sciences University, 6M56+XP8, Police Bhavan Rd, Sector 9, Gandhinagar, Gujarat 382007, India
| | - Raja Kaliyaperumal
- Department of Chemistry, St. Joseph University, Chumoukedima, Nagaland 797115, India
| | - Ellappan Vaishnavi
- Department of Chemistry, Sri GVG Visalakshi College for Women, Udumalpet 642128, Tamil Nadu, India
| | - D Rajaraman
- Humanities and Sciences, St. Peters Engineering College, St Peters College Rd, Opposite TS Forest Academy Dullapally, Maisammaguda, Medchal, Hyderabad, Telangana 500043, India
| | - S M Abhijith
- School of Pharmacy, National Forensic Sciences University, 6M56+XP8, Police Bhavan Rd, Sector 9, Gandhinagar, Gujarat 382007, India
| | - Sayee Kannan Ramaraj
- PG & Research Department of Chemistry, Thiagarajar College, Madurai, Tamil Nadu, India
| | - Isai Mathivanan
- Research Department of Zoology, Seethalakshmi Ramaswami College (Autonomous), Affiliated to Bharathidasan University, Tiruchirapalli, Tamil Nadu, India
| |
Collapse
|
4
|
Chuntakaruk H, Hengphasatporn K, Shigeta Y, Aonbangkhen C, Lee VS, Khotavivattana T, Rungrotmongkol T, Hannongbua S. FMO-guided design of darunavir analogs as HIV-1 protease inhibitors. Sci Rep 2024; 14:3639. [PMID: 38351065 PMCID: PMC10864397 DOI: 10.1038/s41598-024-53940-1] [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: 06/03/2023] [Accepted: 02/06/2024] [Indexed: 02/16/2024] Open
Abstract
The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.
Collapse
Affiliation(s)
- Hathaichanok Chuntakaruk
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Kowit Hengphasatporn
- 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
| | - Chanat Aonbangkhen
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand
| | - Vannajan Sanghiran Lee
- Chemistry Department, Faculty of Science, University Malaya, Kuala Lumpur, 50603, Malaysia
| | - Tanatorn Khotavivattana
- Center of Excellence in Natural Products Chemistry, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Thanyada Rungrotmongkol
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, 10330, Thailand.
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| |
Collapse
|
5
|
Karnchanapandh K, Hanpaibool C, Sanachai K, Rungrotmongkol T. Elucidation of bezlotoxumab binding specificity to toxin B in Clostridioides difficile. J Biomol Struct Dyn 2024; 42:1617-1628. [PMID: 37098802 DOI: 10.1080/07391102.2023.2201360] [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: 12/26/2022] [Accepted: 04/05/2023] [Indexed: 04/27/2023]
Abstract
C. difficile or Clostridioides difficile infection (CDI) is currently one of the major causes of epidemics worldwide. Toxin B from Clostridioides difficile toxin B (TcdB) infection is the main target protein inhibiting CDI recurrence. Clinical research suggested that bezlotoxumab's (Bez) efficiency is significantly reduced in neutralizing the B2 strain compared to the B1 strain. The monoclonal antibody (mAb) functions by binding to the epitope 1 and 2 regions in the combined repetitive oligopeptide (CROP) domain. Some binding residues are distinctively different between B1 and B2 strains. In this work, we aimed to elucidate and compare insights into the interaction of toxins B1 and B2 in complex with Bez by using all-atom molecular dynamics (MD) simulations and binding free energy calculations. The predicted ΔGbinding values suggested that the antibody (Ab) could bind to toxin B1 significantly better than B2, supported by higher salt bridge and hydrogen bonding (H-bonding) interactions, as well as the number of contact residues between the two focused proteins. The toxin B1 residues important for binding with Bez were E1878, T1901, E1902, F1905, N1941, V1946, N2031, T2032, E2033, V2076, V2077, and E2092. The lower susceptibility of Bez towards toxin B2 was primarily due to a change of residue E2033 from glutamate to alanine (A2033) and the loss of E1878 and E1902 contributions, as determined by the intermolecular interaction changes from the dynamic residue interaction network (dRIN) analysis. The obtained data strengthen our understanding of Bez/toxin B binding.
Collapse
Affiliation(s)
- Kun Karnchanapandh
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
| | - Chonnikan Hanpaibool
- Center of Excellence in Biocatalyst and Sustainable Biotechnology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Kamonpan Sanachai
- Department of Biochemistry, Faculty of Science, Khon Kaen University, Khon Kaen, Thailand
| | - Thanyada Rungrotmongkol
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
- Center of Excellence in Biocatalyst and Sustainable Biotechnology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| |
Collapse
|
6
|
Ali S, Saokaew P, Aman A, Todsaporn D, Sanachai K, Krusong K, Hannongbua S, Wolschann P, Mahalapbutr P, Rungrotmongkol T. Enhancing solubility and stability of piperine using β-cyclodextrin derivatives: computational and experimental investigations. J Biomol Struct Dyn 2024:1-14. [PMID: 38260962 DOI: 10.1080/07391102.2024.2305696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 10/20/2023] [Indexed: 01/24/2024]
Abstract
Piperine (PP), a natural alkaloid found in black pepper, possesses significant bioactivities. However, its use in pharmaceutical applications is hindered by low water solubility and susceptibility to UV light degradation. To overcome these challenges, we investigated the potential of β-cyclodextrin (βCD) and its derivatives with dimethyl (DMβCD), hydroxy-propyl (HPβCD) and sulfobutyl-ether (SBEβCD) substitutions to enhance the solubility and stability of PP. This study employed computational and experimental approaches to examine the complexation between PP and βCDs. The results revealed the formation of two types of inclusion complexes: the P-form and M-form involving the insertion of piperidine moiety and the methylene-di-oxy-phenyl moiety, respectively. These complexes primarily rely on van der Waals interactions. Among the three derivatives, the PP/SBEβCD complex exhibited the highest stability followed by HPβCD, as attributed to maximum atom contacts and minimal solvent accessibility. Solubility studies confirmed the formation of inclusion complexes in a 1:1 ratio. Notably, the stability constant of the inclusion complex was approximately two-fold higher with SBEβCD and HPβCD compared to βCD. The DSC thermograms provided confirmation of the formation of the inclusion complex between the host and guest. These findings highlight the potential of βCD derivatives to effectively encapsulate PP, improving its solubility and presenting new opportunities for its pharmaceutical applications.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Saba Ali
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Phattharapawn Saokaew
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Aamir Aman
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
| | - Duangjai Todsaporn
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Kamonpan Sanachai
- Department of Biochemistry, Faculty of Science, Khon Kaen University, Khon Kaen, Thailand
| | - Kuakarun Krusong
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Supot Hannongbua
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Peter Wolschann
- Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Panupong Mahalapbutr
- Department of Biochemistry, Center for Translational Medicine, Faculty of Medicine, Khon Kaen University, Khon Kaen, Thailand
| | - Thanyada Rungrotmongkol
- Center of Excellence in Structural and Computational Biology, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
| |
Collapse
|
7
|
Punia Bangar S, Kajla P, Chaudhary V, Sharma N, Ozogul F. Luteolin: A flavone with myriads of bioactivities and food applications. FOOD BIOSCI 2023. [DOI: 10.1016/j.fbio.2023.102366] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
|
8
|
Wardana AP, Abdjan MI, Aminah NS, Fahmi MZ, Siswanto I, Kristanti AN, Saputra MA, Takaya Y. 3,4,3'-Tri- O-methylellagic acid as an anticancer agent: in vitro and in silico studies. RSC Adv 2022; 12:29884-29891. [PMID: 36321100 PMCID: PMC9580503 DOI: 10.1039/d2ra05246f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2022] [Accepted: 09/23/2022] [Indexed: 11/06/2022] Open
Abstract
We report a natural product compound isolated from Syzygium polycephalum known as 3,4,3'-tri-O-methylellagic acid (T-EA) as a candidate drug for cancer treatment. The characterization of the isolated T-EA compound was carried out using various spectroscopic methods. The in vitro evaluation showcased the inhibition activity of T-EA towards the T47D and HeLa cell lines with EC50 values of 55.35 ± 6.28 μg mL-1 and 12.57 ± 2.22 μg mL-1, respectively. Meanwhile, the in silico evaluation aimed to understand the interaction of T-EA with enzymes responsible for cancer regulation at the molecular level by targeting the hindrance of cyclin-dependent kinase 9 (CDK9) and sirtuin 1 (SIRT1) enzymes. T-EA showed a binding free energy towards the SIRT1 protein of ΔG bind (MM-GBSA): -30.98 ± 0.25 kcal mol-1 and ΔG bind (MM-PBSA): -24.07 ± 0.30 kcal mol-1, while that of CDK9 was ΔG bind (MM-GBSA): -29.50 ± 0.22 kcal mol-1 and ΔG bind (MM-PBSA): -25.87 ± 0.40 kcal mol-1. The obtained results from this research could be considered as important information on 3,4,3'-tri-O-methylellagic acid as a drug to treat cervical and breast cancers.
Collapse
Affiliation(s)
- Andika Pramudya Wardana
- PhD Student of Mathematics and Natural Sciences, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia,Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501
| | - Muhammad Ikhlas Abdjan
- Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501
| | - Nanik Siti Aminah
- Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501,Biotechnology of Tropical Medicinal Plants Research Group, Universitas AirlanggaIndonesia
| | - Mochamad Zakki Fahmi
- Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501
| | - Imam Siswanto
- Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501,Bioinformatic Laboratory, UCoE Research Center for Bio-Molecule Engineering, Universitas AirlanggaSurabayaIndonesia
| | - Alfinda Novi Kristanti
- Department of Chemistry, Faculty of Science and Technology, Universitas AirlanggaSurabaya 60115Indonesia+62-31-5936502+62-31-5936501,Biotechnology of Tropical Medicinal Plants Research Group, Universitas AirlanggaIndonesia
| | - Mirza Ardella Saputra
- Nanotechnology Engineering, Faculty of Advanced Technology and Multidiscipline, Universitas AirlanggaSurabaya 60115Indonesia
| | - Yoshiaki Takaya
- Faculty of Pharmacy, Meijo University150 Yagotoyama, TempakuNagoya468-8503Japan
| |
Collapse
|
9
|
Kogut MM, Grabowska O, Wyrzykowski D, Samsonov SA. Affinity and putative entrance mechanisms of alkyl sulfates into the β-CD cavity. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119978] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
10
|
Aurisin A Complexed with 2,6-Di- O-methyl-β-cyclodextrin Enhances Aqueous Solubility, Thermal Stability, and Antiproliferative Activity against Lung Cancer Cells. Int J Mol Sci 2022; 23:ijms23179776. [PMID: 36077178 PMCID: PMC9456185 DOI: 10.3390/ijms23179776] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Revised: 08/21/2022] [Accepted: 08/22/2022] [Indexed: 11/16/2022] Open
Abstract
Aurisin A (AA), an aristolane dimer sesquiterpene isolated from the luminescent mushroom Neonothopanus nambi, exhibits various biological and pharmacological effects. However, its poor solubility limits its use for further medicinal applications. This study aimed to improve the water solubility of AA via complexation with β-cyclodextrin (βCD) and its derivatives (2,6-di-O-methyl-βCD (DMβCD) and 2-hydroxypropyl-βCD (HPβCD). A phase solubility analysis demonstrated that the solubility of AA linearly enhanced with increasing concentrations of βCDs (ranked in the order of AA/DMβCD > AA/HPβCD > AA/βCD). Notably, βCDs, especially DMβCD, increased the thermal stability of the inclusion complexes. The thermodynamic study indicated that the complexation between AA and βCD(s) was a spontaneous endothermic reaction, and AA/DMβCD possesses the highest binding strength. The complex formation between AA and DMβCD was confirmed by means of FT-IR, DSC, and SEM. Molecular dynamics simulations revealed that the stability and compactness of the AA/DMβCD complex were higher than those of the DMβCD alone. The encapsulation of AA led to increased intramolecular H-bond formations on the wider rim of DMβCD, enhancing the complex stability. The antiproliferative activity of AA against A549 and H1975 lung cancer cells was significantly improved by complexation with DMβCD. Altogether, the satisfactory water solubility, high thermal stability, and enhanced antitumor potential of the AA/DMβCD inclusion complex would be useful for its application as healthcare products or herbal medicines.
Collapse
|
11
|
Mahalapbutr P, Leechaisit R, Thongnum A, Todsaporn D, Prachayasittikul V, Rungrotmongkol T, Prachayasittikul S, Ruchirawat S, Prachayasittikul V, Pingaew R. Discovery of Anilino-1,4-naphthoquinones as Potent EGFR Tyrosine Kinase Inhibitors: Synthesis, Biological Evaluation, and Comprehensive Molecular Modeling. ACS OMEGA 2022; 7:17881-17893. [PMID: 35664590 PMCID: PMC9161259 DOI: 10.1021/acsomega.2c01188] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2022] [Accepted: 04/21/2022] [Indexed: 05/02/2023]
Abstract
Epidermal growth factor receptor (EGFR) has been recognized as one of the attractive targets for anticancer drug development. Herein, a set of anilino-1,4-naphthoquinone derivatives (3-18) was synthesized and investigated for their anticancer and EGFR inhibitory potentials. Among all tested compounds, three derivatives (3, 8, and 10) were selected for studying EGFR inhibitory activity (in vitro and in silico) due to their most potent cytotoxic activities against six tested cancer cell lines (i.e., HuCCA-1, HepG2, A549, MOLT-3, MDA-MB-231, and T47D; IC50 values = 1.75-27.91 μM), high selectivity index (>20), and good predicted drug-like properties. The experimental results showed that these three promising compounds are potent EGFR inhibitors with nanomolar IC50 values (3.96-18.64 nM). Interestingly, the most potent compound 3 bearing 4-methyl substituent on the phenyl ring displayed 4-fold higher potency than the known EGFR inhibitor, erlotinib. Molecular docking, molecular dynamics simulation, and MM/GBSA-based free energy calculation revealed that van der Waals force played a major role in the accommodations of compound 3 within the ATP-binding pocket of EGFR. Additionally, the 4-CH3 moiety of the compound was noted to be a key chemical feature contributing to the highly potent EGFR inhibitory activity via its formations of alkyl interactions with A743, K745, M766, and L788 residues as well as additional interactions with M766 and T790.
Collapse
Affiliation(s)
- Panupong Mahalapbutr
- Department
of Biochemistry, and Center for Translational Medicine, Faculty of
Medicine, Khon Kaen University, Khon Kaen 40002, Thailand
| | - Ronnakorn Leechaisit
- Department
of Chemistry, Faculty of Science, Srinakharinwirot
University, Bangkok 10110, Thailand
| | - Anusit Thongnum
- Department
of Physics, Faculty of Science, Srinakharinwirot
University, Bangkok 10110, Thailand
| | - Duangjai Todsaporn
- Structural
and Computational Biology Research Unit, Department of Biochemistry,
Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
| | - Veda Prachayasittikul
- Center
of Data Mining and Biomedical Informatics, Faculty of Medical Technology, Mahidol University, Bangkok 10700, Thailand
| | - Thanyada Rungrotmongkol
- Structural
and Computational Biology 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
| | - Supaluk Prachayasittikul
- Center
of Data Mining and Biomedical Informatics, Faculty of Medical Technology, Mahidol University, Bangkok 10700, Thailand
| | - Somsak Ruchirawat
- Laboratory
of Medicinal Chemistry and Program in Chemical Sciences, Chulabhorn Research Institute, Chulabhorn Graduate
Institute, Bangkok 10210, Thailand
- Commission
on Higher Education, Ministry of Education, Center of Excellence on Environmental Health and Toxicology (EHT), Bangkok 10400, Thailand
| | - Virapong Prachayasittikul
- Department
of Clinical Microbiology and Applied Technology, Faculty of Medical
Technology, Mahidol University, Bangkok 10700, Thailand
| | - Ratchanok Pingaew
- Department
of Chemistry, Faculty of Science, Srinakharinwirot
University, Bangkok 10110, Thailand
| |
Collapse
|
12
|
Molecular simulation and experimental study on the inclusion of rutin with β-cyclodextrin and its derivative. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.132359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
13
|
Baezzat MR, Jahromi FZ. Differential Pulse Voltammetric Determination of Acetaminophen Using Carbon Paste Electrode Modified with β-Cyclodextrin/Gold/Titanium Dioxide Nanocomposite. JOURNAL OF ANALYTICAL CHEMISTRY 2022. [DOI: 10.1134/s1061934822010026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
14
|
Molecular encapsulation of emodin with various β-cyclodextrin derivatives: A computational study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118002] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
15
|
Kerdpol K, Nutho B, Krusong K, Poo-arporn RP, Rungrotmongkol T, Hannongbua S. Encapsulation of α-tocopherol in large-ring cyclodextrin containing 26 α-D-glucopyranose units: A molecular dynamics study. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116802] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
|
16
|
Mahalapbutr P, Charoenwongpaiboon T, Phongern C, Kongtaworn N, Hannongbua S, Rungrotmongkol T. Molecular encapsulation of a key odor-active 2-acetyl-1-pyrroline in aromatic rice with β-cyclodextrin derivatives. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116394] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
|
17
|
Mazurek AH, Szeleszczuk Ł, Gubica T. Application of Molecular Dynamics Simulations in the Analysis of Cyclodextrin Complexes. Int J Mol Sci 2021; 22:9422. [PMID: 34502331 PMCID: PMC8431145 DOI: 10.3390/ijms22179422] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2021] [Revised: 08/27/2021] [Accepted: 08/27/2021] [Indexed: 11/17/2022] Open
Abstract
Cyclodextrins (CDs) are highly respected for their ability to form inclusion complexes via host-guest noncovalent interactions and, thus, ensofance other molecular properties. Various molecular modeling methods have found their applications in the analysis of those complexes. However, as showed in this review, molecular dynamics (MD) simulations could provide the information unobtainable by any other means. It is therefore not surprising that published works on MD simulations used in this field have rapidly increased since the early 2010s. This review provides an overview of the successful applications of MD simulations in the studies on CD complexes. Information that is crucial for MD simulations, such as application of force fields, the length of the simulation, or solvent treatment method, are thoroughly discussed. Therefore, this work can serve as a guide to properly set up such calculations and analyze their results.
Collapse
Affiliation(s)
- Anna Helena Mazurek
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Doctoral School, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| | - Łukasz Szeleszczuk
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| | - Tomasz Gubica
- Department of Physical Chemistry, Chair of Physical Pharmacy and Bioanalysis, Faculty of Pharmacy, Medical University of Warsaw, Banacha 1 Street, 02-093 Warsaw, Poland;
| |
Collapse
|
18
|
Substrate binding mechanism of glycerophosphodiesterase towards organophosphate pesticides. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
19
|
Aiebchun T, Mahalapbutr P, Auepattanapong A, Khaikate O, Seetaha S, Tabtimmai L, Kuhakarn C, Choowongkomon K, Rungrotmongkol T. Identification of Vinyl Sulfone Derivatives as EGFR Tyrosine Kinase Inhibitor: In Vitro and In Silico Studies. Molecules 2021; 26:molecules26082211. [PMID: 33921332 PMCID: PMC8069501 DOI: 10.3390/molecules26082211] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2020] [Revised: 04/05/2021] [Accepted: 04/07/2021] [Indexed: 01/17/2023] Open
Abstract
Epidermal growth factor receptor (EGFR), overexpressed in many types of cancer, has been proved as a high potential target for targeted cancer therapy due to its role in regulating proliferation and survival of cancer cells. In the present study, a series of designed vinyl sulfone derivatives was screened against EGFR tyrosine kinase (EGFR-TK) using in silico and in vitro studies. The molecular docking results suggested that, among 78 vinyl sulfones, there were eight compounds that could interact well with the EGFR-TK at the ATP-binding site. Afterwards, these screened compounds were tested for the inhibitory activity towards EGFR-TK using ADP-Glo™ kinase assay, and we found that only VF16 compound exhibited promising inhibitory activity against EGFR-TK with the IC50 value of 7.85 ± 0.88 nM. In addition, VF16 showed a high cytotoxicity with IC50 values of 33.52 ± 2.57, 54.63 ± 0.09, and 30.38 ± 1.37 µM against the A431, A549, and H1975 cancer cell lines, respectively. From 500-ns MD simulation, the structural stability of VF16 in complex with EGFR-TK was quite stable, suggesting that this compound could be a novel small molecule inhibitor targeting EGFR-TK.
Collapse
Affiliation(s)
- Thitinan Aiebchun
- Biocatalyst and Environmental Biotechnology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand;
| | - Panupong Mahalapbutr
- Department of Biochemistry, Faculty of Medicine, Khon Kaen University, Khon Kaen 40002, Thailand;
| | - Atima Auepattanapong
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok 10700, Thailand; (A.A.); (O.K.); (C.K.)
| | - Onnicha Khaikate
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok 10700, Thailand; (A.A.); (O.K.); (C.K.)
| | - Supaphorn Seetaha
- Department of Biochemistry, Faculty of Science, Kasetsart University, Chatuchak, Bangkok 10900, Thailand;
| | - Lueacha Tabtimmai
- Department of Biotechnology, Faculty of Applied Science, King Mongkut’s University of Technology of North Bangkok, Bangkok 10800, Thailand;
| | - Chutima Kuhakarn
- Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH-CIC), Faculty of Science, Mahidol University, Bangkok 10700, Thailand; (A.A.); (O.K.); (C.K.)
| | - Kiattawee Choowongkomon
- Department of Biochemistry, Faculty of Science, Kasetsart University, Chatuchak, Bangkok 10900, Thailand;
- Correspondence: (K.C.); (T.R.); Tel.: +66-2218-5426 (T.R.); Fax: +66-2218-5418 (T.R.)
| | - 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, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand
- Correspondence: (K.C.); (T.R.); Tel.: +66-2218-5426 (T.R.); Fax: +66-2218-5418 (T.R.)
| |
Collapse
|
20
|
Sanachai K, Aiebchun T, Mahalapbutr P, Seetaha S, Tabtimmai L, Maitarad P, Xenikakis I, Geronikaki A, Choowongkomon K, Rungrotmongkol T. Discovery of novel JAK2 and EGFR inhibitors from a series of thiazole-based chalcone derivatives. RSC Med Chem 2021; 12:430-438. [PMID: 34046625 PMCID: PMC8130606 DOI: 10.1039/d0md00436g] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2020] [Accepted: 01/25/2021] [Indexed: 01/12/2023] Open
Abstract
The Janus kinase (JAK) and epidermal growth factor receptor (EGFR) have been considered as potential targets for cancer therapy due to their role in regulating proliferation and survival of cancer cells. In the present study, the aromatic alkyl-amino analogs of thiazole-based chalcone were selected to experimentally and theoretically investigate their inhibitory activity against JAK2 and EGFR proteins as well as their anti-cancer effects on human cancer cell lines expressing JAK2 (TF1 and HEL) and EGFR (A549 and A431). In vitro cytotoxicity screening results demonstrated that the HEL erythroleukemia cell line was susceptible to compounds 11 and 12, whereas the A431 lung cancer cell line was vulnerable to compound 25. However, TF1 and A549 cells were not sensitive to our thiazole derivatives. From kinase inhibition assay results, compound 25 was found to be a dual inhibitor against JAK2 and EGFR, whereas compounds 11 and 12 selectively inhibited the JAK2 protein. According to the molecular docking analysis, compounds 11, 12 and 25 formed hydrogen bonds with the hinge region residues Lys857, Leu932 and Glu930 and hydrophobically came into contact with Leu983 at the catalytic site of JAK2, while compound 25 formed a hydrogen bond with Met769 at the hinge region, Lys721 near a glycine loop, and Asp831 at the activation loop of EGFR. Altogether, these potent thiazole derivatives, following Lipinski's rule of five, could likely be developed as a promising JAK2/EGFR targeted drug(s) for cancer therapy.
Collapse
Affiliation(s)
- Kamonpan Sanachai
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University Bangkok 10330 Thailand +662 2185418 +662 2185426
| | - Thitinan Aiebchun
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University Bangkok 10330 Thailand +662 2185418 +662 2185426
| | - Panupong Mahalapbutr
- Department of Biochemistry, Faculty of Medicine, Khon Kaen University Khon Kaen 40002 Thailand
| | - Supaphorn Seetaha
- Department of Biochemistry, Faculty of Science, Kasetsart University Bangkok 10900 Thailand
| | - Lueacha Tabtimmai
- Department of Biotechnology, Faculty of Applied Science, King Mongkut's University of Technology of North Bangkok Bangkok Thailand
| | - Phornphimon Maitarad
- Research Center of Nano Science and Technology, Shanghai University Shanghai 200444 PR China
| | - Iakovos Xenikakis
- Department of Pharmaceutical Chemistry, School of Pharmacy, Aristotle University of Thessaloniki Thessaloniki 54124 Greece
| | - Athina Geronikaki
- Department of Pharmaceutical Chemistry, School of Pharmacy, Aristotle University of Thessaloniki Thessaloniki 54124 Greece
| | | | - Thanyada Rungrotmongkol
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University Bangkok 10330 Thailand +662 2185418 +662 2185426
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University Bangkok 10330 Thailand
| |
Collapse
|
21
|
Karnchanapandh K, Hanpaibool C, Mahalapbutr P, Rungrotmongkol T. Source of oseltamivir resistance due to single E276D, R292K, and double E276D/R292K mutations in H10N4 influenza neuraminidase. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.115294] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
22
|
Computational study on peptidomimetic inhibitors against SARS-CoV-2 main protease. J Mol Liq 2020; 322:114999. [PMID: 33518853 PMCID: PMC7832253 DOI: 10.1016/j.molliq.2020.114999] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Revised: 10/25/2020] [Accepted: 12/06/2020] [Indexed: 12/11/2022]
Abstract
The emergence outbreak caused by a novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has received significant attention on the global risks. Due to itscrucial role in viral replication, the main protease 3CLpro is an important target for drug discovery and development to combat COVID-19. In this work, the structural and dynamic behaviors as well as binding efficiency of the four peptidomimetic inhibitors (N3, 11a, 13b, and 14b) recently co-crystalized with SARS-CoV-2 3CLpro were studied and compared using all-atom molecular dynamics (MD) simulations and solvated interaction energy-based binding free energy calculations. The per-residue decomposition free energy results suggested that the key residues involved in inhibitors binding were H41, M49, L141-C145, H163-E166, P168, and Q189-T190 in the domains I and II. The van der Waals interaction yielded the main energy contribution stabilizing all the focused inhibitors. Besides, their hydrogen bond formations with F140, G143, C145, H164, E166, and Q189 residues in the substrate-binding pocket were also essential for strengthening the molecular complexation. The predicted binding affinity of the four peptidomimetic inhibitors agreed with the reported experimental data, and the 13b showed the most efficient binding to SARS-CoV-2 3CLpro. From rational drug design strategies based on 13b, the polar moieties (e.g., benzamide) and the bulky N-terminal protecting groups (e.g., thiazole) should be introduced to P1' and P4 sites in order to enhance H-bonds and hydrophobic interactions, respectively. We hope that the obtained structural and energetic information could be beneficial for developing novel SARS-CoV-2 3CLpro inhibitors with higher inhibitory potency to combat COVID-19.
Collapse
|
23
|
Abdjan MI, Aminah NS, Siswanto I, Thant TM, Kristanti AN, Takaya Y. In silico approach: biological prediction of nordentatin derivatives as anticancer agent inhibitors in the cAMP pathway. RSC Adv 2020; 10:42733-42743. [PMID: 35514899 PMCID: PMC9058016 DOI: 10.1039/d0ra07838g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2020] [Accepted: 11/06/2020] [Indexed: 11/23/2022] Open
Abstract
A combination of computational techniques has been carried out to predict the binding of nordentatin derivatives based on pyranocoumarin semi-synthesis with the target protein from the expression of the PDE4B gene. The inhibition of the cAMP pathway is the main target of anti-cancer drugs, which is responsible for uncontrolled cell division in cancer. Modeling was done using a combination of semi-empirical methods and the density functional theory (PM3-DFT/6-31G*/B3LYP) to obtain the optimal structure of a small ligand that could be modeled. Studies on the interaction of the ligands and amino acid residues on protein targets were carried out using a combination of molecular docking and molecular dynamic simulation. Molecular docking based on functional grid scores showed a very good native ligand pose with an RMSD of 0.93 Å in determining the initial coordinates of the ligand-receptor interactions. Furthermore, the amino acid residues responsible for interaction through H-bonds were Tyr103, His104, His177, Met217, and Gln313. The binding free energy (kcal mol-1) results of the candidates were PS-1 (-36.84 ± 0.31), PS-2 (-35.34 ± 0.28), PS-3 (-26.65 ± 0.30), PS-5 (-42.66 ± 0.26), PS-7 (-35.33 ± 0.23), and PS-9 (-32.57 ± 0.20), which are smaller than that of the native ligand Z72 (-24.20 ± 0.19), and thus these have good potential as drugs that can inhibit the cAMP pathway. These results provide theoretical information for the efficient inhibition of the cAMP pathway in the future.
Collapse
Affiliation(s)
- Muhammad Ikhlas Abdjan
- Departement of Chemistry, Faculty of Science and Technology, Universitas Airlangga Surabaya 60115 Indonesia
| | - Nanik Siti Aminah
- Departement of Chemistry, Faculty of Science and Technology, Universitas Airlangga Surabaya 60115 Indonesia
- Biotechnology of Tropical Medicinal Plants Research Group, Universitas Airlangga Indonesia
| | - Imam Siswanto
- Departement of Chemistry, Faculty of Science and Technology, Universitas Airlangga Surabaya 60115 Indonesia
| | - Tin Myo Thant
- Department of Chemistry, Faculty of Science and Technology, Universitas Airlangga, Komplek Kampus C Jl. Mulyorejo Surabaya Indonesia 60115
- Department of Chemistry, Mandalar Degree College Mandalay Myanmar
| | - Alfinda Novi Kristanti
- Departement of Chemistry, Faculty of Science and Technology, Universitas Airlangga Surabaya 60115 Indonesia
- Biotechnology of Tropical Medicinal Plants Research Group, Universitas Airlangga Indonesia
| | - Yoshiaki Takaya
- Faculty of Pharmacy, Meijo University 150 Yagotoyama, Tempaku Nagoya 468-8503 Japan
| |
Collapse
|
24
|
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
|
25
|
Kammarabutr J, Mahalapbutr P, Okumura H, Wolschann P, Rungrotmongkol T. Structural dynamics and susceptibility of anti-HIV drugs against HBV reverse transcriptase. J Biomol Struct Dyn 2020; 39:2502-2511. [PMID: 32308149 DOI: 10.1080/07391102.2020.1751715] [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: 02/08/2023]
Abstract
Hepatitis B virus (HBV), a small enveloped DNA virus, attacks the human liver causing both acute and chronic diseases. Current therapeutic drugs use the nucleos(t)ide analogues (NAs) as a competitive inhibitor against HBV reverse transcriptase (HBV-RT), an essential enzyme pivotally involved in viral replication. Unfortunately, this treatment still causes the development of resistant variants of HBV against NAs. As HBV-RT is homologous to the human immunodeficiency virus reverse transcriptase (HIV-RT), it is reasonable to treat HBV-RT with anti-HIV drugs. In the present study, we aimed to investigate the structural dynamics and susceptibility of the known anti-HIV drugs (stavudine [d4T], didanosine [DDI], and zidovudine [ZDV]) against HBV-RT enzyme in comparison to the anti-HBV drug lamivudine (3TC) and deoxythymidine triphosphate (dTTP) substrate using several computational approaches. The ΔGbindresidue calculations revealed that seven polar residues (K32, R41, D83, S85, D205, N236, and K239) and three hydrophobic residues (A86, A87, and F88) of HBV-RT as well as the adjacent DNA strands play an important role in the ligand binding. In addition, the H-bond pattern of d4T is similar to that of 3TC, especially at the residues A86 and A87. Such interactions promote the favorable conformation of ligand in the HBV-RT binding pocket, while the several different conformations of ligand are found in the unbound state. The predicted binding free energy results based on QM/MM-GBSA and MM/GB(PB)SA methods suggested that the susceptibility towards HBV-RT of d4T and ZDV is higher than that of 3TC and dTTP. Altogether, this work sheds light on the potentiality of d4T and ZDV as a promising drug for HBV-infected patients harboring 3TC resistance.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Jirayu Kammarabutr
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok, Thailand
| | - Panupong Mahalapbutr
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok, Thailand
| | - Hisashi Okumura
- Institute for Molecular Science, Research Center for Computational Science, Okazaki, Japan.,Department of Structural Molecular Science, The Graduate University for Advanced Studies (SOKENDAI), Okazaki, Japan
| | - Peter Wolschann
- Department of Pharmaceutical Chemistry, Faculty of Life Sciences, University of Vienna, Vienna, Austria.,Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Thanyada Rungrotmongkol
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Pathumwan, Bangkok, Thailand.,Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| |
Collapse
|
26
|
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
|
27
|
Enhanced Solubility and Anticancer Potential of Mansonone G By β-Cyclodextrin-Based Host-Guest Complexation: A Computational and Experimental Study. Biomolecules 2019; 9:biom9100545. [PMID: 31569832 PMCID: PMC6843486 DOI: 10.3390/biom9100545] [Citation(s) in RCA: 38] [Impact Index Per Article: 7.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2019] [Revised: 09/24/2019] [Accepted: 09/26/2019] [Indexed: 11/17/2022] Open
Abstract
Mansonone G (MG), a plant-derived compound isolated from the heartwood of Mansonia gagei, possesses a potent antitumor effect on several kinds of malignancy. However, its poor solubility limits the use for practical applications. Beta-cyclodextrin (βCD), a cyclic oligosaccharide composed of seven (1→4)-linked α-D-glucopyranose units, is capable of encapsulating a variety of poorly soluble compounds into its hydrophobic interior. In this work, we aimed to enhance the water solubility and the anticancer activity of MG by complexation with βCD and its derivatives (2,6-di-O-methyl-βCD (DMβCD) and hydroxypropyl-βCD). The 90-ns molecular dynamics simulations and MM/GBSA-based binding free energy results suggested that DMβCD was the most preferential host molecule for MG inclusion complexation. The inclusion complex formation between MG and βCD(s) was confirmed by DSC and SEM techniques. Notably, the MG/βCDs inclusion complexes exerted significantly higher cytotoxic effect (~2–7 fold) on A549 lung cancer cells than the uncomplexed MG.
Collapse
|
28
|
Hedayati N, Bemani Naeini M, Mohammadinejad A, Mohajeri SA. Beneficial effects of celery (
Apium graveolens
) on metabolic syndrome: A review of the existing evidences. Phytother Res 2019; 33:3040-3053. [DOI: 10.1002/ptr.6492] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Revised: 08/08/2019] [Accepted: 08/10/2019] [Indexed: 12/31/2022]
Affiliation(s)
- Narges Hedayati
- Pharmaceutical Research Center, Pharmaceutical Technology InstituteMashhad University of Medical Sciences Mashhad Iran
- Department of Pharmacodynamics and Toxicology, School of PharmacyMashhad University of Medical Sciences Mashhad Iran
| | - Mehri Bemani Naeini
- Nanotechnology Research CenterMashhad University of Medical Sciences Mashhad Iran
| | - Arash Mohammadinejad
- Pharmaceutical Research Center, Pharmaceutical Technology InstituteMashhad University of Medical Sciences Mashhad Iran
- Department of ChemistryPayame Noor University Tehran Iran
- Department of Pharmacodynamics and Toxicology, School of PharmacyMashhad University of Medical Sciences Mashhad Iran
| | - Seyed Ahmad Mohajeri
- Pharmaceutical Research Center, Pharmaceutical Technology InstituteMashhad University of Medical Sciences Mashhad Iran
- Department of Pharmacodynamics and Toxicology, School of PharmacyMashhad University of Medical Sciences Mashhad Iran
| |
Collapse
|
29
|
Comparative analysis of phenolic compounds in four taxa of Erigeron acris s. l. (Asteraceae). Biologia (Bratisl) 2019. [DOI: 10.2478/s11756-019-00332-w] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
Abstract
AbstractThe aim of the present work was to investigate and compare the content of phenolic compounds in four taxa of Erigeron acris L. s. l.: E. acris (EAA), E. acris subsp. droebachiensis (O.F. Müll.) Arcang. (EAD), E. acris subsp. serotinus (Weihe) Greuter (EAS) and E. ×huelsenii Vatke (EH), a hybrid between E. acris and E. canadensis L. The total flavonoid content was determined by Christ-Müller method and the total phenolic acid content was determined by the method utilizing Arnov’s reagent. The method using ultra high performance liquid chromatography with photodiode array detection (UHPLC-PDA) was applied for the separation, identification and quantification of nine phenolic compounds (protocatechuic acid, chlorogenic acid, caffeic acid, 6′-O-caffeoylerigeroside, scutellarein-7-O-β-D-glucuronide, quercetin 3-O-glucoside, 4,5-dicaffeoylquinic acid, quercetin and luteolin) in the aerial parts of E. acris s. l. The chromatographic separation was carried out using a BEH C18 column packed with 1.7-μm particles and gradient elution with a mobile phase of water and methanol, both containing 0.02% (v/v) trifluoroacetic acid. The four investigated taxa of E. acris s. l. differed in the composition and the content of phenolic compounds. The main substances determined in the methanolic herbal extracts were: scutellarein-7-O-β-D-glucuronide (EAA, EAS, EAD and EH), 6′-O-caffeoylerigeroside (EAA, EAD and EH) and chlorogenic acid (EAS and EH). Moreover, the results indicated that five of the nine tested compounds were found in all investigated extracts from herbs of E. acris s. l. Two of them (6′-O-caffeoylerigeroside and scutellarein-7-O-β-D-glucuronide) could be selected as potential chemotaxonomic markers of the genus Erigeron L.
Collapse
|
30
|
Butoxy Mansonone G Inhibits STAT3 and Akt Signaling Pathways in Non-Small Cell Lung Cancers: Combined Experimental and Theoretical Investigations. Cancers (Basel) 2019; 11:cancers11040437. [PMID: 30925736 PMCID: PMC6521096 DOI: 10.3390/cancers11040437] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2019] [Revised: 03/22/2019] [Accepted: 03/25/2019] [Indexed: 01/20/2023] Open
Abstract
Epidermal growth factor receptor (EGFR) is the key molecular target for non-small cell lung cancer (NSCLC) due to its major contribution to complex signaling cascades modulating the survival of cancer cells. Targeting EGFR-mediated signaling pathways has been proved as a potential strategy for NSCLC treatment. In the present study, mansonone G (MG), a naturally occurring quinone-containing compound, and its semi-synthetic ether derivatives were subjected to investigate the anticancer effects on human NSCLC cell lines expressing wild-type EGFR (A549) and mutant EGFR (H1975). In vitro cytotoxicity screening results demonstrated that butoxy MG (MG3) exhibits the potent cytotoxic effect on both A549 (IC50 of 8.54 μM) and H1975 (IC50 of 4.21 μM) NSCLC cell lines with low toxicity against PCS201-010 normal fibroblast cells (IC50 of 21.16 μM). Western blotting and flow cytometric analyses revealed that MG3 induces a caspase-dependent apoptosis mechanism through: (i) inhibition of p-STAT3 and p-Akt without affecting upstream p-EGFR and (ii) activation of p-Erk. The 500-ns molecular dynamics simulations and the molecular mechanics combined with generalized Born surface area (MM/GBSA)-based binding free energy calculations suggested that MG3 could possibly interact with STAT3 SH2 domain and ATP-binding pocket of Akt. According to principal component analysis, the binding of MG3 toward STAT3 and Akt dramatically altered the conformation of proteins, especially the residues in the active site, stabilizing MG3 mainly through van der Waals interactions.
Collapse
|
31
|
Kerdpol K, Kicuntod J, Wolschann P, Mori S, Rungnim C, Kunaseth M, Okumura H, Kungwan N, Rungrotmongkol T. Cavity Closure of 2-Hydroxypropyl-β-Cyclodextrin: Replica Exchange Molecular Dynamics Simulations. Polymers (Basel) 2019; 11:polym11010145. [PMID: 30960130 PMCID: PMC6401915 DOI: 10.3390/polym11010145] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2018] [Revised: 01/09/2019] [Accepted: 01/11/2019] [Indexed: 12/15/2022] Open
Abstract
2-Hydroxypropyl-β-cyclodextrin (HPβCD) has unique properties to enhance the stability and the solubility of low water-soluble compounds by inclusion complexation. An understanding of the structural properties of HPβCD and its derivatives, based on the number of 2-hydroxypropyl (HP) substituents at the α-d-glucopyranose subunits is rather important. In this work, replica exchange molecular dynamics simulations were performed to investigate the conformational changes of single- and double-sided HP-substitution, called 6-HPβCDs and 2,6-HPβCDs, respectively. The results show that the glucose subunits in both 6-HPβCDs and 2,6-HPβCDs have a lower chance of flipping than in βCD. Also, HP groups occasionally block the hydrophobic cavity of HPβCDs, thus hindering drug inclusion. We found that HPβCDs with a high number of HP-substitutions are more likely to be blocked, while HPβCDs with double-sided HP-substitutions have an even higher probability of being blocked. Overall, 6-HPβCDs with three and four HP-substitutions are highlighted as the most suitable structures for guest encapsulation, based on our conformational analyses, such as structural distortion, the radius of gyration, circularity, and cavity self-closure of the HPβCDs.
Collapse
Affiliation(s)
- Khanittha Kerdpol
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand.
| | - Jintawee Kicuntod
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand.
| | - Peter Wolschann
- Structural and Computational Biology Research Unit, Department of Biochemistry, 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.
| | - Seiji Mori
- Institute of Quantum Beam Science, Graduate School of Science and Engineering, Ibaraki University, 2-1-1 Bunkyo, Mito, Ibaraki 310-8512, Japan.
| | - Chompoonut Rungnim
- National Nanotechnology Center (NANOTEC), National Science and Technology Development Agency (NSTDA), Pathum Thani 12120, Thailand.
| | - Manaschai Kunaseth
- National Nanotechnology Center (NANOTEC), National Science and Technology Development Agency (NSTDA), Pathum Thani 12120, Thailand.
| | - Hisashi Okumura
- Research Center for Computational Science, Institute for Molecular Science, Okazaki, Aichi 444-8585, 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.
| | - Thanyada Rungrotmongkol
- Structural and Computational Biology 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.
- Molecular Sensory Science Center, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Patumwan, Bangkok 10330, Thailand.
| |
Collapse
|