人类细胞色素P450 2E1酶在不同乙醇浓度下构象变化的分子动力学模拟

细胞色素P450(CYP) 2E1家族酶是一种具有双重功能的单加氧酶, 能够参与市场上6%药物的代谢而具有重要的作用. 这类酶与酒精的消耗、 糖尿病、 肥胖症以及厌食症等密切相关, 引起了广泛的研究兴趣. 目前尚未见从原子水平上对这种酶在不同乙醇浓度下构象行为的研究. 基于此, 本文研究了花生四烯酸(AA)与CYP2E1复合物结构在不同乙醇浓度下构象与能量变化的特点. 对于在不同乙醇浓度下AA与CYP2E1的复合物结构, 采用分子动力学模拟结合自由能计算的方法进行研究. 分子动力学模拟结果表明, His109和Lys243氨基酸残基对AA与CYP2E1的结合起到了至关重要的作用. 当体系的乙醇浓度较高时, AA的结合能力有所下降, 这种结合能力的下降是由于AA与CYP2E1之间氢键相互作用力的减弱所致. 本研究对于AA与CYP2E1复合物结构在不同乙醇浓度下, AA分子与CYP2E1分子结合能力下降以及CYP2E1的构象变化给出了详细的解释. 本研究工作得到的结论对于实验和理论研究均有重要意义, 可为后续细胞色素P450酶类催化活性的研究提供理论支持.     

己烯雌酚与CYP2C9的计算模拟及多光谱法研究

为了研究己烯雌酚与CYP2C9的结合作用机制,结合分子模拟,荧光光谱等光谱实验多角度分析了两者的结合情况。首先从分子对接得到己烯雌酚与CYP2C9相互作用的最佳结合构象,然后通过动力学模拟研究了己烯雌酚与CYP2C9的复合物的稳定性以及构象变化,最后用多光谱法进行实验印证。对接结果表明己烯雌酚与CYP2C9反应可以自发进行,动力学模拟结果表明两者具有较强的结合能力。同时,荧光光谱实验得出两者的结合机制为静态猝灭且形成一个结合位点,热力学参数证明两者结合作用为疏水作用力;紫外光谱实验进一步证明两者结合后的CYP2C9的构象和周围环境发生改变。通过拟合分析酰胺I带,红外定量分析结果表明,与己烯雌酚作用后的CYP2C9蛋白质的二级结构发生改变。综上结果表明CYP2C9与己烯雌酚可以发生相互作用,理论与实验进行相互印证,为进一步研究CYP2C9的催化代谢机制提供参考。     

血管紧张素转换酶与抑制肽结合模式的分子动力学研究

应用分子模拟方法研究了血管紧张素转换酶(Angiotensin-converting enzyme,ACE)C端结构域(C-domain)与两种抑制肽(RIGLF/AHEPVK)的结合机制,预测了两个体系的结合模式,提出在C-domain-RIGLF中His353,Asp377,Asp453,Phe457,His513,Tyr523和Phe527为RIGLF主要结合残基,而在C-domainAHEPVK中Gln281,His353,Ser355,Glu384,Lys511,His513和Tyr523等残基起关键作用.应用结合自由能计算比较了两个体系的结合能力,结果表明,RIGLF和AHEPVK均与C-domain活性位点残基存在较强作用,且AHEPVK对C-domain的结合能力较强,与实验结果一致.     

细胞色素P450酶催化的烷烃选择性羟化

碳氢键选择氧化是合成化学领域的重要课题,其中烷烃选择性羟化反应更是面临着化学选择性、区域选择性和立体选择性等多重挑战.细胞色素P450酶广泛分布于动植物和微生物体内,是公认的多功能生物氧化催化剂. P450酶对惰性C—H键的选择性氧化具有独特优势,在催化烷烃选择性羟化反应方面拥有巨大潜力.本综述简述了P450单加氧酶及其催化烷烃选择性羟化的反应机理,梳理了来自CYP153家族、CYP52家族和其他家族的天然P450酶催化各类烷烃底物的氧化反应和选择性,讨论了理性设计和定向进化策略在开发烷烃羟化P450突变酶过程中的经典案例,介绍了底物工程、诱饵分子、双功能小分子协同催化等几种化学活化P450酶的策略及其在烷烃羟化上的应用,探讨了P450酶在烷烃选择性羟化方面所面临的挑战和解决途径,并展望了其应用前景.     

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的药物分子设计提供了一定的结构信息.     

UPLC-MS/MS结合多探针底物方法研究刺五加叶中黄酮苷类成分对CYP450活性的影响

