Inhibitors of prolyl endopeptidase: Characterization of the pharmacophoric pattern using conformational analysis and 3D-QSAR

Summary A structure-activity study has been carried out on several compounds known as inhibitors of the serine protease prolyl endopeptidase. Conformational analysis has been done using different molecular mechanics methods such as molecular dynamics, or a randomized conformational search method. The conformers obtained were classified using geometric and energetic criteria. A pattern recognition analysis was done in order to divide conformers according to families. The resulting dominant families, for all compounds investigated, showed very similar geometric features. Based on the lowest energy conformers obtained after randomized conformational analysis, a 3D-QSAR model was established using the CoMFA approach. The validity of this model was verified by prediciting correctly the activity of other molecules not used in the construction of this model.     

Impact of Stable Protein-protein Interaction on Protein Conformational Space

Backbone dihedral angle based clustering approach was applied to investigate the effect of protein complexation on backbone conformational space and the effect on protein dynamics. Three representative enzyme-inhibitor complexes and their comprised proteins were used as models for small-and moderate-sized globular proteins to compare available backbone conformational space before and after complexation. Microsecond time scale molecular dynamic simulations were generated to ensure sufficient statistics. The result suggests that stable protein-protein interactions lead to redistribution of protein backbone mobility and restriction of the protein backbone conformational space, especially for short time scale motions. Surprisingly, these effects are found to be uncorrelated with protein-protein interaction surface. Consistent with many experimental and computational observations, our results indicate that both induced-fit and conformational selection models play roles in stable protein complexation process, with the dominant role being different for different protein complexes.     

多肽抑制剂抑制淀粉质多肽42构象转换的分子动力学模拟和结合自由能计算

应用分子动力学模拟和结合自由能计算方法研究了多肽抑制剂KLVFF、VVIA和LPFFD抑制淀粉质多肽42 (Aβ42)构象转换的分子机理. 结果表明, 三种多肽抑制剂均能够有效抑制Aβ42的二级结构由α-螺旋向β-折叠的构象转换. 另外, 多肽抑制剂降低了Aβ42分子内的疏水相互作用, 减少了多肽分子内远距离的接触, 有效抑制了Aβ42的疏水塌缩, 从而起到稳定其初始构象的作用. 这些抑制剂与Aβ42之间的疏水和静电相互作用(包括氢键)均有利于它们抑制Aβ42的构象转换. 此外, 抑制剂中的带电氨基酸残基可以增强其和Aβ42之间的静电相互作用(包括氢键), 并降低抑制剂之间的聚集, 从而大大增强对Aβ42构象转换的抑制能力. 但脯氨酸的引入会破坏多肽的线性结构, 从而大大降低其与Aβ42 之间的作用力. 上述分子模拟的结果揭示了多肽抑制剂KLVFF、VVIA和LPFFD抑制Aβ42构象转换的分子机理, 对于进一步合理设计Aβ的高效短肽抑制剂具有非常重要的理论指导意义.     

An ab initio Molecular Orbital Study of the Conformational Properties of all-(Z)-Cyclododeca-1,4,7,10-tetraene

Summary.  Ab initio HF/6-31G^* and MP2/6-31G^*//HF/6-31G^* methods were used to calculate the structure optimization and conformational interconversion pathways for all-(Z )-cyclododeca-1,4,7,10-tetraene. This compound adopts the symmetrical crown (C _4v) conformation. Ring inversion takes place via symmetrical intermediates, such as boat-chair (BC, C _s) and twist (C _2h) conformers and requires about 22.3 kJ · mol⁻¹. The calculated strain energies for BC and twist conformers are 5.9 and 13.5 kJ · mol⁻¹. The results of semiempirical AM1 calculations for structural parameters and relative energies of the important geometries of the title compound are in good agreement with the results of ab initio methods. Received July 9, 2001. Accepted September 26, 2001     

Conformational analysis of six- and twelve-membered ring compounds by molecular dynamics

