氢键准则在DMSO水溶液中应用的评价

利用分子动力学模拟方法, 分别采用几何准则和能量准则分析了不同浓度下的二甲基亚砜(DMSO)水溶液的氢键统计和动力学等特性. 结果显示, 两种氢键准则可以很好地反映出溶液的氢键性质随浓度的变化趋势. 通过分析比较发现, 由于几何准则不能有效地排除具有弱对势能的分子对, 因此其统计的氢键数量要大于能量准则的结果.此外, 能量准则对于分子间相对取向的区分存在不足, 进而引起氢键寿命的计算结果偏大.因此,为使氢键分析更加准确, 本文建议使用几何-能量混合型氢键准则.     

智能弹药的实时仿真可视化系统设计

针对仿真可视化系统开发时间长、开发费用高的问题,提出了一种利用Unity3D引擎进行实时仿真可视化系统开发的模块式设计方法。首先,对武器仿真可视化系统进行了分析,提出了智能弹药实时仿真可视化系统的五大模块;然后,通过在3ds Max下建立智能弹药模型和Unity3D下的二次开发,构建了智能弹药仿真可视化系统的仿真环境,以C#为编程工具,设计了仿真驱动、串口数据通讯和仿真输出模块;最终,实现了智能弹药实时仿真可视化系统。实验证明,该系统可以有效地在监视器上实时显示出智能弹药的飞行姿态、飞行参数和飞行轨迹。     

Net charge changes in the calculation of relative ligand‐binding free energies via classical atomistic molecular dynamics simulation

The calculation of binding free energies of charged species to a target molecule is a frequently encountered problem in molecular dynamics studies of (bio‐)chemical thermodynamics. Many important endogenous receptor‐binding molecules, enzyme substrates, or drug molecules have a nonzero net charge. Absolute binding free energies, as well as binding free energies relative to another molecule with a different net charge will be affected by artifacts due to the used effective electrostatic interaction function and associated parameters (e.g., size of the computational box). In the present study, charging contributions to binding free energies of small oligoatomic ions to a series of model host cavities functionalized with different chemical groups are calculated with classical atomistic molecular dynamics simulation. Electrostatic interactions are treated using a lattice‐summation scheme or a cutoff‐truncation scheme with Barker–Watts reaction‐field correction, and the simulations are conducted in boxes of different edge lengths. It is illustrated that the charging free energies of the guest molecules in water and in the host strongly depend on the applied methodology and that neglect of correction terms for the artifacts introduced by the finite size of the simulated system and the use of an effective electrostatic interaction function considerably impairs the thermodynamic interpretation of guest‐host interactions. Application of correction terms for the various artifacts yields consistent results for the charging contribution to binding free energies and is thus a prerequisite for the valid interpretation or prediction of experimental data via molecular dynamics simulation. Analysis and correction of electrostatic artifacts according to the scheme proposed in the present study should therefore be considered an integral part of careful free‐energy calculation studies if changes in the net charge are involved. © 2013 The Authors Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Design of a reaction field using a linear‐combination‐based isotropic periodic sum method

In our previous study (Takahashi et al., J. Chem. Theory Comput. 2012, 8, 4503), we developed the linear‐combination‐based isotropic periodic sum (LIPS) method. The LIPS method is based on the extended isotropic periodic sum theory that produces a ubiquitous interaction potential function to estimate homogeneous and heterogeneous systems. The LIPS theory also provides the procedure to design a periodic reaction field. To demonstrate this, in the present work, a novel reaction field of the LIPS method was developed. The novel reaction field was labeled LIPS‐SW, because it provides an interaction potential function with a shape that resembles that of the switch function method. To evaluate the ability of the LIPS‐SW method to describe in homogeneous and heterogeneous systems, we carried out molecular dynamics (MD) simulations of bulk water and water–vapor interfacial systems using the LIPS‐SW method. The results of these simulations show that the LIPS‐SW method gives higher accuracy than the conventional interaction potential function of the LIPS method. The accuracy of simulating water–vapor interfacial systems was greatly improved, while that of bulk water systems was maintained using the LIPS‐SW method. We conclude that the LIPS‐SW method shows great potential for high‐accuracy, high‐performance computing to allow large scale MD simulations. © 2014 Wiley Periodicals, Inc.     

