蛋白质复合物结合自由能的预测

本文发展了一种合理的经验能量函数和结合自由能算法,并应用到２１个蛋白质复合物的结合自由能预测上。与现今发表的其他工作相比,我们的结果与实验的测定值符合得更好,平均预测精度为１．０ｋｃａｌ／ｍｏｌ,与实验值的相关达到９６％。应用本方法预测一个典型的蛋白质与其抑制剂复合物的结合自由能,在ＳＧＩ－ＩＭＰＡＣＴ Ｒ１００００工作站上约需２分钟。本文的结果还证实,与对蛋白折叠过程的认识不同,亲水原子在蛋白质     

预测蛋白质—蛋白质复合物结构的软对接算法

提出了一种有效的软对接算法 ,用于在已知受体和配体三维结构的条件下预测蛋白质 蛋白质复合物的结构。该算法的分子模型基于Janin提出的简化蛋白质模型 ,并在此基础上有所改进。对蛋白质分子表面的柔性氨基酸残基Arg、Lys、Asp、Glu和Met进行了特殊处理 ,通过软化分子表面的方式考虑了它们的侧链柔性。采用双重过滤技术来排除不合理的对接结构 ,此过滤技术是以复合物界面几何互补性和残基成对偏好性为标准提出的。对所得到的构象进行能量优化 ,之后用打分函数对这些结构进行排序 ,挑选出与复合物天然结构接近的构象。该打分函数包括静电、疏水和范德华相互作用能。用此算法对 2 6个复合物进行了结构预测 ,均找到了近天然结构 ,其中有 2 0个复合物的近天然结构排在了前 10位。改进的分子模型可以在一定程度上描述蛋白质表面残基侧链的柔性 ;双重过滤技术使更多的近天然结构保留下来 ,从而提高了算法成功预测的可能性 ;打分函数可以较合理地评价对接结构。总之 ,此种软对接算法能够对蛋白质分子识别的研究提供有益的帮助。     

蛋白质折叠、动力学以及蛋白质-配体结合的物理化学基础

蛋白质是生物体的重要组成部分并参与细胞内几乎所有的生物学过程.随着越来越多物种基因组序列的测定,准确理解基因产物的功能并探索蛋白质功能多样性的原因,已经成为当前的研究热点.为了研究蛋白质的功能,已有大量蛋白质的静态三维结构被测定.但是,蛋白功能最终受其动力学行为所控制,这包括折叠过程、构象波动、分子运动以及蛋白质-配体相互作用等.基于自由能图谱理论,本文深入讨论了蛋白质动力学的底层物理化学机制,并回答了以下问题:蛋白质为什么能够折叠、以及如何折叠成其天然三维结构?为什么蛋白质的动力学特征是固有的?其动力学行为如何控制蛋白质的功能?讨论结果将有助于后基因组时代生命科学研究中蛋白质结构-功能关系的理解.     

蛋白质复合物结构预测的集成分子对接方法

蛋白质分子间相互作用与识别是当前生命科学研究的热点,分子对接方法是研究这一问题的有效手段．为了推进分子对接方法的发展,欧洲生物信息学中心组织了国际蛋白质复合物结构预测（CAPRI）竞赛．通过参加CAPRI竞赛,逐步摸索出了一套用于蛋白质复合物结构预测的集成蛋白质一蛋白质分子对接方法HoDock,它包括结合位点预测、初始复合物结构采集、精细复合物结构采集、结构成簇和打分排序以及最终复合物结构挑选等主要步骤．本文以最近的CAPRI Target 39为例,具体说明该方法的主要步骤和应用．该方法在CAPRI Target 39竞赛中取得了比较好的结果,预测结构Model 10是所有参赛小组提交的366个结构中仅有的3个正确结构之一,其配体均方根偏差（L_Rmsd）为0．25nm．在对接过程中,首先用理论预测和实验信息相结合的方法来寻找蛋白质结合位点残基,确认CAPRI Target 39A链的A31TRP和A191HIS,B链的B512ARG和B531ARG为可能结合位点残基．同时,用ZDock程序做不依赖结合位点的初步全局刚性对接．然后,根据结合位点信息进行初步局部刚性对接,从全局和局部对接中挑出了11个初始对接复合物结构．进而,用改进的Rosetta Dock程序做精细位置约束对接,并对每组对接中打分排序前200的结构进行成簇聚类．最后,综合分析打分、成簇和结合位点三方面的信息,得到10个蛋白质复合物结构．竞赛结果表明,A191HIS,B512ARG和B531ARG三个结合位点残基预测正确,提交的10个蛋白质复合物结构中有5个复合物受体一配体界面残基预测成功率较高．与其他参赛小组的对接结果比较,表明HoDock方法具有一定优势．这些结果说明我们提出的集成分子对接方法有助于提高蛋白质复合物结构预测的准确率．     

