1
|
Mee-udorn P, Nutho B, Chootrakool R, Maenpuen S, Leartsakulpanich U, Chitnumsub P, Rungrotmongkol T. Structural dynamics and in silico design of pyrazolopyran-based inhibitors against Plasmodium serine hydroxymethyltransferases. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
2
|
Titania Nanosheet Generates Peroxynitrite-Dependent S-Nitrosylation and Enhances p53 Function in Lung Cancer Cells. Pharmaceutics 2021; 13:pharmaceutics13081233. [PMID: 34452194 PMCID: PMC8401232 DOI: 10.3390/pharmaceutics13081233] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Revised: 08/02/2021] [Accepted: 08/03/2021] [Indexed: 12/17/2022] Open
Abstract
Metal nanomaterials can enhance the efficacy of current cancer therapies. Here, we show that Ti0.8O2 nanosheets cause cytotoxicity in several lung cancer cells but not in normal cells. The nanosheet-treated cells showed certain apoptosis characteristics. Protein analysis further indicated the activation of the p53-dependent death mechanism. Transmission electron microscopy (TEM) and scanning electron microscopy (SEM) analyses revealed the cellular uptake of the nanosheets and the induction of cell morphological change. The nanosheets also exhibited a substantial apoptosis effect on drug-resistant metastatic primary lung cancer cells, and it was found that the potency of the nanosheets was dramatically higher than standard drugs. Ti0.8O2 nanosheets induce apoptosis through a molecular mechanism involving peroxynitrite (ONOO−) generation. As peroxynitrite is known to be a potent inducer of S-nitrosylation, we further found that the nanosheets mediated the S-nitrosylation of p53 at C182, resulting in higher protein-protein complex stability, and this was likely to induce the surrounding residues, located in the interface region, to bind more strongly to each other. Molecular dynamics analysis revealed that S-nitrosylation stabilized the p53 dimer with a ΔGbindresidue of <−1.5 kcal/mol. These results provide novel insight on the apoptosis induction effect of the nanosheets via a molecular mechanism involving S-nitrosylation of the p53 protein, emphasizing the mechanism of action of nanomaterials for cancer therapy.
Collapse
|
3
|
Boonma T, Nutho B, Darai N, Rungrotmongkol T, Nunthaboot N. Exploring of paritaprevir and glecaprevir resistance due to A156T mutation of HCV NS3/4A protease: molecular dynamics simulation study. J Biomol Struct Dyn 2021; 40:5283-5294. [PMID: 33430709 DOI: 10.1080/07391102.2020.1869587] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
Abstract
Hepatitis C virus (HCV) NS3/4A serine protease is a promising drug target for the discovery of anti-HCV drugs. However, its amino acid mutations, particularly A156T, commonly lead to rapid emergence of drug resistance. Paritaprevir and glecaprevir, the newly FDA-approved HCV drugs, exhibit distinct resistance profiles against the A156T mutation of HCV NS3/4A serine protease. To illustrate their different molecular resistance mechanisms, molecular dynamics simulations and binding free energy calculations were carried out on the two compounds complexed with both wild-type (WT) and A156T variants of HCV NS3/4A protease. QM/MM-GBSA-based binding free energy calculations revealed that the binding affinities of paritaprevir and glecaprevir towards A156T NS3/4A were significantly reduced by ∼4 kcal/mol with respect to their WT complexes, which were in line with the experimental resistance folds. Moreover, the relatively weak intermolecular interactions with amino acids such as H57, R155, and T156 of NS3 protein, the steric effect and the destabilized protein binding surface, which is caused by the loss of salt bridge between R123 and D168, are the main contributions for the higher fold-loss in potency of glecaprevir due to A156T mutation. An insight into the difference of molecular mechanism of drug resistance against the A156T substitution among the two classes of serine protease inhibitors could be useful for further optimization of new generation HCV NS3/4A inhibitors with enhanced inhibitory potency.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Thitiya Boonma
- Supramolecular Chemistry Research Unit and Department of Chemistry, Faculty of Science, Mahasarakham University, Maha Sarakham, Thailand.,Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH‒CIC), Faculty of Science, Mahasarakham University, Maha Sarakham, Thailand
| | - Bodee Nutho
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Nitchakan Darai
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand
| | - Thanyada Rungrotmongkol
- Program in Bioinformatics and Computational Biology, Graduate School, Chulalongkorn University, Bangkok, Thailand.,Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Nadtanet Nunthaboot
- Supramolecular Chemistry Research Unit and Department of Chemistry, Faculty of Science, Mahasarakham University, Maha Sarakham, Thailand.,Department of Chemistry and Center of Excellence for Innovation in Chemistry (PERCH‒CIC), Faculty of Science, Mahasarakham University, Maha Sarakham, Thailand
| |
Collapse
|
4
|
Abstract
Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.
Collapse
|
5
|
Boonma T, Nutho B, Rungrotmongkol T, Nunthaboot N. Understanding of the drug resistance mechanism of hepatitis C virus NS3/4A to paritaprevir due to D168N/Y mutations: A molecular dynamics simulation perspective. Comput Biol Chem 2019; 83:107154. [PMID: 31751885 DOI: 10.1016/j.compbiolchem.2019.107154] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2019] [Revised: 09/26/2019] [Accepted: 10/21/2019] [Indexed: 02/08/2023]
Abstract
Hepatitis C virus (HCV) NS3/4A protease is an attractive target for the development of antiviral therapy. However, the evolution of antiviral drug resistance is a major problem for treatment of HCV infected patients. Understanding of drug-resistance mechanisms at molecular level is therefore very important for the guidance of further design of antiviral drugs with high efficiency and specificity. Paritaprevir is a potent inhibitor against HCV NS3/4A protease genotype 1a. Unfortunately, this compound is highly susceptible to the substitution at D168 in the protease. In this work, molecular dynamics simulations of paritaprevir complexed with wild-type (WT) and two mutated strains (D168 N and D168Y) were carried out. Due to such mutations, the inhibitor-protein hydrogen bonding between them was weakened and the salt-bridge network among residues R123, R155 and D168 responsible for inhibitor binding was disrupted. Moreover, the per-residue free energy decomposition suggested that the main contributions from key residues such as Q80, V132, K136, G137 and R155 were lost in the D168 N/Y mutations. These lead to a lower binding affinity of paritaprevir for D168 N/Y variants of the HCV NS3/4A protease, consistent with the experimental data. This detailed information could be useful for further design of high potency anti-HCV NS3/4A inhibitors.
Collapse
Affiliation(s)
- Thitiya Boonma
- Supramolecular Chemistry Research Unit and Department of Chemistry, Faculty of Science, Mahasarakham University, Maha Sarakham, 44150, Thailand; Center of Excellence for Innovation in Chemistry (PERCH‒CIC), Faculty of Science, Mahasarakham University, Maha Sarakham, 44150, Thailand
| | - Bodee Nutho
- Center of Excellence in Computational Chemistry (CECC), Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok, 10330, 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
| | - Nadtanet Nunthaboot
- Supramolecular Chemistry Research Unit and Department of Chemistry, Faculty of Science, Mahasarakham University, Maha Sarakham, 44150, Thailand; Center of Excellence for Innovation in Chemistry (PERCH‒CIC), Faculty of Science, Mahasarakham University, Maha Sarakham, 44150, Thailand.
| |
Collapse
|
6
|
Kammarabutr J, Mahalapbutr P, Nutho B, Kungwan N, Rungrotmongkol T. Low susceptibility of asunaprevir towards R155K and D168A point mutations in HCV NS3/4A protease: A molecular dynamics simulation. J Mol Graph Model 2019; 89:122-130. [PMID: 30884449 DOI: 10.1016/j.jmgm.2019.03.006] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Revised: 03/06/2019] [Accepted: 03/06/2019] [Indexed: 01/02/2023]
Abstract
Hepatitis C has become an important health problem that requires expensive treatment and leads to liver tumorigenesis. Hepatitis C virus (HCV), which is the main cause of hepatitis C, has a high mutation rate due to the lack of proofreading activity of the RNA polymerase enzyme. The NS3/4A serine protease is an important target for anti-HCV drug discovery and development because of its crucial role in the cleavage of the polypeptides involved in viral replication. In the present study, all-atom molecular dynamics simulation was performed to elucidate the effect of the single point mutations R155K and D168A in the HCV genotype 1 NS3/4A protease on the structural dynamics, molecular interactions and susceptibility of asunaprevir (ASV), a second-generation NS3/4A protease inhibitor. Principal component analysis indicated that these two mutations converted the direction of motion of residues 123, 155 and 168 in the binding pocket to significantly point outwards from ASV, resulting in a loss of the hydrogen bond network of residues R123···R155···D168. The free energy calculations based on different semiempirical QM/MM-GBSA methods revealed that the binding affinity of ASV with the two mutant forms of the NS3/4A protease was significantly decreased in the order of wild-type < R155K < D168A. This work provided useful structural information regarding the atomistic understanding of acquired drug resistance against ASV caused by the R155K and D168A mutations.
Collapse
Affiliation(s)
- Jirayu Kammarabutr
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Pathumwan, Bangkok, 10330, Thailand
| | - Panupong Mahalapbutr
- Structural and Computational Biology Research Unit, Department of Biochemistry, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Pathumwan, Bangkok, 10330, Thailand
| | - Bodee Nutho
- Program in Biotechnology, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Pathumwan, Bangkok, 10330, Thailand
| | - Nawee Kungwan
- Department of Chemistry, Faculty of Science, Chiang Mai University, 239 Huay Kaew Road, Muang District, 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, 254 Phayathai Road, Pathumwan, Bangkok, 10330, Thailand; Ph.D. Program in Bioinformatics and Computational Biology, Faculty of Science, Chulalongkorn University, Bangkok, 10330, Thailand.
| |
Collapse
|
7
|
Wang C, Greene D, Xiao L, Qi R, Luo R. Recent Developments and Applications of the MMPBSA Method. Front Mol Biosci 2018; 4:87. [PMID: 29367919 PMCID: PMC5768160 DOI: 10.3389/fmolb.2017.00087] [Citation(s) in RCA: 332] [Impact Index Per Article: 55.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Accepted: 11/30/2017] [Indexed: 12/23/2022] Open
Abstract
The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.
Collapse
Affiliation(s)
- Changhao Wang
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, CA, United States
| | - D'Artagnan Greene
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Li Xiao
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
| | - Ruxi Qi
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Ray Luo
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
- Department of Chemical Engineering and Materials Science, University of California, Irvine, Irvine, CA, United States
| |
Collapse
|