利用超高效液相色谱-串联质谱联用(UPLC-MS/MS)的多反应监测(MRM)技术结合多探针底物方法, 研究了刺五加叶中的主要黄酮苷类化合物槲皮苷、金丝桃苷及芦丁对肝细胞色素P450酶(CYP450)亚型CYP1A2, CYP2C, CYP2E1, CYP2D和CYP3A活性的影响. 结果表明, 3种化合物对各CYP亚型酶均有抑制作用, 其中金丝桃苷和槲皮苷对CYP1A2催化的非那西丁的O-脱乙基反应抑制的IC₅₀值分别为46.53和49.75 μmol/L, 金丝桃苷和芦丁对CYP2E1催化的氯唑沙宗的6-羟基化反应抑制的IC₅₀值分别为99.87和86.36 μmol/L. 机理性抑制实验结果表明, 3种化合物对2种亚型酶的抑制作用是随着预孵时间延长而增强的机理性抑制.     

用分子对接方法预测HIV-1整合酶与金精三羧酸抑制剂的相互作用

用分子对接方法(Docking)研究了HIV-1整合酶与其抑制剂金精三羧酸的结合过程.为弄清金属离子在结合中所起的作用,选择含有一个Mg+2或不含Mg+2的两种不同的整合酶受体分别与金精三羧酸对接.结果表明, Mg+2对稳定配体与受体的结合起了重要作用. 金精三羧酸配体与含有一个金属Mg+2的整合酶受体对接,最优结合自由能为-45.19 kJ/mol. 当Mg+2失去后,整合酶的活性中心构象将发生变化,使金精三羧酸抑制剂与整合酶的结合自由能(-24.35 kJ/mol)明显增加. 预测了未知的HIV-1整合酶与其抑制剂金精三羧酸的复合物结构, 并可对基于结构的抗HIV-1整合酶的药物设计提供重要信息.     

甲状腺结合前清蛋白的理论研究

甲状腺结合前清蛋白TTR是一种具有重要生理功能的蛋白质,它是约30种与淀粉样疾病相关的非同源蛋白中的一种。与TTR相关的淀粉样疾病主要有：家族淀粉化心肌疾病,家族淀粉化神经系统疾病,老年系统性淀粉样病变,以及中枢神经系统选择性淀粉化疾病等。这些疾病是由TTR四聚体解聚过程中错误折叠形成cross-β-sheet结构形态的淀粉样纤维所导致。本文介绍了TTR的生理功能及结构特征,并综述了到目前为止用分子动力学模拟、分子对接和定量构效关系等方法在研究TTR淀粉样机理及TTR和小分子相互作用过程中的计算化学研究成果,为基于TTR结构的TTR淀粉样抑制剂药物分子的设计和筛选提供有力参考。     

Catalytic Mechanism of Cytochrome P450 2D6 for 4-Hydroxylation of Aripiprazole： A QM/MM Study

Drug metabolism is an important issue in drug discovery. Understanding how a drug is metabolized in the body will provide helpful information for lead optimization. Cytochrome P450 2D6 （CYP2D6） is a key enzyme for drug metabolism and responsible for the metabolism of about one third marketed drugs. Aripiprazole is an atypical an- tipsychotic and metabolized by CYP2D6 to its hydroxylated form. In this study, a series of computational methods were performed to understand how CYP2D6 accomplishes the 4-hydroxylation of aripiprazole. Molecular docking and molecular dynamics simulations were first performed to prepare the initial conformations for QM/MM calcula- tions. The results revealed two possible conformations for the drug-CYP2D6 complex. The ONIOM method for QM/MM calculations was then carried out to show detailed reaction pathways for the CYP2D6-catalyzed aripipra- zole hydroxylation reaction, which demonstrated that the dominant reactive channel was electrophilic and involved an initial attack on the n-system of the dichlorophenyl group of aripiprazole to produce cation δ-complex. Further- more, the product complex for each conformation was thermodynamically stable, which is in good agreement with previous reports.     

抑制剂BD2对PKA与PKC βII的抑制选择性研究

