Computer Simulation of Molecular Dynamics: Methodology,Applications, and Perspectives in Chemistry

During recent decades it has become feasible to simulate the dynamics of molecular systems on a computer. The method of molecular dynamics (MD) solves Newton's equations of motion for a molecular system, which results in trajectories for all atoms in the system. From these atomic trajectories a variety of properties can be calculated. The aim of computer simulations of molecular systems is to compute macroscopic behavior from microscopic interactions. The main contributions a microscopic consideration can offer are (1) the understanding and (2) interpretation of experimental results, (3) semiquantitative estimates of experimental results, and (4) the capability to interpolate or extrapolate experimental data into regions that are only difficultly accessible in the laboratory. One of the two basic problems in the field of molecular modeling and simulation is how to efficiently search the vast configuration space which is spanned by all possible molecular conformations for the global low (free) energy regions which will be populated by a molecular system in thermal equilibrium. The other basic problem is the derivation of a sufficiently accurate interaction energy function or force field for the molecular system of interest. An important part of the art of computer simulation is to choose the unavoidable assumptions, approximations and simplifications of the molecular model and computational procedure such that their contributions to the overall inaccuracy are of comparable size, without affecting significantly the property of interest. Methodology and some practical applications of computer simulation in the field of (bio)chemistry will be reviewed.     

分子力学在有机化学中的应用

本文评述了近年来分子力学在有机化学中的应用,特别是在有机化学反应和核磁共振中的应用。     

Modeling of peptides containing D-amino acids: implications on cyclization

