Influence of Intermolecular Vibrations on the Electronic Coupling in Organic Semiconductors: The Case of Anthracene and Perfluoropentacene

We have performed classical molecular dynamics simulations and quantum‐chemical calculations on molecular crystals of anthracene and perfluoropentacene. Our goal is to characterize the amplitudes of the room‐temperature molecular displacements and the corresponding thermal fluctuations in electronic transfer integrals, which constitute a key parameter for charge transport in organic semiconductors. Our calculations show that the thermal fluctuations lead to Gaussian‐like distributions of the transfer integrals centered around the values obtained for the equilibrium crystal geometry. The calculated distributions have been plugged into Monte‐Carlo simulations of hopping transport, which show that lattice vibrations impact charge transport properties to various degrees depending on the actual crystal structure.     

A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-containing active sites. Furthermore, an energy decomposition analysis was performed to determine the contributions of individual residues to the enzyme-inhibitor interactions on averaged structures obtained from rather extensive conformational sampling. Analysis of the results reveals first that there is a correlation between protein-ligand interaction energy and experimental strand transfer into human chromosomes and secondly that the Asn-155, Lys-156 and Lys-159 residues and the Mg(2+) ion are crucial to anti-HIV IN activity. These results may explain the available experimental data.     

Influence of Aromatic Residues on the Material Characteristics of Aβ Amyloid Protofibrils at the Atomic Scale

Amyloid fibrils, which cause a number of degenerative diseases, are insoluble under physiological conditions and are supported by native contacts. Recently, the effects of the aromatic residues on the Aβ amyloid protofibril were investigated in a ThT fluorescence study. However, the relationship between the material characteristics of the Aβ protofibril and its aromatic residues has not yet been investigated on the atomic scale. Here, we successfully constructed wild‐type (WT) and mutated types of Aβ protofibrils by using molecular dynamics simulations. Through principle component analysis, we established the structural stability and vibrational characteristics of F20L Aβ protofibrils and compared them with WT and other mutated models such as F19L and F19LF20L. In addition, structural stability was assessed by calculating the elastic modulus, which showed that the F20L model has higher values than the other models studied. From our results, it is shown that aromatic residues influence the structural and material characteristics of Aβ protofibrils.     

Theoretical Consideration of the External Donor of Heterogeneous Ziegler–Natta Catalysts Using Molecular Mechanics,Molecular Dynamics,and QSAR Analysis

In the previous paper we calculated the most stable structure of the catalyst and cocatalyst of metallocene catalyst systems using molecular dynamics and molecular mechanics. In this paper we postulate the active site of the heterogeneous Ziegler–Natta catalyst, and analyze it by the same methods to clarify the main factors that control the activity and molecular weight distribution of the heterogeneous catalyst systems. The roles of the external donor were also considered, and we found that the interaction energy between an external donor and the active site of the catalyst, as well as the structural factors of the external donor itself, are strongly related to both activity and molecular weight distribution. These results reveal that molecular dynamics and molecular mechanics calculations are also effective methods to screen the heterogeneous catalyst systems.     

Is Regioselectivity in the Enzyme‐Catalyzed Hydroperoxidation of Arachidonic Acid Necessarily Determined by Hydrogen Abstraction? The Case of Rabbit Leu597Ala/Ile663Ala ALOX15 Mutant

Molecular dynamics simulations and quantum mechanics/molecular mechanics calculations were performed on the in silico Leu597Ala/Ile663Ala double mutant of rabbit ALOX15 (12/15 lipoxygenase). The computational results suggested that subtle steric hindrance by the conserved Leu597 and C‐terminal Ile663 residues disturbed H10 abstractions in wildtype ALOX15 (which abstracts H13), but if these two bulky residues were mutated to smaller ones, H10 abstraction was no longer impeded and the regioselectivity of the initial H‐abstraction step was changed. However, site‐directed mutagenesis with HPLC analysis of the products of the whole oxidation process showed that the regioselectivity of the hydroperoxidation was not altered. This disagreement may be explained by the conformational reorganization of the system needed to rotate the ?OO^. group from an antarafacial to a suprafacial arrangement prior to back‐hydrogen transfer. After H10 abstraction and O₂ insertion, the evolution of the peroxy radical at C12 was sterically impeded, whereas peroxyl group rotation at C15 (after H13 abstraction) could easily evolve to a suprafacial arrangement, which thus led to the final product. For this reason, the global regiospecificity was not affected in the mutant. These findings exemplify that the regioselectivity of initial hydrogen abstraction and the regioselectivity of the final product do not necessarily coincide (in fact, they can be opposite) for the hydroperoxidation of arachidonic acid catalyzed by a lipoxygenase.     

The Catalytic Mechanism of the Marine‐Derived Macrocyclase PatGmac

Cyclic peptides are a class of compounds with high therapeutic potential, possessing bioactivities including antitumor and antiviral (including anti‐HIV). Despite their desirability, efficient design and production of these compounds has not been achieved to date. The catalytic mechanism of patellamide macrocyclization by the PatG macrocyclase domain has been computationally investigated by using quantum mechanics/molecular mechanics methodology, specifically ONIOM(M06/6‐311++G(2d,2p):ff94//B3LYP/6‐31G(d):ff94). The mechanism proposed herein begins with a proton transfer from Ser783 to His 618 and from the latter to Asp548. Nucleophilic attack of Ser783 on the substrate leads to the formation of an acyl–enzyme covalent complex. The leaving group Ala‐Tyr‐Asp‐Gly (AYDG) of the substrate is protonated by the substrate's N terminus, leading to the breakage of the P1?P1′ bond. Finally, the substrate's N terminus attacks the P1 residue, decomposing the acyl–enzyme complex forming the macrocycle. The formation and decomposition of the acyl–enzyme complex have the highest activation free energies (21.1 kcal mol^?1 and 19.8 kcal mol^?1 respectively), typical of serine proteases. Understanding the mechanism behind the macrocyclization of patellamides will be important to the application of the enzymes in the pharmaceutical and biotechnological industries.     

Molecular Dynamics Simulations of Small Peptides: Can One Derive Conformational Preferences from ROESY Spectra?

Folding properties of beta-peptides were investigated by means of NMR experiments and MD simulations of beta-dipeptides, which serve as small test systems to study the influence of stereocenters and side chains on hydrogen-bond and consequently on secondary-structure formation. Two stereoisomers, SR and SS, of a Val-Phe dipeptide, and of the corresponding Ala-Ala dipeptide, and a Gly-Gly dipeptide were simulated in methanol for 40 ns. In agreement with experiment, the isomers of the Val-Phe dipeptide adopt quite different conformers at 298 K, the differences being reduced at 340 K. Interestingly, the SR isomer shows enhanced hydrogen bonding at the higher temperature. The adopted conformations are primarily determined by the R or S side chain substitution, and less by the type of side chain. Back-calculation of (1)H ROESY spectra and (3)J coupling constants from the MD simulations and comparison with the experimental data for the Val-Phe dipeptides shows good agreement between simulation and experiment, and reveals possible problems and pitfalls, when deriving structural properties of a small and extremely flexible molecule from NMR data only. Inclusion of all aspects of internal dynamics is essential to the correct prediction of the NMR spectra of these small molecules. Cross comparison of calculated with experimental spectra for both isomers shows that only a few out of many ROESY peaks reflect the sizeable conformational differences between the isomers at 298 K.     

新的钌环丙烷化催化体系的探索(Ⅱ) 氮膦配体立体化学结构控制顺反比选择性的初步研究

本文研究了一个新的钌催化的苯乙烯环丙烷化反应催化体系。当该催化体系用三氟甲磺酸银预处理后其催化的环丙烷化反应产物的顺反比接近1：1,且顺式产物略微占优。这在环丙烷化反应中是不多见的。我们用分子力学／动力学对该催化体系进行了模拟。结果表明键能并非是决定产物顺反比的主要因素,但是对于非键能来说,由于中间体M1与苯乙烯作用生成顺式产物时的范德华作用非常强,它弥补了由静电能所引起的能量上的不利因素。所以总能量对生成顺式产物较为有利。但是由于生成顺式产物与反式产物的总能量差别不大（2．5kcal.mol^-1),所以环丙烷化产物的顺反比接近1：1。这是氮膦配体立体化学结构控制产物选择性的结果。     

Comparison of properties of Aib-rich peptides in crystal and solution: a molecular dynamics study.

In order to study the differences of the structural properties of Aib-rich peptides in solution and in the crystalline state, molecular dynamics (MD) simulations of the Aib-containing peptide II (pBrBz-(Aib)5-Leu-(Aib)2-OMe) were performed in the crystalline state, starting from two different conformers obtained experimentally by X-ray diffraction. The structural properties as derived from X-ray crystallography (e.g., torsional angles and hydrogen bonds) are well-reproduced in both constant-volume and constant-pressure simulations, although the force-field parameters used result in a too-high density of the crystals. Through comparison with the results from previous MD and nuclear magnetic resonance (NMR) studies of the very similar peptide I (Z-(Aib)s-Leu-(Aib)2-OMe) in dimethylsulfoxide (DMSO) solution, it is found that, in the crystal simulation, the conformational distribution of peptide II is much narrower than that in the solution simulation of peptide. I. This leads to a significant difference in 3 [symbol: see text] (HN, HC alpha) coupling constant values, in agreement with experimental data, whereas the NOE intensities or proton-proton distance bounds appear insensitive to the difference in conformational distribution. For small peptides the differences between their conformational distribution in the crystalline form and in solution may be much larger than for proteins, a fact which should be kept in mind when interpreting molecular properties in the solution state by using X-ray crystallographic data.     

分子力学和分子动力学方法研究不同变质程度烟煤的分子结构

采用分子力学和分子动力学模拟方法 ,研究了不同变质程度烟煤的三维分子构型和能量参数。结果表明 :随变质程度的增加 ,煤分子内平行的芳香片层结构增大。不同变质程度烟煤分子在团聚前后的成键相互作用能仅略有变化 ,其中扭转能Et 的变化相对较为显著 ,并且扭转能随煤阶的增加呈逐渐减小的趋势。非成键作用能 ,特别是超过三个原子的范德华作用能 ,是煤中的重要相互作用能 ,是模型分子团聚的重要驱动力 ,对于烟煤分子聚集状态的形成起着决定性的作用。同时随煤阶的升高 ,超过三个原子的范德华作用能逐渐增加。模型分子的总势能随煤阶的变化呈两头高、中间低的趋势 ,与煤的一些宏观物理性质有一定的相一致性。