蛋白激酶A (PKA)和蛋白激酶C (PKC)的过度表达导致细胞生长分化异常, 是治疗肿瘤的潜在靶点. 抑制剂BD2对PKA和PKC抑制作用存在高选择性. 为了探讨BD2高选择性机制, 本工作以PKA与BD2复合物的晶体结构为模板, 通过同源模建结合分子对接的方法构建PKC βII与BD2复合物的结构, 并对PKA-BD2复合物和PKC-BD2复合物进行了2.5 ns的分子动力学模拟, 运用MM-GBSA方法计算了结合自由能, 通过能量分解的方法考察PKA和PKC的主要残基与BD2之间的相互作用和识别机制. 结合能分析结果很好地描述了BD2对PKA抑制活性比其对PKC抑制活性高这一实验现象. 氢键分析和能量分解结果共同说明了BD2的B环和酰胺链部分与PKA和PKC中相应位点的残基之间的相互作用存在差异, 这是BD2存在选择性的内在因素. BD2高选择性作用机制的阐明为进一步基于结构的balanol类抑制剂的结构设计和优化提供了合理的指导.     

Molecular characterization,modeling and docking of CYP107CB2 from Bacillus lehensis G1, an alkaliphile

Cytochrome P450s are a superfamily of heme monooxygenases which catalyze a wide range of biochemical reactions. The reactions involve the introduction of an oxygen atom into an inactivated carbon of a compound which is essential to produce an intermediate of a hydroxylated product. The diversity of chemical reactions catalyzed by cytochrome P450s has led to their increased demand in numerous industrial and biotechnology applications. A recent study showed that a gene sequence encoding a CYP was found in the genome of Bacillus lehensis G1, and this gene shared structural similarity with the bacterial vitamin D hydroxylase (Vdh) from Pseudonocardia autotrophica. The objectives of present study was to mine, for a novel CYP from a new isolate B. lehensis G1 alkaliphile and determine the biological properties and functionalities of CYP in this bacterium. Our study employed the usage of computational methods to search for the novel CYP from CYP structural databases to identify the conserved pattern, functional domain and sequence properties of the uncharacterized CYP from B. lehensis G1. A computational homology model of the protein’s structure was generated and a docking analysis was performed to provide useful structural knowledge on the enzyme’s possible substrate and their interaction. Sequence analysis indicated that the newly identified CYP, termed CYP107CB2, contained the fingerprint heme binding sequence motif FxxGxxxCxG at position 336-345 as well as other highly conserved motifs characteristic of cytochrome P450 proteins. Using docking studies, we identified Ser-79, Leu-81, Val-231, Val-279, Val-383, Ala-232, Thr-236 and Thr-283 as important active site residues capable of stabilizing interactions with several potential substrates, including vitamin D₃, 25-hydroxyvitamin D₃ and 1α-hydroxyvitamin D₃, in which all substrates docked proximally to the enzyme’s heme center. Biochemical analysis indicated that CYP107CB2 is a biologically active protein to produce 1α,25-dihydroxyvitamin D₃ from 1α-hydroxyvitamin D₃. Based on these results, we conclude that the novel CYP107CB2 identified from B. lehensis G1 is a putative vitamin D hydroxylase which is possibly capable of catalyzing the bioconversion of parental vitamin D₃ to calcitriol, or related metabolic products.     

人类2-氨基3-羧基粘康酸6-半醛脱羧酶(ACMSD)与底物及抑制剂作用模型的理论研究

利用同源模建和分子动力学模拟方法构建了人类2-氨基3-羧基粘康酸6-半醛脱羧酶(hACMSD)的三维结构, 并利用Profile-3D和Procheck等方法评估了模型的可靠性. 在此基础上, 用分子对接程序(Affinity), 将其底物2-氨基3-羧基粘康酸6-半醛(ACMS)和抑制剂喹啉酸(QA)分别与hACMSD进行对接, 获得了复合物结构的理论模型. 通过配体与受体之间相互作用能和结构分析给出了底物和抑制剂的具体结合方式, 明确了hACMSD与底物和抑制剂结合时起重要作用的氨基酸残基.     

Rapid testing of clopidogrel resistance by genotyping of CYP2C19 and CYP2C9 polymorphisms using denaturing on-chip capillary electrophoresis

Antiplatelet therapy is a cornerstone of cardiovascular treatment in patients with coronary artery disease and after myocardial infarction. Clopidogrel has become a popular antiplatelet agent due to its fast action and low frequency of adverse effects. Kinetics of clopidogrel metabolism is driven by enzymatic activity of the Cytochrome P450 system. Genotyping of CYP2C19 and CYP2C9 polymorphisms allows to identify slow metabolizers showing resistance to clopidogrel therapy. Today, a number of PCR-based techniques for single nucleotide polymorphism genotyping directed at clopidogrel resistance polymorphisms are in use. Here, we describe a new alternative genotyping approach combining the separation power of denaturing capillary electrophoresis with the analysis speed and ease of use of Bioanalyzer chipCE platform. Using an upgraded heater control, we present an optimization for allele separation of CYP2C19 I331V, CYP2C9 R144C, and CYP2C9 I359L polymorphisms employing run temperatures of up to 55°C. We demonstrate rapid and accessible approach to reproducible clopidogrel resistance with feasibility and low cost.     