引入显式水分子对接提高蛋白质配体复合物结构的预测精度

在蛋白质复合物界面一般都会存在着一定量的水分子,这些水分子通过空间占据和氢键方式影响蛋白质与配体的位置关系。应用现有的计算机方法研究蛋白质-配体对接时,一般不会显式地考虑水分子的作用。本文显式地将水分子引入蛋白质-配体对接过程,考虑水分子空间占据和氢键能量对复合物对接结构的影响,提出了一种包含水分子的蛋白质-配体对接算法。实验结果表明引入水分子使蛋白质-配体对接质量有明显提高。     

蛋白质-核酸复合物界面氨基酸与核苷酸偏好性分析

蛋白质-核酸相互作用机制到目前还不是很清楚,尤其是蛋白质与RNA的相互作用。目前,可得到的蛋白质-核酸复合物结构数据不断增多,作者收集了Protein Data Bank数据库中所有的蛋白质-核酸复合物结构数据,对复合物中结合残基和结合核苷酸的偏好性进行了统计分析。发现：1）不同功能的蛋白质-核酸复合物间的结合残基数量存在显著差异;2）在蛋白 质-DNA和蛋白质-RNA复合物界面,碱性氨基酸都是最受欢迎的;3）氨基酸的极性大小及方向在决定它是否与RNA分子进行结合时起到重要的作用,同时发现氨基酸侧链形成的空间位阻会影响氨基酸残基与RNA分子的相互作用;4）随着定义结合残基距离阈值的增大,其氨基酸使用的特异性降低,而受欢迎与不受欢迎的氨基酸种类均没有变化。     

光合生物色素-蛋白质复合物的多样性

通过对近年来放氧型光合生物有关研究结果的概括分析,提出放氧型光合生物光系统Ⅰ、光系统Ⅱ及捕光色素-蛋白质复合物具有多样性,并根据其PSⅠ复合物的77 K荧光发射特点,指出放氧型光合生物光合系统中的能量传递方式也可能存在多样性.     

蛋白质- 配体相互作用信息数据库研究进展

研究蛋白质和配体相互作用的结构和亲和力,不仅有助于了解蛋白质的功能,而且对药物研发以及药物作用机制的研究,也具有十 分重要的意义。目前,人们通过人工检索和半自动检索的方式,从文献和蛋白质数据库（Protein Data Bank,PDB）中获得了许多蛋白质- 配体亲和力信息和生物相关配体信息,并构建了许多蛋白质-配体相互作用的信息数据库。对3 个蛋白质-配体亲和力数据库和6 个蛋白质 晶体结构-配体生物相关性数据库进行介绍,并对其主要应用进行简述,希望能为实现高效准确地筛选和设计药物提供一定的帮助。     

顶- 底极性复合物与癌症的关系研究进展

顶-底极性是上皮细胞的一项主要特征,参与细胞形态、迁移、功能维持等多个生物学事件。上皮细胞顶-底极性复合物包括PAR复合物、SCRIB复合物和CRB复合物。丧失极性是细胞癌化的标志之一,并且在人类癌症中也发现了顶-底极性复合物的异常表达。本文将就目前有关顶-底极性复合物在癌症方面的研究进行综述,重点阐述顶-底极性复合物在肿瘤发生、发展过程中的作用及调控机制。     