A molecular dynamics (MD)-based conformational analysis has been performed on a number of cycloalkanes in order to demonstrate the reliability and generality of MD as a tool for conformational analysis. MD simulations on cyclohexane and a series of methyl-substituted cyclohexanes were performed at temperatures between 400 and 1200 K. Depending on the simulation temperature, different types of interconversions (twist-boat–twist-boat, twist- boat–chair and chair–chair) could be observed, and the MD simulations demonstrated the expected correlation between simulation temperature and ring inversion barriers. A series of methyl-substituted 1,3- dioxanes were investigated at 1000 K, and the number of chair–chair interconversions could be quantitatively correlated to the experimentally determined ring inversion barrier. Similarly, the distribution of sampled minimum-energy conformations correlated with the energy-derived Boltzmann distribution. The macrocyclic ring system cyclododecane was subjected to an MD simulation at 1000 K and 71 different conformations could be sampled. These conformations were compared with the results of previously reported conformational analyses using stochastic search methods, and the MD method provided 19 out of the 20 most stable conformations found in the MM2 force field. Finally, the general performance of the MD method for conformational analysis is discussed.     

Molecular dynamics simulations of a nucleoside analogue of 1,4-dihydro-4-oxoquinoline-3-carboxylic acid synthesized as a potential antiviral agent: Conformational studies in vacuum and in water

Conformational analysis of nucleosides may have direct applications to the structure–activity relationship (SAR) studies and in the design of new drug candidates. Although conformational analysis may be accessed in many different ways, in this work it was performed using molecular dynamics (MD) simulation in order to study the dynamic behavior of a nucleoside derivative of 1,4-dihydro-4-oxoquinoline-3-carboxylic acid, synthesized by our group as a potential antiviral agent. The MD simulation was carried out during 10 ns in vacuum and in a box of water at two different temperatures (i.e., 300 and 600 K) using the AMBER force field. The in vacuum MD simulation results are in agreement with the crystallographic structure and with the DFT calculations of the nucleoside, revealing the anti conformer as the more stable one. The simulation in water, however, shows that both conformers may exist at 300 K, the temperature of the in vivo and in vitro assays, revealing that both the syn and anti conformers should be considered in a MD simulation study of the inhibitor–enzyme complex. Simulations are also in agreement with the NOE experiment, which shows that the anti conformer is the preferential one in DMSO-d6 solution at 298 K.     

Structure-activity relationships for apomorphine congeners. Conformational energies vs. biological activities

Summary A series of apomorphine congeners has been studied with respect to their ability to mimic the structural requirements of the dopamine pharmacophore in the potent and stereoselective dopamine receptor agonist (R)-apomorphine. Conformational energies of the mimicking structures calculated by molecular mechanics (MMP2) correlate well with the observed biological activities.     

两类HPPD酶抑制剂的比较分子场分析研究

用比较分子场分析法（CoMFA）研究了环已二酮类及3-烷基酸-2-环已烯酯类化 合物的结构与活性的关系。本研究从蛋白酶与底物动力学模拟的复合物结构出发构 建两类抑制剂化合物分子的构象，并进行了全空间搜索，CoMFA分析得到了较好的 模型（交叉验证回归系数q~2 = 0.779，模型的线性回归系数r~2 = 0.989）。该方 程不仅可以帮助推测抑制剂与受体的结合方式，还可定量地预测结构相近的类似物 活性，为设计合成新的HPPD酶抑制剂提供了理论依据。     

神经氨酸酶与三羟基甲氧基黄酮衍生物相互作用的分子模拟