Determinants of the Inhibition of DprE1 and CYP2C9 by Antitubercular Thiophenes

Mycobacterium tuberculosis (Mtb) DprE1, an essential isomerase for the biosynthesis of the mycobacterial cell wall, is a validated target for tuberculosis (TB) drug development. Here we report the X‐ray crystal structures of DprE1 and the DprE1 resistant mutant (Y314C) in complexes with TCA1 derivatives to elucidate the molecular basis of their inhibitory activities and an unconventional resistance mechanism, which enabled us to optimize the potency of the analogs. The selected lead compound showed excellent in vitro and in vivo activities, and low risk of toxicity profile except for the inhibition of CYP2C9. A crystal structure of CYP2C9 in complex with a TCA1 analog revealed the similar interaction patterns to the DprE1–TCA1 complex. Guided by the structures, an optimized molecule was generated with differential inhibitory activities against DprE1 and CYP2C9, which provides insights for development of a clinical candidate to treat TB.     

细胞色素P450 2f1(CYP2f1)的三维结构模建和与维甲酸的对接研究

利用同源模建和分子动力学模拟方法，模建了细胞色素P450 2f1(CYP2f1)的三维结构．在此基础上，分析了活性位点的组成和结构，并进行了与小分子(维甲酸)的对接研究．研究结果表明，His328，Ser397和Arg417对复合物的结合起重要作用．     

How Mutations Affecting the Ligand-receptor Interactions:a Combined MD and QM/MM Calculation on CYP2E1 and Its Two Mutants

Cytochrome P450(CYP) 2E1 is a dual function monoxygenase with a crucial role in the metabolism of 6% of drugs on the market at present. The enzyme is of tremendous interest for its association with alcohol consumption, diabetes, obesity and fasting. Despite the abundant experimental mutagenesis data, the molecular origin and the structural motifs for the enzymatic activity deficiencies have not been rationalized at the atomic level. In this regard, we have investigated the effects of mutation on the structural and energetic characteristics upon single point mutations in CYP2E1, N219D and S366C. The molecular dynamics(MD) simulation combined with quantum mechanics/molecular mechanics(QM/MM) and noncovalent interaction(NCI) analysis was carried out on CYP2E1 and its two mutants. The results highlight the critical role of Phe207, which is responsible for both structural flexibility and energetic variation, shortening the gap between the theory and the experimentally observed results of enzymatic activity decrease. The underlying molecular mechanism of the enzymatic activity deficiencies for mutants may be attributed to the changes of spatial position of Phe207 in the two mutants. This work provides particular explanations to how mutations affect ligand-receptor interactions based on combined MD and QM/MM calculations. Furthermore, the mutational effects on the activity of CYP2E1 obtained in the present study are beneficial to both the experimental and the computational works of CYPs and may allow researchers to achieve desirable changes in enzymatic activity.     

Mechanism of the decrease in catalytic activity of human cytochrome P450 2C9 polymorphic variants investigated by computational analysis

Cytochrome P450 (CYP) is deeply involved in the metabolism of chemicals including pharmaceuticals. Therefore, polymorphisms of this enzyme have been widely studied to avoid unfavorable side effects of drugs in chemotherapy. In this work, we performed computational analysis of the mechanism of the decrease in enzymatic activity for three typical polymorphisms in CYP 2C9 species: *2, *3, and *5. Based on the equilibrated structure obtained by molecular dynamics simulation, the volume of the binding pocket and the fluctuation of amino residues responsible for substrate holding were compared between the wild type and the three variants. Further docking simulation was carried out to evaluate the appropriateness of the binding pocket to accommodate substrate chemicals. Every polymorphic variant was suggested to be inferior to the wild type in enzymatic ability from the structural viewpoint. F‐G helices were obviously displaced outward in CYP2C9*2. Expansion of the binding pocket, especially the space near F′ helix, was remarkable in CYP2C9*3. Disappearance of the hydrogen bond between K helix and β4 loop was observed in CYP2C9*5. The reduction of catalytic activity of those variants can be explained from the deformation of the binding pocket and the consequent change in binding mode of substrate chemicals. The computational approach is effective for predicting the enzymatic activity of polymorphic variants of CYP. This prediction will be helpful for advanced drug design because calculations forecast unexpected change in drug efficacy for individuals. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010