• International Journal of Aeronautical and Space Sciences
• /
• 제16권3호
• /
• pp.347-359
• /
• 2015

A method of describing the rovibrational energy transitions and coupled chemical reactions in the direct simulation Monte Carlo (DSMC) calculations is constructed for $H(^2S)+H_2(X^1{\Sigma}_g)$ and $He(^1S)+H_2(X^1{\Sigma}_g)$. First, the state-specific total cross sections for each rovibrational states are proposed to describe the state-resolved elastic collisions. The state-resolved method is constructed to describe the rotational-vibrational-translational (RVT) energy transitions and coupled chemical reactions by these state-specific total cross sections and the rovibrational state-to-state transition cross sections of bound-bound and bound-free transitions. The RVT energy transitions and coupled chemical reactions are calculated by the state-resolved method in various heat bath conditions without relying on a macroscopic properties and phenomenological models of the DSMC. In nonequilibrium heat bath calculations, the state-resolved method are validated with those of the master equation calculations and the existing shock-tube experimental data. In bound-free transitions, the parameters of the existing chemical reaction models of the DSMC are proposed through the calibrations in the thermochemical nonequilibrium conditions. When the bound-free transition component of the state-resolved method is replaced by the existing chemical reaction models, the same agreement can be obtained except total collision energy model.

• 약학회지
• /
• 제39권3호
• /
• pp.329-336
• /
• 1995

Conformational free energy calculations using an empirical potential function and a hydration shell model(program CONBIO) were carried out on hypoglycemic agent acetohexamide and tolazamide in the unhydrated and hydrated states. The initial geometry of sulfonylureas was obtained from X-ray crystallographieal data and homologous molecular fragments. In both states, the feasible conformations were obtained from the calculations of conformational energy, conformational entropy, and hydration free energy by varying all the torsion angles of the molecules. From the calculation results, it is known that the conformations] entropy is the major contribution to stabflize the low-free-energy conformations of two sulfonylureas in both states. But, in hydrated state, the hydration does not directly affect each conformations. The intramolecular hydrogen bonding of sulfonylurea hydrogen and 7-membered nitrogen appeared to the conformations of tolazamide in both states. It is thought that the hydrogen bonding decrease steric hindrance on the receptor binding direction. The substitution of alicyclic or N-heterocyclic ring than that of carbons chain of urea moiety may be properly interaction between sulfonylureas and the putative pancreatic receptor.

• Bulletin of the Korean Chemical Society
• /
• 제11권2호
• /
• pp.144-149
• /
• 1990

Conformational free energy calculations using an empirical potential function (ECEPP/2) and the hydration shell model were carried out on the sulfonylurea herbicides of bensulfuron methyl (Londax) and metsulfuron methyl (Ally). The conformational energy was minimized from starting conformations which included possible combinations of torsion angles in the molecule. The conformational entropy of each conformation was computed using a harmonic approximation. To understand the hydration effect on the conformation of the molecule in aqueous solution, the hydration free energy of each group was calculated and compared each other. It was found that the low-free-energy conformations of two molecules in aqueous solution prefer the overall folded structure, in which an interaction between the carbonyl group of ester in aryl ring and the first amido group of urea bridge plays an important role. From the analysis of total free energy, the hydration and conformational entropy are known to be essential in stabilizing low-free-energy conformations of Londax, whereas the conformational energy is proved to be a major contribution to the total free energy of low-free-energy conformations of Ally.

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

• Bulletin of the Korean Chemical Society
• /
• 제35권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.

• Bulletin of the Korean Chemical Society
• /
• 제26권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.

• Archives of Pharmacal Research
• /
• 제13권1호
• /
• pp.43-50
• /
• 1990