以5,7,4'-三羟基-8-甲氧基黄酮(MF)为先导物, 以现有的神经氨酸酶(NA)抑制剂类药物的结构特征为参考, 设计了一系列的三羟基甲氧基黄酮衍生物, 并运用分子对接与分子动力学相结合的方法进行了筛选及作用机制分析. 分子对接结果表明, 功能团(羧基及胍基/氮-乙酰氨基)的引入并未影响衍生物在酶活性腔中的结合位置, 衍生物的结构与相互作用能之间存在一定的联系. 将羧基和胍基作为替代基团引入到MF的C7及C5位上所得的新化合物(9)在所合成的衍生物中具有最好的结合能力(-1172.52 kJ/mol), 远远优于现有先导药物4-(氮-乙酰氨基)-5-胍基-3-(3-戊氧基)安息香酸(BA)和MF与NA的结合能力(-672.12和-347.44 kJ/mol). 进一步的作用机制分析发现, 在神经氨酸酶活性腔中, 化合物9的羧基和胍基的空间取向与现有药物中这两个基团的空间取向一致, 且化合物9与先导药物MF一样, 能与活性腔内保守残基Asp151和Glu227发生较强的相互作用. 因此可认为化合物9是一种具有应用潜质的新型神经氨酸酶抑制剂. 本研究结果为实验研究和设计抗流感药物提供了可行性思路.     

Structure,dynamics, and solvation in a disordered metal–organic coordination polymer: a multiscale study

Metal–organic coordination polymers are a growing class of technologically-important materials in which transition metal ions are connected by multitopic organic chelators to form a 3-D network structure. While the structures of many highly-ordered metal–organic frameworks have been determined, far less structural information is available about the more common disordered materials. Our study combines pair distribution function analysis from total X-ray scattering, ab initio quantum mechanical calculations, and all-atom molecular dynamics to explore the structure and dynamics of a poorly-ordered branched coordination polymer. The polymer structure is highly flexible and dynamic, and is dramatically affected by its solvation state, a finding with far-reaching implications for the incorporation of coordination polymers into nanocomposite materials.     

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     

HIV-1病毒DNA与整合酶结合后的构象变化

用分子动力学(MD)模拟方法优化了HIV-1病毒DNA与整合酶(IN)二聚体(IN2)复合物模型结构, 并分析了HIV-1病毒DNA结合IN2后的构象变化. 结果表明, 按照HIV-1病毒DNA与IN2结合能力的强度, 病毒DNA可分为五个区域: 非结合区、强结合区1、弱结合区、强结合区2和反应区, 并用结合自由能计算验证了该分区的合理性. 与未结合IN2的病毒DNA相比, 复合物模型中病毒DNA除了非结合区碱基外, 其它四个区域的碱基构象变化较大. 复合物模型中病毒DNA主链较大程度地偏离标准B型DNA以及结合部位的小沟变宽都是识别IN的结构基础. 模拟结果与实验数据吻合较好, 为基于HIV-1 IN的药物分子设计提供了一定的结构信息.     

Ab Initio Study of Conformational Properties of ( Z, Z)-, ( E, Z)-, and ( E, E)-Cyclonona-1,5-dienes

Summary. Ab initio calculations at the HF/6-31G^* level of theory for geometry optimization and MP2/6-31G^*//HF/6-31G^* for a single point total energy calculation are reported for the important energy-minimum conformations and transition-state geometries of (Z,Z)-, (E,Z)-, and (E,E)-cyclonona-1,5-dienes. The C₂ symmetric chair conformation of (Z,Z)-cyclonona-1,5-diene is calculated to be the most stable form; the calculated energy barrier for ring inversion of the chair conformation via the C_s symmetric boat-chair geometry is 58.3kJmol^–1. Interconversion between chair and twist-boat-chair (C₁) conformations takes place via the twist (C₁) as intermediate. The unsymmetrical twist conformation of (E,Z)-cyclonona-1,5-diene is the most stable form. Ring inversion of this conformation takes place via the unsymmetrical chair and boat-chair geometries. The calculated strain energy for this process is 63.5kJmol^–1. The interconversion between twist and the boat-chair conformations can take place by swiveling of the trans double bond with respect to the cis double bond and requires 115.6kJmol^–1. The most stable conformation of (E,E)-cyclonona-1,5-diene is the C₂ symmetric twist-boat conformation of the crossed family, which is 5.3kJmol^–1 more stable than the C_s symmetric chair–chair geometry of the parallel family. Interconversion of the crossed and parallel families can take place by swiveling of one of the double bonds and requires 142.0kJmol^–1.     

