Pincus Matthew R, Woo Jannie, Monaco Regina, Lubowsky Jack, Avitable Matthew J, Carty Robert P
Department of Pathology, SUNY Downstate Medical Center, 450 Clarkson Avenue, Brooklyn, NY, 11203, USA,
Protein J. 2014 Dec;33(6):565-74. doi: 10.1007/s10930-014-9589-3.
Using the chain-build-up method based on Empirical Conformational Energies of Peptides Program including solvation, we have computed, the low energy conformations of gonadotrpin-releasing hormone, GnRH, whose sequence is Pyro-Glu(PG)-His-Trp-Ser-Tyr-Gly-Leu-Arg-Pro-Gly-NH2. We have found 5,077 solvated conformations with conformational energies that were within 5 kcal/mole of that of the global minimum. These minima were found to occur in 802 distinct conformational classes of which 25 represented 70 % of the Boltzmann energy-weighted structures. Virtually all of these structures adopted bend conformations from Tyr 5-Leu 8, and 3,861 structures adopted bend conformations at residues 4-7. However, these structures differed significantly from one another, indicating that GnRH does not adopt a well-defined structure in aqueous solution consistent with the absence of a well-defined NMR structure of GnRH in water. A total of 300 of these structures were found to be superimposable on possible NMR structures for GnRH in DMSO with a combined statistical weight of 1.6 %. We found that Gly 6 adopts low energy "starred" states, e.g., C* and D*, that are energetically forbidden to L-amino acids but are low energy for D-amino acids, with a statistical weight of 43 %. This can explain why substitutions of L-amino acids for Gly 6 are known to inactivate GnRH while D-amino acid substitutions enhance its activity. Using these findings, in the accompanying manuscript, we compute the low energy conformations for the substituted GnRHs that enable inference of possible receptor-bound conformations.
我们使用基于肽构象能量经验程序(包括溶剂化)的链构建方法,计算了促性腺激素释放激素(GnRH)的低能构象,其序列为焦谷氨酸(PG)-组氨酸-色氨酸-丝氨酸-酪氨酸-甘氨酸-亮氨酸-精氨酸-脯氨酸-甘氨酸-NH2。我们发现了5077种溶剂化构象,其构象能量在全局最小值的5千卡/摩尔范围内。这些最小值出现在802个不同的构象类别中,其中25个代表了70%的玻尔兹曼能量加权结构。实际上所有这些结构在酪氨酸5-亮氨酸8处都采用了弯曲构象,3861个结构在残基4-7处采用了弯曲构象。然而,这些结构彼此之间存在显著差异,这表明GnRH在水溶液中没有采用明确的结构,这与在水中缺乏明确的GnRH核磁共振结构一致。总共发现其中300个结构可以叠加在GnRH在二甲基亚砜中的可能核磁共振结构上,综合统计权重为1.6%。我们发现甘氨酸6采用低能的“星号化”状态,例如C和D,这些状态对于L-氨基酸在能量上是禁止的,但对于D-氨基酸是低能的,统计权重为43%。这可以解释为什么已知用L-氨基酸取代甘氨酸6会使GnRH失活,而用D-氨基酸取代则会增强其活性。利用这些发现,在随附的手稿中,我们计算了取代GnRH的低能构象,从而能够推断出可能的受体结合构象。