互作蛋白质与复合物亚基的基因表达相关性比较分析

利用在多种应激条件下酵母的基因表达谱数据 ,分别计算互作蛋白质及复合物亚基编码基因的表达相关性。结果发现 ,相对于随机对照组 ,互作蛋白质的编码基因与蛋白质复合物的编码基因表达相关性均显著 (P <0 .0 1) ,即互作蛋白质及复合物亚基有共表达的倾向。通过比较 ,进一步发现蛋白质复合物亚基的基因表达相关性显著高于互作蛋白质的基因表达相关性 (P <0 .0 1) ,这与复合物亚基之间功能联系强于定义不甚确切的互作蛋白之间功能联系现象吻合。     

The Influence of Starting Coordinates in Free Energy Simulations of Ligand Binding to Dihydrofolate Reductase

Abstract

Molecular dynamics (MD) simulation combined with free energy perturbation (FEP) methods have been used to study the key structural differences and relative free energies for the binding of 6-methyl-N5-deazapterin (N8 protonated) and the 8-substituted compound, 6,8-dimethyl-N5-deazapterin (N3 protonated), to dihydrofolate reductase (DHFR). The free energy changes have been calculated using a variety of initial X-ray coordinates derived from bacterial and vertebrate (including human) DHFRs, and both with and without the reduced cofactor nicotinamide adenine dinucleotide (NADPH) bound. Given a sufficiently long simulation time for the FEP calculations (ca. 200 ps), all structures obtained after mutating 6,8-methyl-N5-deazapterin to 6-methyl-N5-deazapterin exhibited hydrogen bond formation between a backbone carbonyl group of DHFR and H(N8) of 6-methyl-N5-deazapterin, analogous to that found in the X-ray crystal structure of N5-deazafolate(N8 protonated) bound to human DHFR. However, both simulation and experiment suggest this additional H-bonding does not greatly enhance thermodynamic stability, with experiment indicating at most a factor of 2 difference in the relative affinities of the two ligand cations for vertebrate DHFR. Moreover, a binding differential of 10 in favour of the protonated 8-substituted compound is found experimentally for bacterial DHFR. The MD/FEP calculations suggest that the relative cost of ligand desolvation may largely cancel the lowering of free energy obtained in the active site, resulting in predicted binding differences within the range indicated by the vertebrate and bacterial DHFR experiments. However, the theoretical free energy changes could not be obtained with the accuracy required for the rationalization of the observed species dependence. While sampling difficulties are known to be inherent in MD simulation methodologies, these studies with several initial coordinate sets have demonstrated the contribution of coordinate choice to this problem. The results indicate that for demanding protein-ligand binding problems such as this one, the accuracy of the method may be no better than ± 2 kcal/mol.     

Free energy surface for rotamers of cis-enol malonaldehyde in aqueous solution studied by molecular dynamics calculations

Two-dimensional free energy surfaces for four rotamers of cis-enol malonaldehyde in water have been investigated by umbrella sampling molecular dynamics (MD) calculations. Biasing potential used in the umbrella sampling calculation was adopted to be the minus of conformational free energy preliminary obtained by the thermodynamic integration MD calculations for the rigid malonaldehyde whose stretching and bending were all fixed. The calculated free energy surface shows that, in water, a rotamer that has an intramolecular hydrogen bond is most stable among the rotamers. This is the same as that in vacuum, while order of relative stability of the other three rotamers is different in water and in vacuum. Inclusion of intramolecular vibrations changed the free energy surface little, i.e. at most 2.6 kJ/mol, which is much smaller than the solvation free energy. Free energy barriers from the most stable intramolecular hydrogen bonded rotamer to the others are lowered by hydration but they are still very high, >50 kJ/mol, such that the malonaldehyde molecule spends most of its time in water taking this conformation. Thus, reaction coordinate for intramolecular proton transfer reaction in water may be constructed assuming this rotamer.     

