Exendin-4中13号残基的分子动力学模拟

Exendin-4作为胰腺GLP-1受体上的一种有效的激活剂, 是一种含有39个氨基酸残基的多肽, 其第13号氨基酸Gln突变为Tyr, 使活性增强. 应用分子动力学模拟方法, 分别优化了突变前后, Exendin-4与蛋白的复合物结构, 并对整体结构的性质、静电势、相互作用模式及能量进行了分析. 阐明了Gln突变为Tyr的活性增强的内在原因, 结果表明, 突变的Exendin-4能够通过改变自身结构的局部柔性调整与蛋白受体相互作用, 从而可以改善Exendin-4与其蛋白受体的结合能力.     

细菌化学趋向性受体聚集体的相互作用

细菌化学趋向性受体的最小结构单元为二聚体,在细胞膜上这些二聚体会聚集成大团簇. X射线晶体结构和低分辨电镜结构测定表明,这些团簇有两类不同的形式,一种是在晶体结构中观察到的倒金字塔式二聚体的三聚体重复形成的聚集,另一种为由二聚体尾部相互盘绕形成的拉链状聚集. 有关拉链状聚集的详细分子模型目前尚不清楚. 本文使用蛋白质-蛋白质对接的方法研究了大肠杆菌丝氨酸化学趋向性受体Tsr 二聚体之间的相互作用. 分子对接计算表明,倒金字塔式聚集和拉链状聚集的基本复合物都是可以出现的,相应复合物的分子动力学模拟表明这些结构都具有一定的稳定性. 对于所获得的拉链状聚集体的基本复合物结构模型进行了详细的二聚体作用界面分析,发现二聚体间主要通过静电和疏水作用形成复合物,其中Arg388、Phe373 和Ile377是形成拉链状聚集的关键作用残基. 所建立的Tsr 拉链状聚集的结构模型有助于揭示细菌化学趋向性受体在细胞膜上聚集的分子机制,为进一步的聚集理论及模拟研究提供了基础.     

残基突变对P53-DNA结合域肽段构象影响的分子动力学模拟

基于分子动力学模拟方法研究了R249S、R248W 和G245S 突变对P53-DNA 结合域肽段(残基230-258)结构的影响. 采用GROMACS 软件包和GROMOS 43A1分子力场, 分别对野生型wtP53肽段、单点突变型P53-R249S肽段、两点突变型P53-R249S/R248W肽段和三点突变型P53-R249S/R248W/G245S肽段进行了4组独立的分子动力学模拟, 每组体系模拟时间为500 ns. 研究结果表明: R249S单残基突变影响肽段残基形成二级结构的情况, 但不改变肽段三维结构的模式, 同时使该肽段结构相对稳定; R249S和R248W两残基同时突变会加剧R249S突变对肽段的影响, 同时导致三维结构发生较大变化, 构象弯曲呈现双turn 结构, 肽段稳定性进一步增大; G245S突变对肽段的影响与R249S和R248W同时突变对其结构的影响相反, 在两残基突变的基础上, G245S突变会使原突变引起的变化减弱甚至消失, 同时使得该肽段结构稳定性减小. 该研究对认识肿瘤致病分子机制和设计新药物有重要意义.     

细菌化学趋向性受体聚集体的相互作用

细菌化学趋向性受体的最小结构单元为二聚体,在细胞膜上这些二聚体会聚集成大团簇.X射线晶体结构和低分辨电镜结构测定表明,这些团簇有两类不同的形式,一种是在晶体结构中观察到的倒金字塔式二聚体的三聚体重复形成的聚集,另一种为由二聚体尾部相互盘绕形成的拉链状聚集.有关拉链状聚集的详细分子模型目前尚不清楚.本文使用蛋白质-蛋白质对接的方法研究了大肠杆菌丝氨酸化学趋向性受体Tsr二聚体之间的相互作用.分子对接计算表明,倒金字塔式聚集和拉链状聚集的基本复合物都是可以出现的,相应复合物的分子动力学模拟表明这些结构都具有一定的稳定性.对于所获得的拉链状聚集体的基本复合物结构模型进行了详细的二聚体作用界面分析,发现二聚体间主要通过静电和疏水作用形成复合物,其中Arg388、Phe373和Ile377是形成拉链状聚集的关键作用残基.所建立的Tsr拉链状聚集的结构模型有助于揭示细菌化学趋向性受体在细胞膜上聚集的分子机制,为进一步的聚集理论及模拟研究提供了基础.     

