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

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

蛋白质-蛋白质分子对接方法中分子柔性处理与近天然结构筛选的研究

蛋白质-蛋白质分子对接方法是研究蛋白质分子间相互作用与识别的重要理论方法。该方法主要涉及复合物结合模式的构象搜索和近天然结构的筛选两个问题。在构象搜索中,分子柔性的处理是重点也是难点,围绕这一问题,近年来提出了许多新的方法。针对近天然结构的筛选问题,目前主要采用三种解决策略:结合位点信息的利用、相似结构的聚类和打分函数对结构的评价。本文围绕以上问题,就国内外研究进展和本研究小组的工作作详细的综述,并对进一步的研究方向进行了展望。     

蛋白质-蛋白质分子对接方法研究进展

蛋白质分子间相互作用与识别是21世纪生命科学研究的前沿和热点.分子对接方法是研究这一课题有效的计算机模拟手段.通常,蛋白质-蛋白质分子对接包括四个阶段:搜索受体与配体分子间的结合模式,过滤对接结构以排除不合理的结合模式,优化结构,用精细的打分函数评价、排序对接模式并挑选近天然构象.结合国内外研究蛋白质-蛋白质分子对接方法进展和本研究小组的工作,对以上四个环节做了详细的综述.另外,还分析了目前存在的主要问题,并提出对未来工作的展望.     

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

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

蛋白质- 核酸对接方法研究进展

分子对接技术作为预测蛋白质-核酸复合物结构的有效方法,为研究在生物学过程中蛋白质-核酸的相互作用提供了重要的工具。本文首先分析了当前蛋白质-核酸对接研究中的主要困难,例如构象变化和核糖磷酸骨架的带电性问题。然后从构象搜索、打分函数、柔性策略三个方面比较和总结了蛋白质-核酸对接中主要的计算方法。最后回顾了蛋白质-核酸对接计算模型的应用,并对未来的工作进行了展望。     

蛋白质-配体柔性对接综述

蛋白质与配体的相互作用在药物发现与筛选、功能食品开发等方面,都具有重要的科学意义。研究蛋白质-配体对接的目的在于通过已知配体和目标蛋白质的三维结构,运用计算手段来预测和评估蛋白质-配体复合物的三维构象,从而更好地理解蛋白质-配体之间的相互作用。随着已解析蛋白质单体结构的数量和计算能力的不断增加,蛋白质-配体对接也越来越现实并可靠。作者对蛋白质-配体柔性对接过程中所涉及的蛋白质的柔性、配体的准备、构象空间的采样方法、打分函数及其后期处理等方面进行了综述。     

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

蛋白质分子间相互作用与识别是当前生命科学研究的热点,分子对接方法是研究这一问题的有效手段．为了推进分子对接方法的发展,欧洲生物信息学中心组织了国际蛋白质复合物结构预测（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方法具有一定优势．这些结果说明我们提出的集成分子对接方法有助于提高蛋白质复合物结构预测的准确率．     

蛋白质与极性配体复合物的绝对结合自由能计算

介绍了用分子动力学模拟与热力学积分法相结合,模拟蛋白质与配体的绝对结合自由能的方法.通过分子转换法,使蛋白质分子（包括水分子）与配体小分子之间的相互作用逐渐减弱 （或增强）至完全消失（或完全出现）. 运用体约束方法,计算了配体与受体结合后平动、转动自由度的丧失即熵效应所引起的自由能变化.以胰蛋白酶双突变体（D189G/G226D）与极性配体苯甲脒为例,研究了蛋白质活性部位与极性配体的相互作用对结合自由能的影响,该复合物绝对结合自由能的模拟结果(－15.5 kJ/mol)与实验值(－10.5 kJ/mol)相近.     

hIL-6·hIL-6Rα·gp130三元复合物的结构预测

通过Ｄｅｌｐｈｉ方法分析人白介素６（ｈｕｍａｎ ｉｎｔｅｒｌｅｕｋｉｎ－６,ｈＩＬ－６）／人白介素６受体α亚基（ｈｕｍａｎ ｉｎｔｅｒｌｅｕｋｉｎ－６ ｒｅｃｅｐｔｏｒ α ｓｕｂｕｎｉｔ,ｈＩＬ－６Ｒα）复合物、ｇｐ１３０（βｓｕｂｕｎｉｔ）的空间构象的表现静电分布,利用分子对接方法研究ｇｐ１３０和ｈＩＬ－６／ｈＩＬ－６Ｒα复合物作用成三元复合物的空间构象,经过分子力学优化、分子动力学常温动态模     

侧链能量在蛋白质对接优化中的应用

一般的蛋白质对接程序能够提供大量的待选构象,但其中仅含有少量的正确构象。现在对接的主要工作在于如何从这些大量构象中挑出正确构象。我们先前的研究工作证明蛋白质界面比非界面表面具有更高的能量。在这里,我们使用由chen等人提出的一个用于检验、设计对接程序的蛋白质复合物标准库中的非抗原-抗体复合物,将侧链能量运用到对接中,并比较了侧链能量和残基配对倾向性、残基组成倾向性、残基保守性在对接中的表现。单独使用这四项的正确构象的平均百排分位排序分别为：38．6±19．6、26．3±20．8、22．7±16．6和37．8±26．1,但是对于个别蛋白,侧链能量的表现要优于其它的三个参数。我们将四个参数综合起来考虑,发展了一个新的打分函数,平均百排分位排序为22．2±7．8,并且提高了筛选效率。     

Discrimination of near-native protein structures from misfolded models by empirical free energy functions

Free energy potentials, combining molecular mechanics with empirical solvation and entropic terms, are used to discriminate native and near-native protein conformations from slightly misfolded decoys. Since the functional forms of these potentials vary within the field, it is of interest to determine the contributions of individual free energy terms and their combinations to the discriminative power of the potential. This is achieved in terms of quantitative measures of discrimination that include the correlation coefficient between RMSD and free energy, and a new measure labeled the minimum discriminatory slope (MDS). In terms of these criteria, the internal energy is shown to be a good discriminator on its own, which implies that even well-constructed decoys are substantially more strained than the native protein structure. The discrimination improves if, in addition to the internal energy, the free energy expression includes the electrostatic energy, calculated by assuming non-ionized side chains, and an empirical solvation term, with the classical atomic solvation parameter model providing slightly better discrimination than a structure-based atomic contact potential. Finally, the inclusion of a term representing the side chain entropy change, and calculated by an established empirical scale, is so inaccurate that it makes the discrimination worse. It is shown that both the correlation coefficient and the MDS value (or its dimensionless form) are needed for an objective assessment of a potential, and that together they provide much more information on the origins of discrimination than simple inspection of the RMSD-free energy plots.     

Empirical free energy calculation: comparison to calorimetric data.

An effective free energy potential, developed originally for binding free energy calculation, is compared to calorimetric data on protein unfolding, described by a linear combination of changes in polar and nonpolar surface areas. The potential consists of a molecular mechanics energy term calculated for a reference medium (vapor or nonpolar liquid), and empirical terms representing solvation and entropic effects. It is shown that, under suitable conditions, the free energy function agrees well with the calorimetric expression. An additional result of the comparison is an independent estimate of the side-chain entropy loss, which is shown to agree with a structure-based entropy scale. These findings confirm that simple functions can be used to estimate the free energy change in complex systems, and that a binding free energy evaluation model can describe the thermodynamics of protein unfolding correctly. Furthermore, it is shown that folding and binding leave the sum of solute-solute and solute-solvent van der Waals interactions nearly invariant and, due to this invariance, it may be advantageous to use a nonpolar liquid rather than vacuum as the reference medium.     

An evaluation of Poisson-Boltzmann electrostatic free energy calculations through comparison with experimental mutagenesis data

For systems involving highly and oppositely charged proteins, electrostatic forces dominate association and contribute to biomolecular complex stability. Using experimental or theoretical alanine-scanning mutagenesis, it is possible to elucidate the contribution of individual ionizable amino acids to protein association. We evaluated our electrostatic free energy calculations by comparing calculated and experimental data for alanine mutants of five protein complexes. We calculated Poisson-Boltzmann electrostatic free energies based on a thermodynamic cycle, which incorporates association in a reference (Coulombic) and solvated (solution) state, as well as solvation effects. We observe that Coulombic and solvation free energy values correlate with experimental data in highly and oppositely charged systems, but not in systems comprised of similarly charged proteins. We also observe that correlation between solution and experimental free energies is dependent on dielectric coefficient selection for the protein interior. Free energy correlations improve as protein dielectric coefficient increases, suggesting that the protein interior experiences moderate dielectric screening, despite being shielded from solvent. We propose that higher dielectric coefficients may be necessary to more accurately predict protein-protein association. Additionally, our data suggest that Coulombic potential calculations alone may be sufficient to predict relative binding of protein mutants.     

Automated computational framework for the analysis of electrostatic similarities of proteins

Charge plays an important role in protein-protein interactions. In the case of excessively charged proteins, their electrostatic potentials contribute to the processes of recognition and binding with other proteins or ligands. We present an automated computational framework for determining the contribution of each charged amino acid to the electrostatic properties of proteins, at atomic resolution level. This framework involves computational alanine scans, calculation of Poisson-Boltzmann electrostatic potentials, calculation of electrostatic similarity distances (ESDs), hierarchical clustering analysis of ESDs, calculation of solvation free energies of association, and visualization of the spatial distributions of electrostatic potentials. The framework is useful to classify families of mutants with similar electrostatic properties and to compare them with the parent proteins in the complex. The alanine scan mutants introduce perturbations in the local electrostatic properties of the proteins and aim in delineating the contribution of each mutated amino acid in the spatial distribution of electrostatic potential, and in biological function when electrostatics is a dominant contributing factor in protein-protein interactions. The framework can be used to design new proteins with tailored electrostatic properties, such as immune system regulators, inhibitors, and vaccines, and in guiding experimental studies. We present an example for the interaction of the immune system protein C3d (the d-fragment of complement protein C3) with its receptor CR2, and we discuss our data in view of a binding site controversy.     

Selecting near-native conformations in homology modeling: the role of molecular mechanics and solvation terms

A free energy function, combining molecular mechanics energy with empirical solvation and entropic terms, is used for ranking near-native conformations that occur in the conformational search steps of homology modeling, i.e., side-chain search and loop closure calculations. Correlations between the free energy and RMS deviation from the X-ray structure are established. It is shown that generally both molecular mechanics and solvation/entropic terms should be included in the potential. The identification of near-native backbone conformations is accomplished primarily by the molecular mechanics term that becomes the dominant contribution to the free energy if the backbone is even slightly strained, as frequently occurs in loop closure calculations. Both terms become equally important if a sufficiently accurate backbone conformation is found. Finally, the selection of the best side-chain positions for a fixed backbone is almost completely governed by the solvation term. The discriminatory power of the combined potential is demonstrated by evaluating the free energies of protein models submitted to the first meeting on Critical Assessment of techniques for protein Structure Prediction (CASP1), and comparing them to the free energies of the native conformations.     

Combination of scoring functions improves discrimination in protein-protein docking

Two structure-based potentials are used for both filtering (i.e., selecting a subset of conformations generated by rigid-body docking), and rescoring and ranking the selected conformations. ACP (atomic contact potential) is an atom-level extension of the Miyazawa-Jernigan potential parameterized on protein structures, whereas RPScore (residue pair potential score) is a residue-level potential, based on interactions in protein-protein complexes. These potentials are combined with other energy terms and applied to 13 sets of protein decoys, as well as to the results of docking 10 pairs of unbound proteins. For both potentials, the ability to discriminate between near-native and non-native docked structures is substantially improved by refining the structures and by adding a van der Waals energy term. It is observed that ACP and RPScore complement each other in a number of ways (e.g., although RPScore yields more hits than ACP, mainly as a result of its better performance for charged complexes, ACP usually ranks the near-native complexes better). As a general solution to the protein-docking problem, we have found that the best discrimination strategies combine either an RPScore filter with an ACP-based scoring function, or an ACP-based filter with an RPScore-based scoring function. Thus, ACP and RPScore capture complementary structural information, and combining them in a multistage postprocessing protocol provides substantially better discrimination than the use of the same potential for both filtering and ranking the docked conformations.     

Identification of protein-protein interaction sites from docking energy landscapes

Protein recognition is one of the most challenging and intriguing problems in structural biology. Despite all the available structural, sequence and biophysical information about protein-protein complexes, the physico-chemical patterns, if any, that make a protein surface likely to be involved in protein-protein interactions, remain elusive. Here, we apply protein docking simulations and analysis of the interaction energy landscapes to identify protein-protein interaction sites. The new protocol for global docking based on multi-start global energy optimization of an all-atom model of the ligand, with detailed receptor potentials and atomic solvation parameters optimized in a training set of 24 complexes, explores the conformational space around the whole receptor without restrictions. The ensembles of the rigid-body docking solutions generated by the simulations were subsequently used to project the docking energy landscapes onto the protein surfaces. We found that highly populated low-energy regions consistently corresponded to actual binding sites. The procedure was validated on a test set of 21 known protein-protein complexes not used in the training set. As much as 81% of the predicted high-propensity patch residues were located correctly in the native interfaces. This approach can guide the design of mutations on the surfaces of proteins, provide geometrical details of a possible interaction, and help to annotate protein surfaces in structural proteomics.     

