• Journal of the Korean Chemical Society
• /
• v.56 no.2
• /
• pp.195-200
• /
• 2012

In this study we present a new approach for computing the partial atomic charge derived from the wavefunctions of molecules. This charge, which we call the "y_charge", was calculated by taking into account the energy level and orbital populations in each molecular orbital (MO). The charge calculations were performed in the software, which was developed by us, developed using the C# programming language. Partial atomic charges cannot be calculated directly from quantum mechanics. According to a partitioning function, the electron density of constituent molecular atoms depends on the electrostatic attraction field of the nucleus. Taking into account the Boltzmann population of each MO as a function of its energy and temperature we obtain a formula of partial charges.

• Bulletin of the Korean Chemical Society
• /
• v.16 no.10
• /
• pp.915-923
• /
• 1995

The parameters for an empirical net atomic charge calculation method, Modified Partial Equalization of Orbital Electronegativity (MPEOE), were determined for the atoms in organosilicon compounds and zeolites. For the organosilicon family, the empirical parameters were determined by introducing both experimental and ab initio observables as constraints, these are the experimental and ab initio dipole moments, and the ab initio electrostatic potential of the organosilicon molecules. The Mulliken population was also introduced though it is not a quantum mechanical observable. For the parameter optimization of the atoms in the aluminosilicates, the dipole moments and the electrostatic potentials which calculated from the 6-31G** ab initio wave function were used as constraints. The empirically calculated atomic charges of the organosilicons could reproduce both the experimental and the ab inito dipole moments well. The empirical atomic charges of the aluminosilicates could reproduce the ab initio electrostatic potentials well also.

• Journal of Integrative Natural Science
• /
• v.3 no.4
• /
• pp.226-230
• /
• 2010

Calculation of partial charge is important in chemistry. However, because there are many methods developed, it is of considerable interest to know how to calculate and apply properly to address various chemical problems. For basis set, usually double zeta quality is acceptable, and double zeta polarization function would be enough for most cases. To describe electronic state more accurately, Many electron configurations would be necessary to describe highly strained or anionic species. The NPA population introduced new concept about amide bonds, i.e., the planar geometry of nitrogen atom may not come from resonance, but from the lowering of p-orbital energy by electronegative carbonyl carbon atom. The issues for hypervalent atomic charges was also addressed by various charge derivation scheme. When the charge schemes were applied to organolithium compounds, the ionic nature of boding was revealed. This comes from the fact that previous Mulliken partial atomic charges overemphasized the covalent character, wihout much justification. The other partial charge derivation schemes such as NPA(natural population analysis), IPP (Integrated Projected Population) showed that much more ionic picture. ESP potential derived charges are generally believed to be suitable to describe intermolecular interactions, therefore they are used for molecular dynamics simulations and CoMFA (comparative molecular field analysis). The charge derivation schemes using multipole polarization was mainly applied to reproduce experimental infrared spectroscopy. In some reports these schemes are also suitable for intermecular electrostatic interactions. Charges derived from electron density gradient have shown the some bonds are not straight, but actually bent. The proper choice of charge-calculation method along with suitable level of theory and basis set are briefly discussed.

• Journal of Integrative Natural Science
• /
• v.3 no.4
• /
• pp.231-236
• /
• 2010

Partial charge is an important and fundamental concept which can explain many aspects of chemistry. Since a molecule can be regarded as neclei surrounded by electron cloud, there is no way to define a partial charge accurately. Nevertheless, there have been many attempts to define these seemingly impossible parameters, since they would facilitate the understanding of molecular properties such as molecular dipole moment, solvation, hydrogen bonding, molecular spectroscopy, chemical reaction, etc. Common methods are based on the charge equalization, orbital occupancy, charge density, and electric multipole moments, and electrostatic potential fitting. Methods based on the charge equalization using electronegativity are very fast, and therefore they have been used to study many compounds. Methods to subdivide orbital occupancy using basis set conversion, relies on the notion that molecular orbitals are composed of atomic orbitals. The main idea is to reduce overlap integral between two nuclei using converted orthogonal basis sets. Using some quantum mechanical observables like electrostatic potential or charge multipole moments. Using potential grids obtained from wavefunction, partial charges can be fitted. these charges are most useful to describe intermolecular electrostatic interactions. Methods to using dipole moment and its derivatives, seems to be sensitive the level of theory, Dividing electron density using density gradient being the most rigorous theoretically among various schemes, bears best potential to describe the charge the most adequately in the future.