亮氨酸拉链型脂肽对脂质体温敏性调节的分子模拟

含亮氨酸拉链型脂肽的温敏性脂质体被认为是抗癌药物的优良载体。亮氨酸拉链型脂肽的主要氨基酸残基序列为[VAQLEVK-VAQLESK-VSKLESK-VSSLESK],嵌入脂质体后可以有效改善脂质体的温敏性。本文首先采用隐式溶剂副本交换分子动力学方法,对N端修饰的亮氨酸拉链单链的折叠状态进行了模拟,得到了亮氨酸拉链单链的转变温度。并对包含该种新型亮氨酸拉链型脂肽的DPPC脂质体进行常规分子动力学模拟,研究了2种不同头基的亮氨酸拉链型脂肽(ALA,C3CO)二聚体嵌入后DPPC脂质体的相转变温度变化,证明了亮氨酸拉链型脂肽对于该脂质体温敏性的控制作用。利用这一规律,可以对亮氨酸拉链型脂肽进行优化改良,得到效果更佳的温敏脂质体,对于抗癌药物载体的开发有着重要的意义。     

胰高血糖素样肽-1受体胞外区域突变体复合物的动力学研究

采用分子动力学模拟方法研究了胰高血糖素样肽-1(GLP-1)与GLP-1受体(GLP-1R)胞外区域的相互作用.结果表明,配体的结合导致受体的构象发生改变,Loop2区域的氨基酸Pro90和Trp91以及C末端的Glu128向配体移动.根据保守位点突变受体(P73A,V81L,Y88A,P90A和W91A)后所得多肽模拟数据,发现Loop2区域在突变体中的结构和柔性均发生了明显变化,Trp91和Tyr88的突变将导致配体亲和力丧失.研究结果证明,P73A突变型受体和野生型受体分别与配体相互作用后,二者数值差别不大,因此Pro73不是关键残基;V81L突变体则会导致配体亲和力的丧失.该结果为GLP-1药物设计提供了重要理论依据.     

分子动力学模拟残基突变对芳香烃受体配体结合区的影响

芳香烃受体(Aryl hydrocarbon receptor,AhR)属于配体依赖性的转录因子蛋白。本文通过对AhR配体结合区域(Ligand binding domain,LBD)的结构功能及物种特异性分析,发现在其结合腔口有一些关键残基可能起到"门控"作用,进一步将野生型(WT)和3个突变模型(Phe289Ala、Tyr316Ala、Ile319Ala)进行分子动力学模拟,从蛋白稳定性、蛋白结构变化、蛋白结合腔变化及蛋白和配体结合能力4个方面分析3个残基的门控作用。研究发现,Phe289、Tyr316、Ile319氨基酸残基通过形成疏水作用为AhR LBD起到"门控"作用;而将这些氨基酸分别突变后,其蛋白稳定性降低,整体运动性增加,配体亲和力减弱,其中Tyr316、Ile319对腔内体积影响较大,Phe289使腔内环境稳定性降低。本研究可为基于芳香烃受体的药物设计提供相关理论指导。     

分子动力学模拟Cu2+对α-突触核蛋白(1-17)肽段构象变化的影响

