RD'nin aktif bileşenleri
TCMSP veri tabanından "Rhizoma Drynariae" anahtar kelimesi kullanılarak toplam 71 ham kimyasal bileşen elde edilmiştir. İkincil tarama ADME kriterlerine (OB ≥ 30%, DL ≥ 0.18) göre gerçekleştirilmiş ve sonuç olarak uygun oral absorbsiyona ve ilaçlaştırılabilirlik özelliğine sahip 18 biyoaktif bileşen belirlenmiştir. Tablo 1, bu 18 aktif bileşene ait ayrıntılı temel bilgileri özetlemektedir. Bu 18 bileşik, RD'nin RCI'ya karşı etkisindeki çok bileşenli materyal temelini oluşturmuştur.
RD'nin ilgili gen hedefleri
RD aktif bileşenleriyle eşleşen hedef genler için TCMSP, ETCM ve SymMap veri tabanları taranmış ve toplamda 327 orijinal hedef girişi toplanmıştır. Verilerin entegre edilmesi ve tekrarlanan gen sembollerinin çıkarılmasının ardından, sonraki kesişim analizi için RD'nin 213 adet tekrarsız potansiyel hedef geni elde edilmiştir.
RCI için ilgili gen hedefleri
Genecards ve OMIM veri tabanları "rotator cuff injury" anahtar kelimesi kullanılarak taranmış ve RCI ile ilgili 989 ham gen girişi elde edilmiştir. Veri entegrasyonu ve mükerrer kayıtların temizlenmesinin ardından, RCI progresyonu ile yakından ilişkili 870 benzersiz insan patojenik geni belirlenmiştir.
PPI ağının oluşturulması
213 RD hedef geni ile 870 RCI patojenik geni arasında yapılan kesişim analizi, RCI tedavisi için RD'nin potansiyel terapötik hedefleri olarak tanımlanan 71 örtüşen gen belirlemiştir. Kesişim ilişkisini gösteren Venn diyagramı Venny 2.0 yazılımı ile oluşturulmuştur (Şekil 2). 71 örtüşen hedef, RCI'yı etkileyen RD'nin PPI ağını kurmak için STRING'e aktarılmıştır (Şekil 3). Tüm PPI verileri Cytoscape'e aktarılmış ve cytoHubba ile en yüksek MCC skorlarına sahip ilk 10 hub gen taranmıştır: TNF, IL6, VEGFA, PTGS2, CASP3, AKT1, CCL2, PPARG, MAPK3, HMOX1. Bunların MCC topolojik bilgileri Tablo 2'de özetlenmiş ve Şekil 4'te görselleştirilmiştir. Kesişen 71 gen, RD'nin RCI'ya karşı terapötik etkisini aracılarla sağlayan temel hedefleri temsil etmektedir.
GO ve KEGG zenginleştirme analizi
RCI tedavisinde RD aday hedeflerinin moleküler mekanizmasını araştırmak için, çakışan terapötik hedefler üzerinde clusterProfiler paketi kullanılarak GO ve KEGG zenginleştirme analizleri gerçekleştirilmiştir. BP, CC ve MF için ayrı GO zenginleştirme istatistikleri yürütülmüştür. Şekil 5A–B, 6A–B ve 7A–B; BP, CC ve MF'nin en anlamlı zenginleşmiş ilk 10 girdisini göstermekte ve her şeklin B alt grafiği her bir fonksiyonel terimin gen dağılımını sergilemektedir. RD'nin RCI üzerindeki müdahalesinde yer alan temel yolakları ortaya çıkarmak için ayrıca KEGG zenginleştirme analizi yapılmış ve zenginleşmiş ilk 30 yolak Şekil 8'de görselleştirilmiştir.
Moleküler kenetlenme simülasyonu
En yüksek 10 hub gen protein reseptörleri olarak belirlenmiş ve tüm 18 RD aktif bileşeni, AutoDock Vina ve Open Babel aracılığıyla yapılacak docking hesaplamaları için ligand olarak benimsenmiştir. Ligand-reseptör bağlanma konformasyonları PyMOL ile görselleştirilmiştir. Docking sonuçları, çoğu hub hedef proteinin RD aktif bileşenleri ile stabil bağlanma yapıları oluşturduğunu göstermiştir. Tablo 3 ve Şekil 9, optimal bağlanma afinitesine sahip en iyi 20 kompleksin tüm docking parametrelerini ve üç boyutlu bağlanma diyagramlarını kaydetmektedir.
Mevcut ağ farmakolojisi analizinde, RD'nin 18 biyoaktif bileşeni ve 213 tekrarlayan olmayan ilaçla ilişkili hedef geri çağrılmıştır. 870 RCI ilişkili genle kesiştirme yapıldıktan sonra, 71 ortak terapötik aday hedef elde edilmiştir. PPI ağı oluşturma ve cytoHubba-MCC taraması sonucunda on adet merkez gen belirlenmiştir: TNF, IL6, VEGFA, PTGS2, CASP3, AKT1, CCL2, PPARG, MAPK3 ve HMOX1. GO zenginleştirme analizi BP, CC ve MF'yi kapsarken, KEGG zenginleştirmesi birden fazla zenginleştirilmiş sinyal yolunu tanımlamıştır. Moleküler kenetleme simülasyonu, RD'nin ana aktif bileşenleri ile bu merkez hedef proteinler arasındaki uygun bağlanma afinitesini doğrulamıştır.
VERİERİLEBİLİRLİK:
Hedef listeleri, PPI ağı dosyaları, GO ve KEGG zenginleştirme sonuçları ve moleküler kenetlenme matrisi dahil olmak üzere bu çalışmada üretilen tüm ham ve işlenmiş veriler Science Data Bank'e kaydedilmiştir. Ham veri erişim bağlantısı https://www.scidb.cn/anonymous/alVOVlJu şeklindedir.

