Dataset Information

Development of a polarizable intermolecular potential function (PIPF) for liquid amides and alkanes.

ABSTRACT: A polarizable intermolecular potential function (PIPF) employing the Thole interacting dipole (TID) polarization model has been developed for liquid alkanes and amides. In connection with the internal bonding terms of the CHARMM22 force field, the present PIPF-CHARMM potential provides an adequate description of structural and thermodynamic properties for liquid alkanes and for liquid amides through molecular dynamics simulations. The computed heats of vaporization and liquid density are within 1.4% of experimental values. Polarization effects play a major role in liquid amides, which are reflected by an increase of 1.5 to 1.8 D in molecular dipole moment for primary and secondary amides. Furthermore, the computed polarization energies contribute to the total intermolecular interaction energy by 6 to 24%. The ability of the PIPF-CHARMM force field to treat protein backbone structures is tested by examining the potential energy surface of the amide bond rotation in N-methylacetamide and the Ramachandran surface for alanine dipeptide. The agreement with ab initio MP2 results and with the original CHARMM22 force field is encouraging, suggesting the PIPF-CHARMM potential can be used as a starting point to construct a complete polarizable force field for proteins.

SUBMITTER: Xie W

PROVIDER: S-EPMC2572772 | biostudies-other | 2007

REPOSITORIES: biostudies-other

ACCESS DATA

Publications

Development of a polarizable intermolecular potential function (PIPF) for liquid amides and alkanes.

Xie Wangshen W Pu Jingzhi J Mackerell Alexander D AD Gao Jiali J

Journal of chemical theory and computation 20070101 6

A polarizable intermolecular potential function (PIPF) employing the Thole interacting dipole (TID) polarization model has been developed for liquid alkanes and amides. In connection with the internal bonding terms of the CHARMM22 force field, the present PIPF-CHARMM potential provides an adequate description of structural and thermodynamic properties for liquid alkanes and for liquid amides through molecular dynamics simulations. The computed heats of vaporization and liquid density are within ...[more]

PMID: 18958290

Similar Datasets

Project description:We report a set of rare copper-catalyzed reactions of alkanes with simple amides, sulfonamides, and imides (i.e., benzamides, tosylamides, carbamates, and phthalimide) to form the corresponding N-alkyl products. The reactions lead to functionalization at secondary C-H bonds over tertiary C-H bonds and even occur at primary C-H bonds. [(phen)Cu(phth)] (1-phth) and [(phen)Cu(phth)2] (1-phth2), which are potential intermediates in the reaction, have been isolated and fully characterized. The stoichiometric reactions of 1-phth and 1-phth2 with alkanes, alkyl radicals, and radical probes were investigated to elucidate the mechanism of the amidation. The catalytic and stoichiometric reactions require both copper and tBuOOtBu for the generation of N-alkyl product. Neither 1-phth nor 1-phth2 reacted with excess cyclohexane at 100 °C without tBuOOtBu. However, the reactions of 1-phth and 1-phth2 with tBuOOtBu afforded N-cyclohexylphthalimide (Cy-phth), N-methylphthalimide, and tert-butoxycyclohexane (Cy-OtBu) in approximate ratios of 70:20:30, respectively. Reactions with radical traps support the intermediacy of a tert-butoxy radical, which forms an alkyl radical intermediate. The intermediacy of an alkyl radical was evidenced by the catalytic reaction of cyclohexane with benzamide in the presence of CBr4, which formed exclusively bromocyclohexane. Furthermore, stoichiometric reactions of [(phen)Cu(phth)2] with tBuOOtBu and (Ph(Me)2CO)2 at 100 °C without cyclohexane afforded N-methylphthalimide (Me-phth) from ?-Me scission of the alkoxy radicals to form a methyl radical. Separate reactions of cyclohexane and d12-cyclohexane with benzamide showed that the turnover-limiting step in the catalytic reaction is the C-H cleavage of cyclohexane by a tert-butoxy radical. These mechanistic data imply that the tert-butoxy radical reacts with the C-H bonds of alkanes, and the subsequent alkyl radical combines with 1-phth2 to form the corresponding N-alkyl imide product.

Project description:In the companion paper, we presented a set of induced dipole interaction models using four types of screening functions, which include the Applequist (no screening), the Thole linear, the Thole exponential model, and the Thole Tinker-like (another form of exponential screening function) functions. In this work, we evaluate the performance of polarizability models using a large set of amino acid analog pairs in conformations that are frequently observed in protein structures as a benchmark. For each amino acid pair, we calculated quantum mechanical interaction energies at the MP2/aug-cc-pVTZ//MP2/6-311++G(d,p) level with the basis set superposition error (BSSE) correction and compared them with molecular mechanics results. Encouragingly, all polarizable models significantly outperform the additive F94 and F03 models (mimicking AMBER ff94/ff99 and ff03 force fields, respectively) in reproducing the BSSE-corrected quantum mechanical interaction energies. In particular, the root-mean-square errors (RMSEs) for three Thole models in Set A (where the 1-2 and 1-3 interactions are turned off and all 1-4 interactions are included) are 1.456, 1.417, and 1.406 kcal/mol for model AL (Thole Linear), model AE (Thole exponential), and model AT (Thole Tinker-like), respectively. In contrast, the RMSEs are 3.729 and 3.433 kcal/mol for F94 and F03 models, respectively. A similar trend was observed for the average unsigned errors (AUEs), which are 1.057, 1.025, 1.011, 2.219, and 2.070 kcal/mol for AL, AE, AT, F94/ff99, and F03, respectively. Analyses based on the trend line slopes indicate that the two fixed charge models substantially underestimate the relative strengths of noncharge-charge interactions by 24 (F03) and 35% (F94), respectively, whereas the four polarizable models overestimate the relative strengths by 5 (AT), 3 (AL, AE), and 13% (AA), respectively. Agreement was further improved by adjusting the van der Waals parameters. Judging from the notably improved accuracy in comparison with the fixed charge models, the polarizable models are expected to form the foundation for the development of high quality polarizable force fields for protein and nucleic acid simulations.

Dataset Information

Development of a polarizable intermolecular potential function (PIPF) for liquid amides and alkanes.

Publications

Development of a polarizable intermolecular potential function (PIPF) for liquid amides and alkanes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets