• Bulletin of the Korean Chemical Society
• /
• v.33 no.3
• /
• pp.828-832
• /
• 2012

We performed replica exchange molecular dynamics (REMD) simulations of the tripzip2 peptide (betahairpin) using the GB implicit and TI3P explicit solvation models. By comparing the resulting free energy surfaces of these two solvation model, we found that the GB solvation model produced a distorted free energy map, but the explicit solvation model yielded a reasonable free energy landscape with a precise location of the native structure in its global free energy minimum state. Our result showed that in particular, the GB solvation model failed to describe the tryptophan packing of trpzip2, leading to a distorted free energy landscape. When the GB solvation model is replaced with the explicit solvation model, the distortion of free energy shape disappears with the native-like structure in the lowest free energy minimum state and the experimentally observed tryptophan packing is precisely recovered. This finding indicates that the main source of this problem is due to artifact of the GB solvation model. Therefore, further efforts to refine this model are needed for better predictions of various aromatic side chain packing forms in proteins.

• Bulletin of the Korean Chemical Society
• /
• v.31 no.5
• /
• pp.1247-1251
• /
• 2010

We propose an improved solvent contact model to estimate the solvation free energies of amino acids from individual atomic contributions. The modification of the solvation model involves the optimization of three kinds of parameters in the solvation free energy function: atomic fragmental volume, maximum atomic occupancy, and atomic solvation parameters. All of these atomic parameters for 17 atom types are developed by the operation of a standard genetic algorithm in such a way to minimize the difference between experimental and calculated solvation free energies. The present solvation model is able to predict the experimental solvation free energies of amino acids with the squared correlation coefficients of 0.94 and 0.93 for the parameterization with Gaussian and screened Coulomb potential as the envelope functions, respectively. This result indicates that the improved solvent contact model with the newly developed atomic parameters would be a useful tool for the estimation of the molecular solvation free energy of a protein in aqueous solution.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.12
• /
• pp.1742-1750
• /
• 2003

In order to characterize the hydrophobic parameters of N-acetyl amino acid amides in 1-octanol/water, a theoretical calculation was carried out using a solvation free energy density model. The hydrophobicity parameters of the molecules are obtained with the consideration of the solvation free energy over the solvent volume surrounding the solute, using a grid model. Our method can account for the solvent accessible surface area of the molecules according to conformational variations. Through a comparison of the hydrophobicity of our calculation and that of other experimental/theoretical works, the solvation free energy density model is proven to be a useful tool for the evaluation of the hydrophobicity of amino acids and peptides. In order to evaluate the solvation free energy density model as a method of calculating the activity of drugs using the hydrophobicity of its building blocks, the contracture of Bradykinin potentiating pentapeptide was also predicted from the hydrophobicity of each residue. The solvation free energy density model can be used to employ descriptors for the prediction of peptide activities in drug discovery, as well as to calculate the hydrophobicity of amino acids.

• Bulletin of the Korean Chemical Society
• /
• v.26 no.4
• /
• pp.589-593
• /
• 2005

The solvation free energy of proton in methanol was calculated by B3LYP flavor of density functional calculations in combination with the Poisson-Boltzmann continuum solvation model. In order to check the adequacy of the computation level, the free energies of clustering in the gas phase were compared with the experimental data. The solvents were taken into account in a hybrid manner, i.e. one to five molecules of methanol were explicitly considered while other solvent molecules were represented with an implicit solvation model.

• Journal of the Korean Chemical Society
• /
• v.57 no.5
• /
• pp.554-559
• /
• 2013

In this study, we have investigated potential energy of ${\beta}$-D-glucopyranose in vacuum and implicit water condition. By Comparing two conditions we find that how solvation energy influence ${\beta}$-D-glucopyranose structure. We use AMBER package program and GLYCAM_06 force field. Solvation model was used for the generalized Born model with Hawkins, Cramer, Truhlar has been proposed. We conclude that difference of contour map of two conditions is caused by solvation effect by reducing hydrogen bonding interaction.