Şekil 1: Bu çalışmada uygulanan biyoinformatik işlem hattını özetleyen şematik iş akışı. Soldan başlanarak, aday biyoaktif bileşikler ve bunlara karşılık gelen RD hedefleri TCMSP, ETCM ve SymMap veri tabanlarından elde edilmiştir. Sağ tarafta, RCI ile ilişkili genler GeneCards ve OMIM'den alınmıştır. Bu iki veri setinden gelen örtüşen gen hedefleri bir Venn diyagramı kullanılarak belirlenmiş ve potansiyel terapötik hedefler olarak tanımlanmıştır. Ortak hedefler daha sonra iki ayrı kolda işlenmiştir: bir kol STRING ve Cytoscape kullanarak bir protein-protein etkileşim ağı oluşturmuş; diğer kol ise R paketi clusterProfiler kullanarak GO ve KEGG zenginleştirme analizleri gerçekleştirmiştir. PPI ağından çıkarılan hub genler; yapısal hazırlık ve simülasyonları PubChem, PDB, AutoDock Vina ve PyMOL kullanılarak tamamlanan moleküler kenetleme ile daha fazla doğrulanmıştır. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, TCMSP = The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, ETCM = The Encyclopedia of Traditional Chinese Medicine, OMIM = The Online Mendelian Inheritance in Man, GO = gen ontolojisi, KEGG = Kyoto encyclopedia of genes and genomes. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 2: RD ve RCI arasındaki ortak hedeflerin Venn diyagramı. Mavi daire, RCI ile ilişkili hedefleri temsil eder. Sarı daire, RD'nin potansiyel terapötik hedeflerini temsil eder. Kesişen zeytin yeşili bölge, RD ve RCI arasındaki ortak hedefleri belirtir. Her bölgenin içindeki sayılar ve yüzdeler, benzersiz veya ortak hedeflerin sayısını ve oranını ifade eder. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayınız.