利用分子动力学模拟研究铜离子(Cu2+)对α-突触核蛋白1-17号氨基酸肽段(α-synuclein(1-17))构象变化的影响,采用GROMOS 43A1力场对Cu2+-α-synuclein(1-17)复合体和α-synuclein(1-17)肽段单体分别进行了6组独立的分子动力学模拟,每组模拟时间为500ns,总模拟时间为3μs.研究结果表明:Cu2+与α-synuclein(1-17)肽段结合使其更易向β折叠片结构折叠,促进了其二级结构的形成,增强了构象的稳定性;Cu2+增大了α-synuclein肽段疏水残基的溶剂可及表面积,增强了其疏水残基的暴露程度.自由能分析指出,Cu2+-α-synuclein(1-17)复合体的自由能比α-synuclein(1-17)肽段低,构象稳定,采样空间紧密,其自由能极小构象为β折叠片结构.构象聚类分析进一步表明,Cu2+使得α-synuclein(1-17)肽段构象趋于稳定.总之,Cu2+诱导固有无序蛋白α-synuclein(1-17)肽段由无序向有序转变,降低了构象的自由能,同时Cu2+增强了α-synuclein(1-17)肽段的疏水性,使得α-synuclein肽段因疏水作用更倾向于形成β折叠片结构,加速其疏水性聚集.     

分子动力学模拟研究点突变(Met108→Leu108)对树胶醛糖结合蛋白与配体作用的影响

为了研究点突变(Met108→Leu108)对树胶醛糖结合蛋白(ABP)与配体结合能力的影响,对ABP、ABP结合树胶醛糖复合物及ABP结合半乳糖复合物以及它们各自的突变体分别进行60 ns的分子动力学模拟.模拟结果表明,108号残基突变前后,电子等排体的两个氨基酸残基,使蛋白与配体间的范德华相互作用发生明显变化,同时导致蛋白的内部运动也发生变化,进而影响蛋白与配体的相互作用.进一步分析表明,突变前后的蛋白构象变化都趋向于两个结构域张开,而与配体的结合可减缓张开程度.     

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

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

Effect of cobratoxin binding on the normal mode vibration within acetylcholine binding protein

Recent crystal structures of the acetylcholine binding protein (AChBP) have revealed surprisingly small structural alterations upon ligand binding. Here we investigate the extent to which ligand binding may affect receptor dynamics. AChBP is a homologue of the extracellular component of ligand-gated ion channels (LGICs). We have previously used an elastic network normal-mode analysis to propose a gating mechanism for the LGICs and to suggest the effects of various ligands on such motions. However, the difficulties with elastic network methods lie in their inability to account for the modest effects of a small ligand or mutation on ion channel motion. Here, we report the successful application of an elastic network normal mode technique to measure the effects of large ligand binding on receptor dynamics. The present calculations demonstrate a clear alteration in the native symmetric motions of a protein due to the presence of large protein cobratoxin ligands. In particular, normal-mode analysis revealed that cobratoxin binding to this protein significantly dampened the axially symmetric motion of the AChBP that may be associated with channel gating in the full nAChR. The results suggest that alterations in receptor dynamics could be a general feature of ligand binding.     

Top leads for swine influenza A/H1N1 virus revealed by steered molecular dynamics approach

Since March 2009, the rapid spread of infection during the recent A/H1N1 swine flu pandemic has raised concerns of a far more dangerous outcome should this virus become resistant to current drug therapies. Currently oseltamivir (tamiflu) is intensively used for the treatment of influenza and is reported effective for 2009 A/H1N1 virus. However, as this virus is evolving fast, some drug-resistant strains are emerging. Therefore, it is critical to seek alternative treatments and identify roots of the drug resistance. In this paper, we use the steered molecular dynamics (SMD) approach to estimate the binding affinity of ligands to the glycoprotein neuraminidase. Our idea is based on the hypothesis that the larger is the force needed to unbind a ligand from a receptor the higher its binding affinity. Using all-atom models with Gromos force field 43a1 and explicit water, we have studied the binding ability of 32 ligands to glycoprotein neuraminidase from swine flu virus A/H1N1. The electrostatic interaction is shown to play a more important role in binding affinity than the van der Waals one. We have found that four ligands 141562, 5069, 46080, and 117079 from the NSC set are the most promising candidates to cope with this virus, while peramivir, oseltamivir, and zanamivir are ranked 8, 11, and 20. The observation that these four ligands are better than existing commercial drugs has been also confirmed by our results on the binding free energies obtained by the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method. Our prediction may be useful for the therapeutic application.     