Acceleration of self‐consistent field convergence in ab initio molecular dynamics simulation with multiconfigurational wave function

The Lagrange interpolation of molecular orbital (LIMO) method, which reduces the number of self‐consistent field iterations in ab initio molecular dynamics simulations with the Hartree–Fock method and the Kohn–Sham density functional theories, is extended to the theory of multiconfigurational wave functions. We examine two types of treatments for the active orbitals that are partially occupied. The first treatment, as denoted by LIMO(C), is a simple application of the conventional LIMO method to the union of the inactive core and the active orbitals. The second, as denoted by LIMO(S), separately treats the inactive core and the active orbitals. Numerical tests to compare the two treatments clarify that LIMO(S) is superior to LIMO(C). Further applications of LIMO(S) to various systems demonstrate its effectiveness and robustness. © 2014 Wiley Periodicals, Inc.     

CYP2C9酶与Warfarin结合模型的立体选择性理论研究

对CYP2C9酶与S-Warfarin复合物的晶体结构进行分子对接、分子动力学模拟、通道分析及结合自由能计算,发现原晶体结构中的结合模式为"亚稳态",提出了CYP2C9与S-Warfarin结合的可催化模式;比较了CYP2C9与S-和R-Warfarin结合的异同,确定了在结合过程中起重要作用的锚定氨基酸残基,尤其是位于活性位点区域的苯丙氨酸簇.在结合过程中这些残基通过芳香环的移动对稳定底物的结合模式起到至关重要的作用,阐明了该酶呈现相关底物选择性的原因.对于CYP2C9与底物对接模式及立体选择性的研究有助于在分子层面上理解特异性底物与酶的结合特点,为潜在的药物设计提供了合理可信的理论依据.     

分子动力学模拟研究点突变(Met108→Leu108)对树胶醛糖结合蛋白与配体作用的影响

为了研究点突变(Met108→Leu108)对树胶醛糖结合蛋白(ABP)与配体结合能力的影响,对ABP、ABP结合树胶醛糖复合物及ABP结合半乳糖复合物以及它们各自的突变体分别进行60 ns的分子动力学模拟.模拟结果表明,108号残基突变前后,电子等排体的两个氨基酸残基,使蛋白与配体间的范德华相互作用发生明显变化,同时导致蛋白的内部运动也发生变化,进而影响蛋白与配体的相互作用.进一步分析表明,突变前后的蛋白构象变化都趋向于两个结构域张开,而与配体的结合可减缓张开程度.     

Finite Element Based Micro-mechanical Modeling of Interphase in Filler Reinforced Elastomers

Characteristic properties of elastomers can be tailored by embedding them with filler particles. Along with enhancing the overall properties of the composite, filler particles also induce some inelastic effects. In this paper, a finite element computational model is used to study the effect of microstructure morphology in filled elastomers, on its macroscopic large deformation behavior. A multiphase material model that accounts for the hypothesis of shift in glass transition temperature in the vicinity of the filler particle is developed to simulate the interphase between the fillers and the matrix. It also accounts for the breakdown and re-aggregation of filler networks under cyclic loading. Examples at the microstructural level, demonstrating the dynamics of the interphase using the developed multiphase model have been successfully simulated. The obtained results are in good qualitative agreement with the Mullins effect. Therefore, computational experiments using this methodology enable the prediction of the experimentally observed softening behavior in filled elastomers based on its microstructure evolution.     

Development of Complex Mathematical Model of Light Naphtha Isomerization and Rectification Processes

The technique of developing a mathematical model of catalytic isomerization of light naphtha is stated Using experimental data from an industrial isomerization unit shows adequacy of the mathematical model to the real process. The paper presents a method for optimizing the operation of the plant together with catalytic isomerization unit and separation columns. Selection of optimal modes of separation columns allows achieving the desired flow separation between units, as well as extension of the life of the catalyst SI-2.