Study of protein-protein interaction using conformational space annealing

We apply conformational space annealing (CSA), an efficient global optimization method, to the study of protein-protein interaction. The CSA is incorporated into the Tinker molecular modeling package along with a B-spline method for CAPRI Round 5 experiments. We have used an energy function for the protein-protein interaction that consists of electrostatic interaction, van der Waals interaction, and solvation energy terms represented by the occupancy desolvation method. The parameters of the AMBER94 all-atom empirical force field are used. Each energy term is calculated by precalculated grid potentials and B-spline method approximation. The ligand protein is placed inside a sphere of 50 A radius centered at an appropriate location, and the CSA rigid docking studies are carried out to find stable complexes. Up to 10 complexes are selected using the K-mean clustering method and biological information when available. These complexes are energy-minimized for further refinement by considering the flexibility of interacting proteins. The results show that the CSA method has a potential for the study of protein-protein interaction.     

Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes

Absolute binding free energy calculations and free energy decompositions are presented for the protein-protein complexes H-Ras/C-Raf1 and H-Ras/RalGDS. Ras is a central switch in the regulation of cell proliferation and differentiation. In our study, we investigate the capability of the molecular mechanics (MM)-generalized Born surface area (GBSA) approach to estimate absolute binding free energies for the protein-protein complexes. Averaging gas-phase energies, solvation free energies, and entropic contributions over snapshots extracted from trajectories of the unbound proteins and the complexes, calculated binding free energies (Ras-Raf: -15.0(+/-6.3)kcal mol(-1); Ras-RalGDS: -19.5(+/-5.9)kcal mol(-1)) are in fair agreement with experimentally determined values (-9.6 kcal mol(-1); -8.4 kcal mol(-1)), if appropriate ionic strength is taken into account. Structural determinants of the binding affinity of Ras-Raf and Ras-RalGDS are identified by means of free energy decomposition. For the first time, computationally inexpensive generalized Born (GB) calculations are applied in this context to partition solvation free energies along with gas-phase energies between residues of both binding partners. For selected residues, in addition, entropic contributions are estimated by classical statistical mechanics. Comparison of the decomposition results with experimentally determined binding free energy differences for alanine mutants of interface residues yielded correlations with r(2)=0.55 and 0.46 for Ras-Raf and Ras-RalGDS, respectively. Extension of the decomposition reveals residues as far apart as 25A from the binding epitope that can contribute significantly to binding free energy. These "hotspots" are found to show large atomic fluctuations in the unbound proteins, indicating that they reside in structurally less stable regions. Furthermore, hotspot residues experience a significantly larger-than-average decrease in local fluctuations upon complex formation. Finally, by calculating a pair-wise decomposition of interactions, interaction pathways originating in the binding epitope of Raf are found that protrude through the protein structure towards the loop L1. This explains the finding of a conformational change in this region upon complex formation with Ras, and it may trigger a larger structural change in Raf, which is considered to be necessary for activation of the effector by Ras.     

An explorative study on Staphylococcus aureus MurE inhibitor: induced fit docking,binding free energy calculation,and molecular dynamics

Staphylococcus aureus MurE enzyme catalyzes the addition of l-lysine as third residue of the peptidoglycan peptide moiety. Due to the high substrate specificity and its ubiquitous nature among bacteria, MurE enzyme is considered as one of the potential target for the development of new therapeutic agents. In the present work, induced fit docking (IFD), binding free energy calculation, and molecular dynamics (MD) simulation were carried out to elucidate the inhibition potential of 2-thioxothiazolidin-4-one based inhibitor 1 against S. aureus MurE enzyme. The inhibitor 1 formed majority of hydrogen bonds with the central domain residues Asn151, Thr152, Ser180, Arg187, and Lys219. Binding free-energy calculation by MM-GBSA approach showed that van der Waals (ΔG_vdW, ?57.30?kcal/mol) and electrostatic solvation (ΔG_solv, ?36.86?kcal/mol) energy terms are major contributors for the inhibitor binding. Further, 30-ns MD simulation was performed to validate the stability of ligand–protein complex and also to get structural insight into mode of binding. Based on the IFD and MD simulation results, we designed four new compounds D1–D4 with promising binding affinity for the S. aureus MurE enzyme. The designed compounds were subjected to the extra-precision docking and binding free energy was calculated for complexes. Further, a 30-ns MD simulation was performed for D1/4C13 complex.