CN112233730B - 一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 - Google Patents
一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 Download PDFInfo
- Publication number
- CN112233730B CN112233730B CN202011111077.8A CN202011111077A CN112233730B CN 112233730 B CN112233730 B CN 112233730B CN 202011111077 A CN202011111077 A CN 202011111077A CN 112233730 B CN112233730 B CN 112233730B
- Authority
- CN
- China
- Prior art keywords
- receptor
- ligand
- pbdes
- enoyl
- acp reductase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000000694 effects Effects 0.000 title claims abstract description 39
- 101710198510 Enoyl-[acyl-carrier-protein] reductase [NADH] Proteins 0.000 title claims abstract description 30
- 238000010276 construction Methods 0.000 title claims abstract description 10
- 239000003446 ligand Substances 0.000 claims abstract description 69
- 230000027455 binding Effects 0.000 claims abstract description 33
- 238000000329 molecular dynamics simulation Methods 0.000 claims abstract description 31
- 238000003032 molecular docking Methods 0.000 claims abstract description 27
- 102000005962 receptors Human genes 0.000 claims description 77
- 108020003175 receptors Proteins 0.000 claims description 77
- 238000000034 method Methods 0.000 claims description 24
- BAWFJGJZGIEFAR-NNYOXOHSSA-O NAD(+) Chemical compound NC(=O)C1=CC=C[N+]([C@H]2[C@@H]([C@H](O)[C@@H](COP(O)(=O)OP(O)(=O)OC[C@@H]3[C@H]([C@@H](O)[C@@H](O3)N3C4=NC=NC(N)=C4N=C3)O)O2)O)=C1 BAWFJGJZGIEFAR-NNYOXOHSSA-O 0.000 claims description 19
- 229910052739 hydrogen Inorganic materials 0.000 claims description 16
- 239000001257 hydrogen Substances 0.000 claims description 16
- 230000003993 interaction Effects 0.000 claims description 14
- 125000000539 amino acid group Chemical group 0.000 claims description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 9
- 229910052736 halogen Inorganic materials 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 125000004429 atom Chemical group 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims description 4
- 230000000975 bioactive effect Effects 0.000 claims description 3
- 238000005421 electrostatic potential Methods 0.000 claims description 3
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims 2
- 238000000338 in vitro Methods 0.000 abstract description 8
- 238000005259 measurement Methods 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 abstract description 2
- XEFQLINVKFYRCS-UHFFFAOYSA-N Triclosan Chemical compound OC1=CC(Cl)=CC=C1OC1=CC=C(Cl)C=C1Cl XEFQLINVKFYRCS-UHFFFAOYSA-N 0.000 description 16
- 229960003500 triclosan Drugs 0.000 description 16
- 230000007613 environmental effect Effects 0.000 description 7
- 230000002401 inhibitory effect Effects 0.000 description 7
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 6
- RAXXELZNTBOGNW-UHFFFAOYSA-N imidazole Natural products C1=CNC=N1 RAXXELZNTBOGNW-UHFFFAOYSA-N 0.000 description 6
- 102000004169 proteins and genes Human genes 0.000 description 6
- 108090000623 proteins and genes Proteins 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- USIUVYZYUHIAEV-UHFFFAOYSA-N diphenyl ether Chemical class C=1C=CC=CC=1OC1=CC=CC=C1 USIUVYZYUHIAEV-UHFFFAOYSA-N 0.000 description 5
- 230000005764 inhibitory process Effects 0.000 description 5
- 102000004190 Enzymes Human genes 0.000 description 4
- 108090000790 Enzymes Proteins 0.000 description 4
- 235000014113 dietary fatty acids Nutrition 0.000 description 4
- 229930195729 fatty acid Natural products 0.000 description 4
- 239000000194 fatty acid Substances 0.000 description 4
- 150000004665 fatty acids Chemical class 0.000 description 4
- 229930027945 nicotinamide-adenine dinucleotide Natural products 0.000 description 4
- BOPGDPNILDQYTO-NNYOXOHSSA-N nicotinamide-adenine dinucleotide Chemical compound C1=CCC(C(=O)N)=CN1[C@H]1[C@H](O)[C@H](O)[C@@H](COP(O)(=O)OP(O)(=O)OC[C@@H]2[C@H]([C@@H](O)[C@@H](O2)N2C3=NC=NC(N)=C3N=C2)O)O1 BOPGDPNILDQYTO-NNYOXOHSSA-N 0.000 description 4
- 150000003384 small molecules Chemical class 0.000 description 4
- 239000000243 solution Substances 0.000 description 4
- 108010039731 Fatty Acid Synthases Proteins 0.000 description 3
- 238000003556 assay Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 239000011780 sodium chloride Substances 0.000 description 3
- RNFJDJUURJAICM-UHFFFAOYSA-N 2,2,4,4,6,6-hexaphenoxy-1,3,5-triaza-2$l^{5},4$l^{5},6$l^{5}-triphosphacyclohexa-1,3,5-triene Chemical compound N=1P(OC=2C=CC=CC=2)(OC=2C=CC=CC=2)=NP(OC=2C=CC=CC=2)(OC=2C=CC=CC=2)=NP=1(OC=1C=CC=CC=1)OC1=CC=CC=C1 RNFJDJUURJAICM-UHFFFAOYSA-N 0.000 description 2
- 101710146995 Acyl carrier protein Proteins 0.000 description 2
- 102100022089 Acyl-[acyl-carrier-protein] hydrolase Human genes 0.000 description 2
- 241000894006 Bacteria Species 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 125000002252 acyl group Chemical group 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000003063 flame retardant Substances 0.000 description 2
- 239000006166 lysate Substances 0.000 description 2
- 230000003647 oxidation Effects 0.000 description 2
- 238000007254 oxidation reaction Methods 0.000 description 2
- 239000011347 resin Substances 0.000 description 2
- 229920005989 resin Polymers 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 239000011534 wash buffer Substances 0.000 description 2
- QKNYBSVHEMOAJP-UHFFFAOYSA-N 2-amino-2-(hydroxymethyl)propane-1,3-diol;hydron;chloride Chemical compound Cl.OCC(N)(CO)CO QKNYBSVHEMOAJP-UHFFFAOYSA-N 0.000 description 1
- 229920000936 Agarose Polymers 0.000 description 1
- 102000015303 Fatty Acid Synthases Human genes 0.000 description 1
- -1 NAD+ small molecules Chemical class 0.000 description 1
- ZSLZBFCDCINBPY-ZSJPKINUSA-N acetyl-CoA Chemical compound O[C@@H]1[C@H](OP(O)(O)=O)[C@@H](COP(O)(=O)OP(O)(=O)OCC(C)(C)[C@@H](O)C(=O)NCCC(=O)NCCSC(=O)C)O[C@H]1N1C2=NC=NC(N)=C2N=C1 ZSLZBFCDCINBPY-ZSJPKINUSA-N 0.000 description 1
- 150000001413 amino acids Chemical class 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 231100000693 bioaccumulation Toxicity 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 125000000484 butyl group Chemical group [H]C([*])([H])C([H])([H])C([H])([H])C([H])([H])[H] 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011067 equilibration Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000002209 hydrophobic effect Effects 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- BPHPUYQFMNQIOC-NXRLNHOXSA-N isopropyl beta-D-thiogalactopyranoside Chemical compound CC(C)S[C@@H]1O[C@H](CO)[C@H](O)[C@H](O)[C@H]1O BPHPUYQFMNQIOC-NXRLNHOXSA-N 0.000 description 1
- 239000012139 lysis buffer Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 229950006238 nadide Drugs 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002957 persistent organic pollutant Substances 0.000 description 1
- 239000013641 positive control Substances 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- QAQREVBBADEHPA-IEXPHMLFSA-N propionyl-CoA Chemical compound O[C@@H]1[C@H](OP(O)(O)=O)[C@@H](COP(O)(=O)OP(O)(=O)OCC(C)(C)[C@@H](O)C(=O)NCCC(=O)NCCSC(=O)CC)O[C@H]1N1C2=NC=NC(N)=C2N=C1 QAQREVBBADEHPA-IEXPHMLFSA-N 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 239000006228 supernatant Substances 0.000 description 1
- 231100000331 toxic Toxicity 0.000 description 1
- 230000002588 toxic effect Effects 0.000 description 1
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
- C12N9/00—Enzymes; Proenzymes; Compositions thereof; Processes for preparing, activating, inhibiting, separating or purifying enzymes
- C12N9/0004—Oxidoreductases (1.)
- C12N9/001—Oxidoreductases (1.) acting on the CH-CH group of donors (1.3)
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Y—ENZYMES
- C12Y301/00—Hydrolases acting on ester bonds (3.1)
- C12Y301/01—Carboxylic ester hydrolases (3.1.1)
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B15/00—ICT specially adapted for analysing two-dimensional or three-dimensional molecular structures, e.g. structural or functional relations or structure alignment
- G16B15/30—Drug targeting using structural data; Docking or binding prediction
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/30—Detection of binding sites or motifs
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Health & Medical Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Biotechnology (AREA)
- Organic Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Genetics & Genomics (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Biophysics (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- General Engineering & Computer Science (AREA)
- Biochemistry (AREA)
- Pharmacology & Pharmacy (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Microbiology (AREA)
- Biomedical Technology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- General Chemical & Material Sciences (AREA)
- Enzymes And Modification Thereof (AREA)
Abstract
本发明公开了一种区分PBDEs衍生物对烯酰‑ACP还原酶活性效应模型的构建方法，包括受体与配体建模，分子对接，分子动力学模拟；本发明模型计算结果与体外结合活性测定实验结果一致，表明本发明模型的构建方法能够准确计算PBDEs衍生物对烯酰‑ACP还原酶的活性效应。
Description
技术领域
本发明属于有机污染物识别技术领域，具体涉及一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法。
背景技术
脂肪酸的合成由脂肪酸合酶催化是有机体必须的。脂肪酸合酶是一个复合酶系统，在自然界中能够催化乙酰辅酶A和丙酰辅酶A从头合成脂肪酸，在脂肪酸合酶系统中脂肪酸的延长分别由四种酶催化，其中烯酰-ACP还原酶位于脂肪酸延长的最后一步，是催化烯酰-ACP还原成酰基-ACP(ACP，acyl carrier protein,酰基载体蛋白)必须的酶。
多溴二苯醚(PBDEs，Poly Brominated Diphenyl Ethers)作为一种溴系阻燃剂的一大类阻燃物质，由于其优异的阻燃性能，已经越来越广泛的应用于各种消费产品当中。但是随着在环境样品中不断报道PBDEs的检出，该类化合物所造成的环境问题也越来越受到大家特别是环境科学的关注。PBDEs具有环境持久性，远距离传输，生物可累积性及对生物和人体具有毒害效应等特性，对其环境问题的研究已成为当前环境科学的一大热点。
多溴二苯醚(PBDEs)的衍生物包括羟基化多溴二苯醚(HO-PBDEs)和甲氧基化多溴二苯醚(MeO-PBDEs),在各种环境介质中被大量检出。
目前，PBDEs衍生物与烯酰-ACP还原酶的结合活性需要通过实验来进行测定，如何提供一种能够区分不同PBDEs衍生物对烯酰-ACP还原酶结合活性的预测模型是有待解决的技术问题。
发明内容
本发明提供一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，包括，受体与配体建模：删除烯酰-ACP还原酶的B链以及B链中结合的NAD+，删除受体中所有水分子，将结合有NAD+的A链作为下一步的分子对接使用的受体，并将NAD+与TCS加入氢原子，配体分子用sybyl7.3中的Sketch Molecule模块构建，并用sybyl7.3内Minimize模块进行结构优化；
分子对接：使用sybyl7.3中的Surflex-Dock模块将配体分子对接到受体中，对接时，使用Ligand模式，产生受体结合口袋，将配体与受体进行对接；每个配体化合物生成20个结合构象，以得分最高的构象作为最有可能的生物活性构象进行下一步的分子动力学模拟；
分子动力学模拟：在sybyl的Tripos力场下将配体分子分别用分子力学程序minimize进行优化，用powell能量梯度法，原子载电荷用Gasteiger-Huckel电荷，优化后的配体放入database中；在进行分子动力学模拟前，系统先经受5,000步的能量最小化，然后用NVT和NPT来平衡系统；每个系统在NVT系统中被模拟为500ps，在能量最小化的情况下从0到300K逐渐加热，并且在NPT系综中在300K下500ps平衡，然后采用Gromacs5.1.2软件包进行20ns的分子动力学模拟，在分子动力学模拟20ns时分析配体与受体之间的氢键、π-π堆积、卤键相互作用，并每2ps记录一次轨迹。
优选地，还包括在分析配体与受体之间的氢键、π-π堆积、卤键相互作用时，判断配体是否与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积；若配体未与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积，则分析构象变化，判断PRO154-ASN157helix-loop区、PHE251-ASN257 loop-helix区以及THR206-MET216 helix是否出现相互靠近的现象造成受体口袋入口闭合。
优选地，所述分子对接过程中，阈值(Threshold)设为0.5，膨胀系数(Bloat)为0。
优选地，所述分子对接，其中，所述产生受体结合口袋，为将TCS作为阳性物质，以TCS结合位点为中心产生受体结合口袋，预测构象与真实构象之间的RMSD为
优选地，所述在sybyl的Tripos力场下将配体分子分别用分子力学程序minimize进行优化，其中，最大优化次数为100步，0.05KJ/mol。
优选地，在分子动力学模拟过程中，配体分子以及受体分子中的NAD+分子使用Amber的GAFF力场，具体过程为先使用Gaussian09 D01优化分子的结构，计算静电势，拟合RESP电荷，然后用Ambertools中的模块antechamber生成resp文件，使用acpype.py脚本生成Gromacs可用的分子拓扑文件，受体蛋白质的拓扑文件用pdb2gmx生成，使用Amber99sb.ff力场，以距离受体蛋白质-配体分子表面至少2nm构建盒子，并填充TIPP3P水分子，即受体蛋白质与边缘至少有2nm水分子。
优选地，所述分子动力学模拟，包括以TCS结合位点为中心作为受体结合口袋。
本发明的有益效果：本发明模型计算结果与体外结合活性测定实验结果一致，表明本发明模型的构建方法能够准确计算PBDEs衍生物对烯酰-ACP还原酶的活性效应。
附图说明
为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，其中：
图1为TCS作为配体，分子动力学模拟20ns时的构象图。
图2为TCS作为配体，分子动力学模拟0ns时的构象图。
图3为6MeOBDE17在分子动力学模拟20ns时的构象图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合具体实施例对本发明的具体实施方式做详细的说明。
PBDE衍生物对烯酰-ACP还原酶抑制活性的预测方法：
(1)受体与配体建模：烯酰-ACP还原酶(PDB ID:1C14，FabI)的蛋白质分子结构来源于RCSB蛋白质数据库(http://www.rcsb.org/pdb)；在进行对接前删除烯酰-ACP还原酶的B链以及B链中结合的NAD+(烟酰胺腺嘌呤二核苷酸)，删除受体中所有水分子，将结合有NAD+的A链作为下一步的分子对接使用的受体，并将NAD+与TCS(Triclosan，三氯生)加入氢原子；所有配体小分子(即PBDEs衍生物及阳性对照物TCS)均用sybyl7.3(Tripos Inc.,St.Louis,MO,USA)中的Sketch Molecule模块构建，并用sybyl7.3内Minimize模块进行结构优化。
(2)分子对接：使用sybyl7.3中的Surflex-Dock模块将配体小分子对接到受体蛋白当中；对接时，使用Ligand模式，使用TCS作为阳性物质，以TCS结合位点为中心产生受体结合口袋，与烯酰-ACP还原酶结构进行对接，阈值(Threshold)为0.5，膨胀系数(Bloat)为0。
对接后得到的结构和结合有TCS的蛋白质晶体结构进行对比，结果两者之间的差异较小，预测构象与真实构象之间的RMSD(均方根偏差)仅为该值小于X射线晶体学分辨率/>这说明了对接过程中参数设置以及对接结果的可靠性。
对接过程中考虑环的柔性，每个配体化合物生成20个结合构象，以得分最高的构象作为最有可能的生物活性构象进行下一步的分子动力学模拟。结构-功能关系受结合后的构象影响，而结合后的构象与配体和受体之间的相互作用有关，在此我们选择分析分子对接后配体与受体之间的氢键、π-π堆积、卤键相互作用，并与实测的体外配体(PBDEs衍生物)与受体(烯酰-ACP还原酶)结合活性测定实验结果进行比较，如表1所示：
表1配体与受体之间的氢键、π-π堆积、卤键相互作用及体外结合活性测定
从表1可以看出，各个配体(PBDEs衍生物)对接的构象虽然与TCS相近，但是它们与蛋白质受体的相互作用却各不相同，配体大多与受体的PHE94残基形成相互作用，但是分子对接未能区分配体对受体抑制能力的有无。
(3)分子动力学(MD)模拟：①在sybyl的Tripos力场下将表1中16个配体分子分别用分子力学程序minimize进行优化，用powell能量梯度法，原子载电荷用Gasteiger-Huckel电荷，最大优化次数100步，0.05KJ/mol，优化后的16个配体放入database中，以TCS结合位点为中心作为受体结合口袋。
②采用Gromacs5.1.2软件包进行MD模拟；在进行动力学模拟前，系统先经受5,000步的能量最小化，然后进行两个阶段平衡：NVT和NPT来平衡系统；每个系统在NVT系统中被模拟为500ps，在能量最小化的情况下从0到300K逐渐加热，并且在NPT系综中在300K下500ps平衡，最后进行20ns的MD模拟，并每2ps记录一次轨迹。
为了评估受体-配体系统在MD模拟中的稳定性，计算系统骨架的均方根偏差(RMSD)，从RMSD的变化可以得出所有体系都在20ns之前达到了平衡。且每个系统RMSD值变化都小于0.25nm；分析受体蛋白质在20ns内的各个氨基酸残基的波动情况可以发现，在MD模拟过程中，残基波动主要集中在ASP98-ALA197loop区，PRO154-ASN157loop区以及结合口袋中THR194-MET206helix-loop区之间，这说明这些氨基酸残基在MD模拟过程中变动较大，比较灵活。
在模拟过程中配体小分子以及受体分子中的cofactorNAD+小分子使用Amber的GAFF力场，具体过程为先使用Gaussian09D01先优化小分子的结构，计算静电势，拟合RESP电荷，然后用Ambertools中的模块antechamber生成resp文件，使用acpype.py脚本生成Gromacs可用的分子拓扑文件，蛋白质的拓扑文件用pdb2gmx生成，使用Amber99sb.ff力场，以距离蛋白质-小分子配体表面至少2nm构建盒子，并填充TIPP3P水分子，即蛋白质与边缘至少有2nm水分子。
在模拟20ns时分析配体与受体之间的氢键、π-π堆积、卤键相互作用，并与实测的体外配体(PBDEs衍生物)与受体(烯酰-ACP还原酶)结合活性测定实验结果进行比较，如表2所示：
表2配体与受体之间的氢键、π-π堆积、卤键、疏水相互作用及体外结合活性测定实验结果
从表2可知，对于受体(烯酰-ACP还原酶)具有抑制效应的配体大多会与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积，这是预测配体与受体结合的关键因素，但仍有一个例外：配体4OHBDE49对于受体烯酰-ACP还原酶具有抑制作用，但其与NAD+或ALA196没有明显的相互作用。
③分析构象变化：受体与对其有抑制能力的配体结合后，ASP98-ALA197 helix-loop区，PRO154-ASN157helix-loop区以及结合口袋中THR194-MET206 helix-loop区从螺旋结构变为规则卷曲，或者螺旋缩短。而PRO154-ASN157 helix-loop区、PHE251-ASN257loop-helix区以及THR206-MET216helix出现相互靠近的现象，由于该段氨基酸是底物进入活性位点的通道，这几个区域的相互靠近使得受体口袋入口闭合，阻碍了酰基进入活性位点，从而产生抑制作用。而无抑制能力的配体大多无这种构象变化，如图1～3所示。图1为TCS作为配体，分子动力学模拟20ns时，由于上述几个区域相互靠近，口袋入口闭合，因而可预测出TCS对受体产生抑制作用；图2为分子动力学模拟在0ns时(初始状态)的构象，口袋入口是打开状态。图3为6MeOBDE17在分子动力学模拟20ns后的构象，模拟20ns后，口袋入口是打开状态，表明6MeOBDE17对受体无抑制效应。配体4OHBDE49由于出现了PRO154-ASN157helix-loop区、PHE251-ASN257loop-helix区以及THR206-MET216helix相互靠近的现象，使得口袋入口闭合，阻碍了酰基进入活性位点，因此，虽然4OHBDE49未与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积，但仍然对受体产生了抑制作用。
体外结合活性测定实验方法为：制备过表达His标记的FabI的大肠杆菌裂解液，当细菌光密度(OD600)达到0.6以后，用1mMIPTG，16℃诱导表达过夜，之后收集细菌并用PBS洗涤一次，在含有500mM NaCl、10mM Tris HCl和0.1％NP-40(pH＝8.5)的裂解缓冲液中裂解细菌，14000g离心30min，用Ni-NTA琼脂糖(Qiagen)从上清液中纯化FabI。将含Ni-NTA的树脂转移到Poly-Prep色谱柱(Bio-Red)，用洗涤缓冲液(50mM PBS，ph8.2，300mM NaCl和10mm咪唑)平衡两次。在4℃下用FabI裂解物孵育30分钟后，用洗涤缓冲液洗涤树脂两次，用含有50mM PBS、pH 8.2、300mM NaCl和250mM咪唑的缓冲液洗脱FabI蛋白。添加300μM丁基辅酶A，加入100nM重组FabI、400μM NADH和40μM NAD+混合物进行反应。NADH氧化为NAD+的速率用微板阅读器在340nm下每2分钟测量一次。为了测定PBDEs衍生物对FabI酶活性的抑制作用，将不同浓度的PBDEs(0.014-10μM)与100nM重组FabI在冰上孵育60分钟进行反应，之后测定NADH氧化为NAD+的速率。
综上，本发明提供了一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，本发明通过受体与配体建模、分子对接、分子动力学模拟，分子动力学模拟包括配体优化、模拟20ns时分析配体与受体之间的氢键、π-π堆积相互作用，判断配体是否与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积，若形成上述相互作用则表明配体对受体具有抑制作用，若未形成上述相互作用则进一步分析构象变化，判断PRO154-ASN157helix-loop区、PHE251-ASN257 loop-helix区以及THR206-MET216 helix是否出现相互靠近的现象造成受体口袋入口闭合，若出现则表明配体对受体产生抑制作用，否则说明配体未对受体产生抑制作用，本发明配体PBDEs衍生物包括2OHBDE123、4MeOBDE17、4OHBDE49、5MeOBDE47、6MeOBDE17、6MeOBDE47、6MeOBDE85、6MeOBDE90、6OHBDE17、6OHBDE47、6OHBDE90、6OHBDE137、6OHClBDE17、6OHClBDE68，本发明模型计算结果与体外结合活性测定实验结果一致，表明本发明模型的构建方法能够准确计算PBDEs衍生物对烯酰-ACP还原酶的活性效应。
应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
Claims (8)
1.一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：包括，受体与配体建模：删除烯酰-ACP还原酶的B链以及B链中结合的NAD+，删除受体中所有水分子，将结合有NAD+的A链作为下一步的分子对接使用的受体，并将NAD+与TCS加入氢原子，配体分子用sybyl7.3中的Sketch Molecule模块构建，并用sybyl7.3内Minimize模块进行结构优化；
分子对接：使用sybyl7.3中的Surflex-Dock模块将配体分子对接到受体中，对接时，使用Ligand模式，产生受体结合口袋，将配体与受体进行对接；每个配体化合物生成20个结合构象，以得分最高的构象作为最有可能的生物活性构象进行下一步的分子动力学模拟；
分子动力学模拟：在sybyl的Tripos力场下将配体分子分别用分子力学程序minimize进行优化，用powell能量梯度法，原子载电荷用Gasteiger-Huckel电荷，优化后的配体放入database中；在进行分子动力学模拟前，系统先经受5,000步的能量最小化，然后用NVT和NPT来平衡系统；每个系统在NVT系统中被模拟为500ps，在能量最小化的情况下从0到300K逐渐加热，并且在NPT系综中在300K下500ps平衡，然后采用Gromacs5.1.2软件包进行20ns的分子动力学模拟，在分子动力学模拟20ns时分析配体与受体之间的氢键、π-π堆积、卤键相互作用，并每2ps记录一次轨迹；
在分析配体与受体之间的氢键、π-π堆积、卤键相互作用时，判断配体是否与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积；若配体未与受体上结合的NAD+或受体的氨基酸残基ALA196形成氢键或π-π堆积，则分析构象变化，判断PRO154-ASN157helix-loop区、PHE251-ASN257 loop-helix区以及THR206-MET216 helix是否出现相互靠近的现象造成受体口袋入口闭合。
2.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述分子对接过程中，阈值(Threshold)设为0.5，膨胀系数(Bloat)为0。
3.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述分子对接，其中，所述产生受体结合口袋，为将TCS作为阳性物质，以TCS结合位点为中心产生受体结合口袋，预测构象与真实构象之间的RMSD为
4.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述在sybyl的Tripos力场下将配体分子分别用分子力学程序minimize进行优化，其中，最大优化次数为100步，0.05KJ/mol。
5.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：在分子动力学模拟过程中，配体分子以及受体分子中的NAD+分子使用Amber的GAFF力场，具体过程为先使用Gaussian09 D01优化分子的结构，计算静电势，拟合RESP电荷，然后用Ambertools中的模块antechamber生成resp文件，使用acpype.py脚本生成Gromacs可用的分子拓扑文件，受体蛋白质的拓扑文件用pdb2gmx生成，使用Amber99sb.ff力场，以距离受体蛋白质-配体分子表面至少2nm构建盒子，并填充TIPP3P水分子，即受体蛋白质与边缘至少有2nm水分子。
6.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述分子动力学模拟，包括以TCS结合位点为中心作为受体结合口袋。
7.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述配体，为PBDEs衍生物，包括2OHBDE123、4MeOBDE17、4OHBDE49、5MeOBDE47、6MeOBDE17、6MeOBDE47、6MeOBDE85、6MeOBDE90、6OHBDE17、6OHBDE47、6OHBDE90、6OHBDE137、6OHClBDE17、6OHClBDE68。
8.如权利要求1所述的区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法，其特征在于：所述受体，为对烯酰-ACP还原酶，包括FabI。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011111077.8A CN112233730B (zh) | 2020-10-16 | 2020-10-16 | 一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 |
| US17/093,580 US11908550B2 (en) | 2020-10-16 | 2020-11-09 | Method for determining different PBDEs derivatives from their response to activities of enoyl-ACP reductase |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011111077.8A CN112233730B (zh) | 2020-10-16 | 2020-10-16 | 一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112233730A CN112233730A (zh) | 2021-01-15 |
| CN112233730B true CN112233730B (zh) | 2023-11-28 |
Family
ID=74118868
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011111077.8A Active CN112233730B (zh) | 2020-10-16 | 2020-10-16 | 一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US11908550B2 (zh) |
| CN (1) | CN112233730B (zh) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003065978A2 (en) * | 2001-11-28 | 2003-08-14 | Structural Genomix, Inc. | Crystals and structures of a flavin mononucleotide binding protein (fmnbp) |
| CN101381894A (zh) * | 2008-05-30 | 2009-03-11 | 南京大学 | 一种有机物雌激素受体激动和拮抗作用的识别方法 |
| KR20090101652A (ko) * | 2008-03-24 | 2009-09-29 | 한국과학기술원 | 시스템 기법을 이용한 비브리오속 미생물의 약물 표적 예측 |
| EP2329018A1 (en) * | 2008-04-25 | 2011-06-08 | Novalis Ltd. | Method for the annotation of natural product gene-clusters and for the generation of novel biologically active chemical entities from dna sequences in silico |
| CN108351350A (zh) * | 2015-08-25 | 2018-07-31 | 杜克大学 | 使用rna指导型内切核酸酶改善基因组工程特异性的组合物和方法 |
| CN109285584A (zh) * | 2018-09-03 | 2019-01-29 | 南京大学 | 区分物质雄激素与抗雄激素效应模型的构建和应用 |
| CN110426512A (zh) * | 2019-05-21 | 2019-11-08 | 南京大学 | 区分过氧化酶激活增殖受体γ完全激动剂、部分激动剂、拮抗剂活性的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6684162B2 (en) * | 2000-05-08 | 2004-01-27 | Wyeth | Methods for identifying agents that interact with an active site of acyl carrier protein synthase-acyl carrier protein complex |
| US20050004225A1 (en) * | 2003-04-16 | 2005-01-06 | Balendiran Ganesaratnam K. | Oxidoreductase inhibitors and methods of screening and using thereof |
-
2020
- 2020-10-16 CN CN202011111077.8A patent/CN112233730B/zh active Active
- 2020-11-09 US US17/093,580 patent/US11908550B2/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2003065978A2 (en) * | 2001-11-28 | 2003-08-14 | Structural Genomix, Inc. | Crystals and structures of a flavin mononucleotide binding protein (fmnbp) |
| KR20090101652A (ko) * | 2008-03-24 | 2009-09-29 | 한국과학기술원 | 시스템 기법을 이용한 비브리오속 미생물의 약물 표적 예측 |
| EP2329018A1 (en) * | 2008-04-25 | 2011-06-08 | Novalis Ltd. | Method for the annotation of natural product gene-clusters and for the generation of novel biologically active chemical entities from dna sequences in silico |
| CN101381894A (zh) * | 2008-05-30 | 2009-03-11 | 南京大学 | 一种有机物雌激素受体激动和拮抗作用的识别方法 |
| CN108351350A (zh) * | 2015-08-25 | 2018-07-31 | 杜克大学 | 使用rna指导型内切核酸酶改善基因组工程特异性的组合物和方法 |
| CN109285584A (zh) * | 2018-09-03 | 2019-01-29 | 南京大学 | 区分物质雄激素与抗雄激素效应模型的构建和应用 |
| CN110426512A (zh) * | 2019-05-21 | 2019-11-08 | 南京大学 | 区分过氧化酶激活增殖受体γ完全激动剂、部分激动剂、拮抗剂活性的方法 |
Non-Patent Citations (8)
| Title |
|---|
| 5-芳亚甲基-N-(4-羧基苯基)罗丹宁衍生物的微波辅助合成及分子模拟研究;蒋军荣;徐峰;柯中炉;;化学试剂(第05期);481-484、487 * |
| Biochemical and molecular biomarkers in integument biopsies of free-ranging coastal bottlenose dolphins from southern Brazil;Righetti, BPH等;《CHEMOSPHERE》;139-149 * |
| Combined 3D-QSAR, molecular docking and molecular dynamics study on thyroid hormone activity of hydroxylated polybrominated diphenyl ethers to thyroid receptors β;Li, XL等;《TOXICOLOGY AND APPLIED PHARMACOLOGY》;300-307 * |
| Shrinivas D. Joshi等.3D-QSAR and molecular docking studies of 1,3,4-oxadiazoles containing substituted phenoxy fragment as inhibitors of enoyl-acyl carrier protein reductase from Escherichia coli.《Medicinal Chemistry Research 》.2014,4542–4558. * |
| 基于Metadynamic的分子动力学模拟建立受体活性区分模型;石来昊;《中国优秀硕士学位论文全文数据库》;B027-422 * |
| 多氯联苯鸡胚肝细胞生物活性非单调剂量-响应关系的QSARs研究;穆云松;张爱茜;高常安;彭素芬;王连生;;中国科学(B辑:化学)(第12期);1105-1112 * |
| 环境中多溴联苯醚(PBDEs)的代谢转化研究现状;温泉;张俊江;关淼;刘红玲;苏冠勇;于红霞;;环境监控与预警(第04期);34-41 * |
| 选择性磷酸二酯酶4D（PDE4D）抑制剂 的计算模型建立;刘景陶等;《分子科学学报》;249-256 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US20220119850A1 (en) | 2022-04-21 |
| CN112233730A (zh) | 2021-01-15 |
| US11908550B2 (en) | 2024-02-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ross et al. | Enhancing water sampling in free energy calculations with grand canonical Monte Carlo | |
| Gao et al. | A generalized hybrid orbital (GHO) method for the treatment of boundary atoms in combined QM/MM calculations | |
| Shi et al. | The 2.0 Å structure of malarial purine phosphoribosyltransferase in complex with a transition-state analogue inhibitor | |
| Erlanson et al. | Fragment-based drug discovery | |
| Cross et al. | FLAP: GRID molecular interaction fields in virtual screening. validation using the DUD data set | |
| Bjelic et al. | Cold adaptation of enzyme reaction rates | |
| Seifert et al. | Multiple molecular dynamics simulations of human p450 monooxygenase CYP2C9: the molecular basis of substrate binding and regioselectivity toward warfarin | |
| Urich et al. | De novo design of protein kinase inhibitors by in silico identification of hinge region-binding fragments | |
| Al-Sha’er et al. | Discovery of novel CDK1 inhibitors by combining pharmacophore modeling, QSAR analysis and in silico screening followed by in vitro bioassay | |
| Burger et al. | Quantum mechanics/molecular mechanics strategies for docking pose refinement: distinguishing between binders and decoys in cytochrome C peroxidase | |
| Nair et al. | Missing fragments: detecting cooperative binding in fragment-based drug design | |
| Mizutani et al. | Effective handling of induced‐fit motion in flexible docking | |
| Hou et al. | Protein flexibility in docking-based virtual screening: discovery of novel lymphoid-specific tyrosine phosphatase inhibitors using multiple crystal structures | |
| Srivastava et al. | Molecular docking studies on DMDP derivatives as human DHFR inhibitors | |
| Kacar et al. | Resurrecting ancestral genes in bacteria to interpret ancient biosignatures | |
| Wassermann et al. | Comprehensive analysis of single‐and multi‐target activity cliffs formed by currently available bioactive compounds | |
| Fraczek et al. | Assessing molecular docking tools for relative biological activity prediction: a case study of triazole HIV-1 NNRTIs | |
| Tan et al. | Computational methodologies for compound database searching that utilize experimental protein–ligand interaction information | |
| Kamat et al. | Enzymatic deamination of the epigenetic base N-6-methyladenine | |
| Ota et al. | Non‐Boltzmann thermodynamic integration (NBTI) for macromolecular systems: Relative free energy of binding of trypsin to benzamidine and benzylamine | |
| Kamath et al. | The tail wagging the dog: Insights into catalysis in R67 dihydrofolate reductase | |
| Shek et al. | Structural determinants for substrate selectivity in guanine deaminase enzymes of the amidohydrolase superfamily | |
| CN112233730B (zh) | 一种区分PBDEs衍生物对烯酰-ACP还原酶活性效应模型的构建方法 | |
| Feenstra et al. | Combining substrate dynamics, binding statistics, and energy barriers to rationalize regioselective hydroxylation of octane and lauric acid by CYP102A1 and mutants | |
| Vazquez et al. | Assessing the performance of mixed strategies to combine lipophilic molecular similarity and docking in virtual screening |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |