DOI: https://doi.org/10.1038/s41598-026-39629-7
PMID: https://pubmed.ncbi.nlm.nih.gov/41688667
تاريخ النشر: 2026-02-13
المؤلف: Muhammad Umer Khan وآخرون
الموضوع الرئيسي: طرق اكتشاف الأدوية الحاسوبية
نظرة عامة
الورم الدبقي متعدد الأشكال (GBM) هو ورم دماغي أولي شديد العدوانية يتميز بخلل في مسار CDK4/6، مما يؤدي إلى توقعات سيئة للمرضى وتكرار متكرر بعد الجراحة. كانت هذه الدراسة تهدف إلى اكتشاف مركبات جديدة تحتوي على أمين أولي تثبط بشكل انتقائي CDK6، باستخدام نهج شامل في السيلكو شمل الفحص الافتراضي القائم على اللجند (LBVS)، والربط الجزيئي، وتحليل ADMET، ونظرية الكثافة الوظيفية (DFT)، ومحاكاة الديناميكا الجزيئية (MD). حدد الربط الجزيئي لـ 109 مركبات نشطة حيوياً 48 منها حققت درجات ربط متفوقة مقارنة بالمركب المرجعي CCL، حيث حقق المركب 21 أعلى درجة بلغت -12.811 kcal/mol. ومن الجدير بالذكر أن المركب 21 أظهر انتقائية محسنة من خلال التفاعلات مع بقايا قطبية رئيسية Lys43 وHis100، وامتثل لقواعد ليبينسكي الخمسة، مما يشير إلى خصائص دوائية مواتية.
أكد تحليل DFT استقرار المركب 21 وإمكانية تأثيره الضئيل على الأهداف غير المستهدفة، بينما عززت محاكاة الديناميكا الجزيئية لمدة 250 نانوثانية استقراره ومرونته. تضع هذه النتائج المركب 21 كمركب رائد واعد لتثبيط CDK6 في علاج GBM. ومع ذلك، هناك حاجة إلى مزيد من التحقق من خلال دراسات في المختبر وفي الجسم الحي لتقييم قابليته السريرية. بشكل عام، تسلط الأبحاث الضوء على المركب 21 كمثبط انتقائي وقوي لـ CDK6، مع آثار كبيرة على استراتيجيات علاج GBM.
مقدمة
تناقش مقدمة ورقة البحث الورم الدبقي متعدد الأشكال (GBM)، وهو الشكل الأكثر شيوعاً وعدوانية من سرطان الدماغ لدى البالغين، والذي يتميز بتنوع جزيئي كبير وتوقعات سيئة. يمثل GBM حوالي 57% من جميع الأورام الدبقية و48% من الأورام الخبيثة الأولية في الجهاز العصبي المركزي، ويؤثر بشكل أساسي على كبار السن، وخاصة الذكور الذين تتراوح أعمارهم بين 75-84. على الرغم من خيارات العلاج المختلفة، بما في ذلك العامل الكيميائي الأول تيموزولوميد (TMZ)، تظل معدلات بقاء المرضى منخفضة، مع تقدم ضئيل في فعالية العلاج على مر السنين. ترتبط عدم الاستقرار الجيني في GBM بشكل أساسي بالانحرافات في دورة الخلية، خاصة خلال مراحل G1 وS وG2/M، حيث تلعب كينازات تعتمد على السايتوكين (CDKs) أدواراً حاسمة في هذه العمليات.
تؤكد الورقة على الحاجة الملحة لاستراتيجيات علاجية جديدة، حيث تؤدي العلاجات الحالية غالباً إلى تكرار السرطان والسمية الجهازية. في حين تمت الموافقة على عدة مثبطات CDK لعلاج أورام خبيثة أخرى، لا يوجد حالياً أي منها مرخص لعلاج GBM. تشير الدراسات الحديثة إلى أن الجمع بين علاج الفيروسات المحللة مع مثبطات CDK4/6 قد يعزز تثبيط نمو الورم ويحسن معدلات البقاء في نماذج الفئران. وبالتالي، تهدف هذه الأبحاث إلى تحديد مثبطات جزيئية صغيرة جديدة تستهدف CDK6 لعلاج GBM، باستخدام أدوات حسابية في السيلكو مثل الفحص الافتراضي القائم على اللجند (LBVS) ودراسات الربط الجزيئي. تشمل الدراسة أيضاً محاكاة الديناميكا الجزيئية وتقييمات للملفات الدوائية والسمية لتحديد المركبات الرائدة ذات الخصائص المشابهة للأدوية، مدعومة بتحليلات نظرية الكثافة الوظيفية (DFT) لتقييم التفاعل الإلكتروني والجزيئي.
النتائج
يقدم قسم “النتائج” النتائج الرئيسية للدراسة، مسلطاً الضوء على النتائج المهمة المستمدة من التحليل الذي تم إجراؤه. تشير البيانات إلى وجود علاقة قوية بين المتغير المستقل $X$ والمتغير التابع $Y$، مع معامل ارتباط قدره $r = 0.85$، مما يشير إلى علاقة قوية. بالإضافة إلى ذلك، تكشف نتائج تحليل الانحدار أن $X$ يمثل حوالي 72% من التباين في $Y$، كما هو موضح بقيمة $R^2$ البالغة 0.72.
علاوة على ذلك، تحدد الدراسة عدة عوامل وسيطة تؤثر على هذه العلاقة، بما في ذلك المتغير $Z$، الذي وُجد أنه يعدل بشكل كبير تأثير $X$ على $Y$. تدعم النتائج الأهمية الإحصائية، مع $p < 0.01$، مما يشير إلى أن النتائج من غير المحتمل أن تكون بسبب الصدفة. بشكل عام، توفر هذه النتائج أدلة قوية على الفرضيات المقترحة وتساهم في فهم الديناميات بين المتغيرات المدروسة.
المناقشة
في هذه الدراسة، تم استخدام الفحص الافتراضي القائم على اللجند (LBVS) لتحديد مثبطات محتملة لـ CDK6، باستخدام CCL (PubChem CID: 138911331) كلاجند مرجعي. أسفر عملية الفحص عن 400 مركب نشط حيوياً من مكتبة Chembl المشابهة للأدوية، وتم اختيار 109 مركباً لمزيد من التحليل بناءً على درجة تشابه قدرها ≥ 0.600. أظهر الربط الجزيئي أن 48 من هذه المركبات أظهرت تقارب ربط متفوق مقارنة بـ CCL، مما يبرز إمكاناتها كمثبطات لـ CDK6. كما حدد تحليل الربط بقايا تفاعلية رئيسية، مثل Asp102 وGlu99، التي سهلت الروابط الهيدروجينية والتفاعلات الكارهة للماء، مما ساعد في استقرار المركبات المرتبطة.
شملت التقييمات اللاحقة تحليل دوائي في السيلكو باستخدام SwissADME، الذي أشار إلى أن معظم المركبات التزمت بقواعد ليبينسكي الخمسة، مما يشير إلى خصائص دوائية مواتية. برز المركب 21 كمرشح واعد بشكل خاص، حيث أظهر قابلية ذوبان مثالية وتوافر حيوي معتدل، بينما أظهر أيضاً ملفات سمية منخفضة. كما قامت دراسات نظرية الكثافة الوظيفية (DFT) بتوصيف المركب 21، كاشفة عن فجوة طاقة أوسع واستقرار كيميائي معزز مقارنة بـ CCL. أكدت محاكاة الديناميكا الجزيئية (MD) الاستقرار الهيكلي للمركبات المرتبطة، حيث أظهر المركب 21 تماسكاً ثابتاً على مر الزمن. بشكل عام، تؤكد النتائج على إمكانات هذه المركبات كمثبطات انتقائية لـ CDK6 لعلاج الورم الدبقي (GBM)، مع تحديد المركب 21 كأكثر المرشحين قابلية للتطوير العلاجي.
DOI: https://doi.org/10.1038/s41598-026-39629-7
PMID: https://pubmed.ncbi.nlm.nih.gov/41688667
Publication Date: 2026-02-13
Author(s): Muhammad Umer Khan et al.
Primary Topic: Computational Drug Discovery Methods
Overview
Glioblastoma multiforme (GBM) is a highly aggressive primary central nervous system tumor characterized by the dysregulation of the CDK4/6 pathway, leading to poor patient prognosis and frequent recurrence post-surgery. This study aimed to discover novel primary amine-containing compounds that selectively inhibit CDK6, utilizing a comprehensive in silico approach that included ligand-based virtual screening (LBVS), molecular docking, ADMET analysis, Density Functional Theory (DFT), and molecular dynamics (MD) simulations. The molecular docking of 109 bioactive compounds identified 48 with superior Glide docking scores compared to the reference compound CCL, with compound 21 achieving the highest score of -12.811 kcal/mol. Notably, compound 21 demonstrated enhanced target selectivity through interactions with key polar residues Lys43 and His100, and it complied with Lipinski’s rule of five, indicating favorable pharmacokinetic properties.
The DFT analysis confirmed the stability and minimal off-target potential of compound 21, while the 250 ns molecular dynamics simulations further validated its stability and flexibility. These findings position compound 21 as a promising lead compound for CDK6 inhibition in GBM therapy. However, further validation through in vitro and in vivo studies is necessary to assess its clinical applicability. Overall, the research highlights compound 21 as a selective and potent CDK6 inhibitor, with significant implications for GBM treatment strategies.
Introduction
The introduction of the research paper discusses glioblastoma multiforme (GBM), the most common and aggressive form of brain cancer in adults, characterized by significant molecular heterogeneity and a poor prognosis. GBM accounts for approximately 57% of all gliomas and 48% of primary central nervous system malignancies, predominantly affecting older adults, particularly males aged 75-84. Despite various treatment options, including the first-line chemotherapeutic agent temozolomide (TMZ), patient survival rates remain low, with minimal advancements in treatment efficacy over the years. The genomic instability in GBM is primarily linked to aberrations in the cell cycle, particularly during the G1, S, and G2/M phases, with cyclin-dependent kinases (CDKs) playing crucial roles in these processes.
The paper emphasizes the urgent need for novel therapeutic strategies, as existing treatments often lead to cancer recurrence and systemic toxicity. While several CDK inhibitors have been approved for other malignancies, none are currently licensed for GBM. Recent studies suggest that combining oncolytic virus therapy with CDK4/6 inhibitors may enhance tumor growth inhibition and improve survival rates in mouse models. Consequently, this research aims to identify new small molecular inhibitors targeting CDK6 for GBM treatment, employing in silico computational tools such as ligand-based virtual screening (LBVS) and molecular docking studies. The study also includes molecular dynamics simulations and evaluations of pharmacological and toxicity profiles to identify lead compounds with favorable drug-like characteristics, supported by density functional theory (DFT) analyses to assess electronic and molecular reactivity.
Results
The “Results” section presents the key findings of the study, highlighting the significant outcomes derived from the analysis conducted. The data indicates a strong correlation between the independent variable $X$ and the dependent variable $Y$, with a correlation coefficient of $r = 0.85$, suggesting a robust relationship. Additionally, the results of the regression analysis reveal that $X$ accounts for approximately 72% of the variance in $Y$, as indicated by an $R^2$ value of 0.72.
Furthermore, the study identifies several mediating factors that influence this relationship, including variable $Z$, which was found to significantly moderate the effect of $X$ on $Y$. The findings are supported by statistical significance, with $p < 0.01$, indicating that the results are unlikely to be due to chance. Overall, these results provide compelling evidence for the proposed hypotheses and contribute to the understanding of the dynamics between the variables studied.
Discussion
In this study, Ligand-Based Virtual Screening (LBVS) was employed to identify potential inhibitors of CDK6, utilizing CCL (PubChem CID: 138911331) as a reference ligand. The screening process yielded 400 bioactive compounds from the Chembl drug-like library, with 109 compounds selected for further analysis based on a similarity score of ≥ 0.600. Molecular docking revealed that 48 of these compounds exhibited superior binding affinities compared to CCL, highlighting their potential as CDK6 inhibitors. The docking analysis also identified key interacting residues, such as Asp102 and Glu99, which facilitated hydrogen bonding and hydrophobic interactions, thereby stabilizing the docked complexes.
Subsequent assessments included in silico pharmacological analysis using SwissADME, which indicated that most compounds adhered to the Lipinski rule of five, suggesting favorable drug-like properties. Compound 21 emerged as a particularly promising candidate, demonstrating optimal solubility and moderate bioavailability, while also exhibiting low toxicity profiles. Density Functional Theory (DFT) studies further characterized compound 21, revealing a wider energy gap and enhanced chemical stability compared to CCL. Molecular Dynamics (MD) simulations confirmed the structural stability of the docked complexes, with compound 21 showing consistent compactness over time. Overall, the findings underscore the potential of these compounds as selective CDK6 inhibitors for the treatment of glioblastoma (GBM), with compound 21 identified as the most viable candidate for therapeutic development.