类ACE抑制肽的合成及其活性分析

血管紧张素转换酶(angiotensin converting enzyme,ACE)通过作用于维持血压正常的肾素-血管紧张系统(rennin-angiotensin system, RAS)和激肽释放酶 激肽系统(kallikrein-kinin system, KKS),使其失衡导致血压升高．而ACE活性抑制肽可以竞争性地与ACE的活性中心结合,从而抑制ACE的活性,使血压降低．天然来源的ACE抑制肽与传统的降压药物相比效果较好,无毒副作用,对正常血压没有影响,对于高血压的治疗和人类健康具有重要意义. 本文以酪蛋白中提取的ACE活性抑制肽KVLPVP为先导肽,根据ACE抑制肽的结构特点,设计合成一系列的类ACE肽(similar ACE-like peptides). 利用反相高效液相色谱法(RP-HPLC)直接测定其体外ACE抑制活性. 结果表明,当芳香性的氨基酸残基Phe、Tyr、His和疏水性Val残基位于C-端时会提高多肽的ACE抑制活性,尤其是His位于C 端时,ACE抑制活性更强. 通过对比先导肽与所合成的类ACE肽的ACE活性抑制率,可以发现,类ACE肽的ACE活性抑制率均高于先导肽．基于不同氨基酸残基位于C-端时对多肽的ACE抑制活性的研究,可以为降血压药物分子设计和筛选提供基础.     

应用定点突变与分子模建方法研究蛋白酶抑制剂eglin C突变体对蛋白酶furin和kexin的抑制活力

蛋白质前体加工酶参与许多重要蛋白质闪体的加工成熟过程,哺乳动物来源的furin和酵母中的kexin是该家族的重要成员。首先人工合成了编码枯草杆菌蛋白酶抑制剂eglin C的基因片段,组装后在大肠杆菌中得到表达。以定点突变方法在野生型eglin C抑制活性中心的P1、P2和P4位引入碱性氨基酸残基可以将其改造为很强的furin抑制剂（Ki约10^-9mol/L）,和kexin抑制剂（Ki约10^-11mol/L）。同时根据枯草杆菌蛋白酶和eglin C复合物的晶体结构,计算机同源模建了前体加工酶与eglin C突变体结构之间的相互作用,并结合实验数据得到以下结果：（1）P1位引入的碱性残基是该抑制剂活力的前提;（2）P4位碱性残基的引入可以极大地提高抑制剂活力约两个数量级;（3）P2 的碱性残基将有效提高抑制剂的活力。然而同时可以破坏抑制剂本身的稳定性。（4）野生型P3位的疏水性残基参与抑制剂活性环附近疏水核心的构成。     

牛心线粒体ATP酶与其抑制蛋白结合特性的研究

利用荧光偏振光谱和生物分子相互作用分析技术（ＢＩＡ）对线粒体ＡＴＰ酶（Ｆ１）与其抑制蛋白（ＩＦ１）的结合特性进行了研究。结果表明,在弱酸性条件下Ｍｇ－ＡＴＰ能够大大促进ＩＦ１与Ｆ１的结合,结合的ＩＦ１对Ｆ１－ＡＴＰ酶表现出最大的抑制活性。在弱酸性没有Ｍｇ－ＡＴＰ存在的条件下,ＩＦ１也可与Ｆ１发生一定量的结合（２０～３０％）,Ｆ１的水解活性也受到同等程度的抑制,ＩＦ１的抑制活性与ＩＦ１、Ｆ１的结合其间显示了密切的量效关系,提示在ＩＦ１与Ｆ１的作用过程中不存在无抑制活性的结合过渡态。Ｚｎ２＋能够抑制ＩＦ１组氨酸残基的质子化,从而阻断ＩＦ１与Ｆ１之间的结合反应     

基于分子对接及酶抑制动力学探究红景天提取物体外α-葡萄糖苷酶抑制活性CSCD

为探究红景天提取物体外α-葡萄糖苷酶抑制活性及具有抑制作用的主要活性物质。该研究首先构建体外α-葡萄糖苷酶抑制体系,测定其抑制活性,通过酶抑制动力学判断抑制类型,然后采用UHPLC-QE-MS、分子对接进一步探索提取物中抑制α-葡萄糖苷酶的主要活性物质。结果表明,红景天提取物对α-葡萄糖苷酶具有较好的抑制效果,IC_(50)为1.538 mg/mL,抑制类型为竞争与非竞争性混合可逆抑制;UHPLC-QE-MS共检测出1245种化合物,其中脂肪酸类、萜类及其衍生物、黄酮及类黄酮类为化合物最多的3类,分别有107种、85种、66种,其次还鉴定出酚类、氨基酸类、糖类等多类物质;分子对接显示,20种相对含量较高的化合物中11种可与α-葡萄糖苷酶结合,(+)-表儿茶素结合能(-17.08 kJ/mol)最低、结合活性最佳,咖啡酸形成氢键最多为5个,分别与His-515、Arg-437、Glu-432、His-348残基相连,咖啡酸、L-苹果酸、槲皮素和酪醇具有相同结合位点Arg-437。研究旨为天然α-葡萄糖苷酶抑制剂的开发以及红景天资源利用提供了基础研究。     

大肠杆菌PGDH末端缺失突变体的构建及抗反馈抑制效应分析

为了获得具有抗反馈抑制性质的大肠杆菌磷酸甘油酸脱氢酶（PGDH, d-3-phosphoglycerate dehydrogenase, EC 1.1.1.95）,通过对其碱基序列和蛋白质结构分析,用PCR突变法构建突变酶M1（缺失第410位氨基酸）、M2（缺失407~410位氨基酸）、M3（缺失337~410位氨基酸）。M0(野生型)及各突变型基因与pET22b(+)载体连接后,表达融合蛋白。在非变性条件下,由NTA-Ni镍离子螯合亲和层析柱纯化野生型和突变体的酶蛋白。酶活性测定结果表明,M1、M2蛋白酶均保持了原有的野生型磷酸甘油酸脱氢酶活性,且部分解除了终产物L-丝氨酸的反馈抑制作用;M3蛋白酶完全解除了终产物的反馈抑制作用,但酶本身的催化活性略有降低（为野生型的83%）。M0、M1、M2菌株PGDH与L-丝氨酸结合的Ki值分别约为7 μmol/L、20 μmol/L、50 μmol/L,说明该酶C-末端1~4个氨基酸残基对L-丝氨酸和调控区的结合有重要影响。     

HCV NS3/4A蛋白酶与Faldaprevir类似物的分子识别机制及其复合物运动模式

Faldaprevir类似物(Faldaprevir analogue molecule,FAM)能有效抑制HCV NS3/4A蛋白酶的催化活性,是一种潜在抗HCV先导化合物。通过生物信息学统计分析了已报道的HCV NS3/4A蛋白酶晶体结构,得到了FAM-HCV NS3/4A蛋白酶晶体结构。对FAM-HCV NS3/4A蛋白酶复合物进行了20.4 ns的分子动力学模拟,重点从氢键和结合自由能两个角度分析了二者分子识别中的关键残基及结合驱动力。氢键和范德华力是促使FAM特异性结合到蛋白V132?S139、F154?D168、D79?D81和V55的活性口袋中的主要驱动力,这与实验数据较为吻合。耐药性突变实验分析了R155K、D168E/V和V170T定点突变对FAM分子识别的影响,为可能存在的FAM耐药性提供了分子依据。最后,用自由能曲面和构象聚类两个方法探讨了体系的构象变化,给出体系的4种优势构象,为后续的基于HCV NS3/4A蛋白酶结构的Faldaprevir类似物抑制剂分子设计提供一定的理论帮助。     

Kazal型蛋白酶抑制剂结构与功能研究进展

蛋白酶抑制剂广泛存在于生物体内,在许多生命活动过程中发挥必不可少的作用,特别是对蛋白酶活性进行精确调控。其中Kazal型蛋白酶抑制剂是最重要的、研究最为广泛的酶抑制剂之一,该类抑制剂一般由一个或几个结构域组成,每一个结构域具有保守的序列和分子构象,同时发现该类抑制剂与蛋白酶作用的结合部位高度易变,它们大多数暴露于与溶剂接触的环上,其中P1部位是抑制作用的关键部位,抑制剂的专一性由P1部位氨基酸残基的性质决定,其它残基取代结合部位残基对抑制剂-酶的结合常数有显著的影响。Laskowski算法可直接从Kazal型丝氨酸蛋白酶抑制剂的序列推测其与6种丝氨酸蛋白酶之间的抑制常数(Ki)。目前在生物体内发现大量的Kazal型蛋白酶抑制剂,并证实其有重要的生物学功能。     

家蝇半胱氨酸蛋白酶抑制剂基因MdCPI的克隆、表达及活性检测

半胱氨酸蛋白酶抑制剂是具有抑制半胱氨酸蛋白酶活性的一类蛋白超家族。本研究根据EST序列信息,通过RACE技术克隆得到1条家蝇Musca domestica半胱氨酸蛋白酶抑制剂基因MdCPI,该基因含有1个357 bp开放阅读框,编码118个氨基酸残基,推导的多肽N端17个氨基酸残基为信号肽序列。同源分析表明,MdCPI 氨基酸序列与红尾肉蝇Sarcophaga crassipalpis的CPI相似性最高（identity=51%）。以邻接法（NJ）构建的系统树表明,家蝇与其他双翅目昆虫CPI起源于共同的祖先,属于I25A型蛋白家族。为了解家蝇CPI对半胱氨酸蛋白酶的抑制活性,构建pET-17b-MdCPI表达载体,并转入大肠杆菌Escherichia coli BL21（DE3）进行重组表达。研究发现1 μg重组家蝇CPI能够抑制约14 μg木瓜蛋白酶的水解活性。结果表明MdCPI确属CPI家族,可能同其他家族成员具有相似的功能,参与免疫及生理调控。这些结果为研究MdCPI在家蝇体内作用机制奠定了基础。     

从大豆中纯化棉铃虫组织蛋白酶B抑制因子HCB-SoyI

棉铃虫组织蛋白酶B（Helicoverpa armigera cathepsin B,HCB） 在胚胎发育过程中降解卵黄蛋白为氨基酸,供给胚胎发育必需的养料,是棉铃虫胚胎正常发育的重要因素。许多植物种子中存在蛋白酶抑制因子,如大豆和向日葵种子。用棉铃虫组织蛋白酶B为酶源,通过硫酸铵沉淀、排阻层析和离子交换层析等方法,从大豆中分离纯化到了一种对HCB有抑制活性的蛋白酶抑制因子,命名为HCB-SoyI。用牛血清白蛋白为底物进一步证明该抑制因子对HCB具有抑制作用。该抑制因子的纯化为进一步克隆其基因奠定了基础,为转基因抗虫作物研究提供了新的候选靶标。     

Api6/AIM/Spα调节免疫和脂质代谢的生物学功能研究

凋亡抑制因子6（apoptosis inhibitor 6,Api6）,又称作AIM/Spα,是清道夫受体富含半胱氨酸残基超家族新成员。Api6/AIM/Spα由巨噬细胞特异性表达,具有抑制CD4＋/CD8＋双阳性胸腺细胞、T淋巴细胞、NKT淋巴细胞和巨噬细胞凋亡的作用。作为模式识别受体,Api6/AIM/Spα直接与病原体相关分子模式LPS/LTA结合,在机体固有免疫和适应性免疫中发挥重要的作用。近年研究发现,Api6/AIM/Spα可以通过抑制动脉粥样硬化斑块部位巨噬细胞凋亡加重动脉粥样硬化早期斑块的进展,也可以通过抑制脂肪酸合成酶（FAS）的生物学活性提高脂肪细胞的脂解作用,在肥胖的进展中发挥重要作用。重点综述了Api6/AIM/Spα调节免疫和脂质代谢等生物学功能的研究进展。     

Insight into the inhibition mechanism and structure–activity relationship of 2,6-dipicolinic acid and its analogue to New Delhi metallo-β-lactamase-1

