Calculating Ligand Binding and Charge Stabilization in Proteins

计算蛋白质中的配体结合和电荷稳定性

• 批准号: 1022208
• 负责人: Marilyn Gunner
• 金额: $ 108.13万
• 依托单位: CUNY City College
• 依托单位国家: 美国
• 项目类别: Continuing Grant
• 财政年份: 2010
• 资助国家: 美国
• 起止时间: 2010-08-01 至 2017-07-31
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=1022208&HistoricalAwards=false
• 关键词: Calculating; Ligand; Binding; Charge; Stabilization

Charged groups play important roles in protein structure and function. For example, acidic and basic amino acids constitute about 25% of the residues in an average protein and almost 50% of the residues in the active site of proteins where reactions are carried out. Key processes in biology involve electron and proton transfer reactions and binding or transport of charged ligands. This project focuses on the analysis of the change in free energy when charged groups associate with proteins, showing how these reaction equilibria are modified in biology. MCCE (Multi Conformation Continuum Electrostatics), developed at CCNY, which combines continuum electrostatics and molecular mechanics, will be the primary method used. The novel aspect of this work is that the protonation states of protein amino acids are allowed to remain at equilibrium with the bound and free states of a charged ligand. The affinities will be calculated with MCCE, with assists from Gaussian and GROMACS. In MCCE side chain positions and protonation states, ligand positions and site occupancy and where appropriate ligand redox and protonation states are all sampled together in one Monte Carlo (MC) simulation. This project will carry out calculations of: (A) the pKas of buried charged groups, comparing the behavior of wild-type and introduced residues. pKas will be calculated for a new dataset of 100 buried charge mutations in Staphylococcal nuclease that have been determined by Garcia-Moreno (Johns Hopkins). The calculated dielectric relaxation in the mutants and wild type residues will be compared. In addition, this aspect of the project will be used to optimize methodology for calculating the changes caused by mutations. (B) The chloride affinity at different binding sites will be calculated in halorhodopsin (HR). There are structures of the homologous, proton-pumping bacteriorhodopsin (BR) trapped in different intermediate pumping states. These structures will be used as a basis to model the conformational changes that move chloride through HR. The comparison of HR and BR will focus on the nature of the gate allowing ions or protons to be pumped against a concentration gradient. (C) Kds, Ems and pKas will be calculated for quinones in the two quinone binding sites (QA and QB) of photosynthetic reaction centers (RCs). The relationship of the changes in binding affinity (Kd) and shifts in reaction free energy (ÄG°) for protonation (ÄpKa) or redox reactions (ÄEm) for bound substrates or cofactors will be investigated. The thermodynamic relationship ÄG°(reaction in protein) - ÄG°(reaction in solution) = RTln(Kd(product) -Kd(reactant)) will be used to understand previously measured Em and pKa shifts for different quinones in RCs. The affinity of a series of neutral quinones in the QA and QB sites will be calculated and compared with measured values. The aspects of the structure that favor binding different quinones to each site will be identified. The affinity of quinone and semiquinone (SQ) will be compared to find the basis of Em shifts for different quinones in the QA site or for the same quinone in the QA and QB sites.This project will be carried out at CCNY, a public, city university, with an extremely diverse, multi-ethnic student body. The students to be mentored reflect this diversity. The MCCE program, which is publicly available, will be developed, distributed and maintained. Methods of calculating ligand binding and best practices for calculating pKas of mutated residues developed to carry out the work in this project will be incorporated into the distributed program. A web-based pKa databank that currently has 35,000 predicted pKas will be revamped. The physics based analysis of proteins pKas and Ems will be incorporated into an advanced, interdisciplinary biophysics class.

带电基团在蛋白质结构和功能中起着重要作用。例如，酸性和碱性氨基酸构成平均蛋白质中约25%的残基和进行反应的蛋白质活性位点中几乎50%的残基。生物学中的关键过程涉及电子和质子转移反应以及带电配体的结合或运输。 该项目的重点是分析当带电基团与蛋白质结合时自由能的变化，显示这些反应平衡在生物学中是如何被修改的。MCCE（多构象连续静电学），在CCNY开发，它结合了连续静电学和分子力学，将是主要的方法。这项工作的新颖之处在于，蛋白质氨基酸的质子化状态被允许与带电配体的结合和自由状态保持平衡。亲和力将使用MCCE计算，并得到Gaussian和GROMACS的协助。在MCCE侧链的位置和质子化状态，配体的位置和网站占有率，并在适当的配体氧化还原和质子化状态都在一个Monte Carlo（MC）模拟采样在一起。 本项目将进行以下计算：（A）掩埋带电基团的pKa，比较野生型和引入残基的行为。将计算Garcia-Moreno（Johns霍普金斯）测定的葡萄球菌核酸酶中100个埋藏电荷突变的新数据集的pKas。将比较突变体和野生型残基中计算的介电弛豫。此外，该项目的这一方面将用于优化计算突变引起的变化的方法。(B)将计算盐视紫红质（HR）中不同结合位点的氯离子亲和力。存在被捕获在不同中间泵浦态的同源质子泵浦细菌视紫红质（BR）的结构。这些结构将被用作基础来模拟通过HR移动氯的构象变化。HR和BR的比较将集中在允许离子或质子被泵送对抗浓度梯度的门的性质上。(C)将计算光合反应中心（RC）的两个醌结合位点（QA和QB）中醌的Kds、Ems和pKas。将研究结合底物或辅因子的结合亲和力（Kd）变化与质子化（Δ pKa）或氧化还原反应（Δ Em）的反应自由能（Δ G °）变化之间的关系。热力学关系Δ G °（蛋白质中的反应）-Δ G °（溶液中的反应）= RTln（Kd（产物）-Kd（反应物））将用于理解先前测量的RC中不同醌的Em和pKa位移。将计算QA和QB位点中一系列中性醌的亲和力，并与测量值进行比较。将确定有利于将不同醌结合到每个位点的结构方面。醌和半醌（SQ）的亲和力将进行比较，以找到不同的醌在QA网站或相同的醌在QA和QB网站的Em位移的基础。该项目将在CCNY，一个公共的，城市大学，具有非常多样化，多民族的学生团体进行。接受指导的学生反映了这种多样性。MCCE计划是公开的，将被开发，分发和维护。计算配体结合的方法和计算突变残基pKa的最佳实践将被纳入分布式程序中，以开展本项目中的工作。一个基于网络的pKa数据库，目前有35,000预测的pKa将被改造。蛋白质pKas和Ems的物理分析将被纳入先进的跨学科生物物理学课程。