Computation of entropy contribution to protein-ligand binding free energy

The entropy contribution ΔS to protein-ligand binding free energy is studied for nine protein-lipid complexes. The entropy effect from the loss of the translational/rotational degrees of freedom (ΔS _tr) is calculated using the ideal gas approach. The change in the vibrational entropy (ΔS _vib) is calculated using the effective quantum oscillator approach with frequencies derived from the coordinate covariance matrix, so the inharmonic effects are taken into account. The change in the entropy of solvation (ΔS _solv) is considered using the binomial cell model (developed by the authors) for the hydrophobic effect. The entropy contribution from loss of conformations that are available for the free ligand (ΔS _conf) is also estimated. It is revealed that the negative in view of binding term ΔS _tr is only partly compensated by increasing of ΔS _vib, so T(ΔS _tr + ΔS _vib + ΔS _conf) < 0 for all complexes under investigation, but taking into account ΔS _solv leads to significantly increased ΔS. For all complexes except biotin-streptavidin, the results are found to be in reasonable agreement with experimental data. Published in Russian in Biokhimiya, 2007, Vol. 72, No. 7, pp. 963–973.     

Free Energy Perturbations in Ribonuclease T1 Substrate Binding. A Study of the Influence of Simulation Length,Internal Degrees of Freedom and Structure in Free Energy Perturbations

Abstract

We have studied the reliability of free energy perturbation calculations with respect to simulation protocol and simulation length in a real biological system, the binding of two different ligands to wildtype Ribonuclease T ₁ (RNT1) and to a mutant of RNT1 with Glu-46 replaced by Gln (RNT1-Gln46). The binding of the natural substrate 3′ GMP has been compared with the binding of a fluorescent probe, 2-aminopurine 3′ mono phosphate (2AP3′MP). These simulations predict that the mutant binds 2AP3′MP better than 3′GMP. Four complete free energy perturbations were performed that form a closed loop of four free energy differences, which should sum up to zero. This could be used as a tool for searching for systematic errors that are not detected by standard forward ? backward perturbations. The perturbation between 2AP3′MP and 3′GMP is quite straightforward and similar to what has been done by other groups. The perturbation between Glu46 and Gln46 is much more complex, involving as many as twelve atoms and a change of charge. This perturbation needs much longer simulation time, 500-600 ps, than used in free energy perturbations before. The increased simulation time is needed both to reach an equilibrium and to include several phases of fluctuations of the observed parameters in the production run. The extremely long simulation time is not such a severe problem as much of the work might be done on several different machines in parallel and cheap workstations are excellent for these calculations. Problems may also occur with values of the coupling parameter Λ close to 0 or 1, due to the high mobility of atoms as well as insertion/deletion in a previously unoccupied space involved in the perturbation.     

Free energy simulation to investigate the effect of amino acid sequence environment on the severity of osteogenesis imperfecta by glycine mutations in collagen

Molecular dynamics simulations were carried out to calculate the free energy change difference of two collagen-like peptide models for Gly --> Ser mutations causing two different osteogenesis imperfecta phenotypes. These simulations were performed to investigate the impact of local amino acid sequence environment adjacent to a mutation site on the stability of the collagen. The average free energy differences for a Gly --> Ser mutant relative to a wild type are 3.4 kcal/mol and 8.2 kcal/mol for a nonlethal site and a lethal site, respectively. The free energy change differences of mutant containing two Ser residues relative to the wild type at the nonlethal and lethal mutation sites are 4.6 and 9.8 kcal/mol, respectively. Although electrostatic interactions stabilize mutants containing one or two Ser residues at both mutation sites, van der Waals interactions are of sufficient magnitude to cause a net destabilization. The presence of Gln and Arg near the mutation site, which contain large and polar side chains, provide more destabilization than amino acids containing small and nonpolar side chains.     

Ligand Binding Site Structure Influences the Evolution of Protein Complex Function and Topology

