• Bulletin of the Korean Chemical Society
• /
• v.33 no.11
• /
• pp.3719-3726
• /
• 2012

Formation of Z-DNA, a left-handed double helix, from B-DNA, the canonical right-handed double helix, occurs during important biological processes such as gene expression and DNA transcription. Such B-Z transitions can also be induced by high salt concentration in vitro, but the changes in the relative stability of B-DNA and Z-DNA with salt concentration have not been fully explained despite numerous attempts. For example, electrostatic effects alone could not account for salt-induced B-Z transitions in previous studies. In this paper, we propose that the B-Z transition can be explained if counterion entropy is considered along with the electrostatic interactions. This can be achieved by conducting all-atom, explicit-solvent MD simulations followed by MM-PBSA and molecular DFT calculations. Our MD simulations show that counterions tend to bind at specific sites in B-DNA and Z-DNA, and that more ions cluster near Z-DNA than near B-DNA. Moreover, the difference in counterion ordering near B-DNA and Z-DNA is larger at a low salt concentration than at a high concentration. The results imply that the exclusion of counterions by Z-DNA-binding proteins may facilitate Z-DNA formation under physiological conditions.

• Bulletin of the Korean Chemical Society
• /
• v.35 no.6
• /
• pp.1769-1776
• /
• 2014

Binding modes of a series of catechol derivatives such as protein tyrosine phosphatase 1B (PTP1B) inhibitors were identified by molecular modeling techniques. Docking, molecular dynamics simulations and free energy calculations were employed to determine the modes of these new inhibitors. Binding free energies were calculated by involving different energy components using the Molecular Mechanics-Poisson-Boltzmann Surface Area and Generalized Born Surface Area methods. Relatively larger binding energies were obtained for the catechol derivatives compared to one of the PTP1B inhibitors already in use. The Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) free energy decomposition analysis indicated that the hydroxyl functional groups and biphenyl ring system had favorable interactions with Met258, Tyr46, Gln262 and Phe182 residues of PTP1B. The results of hydrogen bound analysis indicated that catechol derivatives, in addition to hydrogen bonding interactions, Val49, Ile219, Gln266, Asp181 and amino acid residues of PTP1B are responsible for governing the inhibitor potency of the compounds. The information generated from the present study should be useful for the design of more potent PTP1B inhibitors as anti-diabetic agents.

• Journal of Integrative Natural Science
• /
• v.5 no.4
• /
• pp.216-219
• /
• 2012

Structure-based drug design possibly benefit from in silico methods that precisely predict the binding affinity of small molecules to target macromolecules. There are many limitations arise from the difficulty of predicting the binding affinity of a small molecule to a biological target with the current scoring functions. There is thus a strong interest in novel methodologies based on MD simulations that claim predictions of greater accuracy than current scoring functions, helpful for a regular use designed for drug discovery in the pharmaceutical industry. Herein, we report a short review on free energy calculations using MMPBSA method a useful method in structure based drug discovery.

• Journal of Integrative Natural Science
• /
• v.13 no.4
• /
• pp.170-180
• /
• 2020

Abnormal regulation, hyperphosphorylation, and aggregation of the tau protein are the hallmark of several types of dementia, including Alzheimer's Disease. Increased activity of Glycogen Synthase Kinase-3β (GSK-3β) in the Central Nervous System (CNS), increased the tau hyperphosphorylation and caused the neurofibrillary tangles (NFTs) formation in the brain cells. Over the last two decades, numerous adenosine triphosphate (ATP) competitive inhibitors have been discovered that show inhibitory activity against GSK-3β. But these compounds exhibited off-target effects which motivated researchers to find new GSK-3β inhibitors. In the present study, we have collected the dataset of 31 C-Glycosylflavones derivatives that showed inhibitory activity against GSK-3β. Among the dataset, the most active compound was docked with the GSK-3β and molecular dynamics (MD) simulation was performed for 50 ns. Based on the 50 ns MD pose of the most active compound, the other dataset compounds were sketched, minimized, and aligned. The 3D-QSAR based Comparative Molecular Field Analysis (CoMFA) model was developed, which showed a reasonable value of q2=0.664 and r2=0.920. The contour maps generated based on the CoMFA model elaborated on the favorable substitutions at the R2 position. This study could assist in the future development of new GSK-3β inhibitors.