Şekil 3: RD ve RCI'yi ilişkilendiren 71 örtüşen gen hedefinden oluşturulan PPI ağı. Ağdaki her daire, resmi gen sembolüyle etiketlenmiş bir hedef geni temsil eder. İki daireyi birbirine bağlayan çizgiler, halka açık biyoinformatik veri tabanlarından elde edilen protein-protein ilişkilerini göstermektedir. Düğüm boyutu ve renk doygunluğu, ağ içindeki bireysel genlerin topolojik önemini yansıtır. Daha büyük dairelere ve daha koyu yeşil tonlara sahip genler, daha yüksek topolojik skorlara ulaşmış ve daha merkezi konumlarda yer almaktadır; daha küçük, soluk sarı daireler ise topolojik önemi daha düşük olan genleri belirtmektedir. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator cuff yaralanması, PPI = protein-protein etkileşimi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 4: cytoHubba eklentisi kullanılarak PPI ağından belirlenen ilk 10 hub gen. Dikdörtgen kutular, her biri resmi gen sembolüyle etiketlenmiş hub genleri temsil eder. Kutuları birbirine bağlayan siyah çizgiler, bu temel hedefler arasındaki protein etkileşimlerini yansıtmaktadır. Parlak kırmızıdan turuncuya ve sarıya doğru giden renk gradyanı, hub genlerin cytoHubba eklentisi tarafından hesaplanan MCC skorlarına göre sıralamasını görselleştirmektedir; parlak kırmızı kutular en yüksek MCC skorlarına ve en büyük ağ önemine sahip genleri, sarı kutular ise nispeten daha düşük skorlu genleri belirtir. Kısaltmalar: PPI = protein-protein etkileşimi, MCC = Maksimal Clique Centrality.Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 5: RD'nin RCI'ya karşı 71 ortak hedef geninin BP zenginleştirme analizi. (A) BP için potansiyel terapötik hedeflerin GO zenginleştirme analizi. En anlamlı olarak zenginleşmiş ilk 10 GO biyolojik süreç terimini gösteren yatay bir çubuk grafiktir. Yatay X-ekseni, her bir GO terimine atanan hedef gen sayısına karşılık gelen gen sayısını temsil eder. Dikey Y-ekseni, tekil zenginleşmiş GO terimlerini listeler. Sağdaki sürekli renk çubuğu -log10(p-değeri) değerini kodlar; daha sıcak kırmızı tonlar daha yüksek -log10(p-değeri) ve daha güçlü istatistiksel anlamlılığı, daha soğuk mor ve mavi tonlar ise daha düşük anlamlılığı belirtir. (B) Potansiyel terapötik hedefler ve BP arasındaki ilişki. Panel B, zenginleşmiş GO terimleri ve bunlara karşılık gelen hedef genler arasındaki ikili ilişkileri gösteren bir akor diyagramıdır. Dış dairesel segmentler iki kategoriye ayrılmıştır: renkli bloklar farklı GO biyolojik süreç terimlerini temsil eder ve gri dikdörtgen işaretler tekil hedef gen sembollerini etiketler. Daire içindeki kavisli şeritler, GO terimlerini eşleşen hedef genlere bağlar ve her şerit bir geni katıldığı GO terimine eşler. Her şeridin rengi, dış halkadaki karşılık gelen GO teriminin tonundan miras alınmıştır. Alttaki lejand, her rengi ilgili GO terimiyle ilişkilendirir. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, BP = biyolojik süreç, GO = gen ontolojisi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 6: RCI tedavisinde RD'nin 71 aday hedefinden türetilen GO-CC terimleri. (A) CC için potansiyel terapötik hedeflerin GO zenginleştirme analizi. (B) Potansiyel terapötik hedefler ile CC arasındaki ilişki. Panel A, en anlamlı şekilde zenginleşmiş ilk 10 GO hücresel bileşen girdisini gösteren yatay bir çubuk grafiktir. Yatay X ekseni, her bir GO terimine atanan hedef gen sayısını belirler ve Counts olarak etiketlenmiştir. Zenginleştirilmiş GO hücresel bileşen tanımları dikey Y ekseni boyunca düzenlenmiştir. Sağdaki renk çubuğu, renk değişimini -log10(pvalue) değerleriyle eşleştirir; parlak kırmızı daha yüksek -log10(pvalue) değerlerini ve daha güçlü istatistiksel anlamlılığı, mavi tonlar ise daha düşük anlamlılık seviyelerini temsil eder. Panel B, zenginleştirilmiş GO hücresel bileşen terimleri ile bunlarla ilişkili hedef genler arasındaki bağlantıları görselleştiren bir akor diyagramı sunar. Dış halkadaki segmentler iki gruba ayrılmıştır: renkli bloklar farklı GO hücresel bileşen terimlerini, gri işaretler ise bireysel hedef gen sembollerini belirtir. Dairenin içindeki kavisli şeritler her geni eşleştiği GO terimine bağlar ve her şerit, dış halkadaki karşılık gelen GO teriminin rengini alır. Alttaki lejant, tüm GO hücresel bileşen kategorileri için renk kodlamasını tanımlar. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, CC = hücresel bileşen, GO = gen ontolojisi Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 7: RCI müdahalesi için RD'nin örtüşen 71 hedef geni üzerinde gerçekleştirilen MF zenginleştirme analizi. (A) MF için potansiyel terapötik hedeflerin GO zenginleştirme analizi. (B) Potansiyel terapötik hedefler ve MF arasındaki ilişki. Panel A, en anlamlı şekilde zenginleşmiş ilk 10 GO moleküler fonksiyon terimini sıralayan yatay bir çubuk grafiğidir. Yatay X eksenindeki değerler, her bir moleküler fonksiyon girdisine atanan toplam hedef gen sayısına karşılık gelen gen sayılarını temsil eder. Dikey Y ekseni, bireysel GO moleküler fonksiyon tanımlarını listeler. Sağdaki renk gradyan çubuğu, renk varyasyonunu -log10(pvalue) metriklerine dönüştürür; parlak kırmızı daha yüksek -log10(pvalue) değerlerine ve daha güçlü istatistiksel anlamlılığa karşılık gelirken, mavi tonlar daha zayıf anlamlılığı gösterir. Panel B, zenginleşmiş GO moleküler fonksiyon kategorileri ile bunlarla ilişkili hedef genler arasındaki ilişkileri eşleyen bir kord diyagramıdır. Dış halkadaki segmentler iki etiket sınıfı içerir: renkli bloklar farklı GO moleküler fonksiyon terimlerini belirtir ve küçük gri dikdörtgen işaretler bireysel hedef gen sembollerini işaretler. Dış halkadan içeriye doğru uzanan kavisli şeritler, her geni eşleştiği fonksiyonel terime bağlar ve her şerit, karşılık gelen GO terimine atanan rengi alır. Alttaki lejant, grafikte gösterilen tüm GO moleküler fonksiyon kategorileri için renk kodlamasını tanımlar. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, MF = moleküler fonksiyon, GO = gen ontolojisi Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayınız.