• Bulletin of the Korean Chemical Society
• /
• v.24 no.3
• /
• pp.369-376
• /
• 2003

Coulson (ZINDO), Mulliken $(MP2/6-31G^*)$ and Natural $(MP2/6-31G^*)$ population analyses of several large molecules were performed by the Fragment Reassociation (FR) method. The agreement between the conventional ZINDO (or conventional MP2) and FR-ZINDO (or FR-MP2) charges of these molecules was excellent. The standard deviations of the FR-ZINDO net atomic charges from the conventional ZINDO net atomic charges were 0.0008 for $C_{10}H_{22}$ (32 atoms), 0.0012 for $NH_2-C_{16}O_2H_{28}-COOH$ (53 atoms), 0.0014 for $NH_3^+-C_{16}O_2H_{28}-COOH$ (54 atoms), 0.0017 for $NH_2-C_{16}O_2H_{28}-COO^-$ (52 atoms), 0.0019 for $NH_3^+-C_{16}O_2H_{28}-COO^-$ (53 atoms), 0.0024 for a conjugated model $(O=CH-(CH=CH)_{15}-C=O-(CH=CH)_{12}-CH=CH_2)$, 118 atoms), 0.0038 for aglycoristocetin $(C_{60}N_7O_{19}H_{52}^+$, 138 atoms), 0.0023 for a polypropylene model complexed with a zirconocene catalyst $(C_{68}H-{121}Zr^+$, 190 atoms) and 0.0013 for magainin $(C_{112}N_{29}O_{28}SH_{177}$, 347 atoms), respectively. The standard deviations of the FR-MP2 Mulliken (or Natural) partial atomic charges from the conventional ones were 0.0016 (or 0.0016) for $C_{10}H_{22}$, 0.0019 (or 0.0018) for $NH_2-C_{16}O_2H_{28}-COOH$ and 0.0033 (or 0.0023) for $NH_3^+-C_{16}O_2H_{28}-COO^-$, respectively. These errors were attributed to the shape of molecules, the choice of fragments and the degree of ionic characters of molecules as well as the choice of methods. The CPU time of aglycoristocetin, conjugated model, polypropylene model complexed with zirconocene and magainin computed by the FR-ZINDO method was respectively 2, 4, 6 and 21 times faster than that by the normal ZINDO method. The CPU time of $NH_2-C_{16}O_2H_{28}-COOH\;and\;NH_3^+-C_{16}O_2H_{28}-COO^-$ computed by the FR-MP2 method was, respectively, 6 and 20 times faster than that by the normal MP2 method. The largest molecule calculated by the FR-ZINDO method was B-DNA (766 atoms). These results will enable us to compute atomic charges of huge molecules near future.

• Journal of Pharmaceutical Investigation
• /
• v.26 no.3
• /
• pp.155-168
• /
• 1996

The active site of cyclooxygenase was modeled by complementary receptor-cavity mapping procedure using 3D structures of the non-steroidal antiinflammatory drugs (NSAIDs). A total of 50 NSAIDs were chosen as data ligands which compete the same site on the enzyme. Partial atomic charges were estimated, and the energetic differences for various conformations were calculated so as to meet the need for a most efficient overlapping of the probably-equivalent functional groups of the ligand molecules. The structure activity relationships of the NSAIDs, if available, were fully considered throughout the modeling. The overall shape of the model obtained is similar to a boot-without-bottom. Most of inner surface of the cavity appeared as hydrophobic; two polar counterparts except the carboxyl-binding position were found. By this model, some clear explanations could be given on the experimental observations which were not satisfiably understood yet.

• Journal of Integrative Natural Science
• /
• v.4 no.2
• /
• pp.99-102
• /
• 2011

