• 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.

• 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.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 Magnetic Resonance Society
• /
• v.19 no.1
• /
• pp.1-10
• /
• 2015

Refinements with atomistic molecular dynamics (MD) simulation have contributed to improving the qualities of NMR structures. In most cases, the calculations with atomistic MD simulation for NMR structures employ generalized-Born implicit solvent model (GBIS) to take into accounts solvation effects. Developments in algorithms and computational capacities have ameliorated GBIS to approximate solvation effects that explicit solvents bring about. However, the quantitative comparison of NMR structures in the latest GBIS and explicit solvents is lacking. In this study, we report the direct comparison of NMR structures that atomistic MD simulation coupled with GBIS and water molecules refined. Two model proteins, GB1 and ubiquitin, were recalculated with experimental distance and torsion angle restraints, under a series of simulated annealing time steps. Whereas the root mean square deviations of the resulting structures were apparently similar, AMBER energies, the most favored regions in Ramachandran plot, and MolProbity clash scores witnessed that GBIS-refined structures had the better geometries. The outperformance by GBIS was distinct in the structure calculations with sparse experimental restraints. We show that the superiority stemmed, at least in parts, from the inclusion of all the pairs of non-bonded interactions. The shorter computational times with GBIS than those for explicit solvents makes GBIS a powerful method for improving structural qualities particularly under the conditions that experimental restraints are insufficient. We also propose a method to separate the native-like folds from non-violating diverged structures.