• 통합자연과학논문집
• /
• 제5권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.

• EDISON SW 활용 경진대회 논문집
• /
• 제3회(2014년)
• /
• pp.63-74
• /
• 2014

Predicting protein loop structures is an important modeling problem since protein loops are often involved in diverse biological functions by participating in enzyme active sites, ligand binding sites, etc. However, loop structure prediction is difficult even when structures of homologous proteins are known due to large sequence and structure variability among loops of homologous proteins. Therefore, an ab initio approach is necessary to solve loop modeling problems. One of the difficulties in the development of ab initio loop modeling method is to derive an accurate scoring function that closely approximates the true free energy function. In particular, entropy as well as energy contribution have to be considered adequately for loops because loops tend to be flexible compared to other parts of protein. In this study, the contribution of conformational entropy is considered in scoring loop conformations by employing "colony energy" which was previously proposed to estimate the free energy for an ensemble of conformations. Loop conformations were generated by using two EDISON_Chem programs GalaxyFill and GalaxySC, and colony energy was designed for this sampling by tuning relevant parameters. On a test set of 40 loops, the accuracy of predicted loop structure improved on average by scoring with the colony energy compared to scoring by energy alone. In addition, high correlation between colony energy and deviation from the native structure suggested that more extensive sampling can further improve the prediction accuracy. In another test on 6 ligand-binding loops that show conformational changes by ligand binding, both ligand-free and ligand-bound states could be identified by using colony energy when no information on the ligand-bound conformation is used.

• EDISON SW 활용 경진대회 논문집
• /
• 제6회(2017년)
• /
• pp.32-38
• /
• 2017

Prion is a group of the proteins known for its infection mechanisms of Creutzfeldt-Jakob disease (CJD) and other diseases. Solved structures and proven biological roles of fungal prions add tremendous potential to conducting computational simulations. Our research focuses on the binding dynamics of HET-s(218-289), one of the heterokaryon fungal prion originated from Podospora anserina, by calculating the binding free energy using umbrella sampling at 300 K. The binding free energy calculated was $-54.5kcal\;mol^{-1}$, relatively similar to the binding energy of other amyloid fibrils. The simulation result suggests the thermodynamic properties of ${\beta}$-solenoid of HET-s prion and its similarity in dissociation pathways compared to amyloids.

• Bulletin of the Korean Chemical Society
• /
• 제35권8호
• /
• pp.2494-2504
• /
• 2014

P38 mitogen activated protein (MAP) kinase is an important anti-inflammatory drug target, which can be activated by responding to various stimuli such as stress and immune response. Based on the conformation of the conserved DFG loop (in or out), binding inhibitors are termed as type-I and II. Type-I inhibitors are ATP competitive, whereas type-II inhibitors bind in DFG-out conformation of allosteric pocket. It remains unclear that how these allosteric inhibitors stabilize the DFG-out conformation and interact. Organosilicon compounds provide unusual opportunity to enhance potency and diversity of drug molecules due to their low toxicity. However, very few examples have been reported to utilize this property. In this regard, we performed docking of an inhibitor (BIRB) and its silicon analog (Si-BIRB) in an allosteric binding pocket of p38. Further, molecular dynamics (MD) simulations were performed to study the dynamic behavior of the simulated complexes. The difference in the biological activity and mechanism of action of the simulated inhibitors could be explained based on the molecular mechanics/generalized Born surface area (MM/GBSA) binding free energy per residue decomposition. MM/GBSA showed that biological activities were related with calculated binding free energy of inhibitors. Analyses of the per-residue decomposed energy indicated that van der Waals and non-polar interactions were predominant in the ligand-protein interactions. Further, crucial residues identified for hydrogen bond, salt bridge and hydrophobic interactions were Tyr35, Lys53, Glu71, Leu74, Leu75, Ile84, Met109, Leu167, Asp168 and Phe169. Our results indicate that stronger hydrophobic interaction of Si-BIRB with the binding site residues could be responsible for its greater binding affinity compared with BIRB.

• Archives of Pharmacal Research
• /
• 제10권1호
• /
• pp.29-35
• /
• 1987

The effect of drug addition on the bovine serum albumin (BSA)-edible dye complex was studied by spectrophotometric method. The edible dyes tested were amranth, erythrosine, tatrazine and sunset yellow. The moles of bound dye per protein mole and free energies for edible dyes bounded were determined at pH 7.4. The values of free energy change by the addition of drughs to BSA-edible dye were ranged fro -6, 260 to 08030 cal/mole. In the wide range of edible dye concentration (0.3-$7{\times}10^{-5}$$^{-5}$ M), acetylsalicylic acid (ASA) showed pattern of displacement different from that of dye. It was assumed that ASA has different binding mechanisms from edible dye.

• Bulletin of the Korean Chemical Society
• /
• 제24권6호
• /
• pp.817-823
• /
• 2003

