• Journal of the Korean Magnetic Resonance Society
• /
• v.26 no.1
• /
• pp.1-9
• /
• 2022

There are two distinct approaches to improving the quality of protein NMR structures during refinement: all-atom force fields and accumulated knowledge-assisted methods that include Rosetta. Mao et al. reported that, for 40 proteins, Rosetta increased the accuracies of their NMR-determined structures with respect to the X-ray crystal structures (Mao et al., J. Am. Chem. Soc. 136, 1893 (2014)). In this study, we calculated 32 structures of those studied by Mao et al. using all-atom force field and implicit solvent model, and we compared the results with those obtained from Rosetta. For a single protein, using only the experimental NOE-derived distances and backbone torsion angle restraints, 20 of the lowest energy structures were extracted as an ensemble from 100 generated structures. Restrained simulated annealing by molecular dynamics simulation searched conformational spaces with a total time step of 1-ns. The use of GPU-accelerated AMBER code allowed the calculations to be completed in hours using a single GPU computer-even for proteins larger than 20 kDa. Remarkably, statistical analyses indicated that the structures determined in this way showed overall higher accuracies to their X-ray structures compared to those refined by Rosetta (p-value < 0.01). Our data demonstrate the capability of sophisticated atomistic force fields in refining NMR structures, particularly when they are coupled with the latest GPU-based calculations. The straightforwardness of the protocol allows its use to be extended to all NMR structures.

• Proceedings of the KIEE Conference
• /
• 2015.07a
• /
• pp.742-743
• /
• 2015

Magnetic couplings are do not require any mechanical contact with the power transmitted to the secondary side according to the primary side. For this reason, well-suited for isolated systems such as vacuums or high pressure. So, this paper presents the force characteristic analysis of the permanent magnet (PM) linear coupling with vertical magnetized using an analytical magnetic field calculations. Based on the definition of governing equations and magnetic vector potential, we obtained the analytical solutions according to the boundary condition for each of the regions. Also, we derived from the force generated in the permanent magnet surface using the Maxwell stress tensor. The analytical results are proved the validity by comparing to the finite element method (FEM).

• Proceedings of the Korean Vacuum Society Conference
• /
• 2010.08a
• /
• pp.93-93
• /
• 2010

In this talk, I will introduce a reactive force field (ReaxFF) molecular dynamics (MD) simulation. In contrast to common MD simulations with empirical FFs, we can predict chemical reactions (bond breaking and formation) in large scale systems with the ReaxFF simulation where all of the ReaxFF parameters are from quantum mechanical calculations such as density functional theory to provide high accuracy. Accordingly, the ReaxFF simulation provides both accuracy of quantum mechanical calculations and description of large scale systems of atomistic simulations at the same time. Here, I will first discuss a theory in the ReaxFF including the differences from other empirical FFs, and then show several applications for studying chemical reactions of SiHx radicals on Si surfaces, which is an important issue in Si process.

• 한국신재생에너지학회:학술대회논문집
• /
• 2010.11a
• /
• pp.113.2-113.2
• /
• 2010

We report that the novel covalent organic frameworks (COFs) are capable of reversibly providing an extremely high uptake capacity of carbon dioxide and hydrogen at room temperature. These COFs are designed based on the multiscale simulations approach via the combination of ab initio calculations and force-field calculations. For this goal, we explore the adsorption sites of carbon dioxide and hydrogen on COFs, their porosity, as well as carbon dioxide adsorption isotherms. We identify the binding sites and energies of $CO_2$ on COFs using ab initio calculations and obtain the carbon dioxide adsorption isotherms using grand canonical ensemble Monte Carlo calculations. Moreover, the calculated adsorption isotherms are compared with the experimental values in order to build the reference model in describing the interactions between the $CO_2/H_2$ and the COFs and in predicting the $CO_2$ and $H_2$ adsorption isotherms of COFs. Finally, we design three new COFs, 2D COF-05, 3D COF-05 (ctn), and 3D COF-05 (bor), for the high capacity $CO_2/H_2$ and $H_2$ storage.

• Bulletin of the Korean Chemical Society
• /
• v.23 no.2
• /
• pp.253-261
• /
• 2002

Vibrational properties of ferrocyanide complex ion, $[Fe(CN)_6]^{4-}$ , have been studied based on the force constants obtained from the density functional calculations at B3LYP/$6-31G^{\ast\ast}$ level by means of the normal mode analysis using new bond angle and linear angle internal coordinates recently developed. Vibrations of ferrocyanide were manipulated by twenty-three symmetry force constants. The angled bending deformations of C-Fe-C, the linear bending deformations of Fe-C${\equiv}$N and the stretching vibrations of Fe-C have been quantitatively assigned to the calculated frequencies. The force constants in the internal coordinates employed in the modified Urey-Bradley type potential were evaluated on the density functional force field applied, and better interaction force constants in the internal coordinates have been proposed. The valence force constants in the general quadratic valence force field were also given. The stretch-stretch interaction and stretch-bending interaction constants are not sensitive to the geometrical displacement in the valence force field.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.6
• /
• pp.824-828
• /
• 2003