1. Download high-res image (279KB)
2. Download full-size image

     

Computer assisted simulations and molecular graphics methods in molecular design

A survey is presented of computer-assisted statistical mechanical methods. The general theoretical background is described and special methods are discussed in detail. Practical procedures allowing for the calculation of binding energies are examined. A recent perturbation-relaxation procedure is summarized.     

Calibration of a Linear Free Energy Estimation Approach for Estimating Stability Constants for Metal-Bacterial Surface Complexes

In this study, we use the measured extent of metal adsorption onto bacterial cells to constrain a linear free energy relationship that allows estimation of unknown stability constants for metal-bacterial surface complexes based on the value of corresponding aqueous metal-acetate stability constants. A previous study (Fein et al., 2001 Fein, J B, Martin, A M and Wightman, P G. 2001. Metal adsorption onto bacterial surface: Development of a predictive approach. Geochim Cosmochim Acta, 65: 4267–4273. [Crossref], [Web of Science ®] , [Google Scholar]) used metal adsorption experiments to constrain a similar relationship, but the experiments were conducted using acid-washed bacteria, and subsequent evidence (Borrok et al., 2004a Borrok, D, Fein, J B, Tischler, M, O'Loughlin, E, Meyer, H, Liss, M and Kemner, K M. 2004a. The effect of acidic solutions and growth conditions on the adsorptive properties of bacterial surfaces. Chem Geol, 209: 107–119. [Crossref], [Web of Science ®] , [Google Scholar]) shows that the acid-washing step affects the extent of adsorption of a number of metals onto bacterial surfaces. We measured the adsorption of Zn, Ni, Co, Sr, and Nd onto Bacillus subtilis in 0.1 M NaClO₄ as a function of pH and metal:bacterial site ratio, using a non-electrostatic discrete four-site model of the bacterial protonation reactions as a basis for the metal adsorption modeling. The adsorption of the divalent cations (Zn, Ni, Co, and Sr) could best be modeled by considering adsorption reactions involving three sites on the bacterial surface; we used a one-site model to account for the Nd data that covered a more restricted pH range. The calculated stability constants for metal-Site 2 bacterial surface complexes are used to re-calibrate the linear free energy relationship previously defined by Fein et al. (2001) Fein, J B, Martin, A M and Wightman, P G. 2001. Metal adsorption onto bacterial surface: Development of a predictive approach. Geochim Cosmochim Acta, 65: 4267–4273. [Crossref], [Web of Science ®] , [Google Scholar]. There is a significant difference between the original and the re-calibrated lines for weakly binding cations such as Sr^{2 +}, but the difference becomes negligible for the stronger-binding cations. Because the linear free energy relationship defined in this study was calibrated from experiments that involved bacteria that were not exposed to acidic conditions, the estimated stability constant values that result from using this relationship are likely to reasonably reflect bacterial adsorption behaviors that occur in realistic geologic settings.     

Protein Structure Prediction with a Combined Solvation Free Energy-Molecular Mechanics Force Field

Abstract

Models of protein structure are frequently used to determine the physical characteristics of a protein when the crystal structure is not available. We developed a procedure to optimize such models, by use of a combined solvation free energy and molecular mechanics force field. Appropriately chosen atomic solvation parameters were defined using the criterion that the resulting protein model should deviate least from the crystal structure upon a forty picosecond molecular dynamics simulation carried out using the combined force field. Several tests were performed to refine the set of atomic solvation parameters which best complement the molecular mechanics forces. Four sets of parameters from the literature were tested and an empirically optimized set is proposed. The parameters are defined on a well characterized small molecule (alanyl dipeptide) and on the highly refined crystal structure of rat trypsin, and then tested on a second highly refined crystal structure of α-lytic protease. The new set of atomic solvation parameters provides a significant improvement over molecular mechanics alone in energy minimization of protein structures. This combined force field also has advantages over the use of explicit solvent as it is possible to take solvent effects into account during energetic conformational searching when modeling a homologous protein structure from a known crystal structure.