Şekil 8: RD'nin RCI'ya karşı 71 ortak hedef geninin KEGG yolak zenginleştirme analizi. Bu nokta grafiği, en çok zenginleşen KEGG yolaklarını görselleştirmektedir. Dikey Y ekseni, zenginleşen münferit yolakların isimlerini listeler. Yatay X ekseni, bir yolakla eşleşen hedef genlerin o yolaktaki toplam gen sayısına oranı olarak tanımlanan GeneRatio değerini temsil eder. Her bir dairesel nokta, zenginleşmiş bir yolağı temsil eder. Nokta boyutu gen sayısına karşılık gelir; daha büyük daireler, yolağa daha fazla hedef genin atandığını gösterir. Nokta rengi p-değerini yansıtır: kırmızı daha düşük p-değerlerini ve daha güçlü istatistiksel anlamlılığı, koyu mavi ise daha yüksek p-değerlerini ve daha zayıf istatistiksel anlamlılığı gösterir. İki yan lejant, gen sayısı ve p-değeri renk eşlemesi ölçeklerini ayrı ayrı tanımlar. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, KEGG = Kyoto gen ve genom ansiklopedisi. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayın.

Şekil 9: En yüksek sıralamaya sahip docking kompleksleri. (A) (PTGS2, eriodictyol, -9.7 kcal/mol). (B) (PTGS2, luteolin, -9.6 kcal/mol). (C) (PTGS2, catechin, -9.5 kcal/mol). (D) (PTGS2, naringenin, -9.5 kcal/mol). (E) (PTGS2, Aureusidin, -9.4 kcal/mol). (F) (PTGS2, kaempferol, -9.2 kcal/mol). (G) (PTGS2, chroman, -9.1 kcal/mol). (H) (TNF, catechin, -9.1 kcal/mol), (I) (PTGS2, digallate, -9 kcal/mol). (J) (TNF, digallate, -9 kcal/mol), (K) (PPARG, chroman, -8.9 kcal/mol); (L) (PPARG, naringenin, -8.9 kcal/mol), (M) (TNF, Aureusidin, -8.9 kcal/mol); (N) (TNF, kaempferol, -8.8 kcal/mol); (O) (MAPK3, chroman, -8.7 kcal/mol); (P) (MAPK3, naringenin, -8.7 kcal/mol); (Q) (PPARG, Aureusidin, -8.7 kcal/mol); (R) (PPARG, kaempferol, -8.7 kcal/mol); (S) (PTGS2, 22_Stigmasten_3_one, -8.7 kcal/mol); (T) (MAPK3, luteolin, -8.6 kcal/mol). A'dan T'ye kadar olan paneller, ilgili protein hedeflerin bağlanma boşluğu içindeki münferit küçük moleküllü bileşiklerin hesaplanan bağlanma geometrisini ayrı ayrı göstermektedir. Siyan renkli çubuk modelleri, protein bağlanma boşluğunu çevreleyen amino asit kalıntılarını betimlemektedir. Küçük moleküllü ligandlar, farklı kimyasal yapıları ayırt etmek için macenta veya yeşil çubuklar şeklinde sunulmuştur. Sarı kesikli çizgiler, ligandlar ve cep kalıntıları arasındaki moleküller arası hidrojen bağlarını belirtmekte olup, sayısal notasyonlar hidrojen bağı mesafelerini angstrom (Å) cinsinden göstermektedir. Yarı saydam beyaz yüzey, her bir protein reseptörünün genel 3D topolojisini ana hatlarıyla belirtmektedir. Tüm bağlanma pozu hesaplamaları, AutoDock Vina yazılımı kullanılarak rijit reseptör moleküler docking yoluyla gerçekleştirilmiştir. Bu şeklin daha büyük bir versiyonunu görüntülemek için lütfen buraya tıklayınız.
| Mol Kimliği | Molekül Adı | oral biyoyararlanım(OB) (%) | ilaç benzerliği (DL) |
| MOL001040 | (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one | 42.36 | 0.21 |
| MOL001978 | Aureusidin | 53.42 | 0.24 |
| MOL002914 | Eriodyctiol (flavanon) | 41.35 | 0.24 |
| MOL000449 | Stigmasterol | 43.83 | 0.76 |
| MOL000358 | beta-sitosterol | 36.91 | 0.75 |
| MOL000422 | kaempferol | 41.88 | 0.24 |
| MOL004328 | naringenin | 59.29 | 0.21 |
| MOL000492 | (+)-kateşin | 54.83 | 0.24 |
| MOL005190 | eriodiktyol | 71.79 | 0.24 |
| MOL000569 | digallat | 61.85 | 0.26 |
| MOL000006 | luteolin | 36.16 | 0.25 |
| MOL009061 | 22-Stigmasten-3-on | 39.25 | 0.76 |
| MOL009063 | Siklolaudenol asetat | 41.66 | 0.79 |
| MOL009075 | sikloartenon | 40.57 | 0.79 |
| MOL009076 | Siklolaudenol | 39.05 | 0.79 |
| MOL009078 | davallioside A_qt | 62.65 | 0.51 |
| MOL009087 | marioside_qt | 70.79 | 0.19 |
| MOL009091 | ksantogalenol | 41.08 | 0.32 |
Tablo 1: RD'den türetilen on sekiz aday biyoaktif bileşenin farmakokinetik profilleri. Bileşikler TCMSP veritabanından alınmış ve iki ADME ile ilgili eşik değere göre filtrelenmiştir: OB ≥ 30% ve DL ≥ 0.18. Sütunlar; her bir taranmış bileşen için TCMSP erişim numarasını (Mol ID), bileşik adını, hesaplanan OB değerini (%) ve DL skorunu listelemektedir. Kısaltmalar: RD = Rhizoma Drynariae, RCI = rotator manşet yaralanması, ADME = absorbsiyon, dağılım, metabolizma, ekskresyon, TCMSP = Geleneksel Çin Tıbbı Sistemleri Farmakoloji Veritabanı ve Analiz Platformu, OB = oral biyoyararlanım, DL = ilaç benzerliği.
| Sıralama | Gen | İsim | Puan |
| 1 | TNF | tümör nekroz faktörü | 328291218366500 |
| 2 | IL6 | interlökin 6 | 328291215913592 |
| 3 | VEGFA | vasküler endotelyal büyüme faktörü A | 328291215896284 |
| 4 | PTGS2 | prostaglandin-endoperoksit sentaz 2 | 328288656326402 |
| 5 | CASP3 | Kaspaz 3 | 328288180323600 |
| 6 | AKT1 | AKT serin/treonin kinaz 1 | 328287745700516 |
| 7 | CCL2 | C-C motif kemokin ligand 2 | 328161521686320 |
| 8 | PPARG | peroksizom proliferasyonla aktive olan reseptör gama | 325302724852296 |
| 9 | MAPK3 | mitojenle aktive olan protein kinaz 3 | 308779166166218 |
| 10 | HMOX1 | hem oksijenaz 1 | 299171945318400 |
Tablo 2: PPI ağından belirlenen en önemli 10 hub gen. Genler, Cytoscape içindeki cytoHubba eklentisi tarafından hesaplanan MCC skorlarına göre sıralanmıştır. MCC, yüksek sayısal değerlerin ağ içindeki daha büyük topolojik önemi yansıttığı bir cytoHubba topolojik metriğidir. Sütunlar; sıralama sırasını, resmi gen sembolünü, tam gen adını ve her bir hub gen için karşılık gelen ham MCC skorunu göstermektedir. Kısaltmalar: PPI = protein-protein etkileşimi, MCC = Maksimal Clique Centrality.
| Reseptörler | Ligandlar | Serbest enerji (kcal/mol) | Şekil 9'daki karşılık gelen seri numaraları |
| PTGS2 | eriodictyol | -9.7 | A |
| PTGS2 | luteolin | -9.6 | B |
| PTGS2 | catechin | -9.5 | C |
| PTGS2 | naringenin | -9.5 | D |
| PTGS2 | Aureusidin | -9.4 | E |
| PTGS2 | kaempferol | -9.2 | F |
| PTGS2 | chroman | -9.1 | G |
| TNF | catechin | -9.1 | H |
| PTGS2 | digallate | -9 | I |
| TNF | digallate | -9 | J |
| PPARG | chroman | -8.9 | K |
| PPARG | naringenin | -8.9 | L |
| TNF | Aureusidin | -8.9 | M |
| TNF | kaempferol | -8.8 | N |
| MAPK3 | chroman | -8.7 | O |
| MAPK3 | naringenin | -8.7 | P |
| PPARG | Aureusidin | -8.7 | Q |
| PPARG | kaempferol | -8.7 | R |
| PTGS2 | 22_Stigmasten_3_one | -8.7 | S |
| MAPK3 | luteolin | -8.6 | T |
Tablo 3: Seçilen 20 reseptör-ligand kompleksi için moleküler kenetlenme bağlanma afinitesi sonuçları. Her satır reseptör proteinini, küçük molekül ligandı, hesaplanan bağlanma serbest enerjisini (kcal/mol) ve Şekil 9'daki ilgili alt grafik etiketini listelemektedir. Daha negatif serbest enerji değerleri, protein cebi ile ligand bileşiği arasındaki daha güçlü bağlanma stabilitesine karşılık gelir. Kısaltmalar: PTGS2 = prostaglandin-endoperoksit sentaz 2, TNF = tümör nekroz faktörü, PPARG = peroksizom proliferatif aktive reseptör gama, MAPK3 = mitojenle aktive olan protein kinaz 3.