Effect of protonation and deprotonation of ionization compounds is an important application in Comparative molecular field analysis (CoMFA). There are enough information's were reported about different CoMFA applications such as Series design and selection of training set, Geometries and optimizations of molecules, Effect of partial atomic charges, bioactive conformations and alignment, Interaction energy fields, Effects of different grid spacing etc. However limited information's are available about the ionization of compounds. This study aimed at the critical review of about the effects of protonation of ionizable molecules and its impact on the predictability of CoMFA models. We also discussed about previous implications and the things needed to be considered to come for a final conclusion about its impact on CoMFA predictability.

• Journal of Integrative Natural Science
• /
• v.7 no.1
• /
• pp.45-49
• /
• 2014

The (c-Jun N-terminal kinase 3) JNK3 is a potential therapeutic target for various neurological disorders. Here, a three dimensional quantitative structure-activity relationship (3D-QSAR) study on phenoxypyridine as JNK3 inhibitors was performed to rationalize the structural requirements responsible for the inhibitory activity of these compounds. The comparative molecular field analysis (CoMFA) using different partial atomic charges, was employed to understand the structural factors affecting JNK3 inhibitory potency. The Gasteiger-Marsili yielded a CoMFA model with cross-validated correlation coefficient ($q^2$) of 0.54 and non-cross-validated correlation coefficient ($r^2$) of 0.93 with five components. Furthermore, contour maps suggested that bulky substitution with oxygen atom in $R^3$ position could enhance the activity considerably. The work suggests that further chemical modifications of the compounds could lead to enhanced activity and could assist in the design of novel JNK3 inhibitors.

• Bulletin of the Korean Chemical Society
• /
• v.32 no.5
• /
• pp.1604-1612
• /
• 2011

We developed three-dimensional quantitative structure activity relationship (3D-QASR) models for 17 adamantyl derivatives as P-glycoprotein (P-gp) inhibitors. Eighteen different partial charge calculation methods were tested to check the feasibility of the 3D-QSAR models. Best predictive comparative molecular field analysis (CoMFA) model was obtained with the Austin Model 1-Bond Charge Correction (AM1-BCC) atomic charge. The 3D-QSAR models were derived with CoMFA and comparative molecular similarity indices analysis (CoMSIA). The final CoMFA model ($q^2$ = 0.764, $r^2$ = 0.988) was calculated with an AM1-BCC charge and electrostatic parameter, whereas the CoMSIA model ($q^2$ = 0.655, $r^2$ = 0.964) was derived with an AM1-BCC charge and combined steric, electrostatic, hydrophobic and HB-acceptor parameters. Leave-five-out (LFO) cross-validation was also performed, which yielded good correlation coefficient for both CoMFA (0.801) and CoMSIA (0.656) models. Robustness of the developed models was checked further with 1000 run bootstrapping analyses, which gave an acceptable correlation coefficient for CoMFA (BS-$r^2$ = 0.997, BS-SD = 0.003) and CoMSIA (BS-$r^2$ = 0.996, BS-SD = 0.018).

• Bulletin of the Korean Chemical Society
• /
• v.31 no.10
• /
• pp.2903-2908
• /
• 2010

Protein S-nitrosation is common in cells under nitrosative stress. In order to model proteins with S-nitrosocysteine (CysSNO) residues, we first developed an Amber force field for S-nitrosoethanethiol (EtSNO) and then transferred it to CysSNO. Partial atomic charges for EtSNO and CysSNO were obtained by a restrained electrostatic potential approach to be compatible with the Amber-99 force field. The force field parameters for bonds and angles in EtSNO were obtained from a generalized Amber force field (GAFF) by running the Antechamber module of the Amber software package. The GAFF parameters for the CC-SN and CS-NO dihedrals were not accurate and thus determined anew. The CC-SN and CS-NO torsional energy profiles of EtSNO were calculated quantum mechanically at the level of B3LYP/cc-pVTZ//HF/6-$31G^*$. Torsional force constants were obtained by fitting the theoretical torsional energies with those obtained from molecular mechanics energy minimization. These parameters for EtSNO reproduced, to a reasonable accuracy, the corresponding torsional energy profiles of the capped tripeptide ACE-CysSNO-NME as well as their structures obtained from quantum mechanical geometry optimization. A molecular dynamics simulation of myoglobin with a CysSNO residue produced a well-behaved trajectory demonstrating that the parameters may be used in modeling other S-nitrosated proteins.