In previous literature, it was found that the activity of New Delhi Metallo-β-lactamase-1 (NDM-1) was inhibited by 2,6-dipicolinic acid (DPA) derivatives. To identify the mechanism of interaction between the inhibitors and NDM-1, molecular dynamics simulations were performed for the complex systems. Via the molecular modelling, inhibitors were found to be able bind to the region of catalytic activity of NDM-1. However, the detailed binding sites of the inhibitors differed with their structures. It was determined that His189, Lys211, Met248, Ser249, His250, and Ser251 are key residues for the binding of inhibitor 36 with NDM-1, and Asp124 is the only critical residue in the NDM-1-DPA complex. Furthermore, because of the interaction of the benzene ring in inhibitor 36 with the side chain of Lys211, inhibitor 36 can form 4 strong hydrogen bonds with protein. For the NDM-1-DPA complex, owing to the absence of the aniline group, DPA can only form a weak interaction with the residues around the binding site of NDM-1, except for Asp124, leading to a weaker inhibitory activity. Therefore, we believe that the strong interaction of the inhibitor with Lys211 results in effective inhibition, and the aniline group is the element required for the inhibitory activity.     

A comparative study of trypsin specificity based on QM/MM molecular dynamics simulation and QM/MM GBSA calculation

Hydrogen bonding and polar interactions play a key role in identification of protein-inhibitor binding specificity. Quantum mechanics/molecular mechanics molecular dynamics (QM/MM MD) simulations combined with DFT and semi-empirical Hamiltonian (AM1d, RM1, PM3, and PM6) methods were performed to study the hydrogen bonding and polar interactions of two inhibitors BEN and BEN1 with trypsin. The results show that the accuracy of treating the hydrogen bonding and polar interactions using QM/MM MD simulation of PM6 can reach the one obtained by the DFT QM/MM MD simulation. Quantum mechanics/molecular mechanics generalized Born surface area (QM/MM-GBSA) method was applied to calculate binding affinities of inhibitors to trypsin and the results suggest that the accuracy of binding affinity prediction can be significantly affected by the accurate treatment of the hydrogen bonding and polar interactions. In addition, the calculated results also reveal the binding specificity of trypsin: (1) the amidinium groups of two inhibitors generate favorable salt bridge interaction with Asp189 and form hydrogen bonding interactions with Ser190 and Gly214, (2) the phenyl of inhibitors can produce favorable van der Waals interactions with the residues His58, Cys191, Gln192, Trp211, Gly212, and Cys215. This systematic and comparative study can provide guidance for the choice of QM/MM MD methods and the designs of new potent inhibitors targeting trypsin.     

Molecular dynamics-solvated interaction energy studies of protein-protein interactions: the MP1-p14 scaffolding complex

Using the MP1-p14 scaffolding complex from the mitogen-activated protein kinase signaling pathway as model system, we explored a structure-based computational protocol to probe and characterize binding affinity hot spots at protein-protein interfaces. Hot spots are located by virtual alanine-scanning consensus predictions over three different energy functions and two different single-structure representations of the complex. Refined binding affinity predictions for select hot-spot mutations are carried out by applying first-principle methods such as the molecular mechanics generalized Born surface area (MM-GBSA) and solvated interaction energy (SIE) to the molecular dynamics (MD) trajectories for mutated and wild-type complexes. Here, predicted hot-spot residues were actually mutated to alanine, and crystal structures of the mutated complexes were determined. Two mutated MP1-p14 complexes were investigated, the p14(Y56A)-mutated complex and the MP1(L63A,L65A)-mutated complex. Alternative ways to generate MD ensembles for mutant complexes, not relying on crystal structures for mutated complexes, were also investigated. The SIE function, fitted on protein-ligand binding affinities, gave absolute binding affinity predictions in excellent agreement with experiment and outperformed standard MM-GBSA predictions when tested on the MD ensembles of Ras-Raf and Ras-RalGDS protein-protein complexes. For wild-type and mutant MP1-p14 complexes, SIE predictions of relative binding affinities were supported by a yeast two-hybrid assay that provided semiquantitative relative interaction strengths. Results on the MP1-mutated complex suggested that SIE predictions deteriorate if mutant MD ensembles are approximated by just mutating the wild-type MD trajectory. The SIE data on the p14-mutated complex indicated feasibility for generating mutant MD ensembles from mutated wild-type crystal structure, despite local structural differences observed upon mutation. For energetic considerations, this would circumvent costly needs to produce and crystallize mutated complexes. The sensitized protein-protein interface afforded by the p14(Y56A) mutation identified here has practical applications in screening-based discovery of first-generation small-molecule hits for further development into specific modulators of the mitogen-activated protein kinase signaling pathway.     