Conformational free energy calculations using an empirical potential function (ECEPP/2) and hydration shell model were carried out on the four-4-aminodiphenyl sulfone analogues of 4, 4'-diamino-2' methyldiphenyl sulfone, 4, 2', 4-triaminodiphenyl sulfone, 4, 4'-diaminodiphenyl sulfone, and 4-aminodiphenyl sulfone as antibacterial agents on Mycobacterium lufu. The conformational energy was minimized from starting conformations which included possible combinations of torsion angles in the molecule. The conformational entropy change of each conformation was computed using a harmonic approximation. To understand the hydration effect on the conformation of the molecule in aqueous solution, the contributions of water-accessible volume and the hydration free energy of each group or atom in the lowest-free-energy conformation was calculated and compared each other. From comparison of the computed lowest-free-energy conformations of four analogues with their antibacterial activities, it is known that the conformation and the hydrophobicity of sulfonyl group and its adjacent carbon atom in each compound are the essential factors to show the strong antibacterial activity.

• Bulletin of the Korean Chemical Society
• /
• 제14권2호
• /
• pp.262-265
• /
• 1993

Conformational free energy calculations using an empirical potential function and a hydration shell model (program CONBIO) were carried out on antifungal agent fluconazole in the unhydrated and hydrated states. The initial geometry of fluconazole was obtained from two minimized fragments of it using a molecular mechanics MMPMI and followed by minimizing with a semiempirical AM1 method. In both states, the feasible conformations were obtained from the calculations of conformational energy, conformational entropy, and hydration free energy by varying all the torsion angles of the molecule. The intramolecular hydrogen bonds of isopropyl hydroxyl hydrogen and triazole nitrogens and the structural flexibility are of significant importance in stabilizing the conformations of fluconazole in both states. Hydration is proved to be one of the essential factors in stabilizing the overall conformation in aqueous solution. Two F atoms of phenyl ring are not identified as an essential key in determining the stable conformations and may be responsible for the interaction with the receptor of fluconazole.

• 한국해양환경ㆍ에너지학회지
• /
• 제8권3호
• /
• pp.116-121
• /
• 2005

본 연구에서는 자유표면을 갖는 대규모의 환경유동의 해를 구하는데 있어 비정수압의 효과를 고려하는 새로운 계산 기법을 제시하였다. Sigma 좌표계에서 시간은 전진차분으로, 공간은 중간차분 및 풍상차분을 혼합하여 2단계 해법을 도입하였다. 계산기법의 타당성과 유효성을 검증하기 위하여, 자유표면을 갖는 정상유동과 비정상유동에 대한 전형적인 예를 설정하여 정수압 계산과 비정수압 계산을 수행하고 이들을 비교하였다. 계산결과 복잡한 해저지형을 갖는 자유표면 유동에 있어, 지정수압효과가 무시할 수 없는 유동영역이 존재함이 입증되었으며, MAC기법과 같은 3차원 유체동역학적 기법을 적용하기에는 비경제적인 문제에 대해 본 계산기법이 매우 유효하게 적용될 수 있다는 가능성이 입증되었다.

• Bulletin of the Korean Chemical Society
• /
• 제33권3호
• /
• pp.911-916
• /
• 2012

Ion selectivity in a simple cyclic peptide nanotube, composed of four cyclo[-(D-Ala-Glu-D-Ala-Gln)$_2-$] units, is investigated by calculating the PMF profiles of $Na^+$, $K^+$, and $Cl^-$ ions permeating through the peptide nanotube in water. The final PMF profiles of the ions obtained from the umbrella sampling (US) method show an excellent agreement with those from the thermodynamic integration (TI) method. The PMF profiles of $Na^+$ and $K^+$ display free energy wells while the PMF curve of $Cl^-$ features free energy barriers, indicating the selectivity of the cyclic peptide nanotube to cations. Decomposition of the total mean force into the contribution from each component in the system is also accomplished by using the TI method. The mean force decomposition profiles of $Na^+$ and $K^+$ demonstrate that the dehydration free energy barriers by water molecules near the channel entrance and inside the channel are completely compensated for by attractive electrostatic interactions between the cations and carbonyl oxygens in the nanotube. In the case of $Cl^-$, the dehydration free energy barriers are not eliminated by an interaction between the anion and the peptide nanotube, leading to the high free energy barriers in the PMF profile. Calculations of the coordination numbers of the ions with oxygen atoms pertaining to either water molecules or carbonyl groups in the peptide nanotube reveal that the stabilization of the cations in the midplane regions of the nanotube arises from the favorable interaction of the cations with the negatively charged carbonyl oxygens.