Molecular Modeling of the Human P-glycoprotein Transporter Protein

像许多跨膜蛋白一样，通过x射线晶体学测定P-gp的结构已被证明是难以捉摸的，尽管付出了很多努力。这源于形成足够质量的晶体以维持蛋白质不同部分的天然物理化学环境的困难。因此，代替直接的实验确定，我们正在努力将所有可用的（通常是间接的）实验数据与基于物理化学的数学方法相结合，以产生一个或多个物理上真实的结构模型。幸运的是，超过三十年的研究已经提供了关于P-gp的丰富信息，从中我们可以获得结构信息。在广泛的范围内，已知该蛋白由两个同源结构域组成，每个结构域都具有六段跨膜成分和核苷酸结合成分。迄今为止，关于P-gp的结构信息最好的两个来源是同源细菌蛋白Sav1866的x射线晶体结构和人类P-gp的低分辨率冷冻电子显微照片。此外，细菌脂质翻转酶MsbA的x射线修正结构有望很快得到，它与P-gp的关系甚至比Sav1866更密切。因此，我们工作的一个主要重点是利用细菌蛋白的晶体结构作为模板开发人类P-gp的同源模型。虽然另一个小组最近发表了这样的模型，但它只是基于P-gp和Sav1866序列的简单比对，不幸的是，这对于跨膜片段并没有很好的定义。相反，我们的努力深入研究了密切相关的MDR蛋白家族和ABC转运蛋白超家族中的残基保存模式。这些信息有助于预测哪些残基暴露在膜的核心和头群层，哪些残基排列在孔中，哪些在两个跨膜结构域的界面上。我们目前正在开发同源家族和超家族的大序列比对。这一结果也将有助于确定相关突变的模式，这有助于识别蛋白质三维结构中近端的残基组。最后，我们将检查所得模型与所有实验数据的一致性，例如位点定向诱变的影响，自然发生的多态性和交联数据。如果基于Sav1866模板的模型无法解释实验结果，我们将寻找使其符合的替代构象。另一个重点是根据电子显微镜获得的密度图开发模型。这样做的好处是，结构数据直接来自人类P-gp，即目标蛋白，而不是来自细菌的同源物，后者可能在结构上不同。这包括人类P-gp的两个跨膜结构域的序列不同，因此在孔的近似双轴周围是不对称的（如显微照片中观察到的），而细菌同源物只包含一个结构域，因此在膜中形成完美对称的同型二聚体。在得到与跨膜电子密度相匹配的一般螺旋坐标后，我们几乎完成了计算机软件来确定氨基酸序列在蛋白质主链结构上的排列。这相当于一个典型的大计算规模的蛋白质线程程序。该程序构建氨基酸片段到跨膜螺旋的每种可能组合，并根据一系列物理化学和实验确定的标准对其进行评分。其中包括每种氨基酸残基的适当环境，保守残基的聚类以及实验交联距离限制。必须克服的一个困难是必须扫描相对大量的排列。通过开发聪明的算法，将所需的计算量减少到当今技术能力的实际范围内，克服了这一问题。随着模型的发展，我们将与细胞生物学实验室的其他小组密切合作，特别是Xia博士、Suresh Ambudkar博士和Michael Gottesman博士的小组，进行实验测试、x射线结构测定和抑制剂开发。

Like many transmembrane proteins, determination of the structure of P-gp by X-ray crystallography has proven elusive despite much effort. This stems from the difficulty of forming sufficient-quality crystals that maintain the native physiochemical environments for the different parts of the protein. Thus, in lieu of direct experimental determination, we are striving to integrate all available (often indirect) experimental data with physiochemically-based mathematical methods to produce one or more physically realistic models of the structure. Fortunately, over three decades of study has provided a wealth of information about P-gp from which we can gleam structural information. On a broad scale, it is known that the protein is composed of two homologous domains, each with a six-segment transmembrane component and a nucleotide-binding component. To date, the best two sources of structural information about P-gp are an X-ray crystal structure of the homologous bacterial protein Sav1866, and low-resolution cryo-electron micrographs of human P-gp. In addition, the corrected X-ray structure of the bacterial lipid flippase MsbA is expected soon, which is even more closely related to P-gp than is Sav1866. Thus, one major focus of our work is to develop a homology model of human P-gp using the crystal structures of the bacterial proteins as templates. While another group has recently published such a model, it was only based on a simple alignment of the P-gp and Sav1866 sequences, and unfortunately, this is not well defined for the transmembrane segments. Rather, our efforts go deeper into examining the patterns of residue conservation within the family of closely related MDR proteins and the superfamily of ABC transporters. This information helps predict which residues are exposed to the core and headgroup layers of the membrane, which residues line the pore, and which are at the interfaces of the two transmembrane domains. We are currently in the process of developing a grand sequence alignment of homologous families and the superfamily. The results of this will also enable the determination of patterns of correlated mutations, which help identify groups of residues that are proximal in the 3-dimensional structure of the protein. Finally, we will examine the resultant models for consistency with all the experimental data, such as the effects of site-directed mutagenesis, naturally occurring polymorphisms, and cross-linking data. Where the model based on the Sav1866 template fails to explain the experimental results, we will search for alternate conformations that bring it into compliance. The other major focus is to develop models from the density maps obtained from electron microscopy. This has the advantage that the structural data is directly from human P-gp, the target protein, and not from a bacterial homolog, which likely differs in structure. This includes the fact that the two transmembrane domains of human P-gp are different in sequence, and thus are asymmetrical around the approximate two-fold axis of the pore (as observed in the micrographs), while the bacterial homologs only contain one domain, and thus form perfectly symmetrical homodimers in the membrane. Having obtained coordinates of generic helices fitted to the transmembrane electron density, we have nearly completed computer software to determine the alignment of the amino acids sequences onto the backbone structure of the protein. This amounts to a typical Protein Threading procedure on a grand computational scale. The program builds every possible combination of amino acid segment to transmembrane helix and scores it according to a series of physiochemical and experimentally-determined criteria. These include the proper environment for each type of amino acid residue, the clustering of conserved residues, and experimental cross-linking distance constraints, One difficulty that had to be surmounted was the relative enormous number of permutations that had to be scanned. This was overcome by developing clever algorithms to reduce the required amount of computations to within the practical scope of present-day technological capabilities. As the models develop, we will be working closely with other groups within the Lab of Cell Biology, especially Dr.s Di Xia's, Suresh Ambudkar's, and Michael Gottesman's sections, to undertake experimental tests, X-ray structure determination, and inhibitor development.