Direct Determination of Protonation States of Histidine Residues in a 2 Å Neutron Structure of Deoxy-Human Normal Adult Hemoglobin and Implications for the Bohr Effect

We have investigated the protonation states of histidine residues (potential Bohr groups) in the deoxy form (T state) of human hemoglobin by direct determination of hydrogen (deuterium) positions with the neutron protein crystallography technique. The reversible binding of protons is key to the allosteric regulation of human hemoglobin. The protonation states of 35 of the 38 His residues were directly determined from neutron scattering omit maps, with 3 of the remaining residues being disordered. Protonation states of 5 equivalent His residues—αHis20, αHis50, αHis89, βHis143, and βHis146—differ between the symmetry-related globin subunits. The distal His residues, αHis58 and βHis63, are protonated in the α1β1 heterodimer and are neutral in α2β2. Buried residue αHis103 is found to be protonated in both subunits. These distal and buried residues have the potential to act as Bohr groups. The observed protonation states of His residues are compared to changes in their pK_a values during the transition from the T to the R state and the results provide some new insights into our understanding of the molecular mechanism of the Bohr effect.     

Intermolecular insights into allosteric inhibition of histone lysine-specific demethylase 1

BackgroundHistone lysine-specific demethylase 1 (LSD1) has become a potential anticancer target for the novel drug discovery. Recent reports have shown that SP2509 and its derivatives strongly inhibit LSD1 as allosteric inhibitors. However, the binding mechanism of these allosteric inhibitors in the allosteric site of LSD1 is not known yet.MethodsThe stability and binding mechanism of allosteric inhibitors in the binding site of LSD1 were evaluated by molecular docking, ligand-based pharmacophore, molecular dynamics (MD) simulations, molecular mechanics generalized born surface area (MM/GBSA) analysis, quantum mechanics/molecular mechanics (QM/MM) calculation and Hirshfeld surface analysis.ResultsThe conformational geometry and the intermolecular interactions of allosteric inhibitors showed high binding affinity towards allosteric site of LSD1 with the neighboring amino acids (Gly358, Cys360, Leu362, Asp375 and Glu379). Meanwhile, MD simulations and MM/GBSA analysis were performed on selected allosteric inhibitors in complex with LSD1 protein, which confirmed the high stability and binding affinity of these inhibitors in the allosteric site of LSD1.ConclusionThe simulation results revealed the crucial factors accounting for allosteric inhibitors of LSD1, including different protein–ligand interactions, the positions and conformations of key residues, and the ligands flexibilities. Meanwhile, a halogen bond interaction between chlorine atom of ligand and key residues Trp531 and His532 was recurrent in our analysis confirming its importance.General significanceOverall, our research analyzed in depth the binding modes of allosteric inhibitors with LSD1 and could provide useful information for the design of novel allosteric inhibitors.     

Intrasteric control of AMPK via the gamma1 subunit AMP allosteric regulatory site