• Bulletin of the Korean Chemical Society
• /
• v.29 no.5
• /
• pp.979-984
• /
• 2008

In this study, 8 solutes (aniline, caffeine, p-cresol, ethyl benzene, methylparaben, phenol, pyridine, and toluene) have been tested in terms of linear solvation energy relationships (LSER). Several micellar liquid chromatography (MLC) systems using cationic surfactant cetyltrimethylammonium bromide (CTAB) and a mixture of water with (methanol, n-propanol, and n-butanol) modifiers were characterized using the LSER solvation parameter model. The effects of the surfactant and modifier concentration on the retention in MLC were discussed. LSER model had demonstrated high potential to predict retention factors with high squared correlation coefficients ($r^2$ > 0.99). A comparison of predicted and experimental retention factors suggests that LSER formalism is able to reproduce adequately the experimental retention factors of the solutes studied in the different experimental conditions investigated. This model is a helpful tool to understand the solute-surfactant interactions and evaluate the retention characteristic of micellar liquid chromatography.

• Bulletin of the Korean Chemical Society
• /
• v.29 no.4
• /
• pp.736-740
• /
• 2008

The interaction of myelin basic protein (MBP) from bovine central nervous system with divalent calcium ion was studied by isothermal titration calorimetry at 27 ${^{\circ}C}$ in aqueous solution. The extended solvation model was used to reproduce the enthalpies of $Co^{2+}$-MBP interaction over the whole $Co^{2+}$ concentrations. The solvation parameters recovered from the solvation model were attributed to the structural change of MBP due to the metal ion interaction. It was found that there is a set of three identical and noninteracting binding sites for $Co^{2+}$ ions. The association equilibrium constant is 0.015 ${\mu}M^{-1}$. The molar enthalpy of binding is $\Delta$H = −14.60 kJ $mol^{-1}$.

• Bulletin of the Korean Chemical Society
• /
• v.32 no.6
• /
• pp.1985-1992
• /
• 2011

The stable conformers of glycine and the inter-conversions between them were studied theoretically at various levels of theory, B3LYP, MP2, CCSD and CCSD(T), in the gas phase and in aqueous solution. In aqueous solution, the structures examined by use of the conductor-like polarizable continuum model (CPCM) with various cavity models, UA0, UAHF, UAKS, UFF, BONDI and PAULING, and by use of a discrete/continuum solvation model with eight water clusters. The Gibbs free energy differences between the neutral (NE) and zwitterionic conformers (ZW), ${\Delta}G_{Z-N}[=G_{ZW}-G_{NE}]$, in aqueous solution were well reproduced by using the BONDI and PAULING cavity models. However the ${\Delta}G_{Z-N}$ values were underestimated in other cavity models, although the ZW conformers existed as stable species in aqueous solution. In the studies of a discrete/continuum solvation model with eight water clusters, gas phase results are still insufficient to reproduce the experimental findings. However the ${\Delta}G_{Z-N}$ values calculated by use of CPCM method in aqueous solution agreed well with the experimental ones.

• Proceeding of EDISON Challenge
• /
• 2015.03a
• /
• pp.70-75
• /
• 2015

Fasciclin 1 (FAS1) is an extracellular protein whose aggregation in cornea leads to visual impairment. While a number of FAS1 mutants have been studied that exhibit enhanced/decreased aggregation propensity, no structural information has been provided so far that is associated with distinct aggregation potential. In this study, we have investigated the structural and thermodynamic characteristics of the wild-type FAS1 and its two mutants, R555Q and R555W, by using molecular dynamics simulations and three-dimensional reference interaction site model (3D-RISM) theory. We find that the hydrophobic solvent accessible surface area increases due to hydrophobic core repacking in the C-terminus caused by the mutation. We also find that the solvation free energy of the mutants increases due to the enhanced non-native H-bonding. These structural and thermodynamic changes upon mutation contribute to understand the aggregation of these mutants.

• Journal of the Korean Chemical Society
• /
• v.52 no.6
• /
• pp.608-613
• /
• 2008