Multi-Target Gene Therapy for Osteoarthritis: A Computational and Structural Analysis Framework

Multi-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Multi-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation Po-Sung(Sinclair)Huang This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-8774255/v3 This work is licensed under a CC BY 4.0 License Status: Posted Version 3 posted You are reading this latest preprint version Show more versions Abstract Background Osteoarthritis (OA) remains therapeutically intractable despite advances in regenerative medicine. Single-target interventions—including mesenchymal stromal cells, platelet-rich plasma, and gene therapies—consistently yield transient symptomatic relief without durable structural modification. The recent Phase 2 failure of GLPG1972, a potent ADAMTS-5 inhibitor (ROCCELLA trial, n = 932), exemplifies this persistent challenge. Objective This study proposes a multi-target gene therapy strategy integrating computational structural analysis, network perturbation modeling, and clinical failure analysis to address the dual-axis nature of OA pathogenesis. Methods This computational study integrated structural analysis of therapeutic proteins (IL-1Ra, SOX9, IGF-1) using PyMOL, molecular docking (AutoDock Vina) to evaluate ADAMTS-5 inhibitor binding, Monte Carlo network perturbation analysis (n = 1,000 iterations) quantifying multi-target synergy via an ECM Recovery Score, Reynolds algorithm shRNA design with BLAST off-target validation, and human-canine sequence homology assessment across all therapeutic targets. Results Structural analysis confirms therapeutic transgenes preserve critical functional domains. Network perturbation demonstrates significantly higher ECM Recovery Scores for multi-target intervention (76.2 ± 8.3) versus single-target control (43.6 ± 12.1; p < 0.001). Analysis of GLPG1972 failure reveals an “Exosite Bypass” mechanism supporting complete enzyme elimination via shRNA rather than catalytic inhibition. Validated shRNA sequences for ADAMTS-5 and MMP-13 achieved Reynolds scores of 8–9/9. Human-canine homology analysis (mean 90.5% identity) supports canine models for translational studies. Conclusions Computational analysis, clinical failure mechanisms, and network modeling converge on the necessity of simultaneous multi-axis targeting. The proposed dual-vector AAV system addressing inflammation (IL-1Ra), anabolism (SOX9, IGF-1), and catabolism (ADAMTS-5/MMP-13 shRNA) provides a framework for experimental validation. This work is hypothesis-generating and requires experimental confirmation in cell-based and animal models. Computational Biology Osteoarthritis gene therapy molecular docking ADAMTS-5 multi-target therapy dual-axis model IL-1Ra SOX9 IGF-1 AAV vectors in silico modeling network perturbation Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Full Text Additional Declarations The authors declare no competing interests. Supplementary Files SUPPLEMENTARYMATERIALS.md Multi-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation replicationpackage.zip Readme_Script_Result VERSIONUPDATENOTES.md Verson_Notes figS3synergyanalysis.png figS2bootstrapci.png figS1expressionheatmap.png Cite Share Download PDF Status: Posted Version 3 posted You are reading this latest preprint version Show more versions Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-8774255","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":589216970,"identity":"ebaf90e5-f891-4e9a-868d-e0ad14bda9b3","order_by":0,"name":"Po-Sung(Sinclair)Huang","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA/klEQVRIiWNgGAWjYDACZiB+wCDBYADmVYBEmBsIa0mAazkDEmEkoIUBrIUBooWxDUzi1yLfznvwQ0KNRb65RPKzh1/n1UbztwO1/KjYhlOLwWG+ZImEYxKWO2ekmRvLbjueO+MwYwNjz5nbuLUw85gxJLBJGBjcSDCTltx2LLcBqIWZsQ23FvlmkJZ/IC3p36Ql5xzLnU9IC8NhoJbENpCWHDPJjw01uRsIaTE4zGMskdgH1HLmTZk0w7EDuRuBWg7i84t8/xnDDx++1RkYHE/fJvmjpi533vnDBx/8qMDjMGTAzMNwGMw4QJx6IGD8wVBHtOJRMApGwSgYOQAAuyhXkRLtQ0oAAAAASUVORK5CYII=","orcid":"","institution":"","correspondingAuthor":true,"prefix":"","firstName":"","middleName":"","lastName":"Po-Sung(Sinclair)Huang","suffix":""}],"badges":[],"createdAt":"2026-02-03 10:26:50","currentVersionCode":3,"declarations":{"humanSubjects":false,"vertebrateSubjects":false,"conflictsOfInterestStatement":false,"humanSubjectEthicalGuidelines":false,"humanSubjectConsent":false,"humanSubjectClinicalTrial":false,"humanSubjectCaseReport":false,"vertebrateSubjectEthicalGuidelines":false},"doi":"10.21203/rs.3.rs-8774255/v3","doiUrl":"https://doi.org/10.21203/rs.3.rs-8774255/v3","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":106117742,"identity":"06272e1e-b502-4cac-b883-378d057ac186","added_by":"auto","created_at":"2026-04-03 16:54:58","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":458094,"visible":true,"origin":"","legend":"\u003cp\u003eThe Dual-Axis Model of OA Pathogenesis. Schematic illustration of two mechanistically distinct but coupled pathological axes in osteoarthritis. Axis I (Inflammatory Signaling, red) includes IL-1β, TNF-α, NF-κB, COX-2, iNOS, PGE₂, and NO, which are pharmacologically suppressible yet insufficient for structural repair. Axis II (Mechano-Structural Degradation, blue) encompasses ADAMTS-5, MMPs, aggrecan loss, collagen degradation, and SOX9 suppression, representing a recalcitrant matrix-destruction axis. Bidirectional arrows indicate that inflammation drives matrix breakdown while DAMPs from degraded cartilage sustain inflammation. Single-target failures occupy the dashed box region (Axis I OR Axis II alone), demonstrating why interventions addressing a single axis cannot break the pathological cycle.\u003c/p\u003e","description":"","filename":"figure5homology.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/055404a0eeca7af696b484d1.png"},{"id":106117803,"identity":"e4e230bf-01ca-403e-9b27-f348a92a6da8","added_by":"auto","created_at":"2026-04-03 16:55:09","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":404042,"visible":true,"origin":"","legend":"\u003cp\u003eNetwork Perturbation Analysis (n = 1,000 Monte Carlo Iterations). (A) ECM Recovery Score distributions for control, IL-1Ra only, SOX9 only, and multi-target interventions. The multi-target strategy yields the highest median score (76.2) with narrower variance; ***p \u0026lt; 0.001 versus all other groups by Welch’s t-test with Bonferroni correction. (B) Component-level changes in key inflammatory and structural nodes. (C) Synergy analysis comparing the expected additive effect (67.4) to the observed multi-target effect (76.2); synergy index 1.13. Expected additive was calculated as: IL-1Ra only + SOX9 only − control (52.3 + 58.7 − 43.6 = 67.4), indicating that combined targeting produces more-than-additive benefit.\u003c/p\u003e","description":"","filename":"figure4shrnadesign.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/8f564249434ab07dd6236e28.png"},{"id":106117740,"identity":"a255d8de-86ea-4b8c-b016-1a497857e4cd","added_by":"auto","created_at":"2026-04-03 16:54:57","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":453296,"visible":true,"origin":"","legend":"\u003cp\u003eGLPG1972 Failure Mechanism – The Exosite Bypass Hypothesis. (A) GLPG1972 binding to the ADAMTS-5 catalytic site blocks zinc coordination (IC₅₀ = 19 nM) but leaves the exosite domain unoccupied. (B) Aggrecan substrate is warehoused at the exosite, pre-positioned for cleavage independent of catalytic-site occupancy. (C) Upon inhibitor dissociation, pre-bound substrate is immediately cleaved. (D) shRNA knockdown eliminates ADAMTS-5 protein, abolishing both catalytic activity and exosite-mediated substrate positioning, explaining the Phase 2 failure of GLPG1972 and supporting complete enzyme elimination as a therapeutic strategy.\u003c/p\u003e","description":"","filename":"figure3glpg1972failure.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/c895c39c9a7b0fbff0be3afa.png"},{"id":106117812,"identity":"f7f53f58-cbac-4a06-9405-189bbcac7113","added_by":"auto","created_at":"2026-04-03 16:55:13","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":247184,"visible":true,"origin":"","legend":"\u003cp\u003eshRNA Design and Validation for ADAMTS-5 and MMP-13. (A) Reynolds algorithm scores for candidate shRNAs; all lead sequences achieve 8–9/9, exceeding the design threshold. (B) Target positions of ADAMTS-5 shRNAs along NM_007038.5 at positions 892, 1456, and 2103. (C) Off-target BLAST analysis showing minimal off-target potential for all selected sequences. (D) Summary of lead candidate sequences, positions, Reynolds scores, and off-target profiles.\u003c/p\u003e","description":"","filename":"figure2networkperturbation.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/3aed959c40ff1e0357bdd9a9.png"},{"id":106117805,"identity":"88005caf-9343-4c02-a9db-83c41348a69c","added_by":"auto","created_at":"2026-04-03 16:55:10","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":250151,"visible":true,"origin":"","legend":"\u003cp\u003eHuman–Canine Homology Analysis of Therapeutic Targets. (A) Overall amino acid sequence identity: SOX9 (98.2%) to IL-1Ra (77.4%), mean 90.5%. (B) Functional domain conservation, with critical regions (SOX9 HMG box, IGF-1 mature peptide, catalytic domains) showing ≥96% identity. (C) Conservation summary supporting functional equivalence between human and canine targets. (D) Advantages of the canine model including spontaneous OA, human-like joint anatomy, optimal scAAV2.5 transduction, and established trial infrastructure.\u003c/p\u003e","description":"","filename":"figure1dualaxismodel.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/e5032248becd50b692e45581.png"},{"id":106402473,"identity":"2154c37d-a2d9-4b7d-aef9-afe643bf25d5","added_by":"auto","created_at":"2026-04-08 09:12:06","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2183165,"visible":true,"origin":"","legend":"","description":"","filename":"MultiTargetGene328RSPreprintV3.pdf","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3_covered_f446c02e-c279-492a-8312-ceba951b67d7.pdf"},{"id":106117688,"identity":"78fc9cc9-c971-4daf-a8fc-cb417c52ea23","added_by":"auto","created_at":"2026-04-03 16:54:45","extension":"md","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":9020,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eMulti-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation\u003c/strong\u003e\u003c/p\u003e","description":"","filename":"SUPPLEMENTARYMATERIALS.md","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/56320306cfd93260f75294cf.md"},{"id":106117715,"identity":"1ddffe24-d602-4629-bf13-bfe38bb9fcaf","added_by":"auto","created_at":"2026-04-03 16:54:52","extension":"zip","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":2091466,"visible":true,"origin":"","legend":"\u003cp\u003eReadme_Script_Result\u003c/p\u003e","description":"","filename":"replicationpackage.zip","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/1ce857294dea6cd7de59f2b3.zip"},{"id":106117741,"identity":"110bb888-4224-4d68-92fc-41318f844a5d","added_by":"auto","created_at":"2026-04-03 16:54:57","extension":"md","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":3554,"visible":true,"origin":"","legend":"\u003cp\u003eVerson_Notes\u003c/p\u003e","description":"","filename":"VERSIONUPDATENOTES.md","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/b26a6ba88e8d906114f3d7c3.md"},{"id":106117729,"identity":"6b4f4115-3046-4d97-ba70-a70abfb3c7e4","added_by":"auto","created_at":"2026-04-03 16:54:55","extension":"png","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":138917,"visible":true,"origin":"","legend":"","description":"","filename":"figS3synergyanalysis.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/3ef2f4ee052f3a4eb13d0350.png"},{"id":106117801,"identity":"e3f837f8-e72a-40ff-9d72-c7b2dd3c0b67","added_by":"auto","created_at":"2026-04-03 16:55:09","extension":"png","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":129096,"visible":true,"origin":"","legend":"","description":"","filename":"figS2bootstrapci.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/96670affc328811d837e69ce.png"},{"id":106117739,"identity":"571e0f98-26ea-4d72-a068-591992786495","added_by":"auto","created_at":"2026-04-03 16:54:57","extension":"png","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":271659,"visible":true,"origin":"","legend":"","description":"","filename":"figS1expressionheatmap.png","url":"https://assets-eu.researchsquare.com/files/rs-8774255/v3/99bf7c2df13eae762ab10687.png"}],"financialInterests":"The authors declare no competing interests.","formattedTitle":"Multi-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation","fulltext":[],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":false,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":true,"highlight":"","institution":"HEC Liege university","isAcceptedByJournal":false,"isAuthorSuppliedPdf":true,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":true,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Osteoarthritis, gene therapy, molecular docking, ADAMTS-5, multi-target therapy, dual-axis model, IL-1Ra, SOX9, IGF-1, AAV vectors, in silico modeling, network perturbation","lastPublishedDoi":"10.21203/rs.3.rs-8774255/v3","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-8774255/v3","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eOsteoarthritis (OA) remains therapeutically intractable despite advances in regenerative medicine. Single-target interventions\u0026mdash;including mesenchymal stromal cells, platelet-rich plasma, and gene therapies\u0026mdash;consistently yield transient symptomatic relief without durable structural modification. The recent Phase 2 failure of GLPG1972, a potent ADAMTS-5 inhibitor (ROCCELLA trial, n\u0026thinsp;=\u0026thinsp;932), exemplifies this persistent challenge.\u003c/p\u003e\u003ch2\u003eObjective\u003c/h2\u003e \u003cp\u003eThis study proposes a multi-target gene therapy strategy integrating computational structural analysis, network perturbation modeling, and clinical failure analysis to address the dual-axis nature of OA pathogenesis.\u003c/p\u003e\u003ch2\u003eMethods\u003c/h2\u003e \u003cp\u003eThis computational study integrated structural analysis of therapeutic proteins (IL-1Ra, SOX9, IGF-1) using PyMOL, molecular docking (AutoDock Vina) to evaluate ADAMTS-5 inhibitor binding, Monte Carlo network perturbation analysis (n\u0026thinsp;=\u0026thinsp;1,000 iterations) quantifying multi-target synergy via an ECM Recovery Score, Reynolds algorithm shRNA design with BLAST off-target validation, and human-canine sequence homology assessment across all therapeutic targets.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eStructural analysis confirms therapeutic transgenes preserve critical functional domains. Network perturbation demonstrates significantly higher ECM Recovery Scores for multi-target intervention (76.2\u0026thinsp;\u0026plusmn;\u0026thinsp;8.3) versus single-target control (43.6\u0026thinsp;\u0026plusmn;\u0026thinsp;12.1; p\u0026thinsp;\u0026lt;\u0026thinsp;0.001). Analysis of GLPG1972 failure reveals an \u0026ldquo;Exosite Bypass\u0026rdquo; mechanism supporting complete enzyme elimination via shRNA rather than catalytic inhibition. Validated shRNA sequences for ADAMTS-5 and MMP-13 achieved Reynolds scores of 8\u0026ndash;9/9. Human-canine homology analysis (mean 90.5% identity) supports canine models for translational studies.\u003c/p\u003e\u003ch2\u003eConclusions\u003c/h2\u003e \u003cp\u003eComputational analysis, clinical failure mechanisms, and network modeling converge on the necessity of simultaneous multi-axis targeting. The proposed dual-vector AAV system addressing inflammation (IL-1Ra), anabolism (SOX9, IGF-1), and catabolism (ADAMTS-5/MMP-13 shRNA) provides a framework for experimental validation. This work is hypothesis-generating and requires experimental confirmation in cell-based and animal models.\u003c/p\u003e","manuscriptTitle":"Multi-Target Gene Therapy for Osteoarthritis: Dual-Axis Modeling and In Silico Validation","msid":"","msnumber":"","nonDraftVersions":[{"code":3,"date":"2026-04-03 16:53:37","doi":"10.21203/rs.3.rs-8774255/v3","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}},{"code":2,"date":"2026-02-10 21:17:53","doi":"10.21203/rs.3.rs-8774255/v2","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}},{"code":1,"date":"2026-02-04 09:11:29","doi":"10.21203/rs.3.rs-8774255/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"0db17a57-92e3-4c6b-813a-c350884bacf0","owner":[],"postedDate":"April 3rd, 2026","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":62694276,"name":"Computational Biology"}],"tags":[],"updatedAt":"2026-02-04T09:11:29+00:00","versionOfRecord":[],"versionCreatedAt":"2026-04-03 16:53:37","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v3","identity":"rs-8774255","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-8774255","identity":"rs-8774255","version":["v3"]},"buildId":"XKTyCvWXoU3ODBz1xrDgd","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}