Cyclic peptides are therapeutically attractive due to their high bioavailability, potential selectivity, and scaffold novelty. Furthermore, the presence of D-residues induces conformational preferences not followed by peptides consisting of naturally abundant L-residues. Therefore, comprehending how amino acids induce turns in peptides, subsequently facilitating cyclization, is significant in peptide design. Here, we performed 20-ns explicit-solvent molecular dynamics simulations for three diastereomeric peptides with stereochemistries: LLLLL, LLLDL, and LDLDL. Experimentally LLLLL and LDLDL readily cyclize, whereas LLLDL cyclizes in low yield. Simulations at 310 K produced conformations with inter-terminal hydrogen bonds that correlated qualitatively with the experimental cyclization trend. Energies obtained for representative structures from quantum chemical (B3LYP/PCM/cc-pVTZ//HF/6-31G*) calculations predicted pseudo-cyclic and extended conformations as the most stable for LLLLL and LLLDL, respectively, in agreement with the experimental data. In contrast, the most stable conformer predicted for peptide LDLDL was not a pseudo-cyclic structure. Moreover, D-residues preferred the experimentally less populated α_L rotamers even when simulations were performed at a higher temperature and with strategically selected starting conformations. Energies calculated with molecular mechanics were consistent only with peptide LLLLL. Thus, the conformational preferences obtained for the all L-amino acid peptide were in agreement with the experimental observations. Moreover, refinement of the force field is expected to provide far-reaching conformational sampling of peptides containing D-residues to further develop force field-based conformational-searching methods. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

离子液体的量化计算及分子动力学模拟研究进展

离子液体是由有机阳离子和无机/有机阴离子构成的盐类,一般在室温或接近于室温下呈液态,因此常被称为室温离子液体(RTIL).依据不同的划分标准,离子液体有多种分类方式:根据年代的不同可将离子液体分为第一代、第二代及第三代离子液体,例如:烷基咪唑和烷基吡啶的金属卤化物盐等[1];根据阳离子的不同可将离子液体分为季鏻     

浅谈结构化学中量子力学公设的讲解

结构化学以量子力学为基础,从微观领域揭示物质结构与性质之间关系。本文剖析了作为结构化学基础的量子力学公设教学过程中遇到的难题,结合教学经验,就量子力学公设的讲解进行探讨,提出一种激发学生学习动力,提高教学效果的方法。     

Molecular Modeling of the Chain Structures of Polybenzoxazines

IntroductionThe benzoxazine polymer could be deemed asan alternative to traditional phenolics because ofthe similar main chain structures.A Novolac typeof phenolic resin has the unit,Ar—CH2 — withsome Ph—CH2 N( R)—CH2 — structure,where Phdenotes the phenolic group.In addition to theadvantages owned by ordinary phenolics,such ashigh temperature resistance,dimensional stability,good electrical properties,flame retardance,andlow smoke generation,this family ofpolybenzoxazines provides t…     

分子模拟在有机聚合物电致发光材料设计中的应用

对聚合物分子模拟理论和一种较为精确的聚合物模拟方法－地简要的介绍，分析模拟计算的结果如能量结构、前线轨道能级、能带结构、跃迁能等与电致发光相关的测定指标如导电性、氧化还原电位、光电子能谱、吸收与发光谱、光氧化性、热稳定性相关联。并讨论了分子模拟辅助电致发光新材料设计上的应用。     

A Conformational Study on Silacyclohexane. Comparison of ab initio (HF,MP2), DFT,and Molecular Mechanics Calculations. Conformational Energy Surface of Silacyclohexane

The structures and relative energies for the basic conformations of silacyclohexane 1 have been calculated using HF, RI‐MP2, RI‐DFT and MM3 methods. All methods predict the chair form to be the dominant conformation and all of them predict structures which are in good agreement with experimental data. The conformational energy surface of 1 has been calculated using MM3. It is found that there are two symmetric lowest energy pathways for the chair‐to‐chair inversion. Each of them consists of two sofa‐like transition states, two twist forms with C₁ symmetry (twist‐C₁), two boat forms with Si in a gunnel position (C₁ symmetry), and one twist form with C₂ symmetry (twist‐C₂). All methods calculate the relative energy to increase in the order chair < twist‐C₂ < twist‐C₁ < boat. At the MP2 level of theory and using TZVP and TZVPP (Si atoms) basis sets the relative energies are calculated to be 3.76, 4.80, and 5.47 kcal mol^–1 for the twist‐C₂, twist‐C₁, and boat conformations, respectively. The energy barrier from the chair to the twisted conformations of 1 is found to be 6.6 and 5.7 kcal mol^–1 from MM3 and RI‐DFT calculations, respectively. The boat form with Si at the prow (C_s symmetry) does not correspond to a local minimum nor a saddle point on the MM3 energy surface, whereas a RI‐DFT optimization under C_s symmetry constraint resulted in a local minimum. In both cases its energy is above that of the chair‐to‐twist‐C₁ transition state, however, and it is clearly not a part of the chair‐to‐chair inversion.     

高岭石-水体系中水分子结构的分子动力学模拟

以Hendricks模型为初始结构, 利用CLAYFF力场对高岭石-水体系进行无晶体学限制的分子动力学模拟. 结果表明, 层间水有三种类型: I型类似于Costanzo提出的“洞水”分子, 其HH矢量(水分子中从一个氢原子位置指向另一个氢原子位置的方向矢量)平行于(001)平面, 而C2轴稍微倾斜于(001)面法线; II型类似于“连接水”, 一个氢氧键指向临近的层间四面体氧形成氢键, 另一个氢氧键与(001)面近似平行; III型水分子在层间近似保持为竖直状, 一个氢与层间四面体氧形成氢键, 而另一个氢与对面层的羟基氧形成氢键. 高岭石羟基氢沿(001)晶面法线的浓度曲线显示一部分羟基指向变为近似平行于(001)面, 羟基氧因此能够暴露出来与层间水分子氢形成氢键. 此外, 模拟中还观察到部分II型水分子氧偏离于层间的平均位置而更靠近四面体层, 这和Costanzo的实验结果一致, 可能是X射线谱图中(002)弱衍射峰出现的原因.     

乙醇醛的分子动态结构

以量子化学计算作为起点, 为最简单的糖类分子——乙醇醛开发了两套分子力学力场参数: 基于肽类的力场和基于醛类的力场. 分子动力学模拟结果表明, 所开发的类醛力场参数能够较好地描述乙醇醛分子在水中的结构以及水分子在其周围的分布. 通过瞬时简正模式分析, 得到了3N-6个模式的瞬时振动频率和振动跃迁偶极矩等振动光谱参数的统计分布及其相关性. 结合量子化学计算和分子动力学模拟对生物分子体系的多元振动光谱参数进行预测和评估, 为从化学键水平出发模拟宽带飞秒二维红外相关光谱提供了一个新方法.     

理论化学在中国的发展

本文回顾1949年以来,我国理论化学的发展历史.全文始于“荒漠渐成绿洲”的总体概括,进入主体全面的论述.按1978年划分为前、后两期.前期为学科的成长成熟期,以创业人及少数从业者的科教活动为中心,后期为学科的繁荣期,以本学科现阶段中青年精英创造佳绩、进入学科主流领域为中心.按时间进程将前后期最有代表性的10个方面研究成...     

Fundamentals of Silicon Chemistry: Molecular States of Silicon-Containing Compounds

Silicon and its compounds have made possible the design of new materials, which, from computers to space travel, have helped to shape the technology of our 20th century. Conversely, the demands of new technology have stimulated the fast development of silicon chemistry as part of the “renaissance” of inorganic chemistry. This article uses selected examples of predominantly organosilicon compounds to discuss in simplified terms the measurement and assignment of suitable spectroscopic “molecular fingerprints” as well as the resulting benefit for the preparative chemist. The comparison of “equivalent” states of “chemically related” molecules is emphasized, based on perturbation arguments and supporting quantum-chemical models. Special attention is given to the relation between structure and energy, which allows us to understand and to predict the connectivity between and the spatial arrangement of silicon “building blocks”, the energy-dependent electron distribution over the effective nuclear potentials of a molecular framework, and, especially, the partly considerable effects of “silicon substituents” on molecular properties. Future-directed extensions and applications include polysilane band structures, Rydberg states of chromophores containing silicon centers, redox reactions and ion-pair formation of silicon-substituted π systems, and molecular dynamic phenomena in solution or on thermal fragmentation in the gas phase. The main objective is a set of clear and practical rules for interpreting measurements and planning experiments.     

Structural Dynamics of Some Large-Ring Cyclodextrins. A Molecular Dynamics Study: An Analysis of Force Field Performance

Molecular dynamics simulations were performed on four large ring cyclodextrins (LR-CDs), CD14, CD21, CD26 and CD28, both in gas phase and in water solution. Four different force fields (parm94, parm99, glycam2000a and MM3*) were used to examine their dependence with the results. The differences were not significant when parm99 and glycam2000a were used. Parm94 and MM3* results differ considerably due to inadequate bending parameters (parm94) or due to a strong stabilization by the intramolecular hydrogen bonds between hydroxyl groups. Simulations in the gas phase show that the CDs are elongated and twisted, and that cavities typical for the native CDs (CDn, n = 6, 7, 8) are not present. Simulations in water solution produced average structures that do not correspond to the conformations in the crystalline state. The LR-CDs are highly flexible and this could favor the formation of inclusion complexes through the intermediation of some other molecules that could fit to the smaller and more specific cavities they have.On sabbatical leave from the Institute of Organic Chemistry with Center of Phytochemistry, Bulgarian Academy of Science, ul. Acad. G. Bonchev, bloc 9, 1113 Sofia, Bulgaria     

分析化学中常用有机试剂的分子结构式及作用原理

Due to the fact that the molecular structure has been removed from most of the analytical chemistry textbook of recent editions, this paper collected more than 90 organic reagents from three classic text books, categorized these molecules according to the application, and analyzed the molecular structure and the mechanism of functions in applications.     

Molecular modeling of intercalation complexes of antitumor active 9-aminoacridine and a [d,e]-anellated isoquinoline derivative with base paired deoxytetranucleotides

Summary Intercalators are molecules capable of sliding between DNA base pairs without breaking up the hydrogen bonds between the DNA bases. On the basis of molecular mechanics calculations structural, models of B-DNA tetranucleotide intercalation complexes of some cytostatic active 9-aminoacridines and of a [d, e]-anellated isoquinoline derivative are presented. The drug complexes are stabilized by energetically favouredvan der Waals interactions and by selective hydrogen bonds between the side chains of the drugs and the DNA bases. Semiempirical quantum chemistry calculations revealed that the chromophoric system of the intercalators is able to form ,-charge-transfer interactions with the purine bases of the base paired deoxytetranucleotides. The theoretical findings are of interest for a more specific drug design of cytostatically active agents.

---

Molecular Modeling von Interkalationskomplexen antitumoraktiver 9-Aminoacridine sowie eines [d, e]-anellierten Isochinolinderivates mit basengepaarten Desoxytetranukleotiden

Zusammenfassung Interkalatoren sind Moleküle, die in der Lage sind, sich zwischen DNA-Basenpaare einzulagern, ohne die Wasserstoffbrücken zwischen den DNA-Basen aufzubrechen. Auf der Basis von molekülmechanischen Rechnungen werden Tetranukleotid-Interkalationskomplexe von verschiedenen zytostatisch aktiven 9-Aminoacridinen und von einem [d, e]-anellierten Isochinolinderivat präsentiert. Die Komplexe werden durch energetisch günstigevan der Waals-Interaktionen sowie durch selektive Wasserstoffbrückenbindungen zwischen den Seitenketten der Wirkstoffe und den DNA-Basen stabilisiert. Semiempirische quantenchemische Rechnungen ergaben, daß der Chromophor der Interkalatoren in der Lage ist, ,-charge-transfer Wechselwirkungen mit den Purinbasen der basengepaarten Desoxytetranukleotide auszubilden. Die theoretischen Ergebnisse sind für ein spezifischeres Wirkstoffdesign zytostatisch aktiver Verbindungen von Interesse.

     

磺酰胺类常咯啉的分子力学参数的确定

应用量子化学ＭＮＤＯ程序及分子力学ＭＭＰ２程序，确定了一些未知磺酰胺类常咯啉的分子力学参数。应用这些参数计算了几个晶体结构已知的磺酰胺类化合物，与真实结构比较，计算结果令人满意。由此说明确定的分子力学参数用于计算磺酰胺类化合物时的正确性和可行性。     

A Study of a Molecular Mechanics Field Used in Simulating Enantiorecognition

Enantiorecognitionisameasureofthediscriminationabilitybetweenenantiomericmolecules.Itisoffundamentalsignificanceinmanyareasofthechemicalandbiologicalsciences,italsoplaysanimportantroleinseparationofenantiomersbychromatography.Althoughtheintermolecularforces(hydrogenbonding,multipolarassociation,dispersionforces,chargetransfercomplexation,hydrophobicassociation.etc.)havebeenthoroughlystudiedandarewelldocumented,buthowtheseforcespreciselyworkinconcerttopromoteintermolecularactionisnotclearyet'.T…     

分子动力学模拟技术在生物分子研究中的进展

Molecular dynamics simulation (MDS) has gained increasing importance in current-day scientific research, as the supplement, guidance, or even replacement of experiments. In this review, we briefly introduce the history of the development of molecular dynamics simulation, focusing on recent progress including new-generation force fields, modern enhanced sampling schemes, and application for the investigation of complex biomolecules.