AMP-activated protein kinase (AMPK) is a alphabetagamma heterotrimer that is activated in response to both hormones and intracellular metabolic stress signals. AMPK is regulated by phosphorylation on the alpha subunit and by AMP allosteric control previously thought to be mediated by both alpha and gamma subunits. Here we present evidence that adjacent gamma subunit pairs of CBS repeat sequences (after Cystathionine Beta Synthase) form an AMP binding site related to, but distinct from the classical AMP binding site in phosphorylase, that can also bind ATP. The AMP binding site of the gamma(1) CBS1/CBS2 pair, modeled on the structures of the CBS sequences present in the inosine monophosphate dehydrogenase crystal structure, contains three arginine residues 70, 152, and 171 and His151. The yeast gamma homolog, snf4 contains a His151Gly substitution, and when this is introduced into gamma(1), AMP allosteric control is substantially lost and explains why the yeast snf1p/snf4p complex is insensitive to AMP. Arg70 in gamma(1) corresponds to the site of mutation in human gamma(2) and pig gamma(3) genes previously identified to cause an unusual cardiac phenotype and glycogen storage disease, respectively. Mutation of any of AMP binding site Arg residues to Gln substantially abolishes AMP allosteric control in expressed AMPK holoenzyme. The Arg/Gln mutations also suppress the previously described inhibitory properties of ATP and render the enzyme constitutively active. We propose that ATP acts as an intrasteric inhibitor by bridging the alpha and gamma subunits and that AMP functions to derepress AMPK activity.     

The effect of fucoidan on tyrosinase: computational molecular dynamics integrating inhibition kinetics

Fucoidan is a complex sulfated polysaccharide extracted from brown seaweed and has a wide variety of biological activities. In this study, we investigated the inhibitory effect of fucoidan on tyrosinase via a combination of inhibition kinetics and computational simulations. Fucoidan reversibly inhibited tyrosinase in a mixed-type manner. Time-interval kinetics showed that the inhibition was processed as first order with biphasic processes. For further insight, we simulated dockings with various sizes of molecular models (monomer to decamer) of fucoidan and showed that the best binding energy change results were obtained from the pentamer (?1.89?kcal/mol) and the hexamer (?1.97?kcal/mol) models of AutoDock Vina. The molecular dynamics simulation confirmed the binding mechanisms between tyrosinase and fucoidan and suggested that fucoidan mostly interacts with several residues including copper ions located in the active site. Our study suggests that fucoidan might be a potential natural antipigment agent.     

Minimal Zn(2+) binding site of amyloid-β

Zinc-induced aggregation of amyloid-β peptide (Aβ) is a hallmark molecular feature of Alzheimer's disease. Here we provide direct thermodynamic evidence that elucidates the role of the Aβ region 6-14 as the minimal Zn(2+) binding site wherein the ion is coordinated by His(6), Glu(11), His(13), and His(14). With the help of isothermal titration calorimetry and quantum mechanics/molecular mechanics simulations, the region 11-14 was determined as the primary zinc recognition site and considered an important drug-target candidate to prevent Zn(2+)-induced aggregation of Aβ.     

Computer simulation of zinc finger motifs from cellular nucleic acid binding protein and their interaction with consensus DNA sequences

We report here a computer simulation of the three-dimensional structures of seven zinc finger motifs from cellular nucleic acid binding protein involved in negative feedback inhibition of cholesterol biosynthesis. The structures are optimised using steric constraints imposed by tetrahedral coordination of the zinc ion with Cys and His residues, by molecular mechanics technique. We have also optimised the structure of a finger-I with GpT sequence. The model for the interaction of seven fingered protein with single-stranded d(GTGCGGTG) from sterol regulatory element (SRE) is given on the basis of these results. We also propose a scheme for recognition of a multifingered regulatory protein with small single-stranded DNA fragments.     

Allosteric inhibitor specificity of Thermotoga maritima 3-deoxy-d-arabino-heptulosonate 7-phosphate synthase

3-Deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAH7PS) catalyses the first step of the shikimate pathway for the biosynthesis of aromatic amino acids. Allosteric regulation of Thermotoga maritima DAH7PS is mediated by l-Tyr binding to a discrete ACT regulatory domain appended to a core catalytic (β/α)₈ barrel. Variants of T. maritima DAH7PS (TmaDAH7PS) were created to probe the role of key residues in inhibitor selection. Substitution Ser31Gly severely reduced inhibition by l-Tyr. In contrast both l-Tyr and l-Phe inhibited the TmaHis29Ala variant, while the variant where Ser31 and His29 were interchanged (His29Ser/Ser31His), was inhibited to a greater extent by l-Phe than l-Tyr. These studies highlight the role and importance of His29 and Ser31 for determining both inhibitory ligand selectivity and the potency of allosteric response by TmaDAH7PS.