Intramolecular charge transfer in π-conjugated oligomers: a theoretical study on the effect of temperature and oxidation state

Intramolecular charge transfer (ICT) of gaseous π-conjugated oligo-phenyleneethynylenes (OPE) induced by a homogeneous applied electric field has been theoretically investigated using a combined approach integrating molecular dynamics (MD) simulations and Perturbed Matrix Method calculations. In line with recent investigations, our results indicate the peculiar role of conformational transitions on OPE electronic properties which reflects on a strong temperature effect on ICT. Electron transfer reactions inducing chemical alteration on OPE, also taken into account in this study, revealed extremely important for explaining non-linear ICT effects and probably plays a central role in the mechanisms underlying molecular transport junctions. Our study further points out the necessity of using MD-based approach for modelling molecular electronics, even when relatively rigid molecular systems are concerned.     

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

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

基于机器学习方法的H1N1神经氨酸苷酶抑制剂的分类预测

流感是一种主要的呼吸道传染病, 在普通人群中有着较高的发病率, 而对于一些年老和高危病人还有较高的死亡率. 研究显示抑制神经氨酸苷酶(NA)可以阻断病毒RNA复制, 因此NA是有效治疗H1N1型流感病毒的重要药物靶标. 通过计算机方法进行虚拟筛选和预测NA抑制剂已经变得越来越重要. 针对酶活性位点进行基于结构的合理药物设计, 开发H1N1 病毒神经氨酸苷酶抑制剂, 已成为药物研究的热点之一. 本文通过多种机器学习方法(支持向量机(SVM)、k-最近相邻法(k-NN)和C4.5决策树(C4.5DT))对已知的神经氨酸苷酶抑制剂(NAIs)与非神经氨酸苷酶抑制剂(non-NAIs)建立分类预测模型. 其中227个结构多样性化合物(72个NAIs与155个non-NAIs)被用于测试分类预测系统, 并用递归变量消除法选择与神经氨酸苷酶抑制剂分类相关的性质描述符以提高预测精度. 本研究对独立验证集的总预测精度为75.9%-92.6%, NA 抑制剂的预测精度为64.3%-78.6%, 非H1N1抑制剂的预测精度为77.5%-97.5%. SVM法给出最好的总预测精度(92.6%). 本研究表明支持向量机等机器学习方法可以有效预测未知数据集中潜在的NA抑制剂, 并有助于发现与其相关的分子描述符.     

人类丝氨酸消旋酶的同源模建及其与多肽类抑制剂的分子对接

利用同源模建和分子动力学模拟方法构建了人类丝氨酸消旋酶(hSR)的三维结构, 并利用profile-3D和procheck方法评估了模型的可靠性. 在此基础上用分子对接程序(affinity)将多肽类抑制剂A和B分别与hSR进行对接, 获得了其复合物结构的理论模型. 通过配体与受体之间相互作用能和结构分析给出了此类抑制剂与hSR的具体结合方式, 明确了hSR与此类抑制剂结合时起重要作用的氨基酸残基, 为基于人类丝氨酸消旋酶三维结构的药物设计提供重要的参考信息.     

<Emphasis Type="Italic">Ab initio</Emphasis> Study of Structural and Conformational Properties of Five- to Nine-Membered Cyclic 1,2,3-Trienes

