• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.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.

• YAKHAK HOEJI
• /
• v.39 no.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
• /
• v.11 no.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.

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

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

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

• Archives of Pharmacal Research
• /
• v.13 no.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
• /
• v.14 no.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.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.3
• /
• pp.116-121
• /
• 2005

In the present paper, a new calculation algorithm far solving large scale environmental or geophysical flows with free surface is proposed where the non-hydrostatic pressure component is taken into consideration. Predictor-corrector fractional step approach with explicit, forward time marching scheme in the sigma coordinate system is employed. In order to validate the present calculation algorithm and to estimate the effects of non-hydrostatic pressure on resultant flow and free surface movements, example calculations are carried out for typical steady and unsteady flow problems. Present method can be applied to the meso-scale free surface flows with complex bottom topography where MAC-like 3-d hydrodynamic calculations are quite ineffective and uneconomic.

• Bulletin of the Korean Chemical Society
• /
• v.33 no.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.