We have performed ab initio calculation on the active site of HIV-1 protease. The FEP method was used to determine the binding free energy of four different of protonated states of HIV-1 protease with inhibitor. The structure of the active site and hole structure was taken from the X-ray crystallographic coordinates of the C₂ symmetric inhibitor A74704 protease bound. The active site was modeled with the fragment molecules of binding pocket, acetic acid/ acetate anion (Asp25, Asp125), formamide (amide bond of Thr26/Gly27, Thr126/ Gly127), and methanol as inhibitor fragment. All possibly protonated states of the active site were considered, which were diprotonated state (0, 0), monoprotonated (-1, 0),(0, -1) and diunprotonated state (-1, -1). Once the binding energy Debind, of each model was calculated, more probabilistic protonated states can be proposed from binding energy. From ab-initio results, the FEP simulations were performed for the three following mutations: Ⅰ) Asp25 … Asp125 → AspH25 … Asp125, ⅱ) Asp25 … Asp125 → Asp25 … AspH125, ⅲ) AspH25 … Asp125 → AspH25 … AspH125. The free energy difference between the four states gives the information of the more realistic protonated state of active site aspartic acid. These results provide a theoretical prediction of the protonation state of the catalytic aspartic residues for A74707 complex, and may be useful for the evaluation of potential therapeutic targets.

• Bulletin of the Korean Chemical Society
• /
• 제20권10호
• /
• pp.1129-1135
• /
• 1999

Free energy perturbation and molecular dynamics simulations were carried out to investigate the relative binding affinities of [17] ketonand (1) toward alkali metal cations in methanol. The binding affinities of 1 toward the alkali metal cations were calculated to be in the order Li+ > Na+ > K+ > Rb+ > Cs+, whereas our recent theoretically predicted and experimentally observed binding affinities for [18]starand (2) were in the order K+ > Rb+ > Cs+ > Na+ > Li+. The extremely different affinities of 1 and 2 toward smaller cations, Li + and Na+ , were explained in terms of the differences in their ability to change the conformation to accommodate cations of different sizes. The carbonyl groups constituting the central cavity of 1 can reorganize to form a cavity with the optimal M+ -O distance, even for the smallest Li+, without imposing serious strain on 1. The highest affinity of 1 for Li+ was predominantly due to the highest Coulombic attraction between the smallest Li+ and the carbonyl oxygens of 1.

• EDISON SW 활용 경진대회 논문집
• /
• 제5회(2016년)
• /
• pp.62-66
• /
• 2016

Tumor suppressor protein p53 loses its function upon binding with the HDM2 protein, and inhibiting the p53-HDM2 interaction is critical to suppress tumor cell growth. Recently, the cyclized helix-loop-helix peptide (cHLH) mimicking the ${\alpha}-helix$ part of the p53 protein has been designed and found to exhibit high binding affinity with HDM2. Here, we report the structural and thermodynamic characteristics of the bound complex of the cHLH peptide with the HDM2 protein. We performed molecular dynamics simulations to investigate the structural features of the cHLH peptide as well as its complex with the HDM2. The binding free energy calculation based on the integral equation theory was also executed to quantify the binding affinity for the cHLH/HDM2 complex and to understand the factors contributing to the binding affinity. We found a variety of factors for the helix stability of the cHLH peptide as well as in the complexation with the HDM2, which may provide an insight into the development of anti-cancer drug designs.

• 센서학회지
• /
• 제14권6호
• /
• pp.369-373
• /
• 2005

High quality CdTe single crystal for the solar cell fabrication was grown by vertical Bridgman method. The etch pits patterns of (111) surfaces of CdTe etched by Nakagawa solution was observed the (111)A compesed of Cd atoms with typical triangle etch pits of pyramid mode. From the photoluminescence measurement on (111)A, we observed free exciton ($E_{x}$) existing only high quality crystal and neutal acceptor bound exciton ($A^{0}$,X) having very strong peak intensity. Then, the full width at half maximum and binding energy of neutral acceptor bound exciton were 7 meV and 5.9 meV, respectively. By Haynes rule, an activation enery of impurity was 59 meV. Therefore, the origins on impurity level acting as a neutral acceptor were associated Ag or Cu elements.

• Journal of Microbiology and Biotechnology
• /
• 제23권12호
• /
• pp.1699-1707
• /
• 2013

During the fermentative production of 1,3-propanediol under high substrate concentrations, accumulation of intracellular 3-hydroxypropionaldehyde will cause premature cessation of cell growth and glycerol consumption. Discovery of oxidoreductases that can convert 3-hydroxypropionaldehyde to 1,3-propanediol using NADPH as cofactor could serve as a solution to this problem. In this paper, the yqhD gene from Klebsiella pneumoniae DSM2026, which was found encoding an aldehyde reductase (KpAR), was cloned and characterized. KpAR showed broad substrate specificity under physiological direction, whereas no catalytic activity was detected in the oxidation direction, and both NADPH and NADH can be utilized as cofactors. The cofactor binding mechanism was then investigated employing homology modeling and molecular dynamics simulations. Hydrogen-bond analysis showed that the hydrogen-bond interactions between KpAR and NADPH are much stronger than that for NADH. Free-energy decomposition dedicated that residues Gly37 to Val41 contribute most to the cofactor preference through polar interactions. In conclusion, this work provides a novel aldehyde reductase that has potential applications in the development of novel genetically engineered strains in the 1,3-propanediol industry, and gives a better understanding of the mechanisms involved in cofactor binding.