Fundamentals of Ligand-Protein Interactions

配体-蛋白质复合物的两个伙伴的构象变化，小分子配体及其蛋白质结合位点（在许多情况下是酶的催化活性位点）是许多药物作用的核心方面，也是药物设计计算方法的关键挑战。在该领域最早的出版物之一中，我们“http://www.ncbi.nlm.nih.gov/sites/entrez?Db=pubmed和Cmd=ShowDetailView和TermToSearch=8581425和ordinalpos=47和itool=EntrezSystem2.PEntrez.Pubmed.Pubmed_ResultsPanel”。Pubmed_RVDocSum”表明，对于出现在蛋白质数据库（PDB）和剑桥结构数据库（CSD）中的一小组配体，柔性化合物通常不以其整体真空能构象与蛋白质结合，甚至通常不以任何局部真空能构象与蛋白质结合。虽然这项研究使用了当时最大的数据集和最好的方法，但实验数据库中的结构数量和可用的软件和硬件资源都呈指数级增长。因此，我们将通过对数量级更多的结构的分析和在计算量子化学理论的高水平上进行的计算来重新审视这个重要的主题。在该项目迄今为止取得的其他里程碑中，我们已经提取了最近在PDB的“http://ligand-expo.rcsb.org/”LigandExpo中可用的所有小分子配体。截至2008年5月，这是一组超过350,000个不同的3D坐标集。我们添加了来自几个不同来源的大量注释。在过滤器链中使用这些注释，我们已经生成了高质量和高可靠性的配体结构的“高质量”子集，根据应用的严格程度，它们的出现次数从大约1000次到大约5000次不等。我们对这些高质量配体的构象能进行了高水平的量子化学计算。在第一轮中，真空能量计算部分在我们自己的Linux集群上运行，部分在美国国立卫生研究院CIT的“http://biowulf.nih.gov/”Biowulf集群上运行。在这个计算量巨大的项目中，同时使用了多达一千个cpu，单个作业占用的cpu时间从几个小时到几周不等。我们获得了大约360次成功完井的结果。这些结果清楚地表明，这些量子化学计算完全证实了高构象能的可能性。他们在2008年10月13日至17日在费城Bryn Mawr举行的“http://echeminfo.com/COMTY_conferencesprog08”eCheminfo 2008互动会议上发表，该会议的主题是“http://echeminfo.com/comty_confprog08pdbligands”PDB配体：分析它们的结构和结合数据，由Marc Nicklaus主持。在这次会议上关于这些问题的讨论的结果是，一个由“关心的科学家”组成的国际小组聚集在一起，称为配体质量工作组，它将尝试通过非正式和更结构化的合作，以及通过信息的自由流动，尝试至少在这种情况下发出更集中的光，如果不能以各种方式改善它。为了探索水环境对配体构象能的可能影响——毕竟，真空并不是药物分子通常运作的地方——我们进行了第二轮量子化学计算，采用了高斯03中的SCI-PCM溶剂模型。这些运行甚至比真空计算对计算机资源的要求更高。为了分析能量不确定性作为位置不确定性的函数，而位置不确定性又是晶体学分辨率的函数，我们在分子力学力场水平上以晶体结构构象为中心进行了具有分辨率相关扭转分布的构象采样。所有这些结果目前正在被合并并撰写在一系列关于这一主题的主要出版物中。

The conformational changes of both partners of a ligand-protein complex, the small-molecule ligand its the protein binding site (in many cases the catalytically active site of an enzyme) are a central aspect many drug actions, as well as a crucial challenge in computational approaches to drug design. In one of the earliest publications in the this field, we "http://www.ncbi.nlm.nih.gov/sites/entrez?Db=pubmed and Cmd=ShowDetailView and TermToSearch=8581425 and ordinalpos=47 and itool=EntrezSystem2.PEntrez.Pubmed.Pubmed_ResultsPanel.Pubmed_RVDocSum" showed for a small set of ligands occurring both in the Protein Data Bank (PDB) and the Cambridge Structural Database (CSD) that flexible compounds are not usually bound to a protein in their global vacuum energy conformation, and oftentime not even in any local vacuum energy conformation. While this study used the largest set of data and best methodology available at that time, both the number of structures in either experimental database and the software and hardware resource available have since grown exponentially. We are thus revisiting this important topic with an analysis of orders of magnitudes more structures, and computations performed at a high level of computational quantum-chemical theory. Among other milestones achieved so far in this project, we have extracted all occurrences of small-molecule ligands recently made available in PDB's "http://ligand-expo.rcsb.org/" LigandExpo. As of May 2008, this is a set of over 350,000 distinct sets of 3D coordinates. We have added extensive annotation coming from several different sources. Using these annotations in a chain of filters, we have generated "high-quality" subsets of ligand structures of high quality and reliability numbering from just about one thousand to about 5,000 occurrences depending on the stringency applied. We have conducted high-level quantum-chemical calculations of conformational energies for these high-quality ligand sets. In the first round, vacuum energy calculations were run partly on our own Linux cluster, partly on the "http://biowulf.nih.gov/" Biowulf cluster of the CIT, NIH. Up to a thousand CPUs were used simultaneously in this computationally massive project, with individual jobs taking from a few hours to several weeks of CPU-time. We obtained results from about 360 runs that completed successfully. These results clearly showed that the possibility for high conformational energies are fully confirmed by these quantum-chemical calculations. They were presented at the "http://echeminfo.com/COMTY_conferencesprog08" eCheminfo 2008 InterAction Meeting at Bryn Mawr, Philadelphia (13-17 October 2008) in the session on "http://echeminfo.com/comty_confprog08pdbligands" PDB Ligands: Analysing their Structure and Binding Data, chaired by Marc Nicklaus. As a result of the discussions about these issues at this session, an international group of "concerned scientists" has come together, called the Ligand Quality Working Group, which will attempt, in collaborations both informal and more structured, and by free flow of information, to attempt to at least shine a more focused light on this situation, if not improve it in various ways. To explore the possible influence of aqueous environment on ligand conformational energies - after all, vacuum is not really where drug molecules typically operate - a second round of quantum chemical calculations was conducted, employing the SCI-PCM solvent model in Gaussian 03. These runs were even more demanding in terms of computer resources than the vacuum calculations. To analyze the energetic uncertainty as a function of the positional uncertainty, which in turn is a function of the crystallographic resolution, we conducted sampling of conformations with a resolution-dependent torsion distribution centered around the crystal structure conformation at the molecular mechanics force field level. All these results are currently being combined and written up in a series of major publications on this topic.