Computer‐aided Discovery of Potential Inhibitors for Transthyretin‐related Amyloidosis

Transthyretin (TTR) is a homotetrameric plasma protein associated with human amyloid diseases. Although Tafamidis has been recently approved for the treatment of TTR familial amyloid polyneuropathy (FAP), there is still a need for more effective drugs in the treatment of TTR amyloidosis diseases. In this study, a computer‐aided approach combining molecular docking, virtual screening and molecular dynamics (MD) simulations was employed to identify potent TTR amyloidosis inhibitors from National Cancer Institute (NCI), Maybridge and Asdi databases. A receptor‐specific scoring function was also developed using comparative binding energy (COMBINE) method to accurately predict the inhibitory activities for the selected compounds during virtual screening. The developed receptor‐specific scoring function demonstrated good predictive ability by yielding strong correlation coefficients between experimental activities and estimated activities for 32 training set and 9 test set compounds, respectively. Moreover, it was successfully applied to rank the candidate compounds from structure‐based virtual screening. Finally, three compounds (NSC220163, MFCD00276817 and SPB06319) were identified as potential leads, which exhibited higher predicted inhibitory activities and higher binding affinities in comparison to the Tafamidis. Our results further suggest that halogen bonding interaction plays a crucial role in stabilizing the TTR‐inhibitor complex. These results indicate that our computational approach could effectively discover more potent TTR amyloidosis inhibitors, which can be further validate by in vitro and in vivo biological tests.     

Binding mode prediction and MD/MMPBSA-based free energy ranking for agonists of REV-ERBα/NCoR

The knowledge of the free energy of binding of small molecules to a macromolecular target is crucial in drug design as is the ability to predict the functional consequences of binding. We highlight how a molecular dynamics (MD)-based approach can be used to predict the free energy of small molecules, and to provide priorities for the synthesis and the validation via in vitro tests. Here, we study the dynamics and energetics of the nuclear receptor REV-ERBα with its co-repressor NCoR and 35 novel agonists. Our in silico approach combines molecular docking, molecular dynamics (MD), solvent-accessible surface area (SASA) and molecular mechanics poisson boltzmann surface area (MMPBSA) calculations. While docking yielded initial hints on the binding modes, their stability was assessed by MD. The SASA calculations revealed that the presence of the ligand led to a higher exposure of hydrophobic REV-ERB residues for NCoR recruitment. MMPBSA was very successful in ranking ligands by potency in a retrospective and prospective manner. Particularly, the prospective MMPBSA ranking-based validations for four compounds, three predicted to be active and one weakly active, were confirmed experimentally.     

Detection of isoforms of recombinant human erythropoietin by various plant lectins after isoelectric focusing

A screening method to determine the binding behavior of lectins toward recombinant human erythropoietin (rHuEPO) was developed. Twenty-three different lectins were tested for this purpose. rHuEPO isoforms were separated by isoelectric focusing using the International Olympic Committee (IOC) and World Anti-Doping Agency (WADA) accredited method for the direct detection of the prohibited doping substance erythropoietin (EPO). For the visualization of the rHuEPO isoforms lectins were used instead of antibodies. Optimization of the screening protocol enabled the detection of a maximum number of rHuEPO isoforms. By means of this protocol information about the binding properties of a lectin toward each individual rHuEPO isoform was accessible. All evaluated lectins showed significant differences in their binding behavior. The most intense response was obtained with WGA, DSL, PHA-E, LEL, PSA, and LCA. While WGA, DSL, PHA-E, and LEL were able to bind all isoforms detected by the standard antibody, LCA and PSA demonstrated a clear preference for rHuEPO isoforms located in the more basic region of the electropherogram. Further lectins tested were ConA, succWGA, PHA-L, RCA, SNA, MAA, STL, ECL, GSL-II, SJA, SBA, UEA-I, Jacalin, PNA, DBA, GSL-I, and VVA. Compared to the lectins mentioned above, they showed reduced sensitivity. Endogenous and recombinant EPO only differ in the composition of their N- and O-glycan moieties. As lectins possess the unique ability to recognize subtle differences in glycan substructures, they represent an interesting approach for their structural characterization. Furthermore, they might be useful for affinity enrichment/purification of rHuEPO in doping control.     