Our previous quantum mechanical calculations using polyalanine model systems showed that the monodipolemacrodipoleinteractions selectively stabilize α-helices and make it possible for α-helices to be formed inhydrophobic environment where the solvent effect is not available. The monodipole-macrodipole interactionsin α-helices were studied molecular mechanically using various point-charge systems available. The resultsshow that all the point-charge systems used in the calculations produce the monodipole-macrodipoleinteractions up to about 60% compared to the results of the quantum mechanical calculations. The results ofmolecular mechanical calculations are explained and discussed compared to the results of the quantummechanical calculations.

• Journal of the Korean Magnetic Resonance Society
• /
• v.17 no.1
• /
• pp.11-18
• /
• 2013

Rapid advances of computational power and method have made it practical to apply the time-consuming calculations with all-atom force fields and sophisticated potential energies into refining NMR structure. Added to the all-atom force field, generalized-Born implicit solvent model (GBIS) contributes substantially to improving the qualities of the resulting NMR structures. GBIS approximates the effects that explicit solvents bring about even with fairly reduced computational times. Although GBIS is employed in the final stage of NMR structure calculation with experimental restraints, the effects by GBIS on structures have been reported notable. However, the detailed effect is little studied in a quantitative way. In this study, we report GBIS refinements of ubiquitin and GB1 structures by six GBIS models of AMBER package with experimental distance and backbone torsion angle restraints. Of GBIS models tested, the calculations with igb=7 option generated the closest structures to those determined by X-ray both in ubiquitin and GB1 from the viewpoints of root-mean-square deviations. Those with igb=5 yielded the second best results. Our data suggest that the degrees of improvements vary under different GBIS models and the proper selection of GBIS model can lead to better results.

• Bioinformatics and Biosystems
• /
• v.1 no.2
• /
• pp.103-108
• /
• 2006

Computational prediction of protein-ligand binding has been widely used as a tool to discover lead compounds fur new drugs. Prediction accuracy is determined in part by the scoring function used in docking calculations. Diverse scoring functions are available, and these can be classified into force-field based, empirical, and knowledge-based functions depending upon the basic assumptions made in development. Among these, force-field based functions consider physical interactions the most in detail. However, the force-field based functions have the drawback of not including the entropic effect while considering only the energy contribution such as dispersion or electrostatic forces. In this article, a method to take into account of the entropic effect using the colony energy is suggested when force-field based scoring functions is used by extracting conformational information obtained from the pre-existing docking program. An improved result for decoy discrimination is illustrated when the method is applied to the DOCK scoring function, and this implies that more accurate docking calculation is possible.

• Interdisciplinary Bio Central
• /
• v.2 no.2
• /
• pp.6.1-6.5
• /
• 2010

Fragment molecular orbital (FMO) method provides a novel tool for ab initio calculations of large biomolecules. This method overcomes the size limitation difficulties in conventional molecular orbital methods and has several advantages compared to classical force field approaches. While there are many features in this method, we here focus on explaining the issues related to protein-ligand binding: FMO method provides useful interaction-analysis tools such as IFIE, CAFI and FILM. FMO calculations can provide not only binding energies, which are well correlated with experimental binding affinity, but also QSAR descriptors. In addition, FMO-derived charges improve the descriptions of electrostatic properties and the correlations between docking scores and experimental binding affinities. These calculations can be performed by the ABINIT-MPX program and the calculation results can be visualized by its proper BioStation Viewer. The acceleration of FMO calculations on various computer facilities is ongoing, and we are also developing methods to deal with cytochrome P450, which belongs to the family of drug metabolic enzymes.

• Bulletin of the Korean Chemical Society
• /
• v.6 no.2
• /
• pp.57-60
• /
• 1985

Normal mode calculations about the extended double six ring (D6R) subunit cluster model of Faujasite-type zeolites have been done by using the valence force field with reasonably adjusted force constants. We have studied for four X, Y zeolites species varying in M$_{f}$(Al/Si + Al) values. The calculated characteristic frequencies of D6R mode (${\nu}_{D_{6}R}$) and the rate of change of ${\nu}_{D_{6}R}$ with the mole fractions of aluminum, M$_{f}$ (Al/Si + Al) values agree well with Flanigen's experimental data; the experimental slope is -79 cm$^{-1}$, while the calculated slope is -82 cm$^{-1}$. Those are the improved results as compared to Blackwell's theoretical study; his predicted slope is -94 cm$^{-1}$.