Summary.  The structures and relative energies of fundamental conformations of cyclopenta-1,2,3-triene, cyclohexa-1,2,3-triene, cylohepta-1,2,3-triene, cycloocta-1,2,3-triene, and cyclonona-1,2,3-triene were calculated by the HF/6-31G^* and MP2/6-31G^*//HF/6-31G^* methods. Only a C _2v symmetric planar conformation is available to cyclopenta-1,2,3-triene and cyclohexa-1,2,3-triene. The calculated energy barrier for ring inversion of the C _S symmetric puckerd conformation of cyclohepta-1,2,3-triene via the planar geometry is 62.2 kJ·mol⁻¹. The C ₂ symmetric twist conformation of cycloocta-1,2,3-triene was calculated to be the most stable one. Conformational racemization of the twist form takes place via the C _S symmetric half-chair geometry, which is by 60.8 kJ·mol⁻¹ less stable than the twist conformer. The C _S symmetric chair and unsymmetrical twist-boat conformations of cyclonona-1,2,3-triene were calculated to have similar energies; their interconversion takes place via an unsymmetrical low-energy (18.4 kJ·mol⁻¹) transition state. The twist (C ₂) and boat (C _S) geometries of cyclonona-1,2,3-triene are higher in energy by 13.2 and 33.9 kJ·mol⁻¹, respectively. Ring inversion in chair and twist-boat conformations takes place via a twist form as intermediate and requires 33.6 kJ·mol⁻¹. Corresponding author. E-mail: isayavar@yahoo.com Received March 25, 2002; accepted April 4, 2002     

Y220C突变体影响p53C蛋白质构象转换的分子动力学模拟

p53是迄今发现突变频率最高的一种肿瘤抑制蛋白质,突变会导致p53抑癌功能丧失并诱导癌症的发生。绝大多数的突变发生在p53的核心DNA结合区域（p53C）,其中Y220C是研究较多的一种突变体。虽然已有研究表明该突变能够降低p53C的结构稳定性,但其影响p53C构象转换的分子机制尚不清晰。本文利用分子动力学（MD）模拟方法研究了p53C突变体Y220C（p53C-Y220C）的结构变化,发现Y220C突变主要影响Y220C cluster区域（包括残基138-164和215-238）,且Y220C突变减少了Y220C cluster的β-折叠含量。进一步分析发现,Y220C突变不仅直接破坏突变氨基酸与周围氨基酸Leu145和Thr155之间的氢键,而且降低了Y220C cluster区域的折叠片S3和S8之间的氢键数量,使Y220C突变所形成的亲水性空腔变大,加速了水分子进入该蛋白质内部,并最终导致了p53C-Y220C变性。MD模拟结果揭示了Y220C突变影响p53C结构转换的分子机制,该研究对p53C-Y220C突变体高效稳定剂的筛选和设计具有重要意义。     

Conformational analysis. Part 16 Conformational free energies in substituted piperidines and piperidinium salts

Summary The conformational free energies (-G^o) of a number of 4-substituted piperidines and piperidinium salts have been determined by the J-value method. For the 4-substituted piperidines (R=Me, Phenyl, CO₂Et, Br, OH, F) the relative conformer energies are almost identical to those of the analogous cyclohexanes.The methyl and phenyl compounds showed no change in the couplings on protonation, implying no change in the conformer energies. In constrast, in the remaining compounds with polar 4-substituents an almost constant stabilisation of the axial conformer of ca. 0.7–0.8 kcal mol^-1 was observed on protonation. In three cases (R=F, OH and Br) the conformational preferences is reversed on protonation and the axial form is favoured.The conformer energies of both the free bases and the piperidinium salts can be quantitatively predicted by molecular mechanics calculations using the COSMIC force-field, in which the electrostatic interactions are calculated by a simple Coulombic model with the partial atomic charges in the molecules given by the CHARGE2 routine, and an effective dielectric constant of five. The precise agreement obtained demonstrates conclusively that the electrostatic interactions between the substituents and the protonated nitrogen are the cause of the conformational changes on protonation, and that these can be modelled successfully using existing force-fields.For Part 15, see Ref. 1.