Insights into the binding of GABA to the insect RDL receptor from atomistic simulations: a comparison of models

The resistance to dieldrin (RDL) receptor is an insect pentameric ligand-gated ion channel (pLGIC). It is activated by the neurotransmitter γ-aminobutyric acid (GABA) binding to its extracellular domain; hence elucidating the atomistic details of this interaction is important for understanding how the RDL receptor functions. As no high resolution structures are currently available, we built homology models of the extracellular domain of the RDL receptor using different templates, including the widely used acetylcholine binding protein and two pLGICs, the Erwinia Chrysanthemi ligand-gated ion channel (ELIC) and the more recently resolved GluCl. We then docked GABA into the selected three dimensional structures, which we used as starting points for classical molecular dynamics simulations. This allowed us to analyze in detail the behavior of GABA in the binding sites, including the hydrogen bond and cation-π interaction networks it formed, the conformers it visited and the possible role of water molecules in mediating the interactions; we also estimated the binding free energies. The models were all stable and showed common features, including interactions consistent with experimental data and similar to other pLGICs; differences could be attributed to the quality of the models, which increases with increasing sequence identity, and the use of a pLGIC template. We supplemented the molecular dynamics information with metadynamics, a rare event method, by exploring the free energy landscape of GABA binding to the RDL receptor. Overall, we show that the GluCl template provided the best models. GABA forming direct salt-bridges with Arg211 and Glu204, and cation-π interactions with an aromatic cage including Tyr109, Phe206 and Tyr254, represents a favorable binding arrangement, and the interaction with Glu204 can also be mediated by a water molecule.     

Mechanism‐Dependent Modulation of Ultrafast Interfacial Water Dynamics in Intrinsically Disordered Protein Complexes

The recognition of intrinsically disordered proteins (IDPs) is highly dependent on dynamics owing to the lack of structure. Here we studied the interplay between dynamics and molecular recognition in IDPs with a combination of time‐resolving tools on timescales ranging from femtoseconds to nanoseconds. We interrogated conformational dynamics and surface water dynamics and its attenuation upon partner binding using two IDPs, IBB and Nup153FG, both of central relevance to the nucleocytoplasmic transport machinery. These proteins bind the same nuclear transport receptor (Importinβ) with drastically different binding mechanisms, coupled folding–binding and fuzzy complex formation, respectively. Solvent fluctuations in the dynamic interface of the Nup153FG‐Importinβ fuzzy complex were largely unperturbed and slightly accelerated relative to the unbound state. In the IBB‐Importinβ complex, on the other hand, substantial relative slowdown of water dynamics was seen in a more rigid interface. These results show a correlation between interfacial water dynamics and the plasticity of IDP complexes, implicating functional relevance for such differential modulation in cellular processes, including nuclear transport.     

Monitoring the Reversibility of GPCR Signaling by Combining Photochromic Ligands with Label-free Impedance Analysis

G protein-coupled cell surface receptors (GPCR) trigger complex intracellular signaling cascades upon agonist binding. Classic pharmacological assays provide information about binding affinities, activation or blockade at different stages of the signaling cascade, but real time dynamics and reversibility of these processes remain often disguised. We show that combining photochromic NPY receptor ligands, which can be toggled in their receptor activation ability by irradiation with light of different wavelengths, with whole cell label-free impedance assays allows observing the cell response to receptor activation and its reversibility over time. The concept demonstrated on NPY receptors may be well applicable to many other GPCRs providing a deeper insight into the time course of intracellular signaling processes.