SEARCH · PubMed Health

Results for “molecular docking”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36Linked to original sources

PubMed2001

Synthesis and biological evaluation of all A-ring stereoisomers of 5,6-trans-2-methyl-1,25-dihydroxyvitamin D(3) and their 20-epimers: possible binding modes of potent A-ring analogues to vitamin D receptor.

BACKGROUND: The secosteroid 1 alpha,25-dihydroxyvitamin D(3) (1) has a wide variety of biological activities, which makes it a promising therapeutic agent for the treatment of cancer, psoriasis and osteoporosis. Insight into the structure-activity relationships of the A-ring of 1 is still needed to assist the development of more potent and selective analogues as candidate chemotherapeutic agents, as well as to define the molecular mode of action. RESULTS: All possible A-ring stereoisomers of 5,6-trans-2-methyl-1,25-dihydroxyvitamin D(3) (6a-h) and their 20-epimers (7a-h) were designed and efficiently synthesized. The dependence of the affinities for vitamin D receptor (VDR) and vitamin D binding protein (DBP), as well as the HL-60 cell differentiation-inducing activity, upon the stereochemistry of the A-ring and at C20 in the side chain was evaluated. CONCLUSIONS: The binding affinities and potency of the 5,6-trans and 5,6-cis analogues were enhanced by a 2-methyl substituent in a certain orientation. Molecular docking studies based upon the X-ray crystal structure of VDR suggested that the axial 2-methyl group would be accommodated in a pocket surrounded by hydrophobic amino acid residues in the ligand binding domain, resulting in enhanced interaction.

Animals↗

PubMed2003

Corruption of host seven-transmembrane proteins by pathogenic microbes: a common theme in animals and plants?

Human diseases like AIDS, malaria, and pneumonia are caused by pathogens that corrupt host chemokine G-protein coupled receptors for molecular docking. Comparatively, little is known about plant host factors that are required for pathogenesis and that may serve as receptors for the entry of pathogenic microbes. Here, we review potential analogies between human chemokine receptors and the plant seven-transmembrane MLO protein, a candidate serving a dual role as docking molecule and defence modulator for the phytopathogenic powdery mildew fungus.

Animals↗

PubMed2026

Chemoproteomics Prioritizes Mitochondrial ADP/ATP Translocase as a Candidate Target Associated with 6PPDQ-Induced Respiratory Toxicity in Rainbow Trout Gill.

6PPD quinone (6PPDQ) is an emerging contaminant that induces acute respiratory toxicity in rainbow trout (Oncorhynchus mykiss), yet its underlying molecular mechanisms remain poorly understood. In the present study, short-term in vivo exposure of rainbow trout to 6PPDQ resulted in substantial accumulation and limited biotransformation of 6PPDQ in the gill, accompanied by pronounced gill structural damage and increased whole-fish oxygen consumption. Taking advantage of the electrophilic reactivity of the quinone moiety of 6PPDQ toward cysteine residues, we applied activity-based protein profiling (ABPP) to gill tissue. ABPP revealed marked alterations in mitochondrial cysteine reactivity and highlighted ADP/ATP translocase (ANT) as a candidate 6PPDQ-interacting mitochondrial protein. A Cys-160-containing ANT peptide within the nucleotide-binding domain of ANT was pinpointed as the covalent binding site through ABPP, Peptide-centric Local Stability Assay (PELSA), and molecular docking. Functional assays using isolated gill mitochondria showed that 6PPDQ elicited an uncoupling-like mitochondrial respiratory response that was partially attenuated by the ANT inhibitor carboxyatractyloside (CATR), supporting the functional involvement of ANT in this gill-based model. Together, these findings nominate ANT as a candidate gill mitochondrial target associated with 6PPDQ-induced acute respiratory toxicity and demonstrate the utility of chemoproteomics for prioritizing mechanistically relevant protein interactions of emerging pollutants.

Animals↗

PubMed2026

Identification of Immune Response-Related Proteomic Biomarkers in Moyamoya Disease Using Serum Olink Proteomics.

Moyamoya disease, a rare chronic cerebrovascular disorder, requires invasive digital subtraction angiography (DSA) for diagnosis. This study employed high-throughput proteomics to identify plasma biomarkers for Moyamoya disease diagnosis. We conducted immunopanel analysis using the Olink platform to evaluate 92 immune-related proteins in plasma samples from 88 Moyamoya disease patients and 88 healthy controls. Key proteins were identified through differential expression analysis, GO, and KEGG enrichment analysis. A diagnostic model was constructed using LASSO regression, Boruta algorithm, and machine learning models including random forest and XGBoost. Validation of these proteins was performed using GEO external data sets, followed by prediction of potential therapeutic drugs and molecular docking validation through pharmacogenomic databases. A total of 44 differentially expressed proteins were identified through the Olink immunopanel, with 12 downregulated and 32 upregulated. GO and KEGG analyses revealed significant enrichment of these proteins in innate immune responses and signaling pathways such as NF-kB and MAPK. Through LASSO, random forest, and protein under-area analysis, four potential biomarkers for Moyamoya disease (MGMT, SIT1, PRDX1, TRAF2) were identified. A diagnostic model using these proteins showed the highest AUC value with the XGBoost model. Additionally, TRAF2 and PRDX1 exhibited significant expression differences in Moyamoya disease patients within the GEO data set. Our study revealed the immune landscape of Moyamoya disease, identified four biomarkers, and established a variety of diagnostic models.

Humans↗

PubMedDate not supplied

Synthesis and in vitro characterization of organometallic rhenium and technetium glucose complexes against Glut 1 and hexokinase.

A series of nine organometallic technetium-99m and rhenium complexes of glucose are presented and characterized in vitro regarding their potential as surrogates of [18F]-2-fluoro-desoxy glucose ([18F]-FDG). The glucose derivatives are functionalized at positions C-1, C-2, C-3, and C-6. Different spacer lengths and chelating systems have been introduced at these sites. For the (radio)labeling, the organometallic precursors [99mTc(H2O)3(CO)3]+ and [ReBr3(CO)3](2-) respectively have been used. The resulting complexes have been characterized chemically and radiochemically. The formation of uniform products has been observed on the macroscopic (Re) and no-carrier-added level (99mTc). The Tc-99m complexes revealed good inertness against ligand exchange (Cys and His) and excellent stability in physiological buffered saline as well as in human plasma over a period of 24 h at 37 degrees C. The rhenium complexes have been tested for competitive inhibition of the (yeast) hexokinase. Only for C-2 derivatized glucose complexes with extended spacer functionalities Ki values in the millimolar and sub-millimolar range have been observed. In silico molecular docking experiments supported these experimental findings. However, the competitive inhibitors are not recognized as a pseudosubstrate of hexokinase. The cellular uptake of all 99mTc-complexes into HT-29 colon carcinoma cells via Glut1 was generally low and unspecific independent of the position at the hexose ring, the chelating systems, or the overall charge of the corresponding metal complexes. The current results seem to preclude the use of these compounds as [18F]-FDG surrogates primarily due to the low cellular uptake via Glut1.

Binding, Competitive↗

PubMed2000

A kinetic and conformational study on the interaction of tetrahydropteridines with tyrosine hydroxylase.

Tetrahydropterins are obligatory cofactors for tyrosine hydroxylase (TH), the rate-limiting enzyme of catecholamine biosynthesis. A series of synthetic analogues of 6(R)-L-erythro-5,6,7, 8-tetrahydrobiopterin (BH(4)) with different substituents in positions C2, N3, C4, N5, C6, C7, and N8 on the ring were used as active site probes for recombinant human TH. The enzyme tolerates rather bulky substituents at C6, as seen by the catalytic efficiency (V(max)/K(m)) and the coupling efficiency (mol of L-DOPA produced/mol of tetrahydropterin oxidized) of the cofactors. Substitutions at C2, C4, N5, and N8 abolish the cofactor activity of the pterin analogues. Molecular docking of BH(4) into the crystal structure of the catalytic domain of ligand-free rat TH results in complexes in which the pteridine ring pi-stacks with Phe300 and the N3 and the amino group at C2 hydrogen bonds with Glu332. The pteridine ring also establishes interactions with Leu294 and Gln310. The distance between C4a in the pteridines and the active site iron was 4.2 +/- 0.5 A for the ensemble of docked conformers. Docking of BH(4) analogues into TH also shows that the most bulky substituents at C6 can be well-accommodated within the large hydrophobic pocket surrounded by Ala297, Ser368, Tyr371, and Trp372, without altering the positioning of the ring. The pterin ring of 7-BH(4) shows proper stacking with Phe300, but the distance between the C4a and the active site iron is 0.6 A longer than for bound BH(4), a finding that may be related to the high degree of uncoupling observed for 7-BH(4).

Animals↗

PubMed2001

Influence of bulky polynuclear carcinogen lesions in a TATA promoter sequence on TATA binding protein-DNA complex formation.

The TATA binding protein (TBP) is an essential component of the transcription initiation complex that recognizes and binds to the minor groove of the TATA DNA duplex consensus sequences. The objective of this study was to determine the effect of a carcinogen-modified adenine residue, positioned site-specifically within a regulatory TATA DNA sequence, on the binding of TBP. Two 25-mer oligonucleotides with stereoisomeric 10S (+)-trans-anti- or 10R (-)-trans-anti-BPDE-N(6)-dA residues at A(1) or A(2) within the TATA sequence element (5'-...TA(1)TAAA...-3')-(5'-...TTTA(2)TA...) were synthesized (anti-BPDE-N(6)-dA denotes an adduct formed from the reaction of r7,t8-dihydroxy-t9,10-epoxy-7,8,9,10-tetrahydobenzo[a]pyrene). The formation of complexes with TBP of these two sequences in the double-stranded forms (1 nM) were studied employing electrophoretic mobility shift assays (EMSA) at different TBP concentrations (0-70 nM). The overall affinity of TBP for the BPDE-modified target DNA sequences was weakly enhanced in the case of the (+)-trans or (-)-trans lesions positioned at site A(1) with K(d) approximately 8 and 6 nM, respectively (K(d) approximately 9 nM for the unmodified TATA DNA). Higher-order TBP-DNA complexes were observed at TBP concentrations in excess of approximately 15 nM. However, the stabilities of the biologically significant monomeric TBP-DNA complexes was dramatically increased or decreased, depending on the position of the lesion (A(1) or A(2)), or on its stereochemical and conformational characteristics. A molecular docking modeling approach was employed to insert the stereoisomeric BPDE residues into the known TATA box-TBP structure [Nikolov, D. B., et al. (1996) Proc. Natl. Acad. Sci. U.S.A. 4862-4867] to rationalize these observations. Native gel electrophoresis experiments with the same duplexes without TBP indicate that none of the modified sequences exhibit unusual bending induced by the lesions, nor that they differ from one another in this respect. These results suggest that the hydrophobic, bulky BPDE residues influence the binding of TBP by mechanisms other than prebending. The efficiency of RNA transcription of TBP-controlled promoters could be strongly influenced by the presence of such bulky lesions that could adversely affect the levels of gene expression.

7,8-Dihydro-7,8-dihydroxybenzo(a)pyrene 9,10-oxide↗

PubMed2005

Interaction of phosphonate analogues of the tetrahedral reaction intermediate with 5-enolpyruvylshikimate-3-phosphate synthase in atomic detail.

The enzyme 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) catalyzes the penultimate step of the shikimate pathway and is the target of the broad-spectrum herbicide glyphosate. Since the functionality of the shikimate pathway is vital not only for plants but also for microorganisms, EPSPS is considered a prospective target for the development of novel antibiotics. We have kinetically analyzed and determined the crystal structures of Escherichia coli EPSPS inhibited by (R)- and (S)-configured phosphonate analogues of the tetrahedral reaction intermediate. Both diastereomers are competitive inhibitors with respect to the substrates of the EPSPS reaction, shikimate-3-phosphate (S3P) and phosphoenolpyruvate (PEP). Remarkably, the (S)-phosphonate (K(iS3P) = 750 nM), whose configuration corresponds to that of the genuine tetrahedral intermediate, is a much weaker inhibitor than the (R)-phosphonate analogue (K(iS3P) = 16 nM). The crystal structures of EPSPS liganded with the (S)- and (R)-phosphonates, at 1.5 and 1.9 A resolution, respectively, revealed that binding of the (R)-phosphonate induces conformational changes of the strictly conserved residues Arg124 and Glu341 within the active site. This appears to give rise to substantial structural alterations in the amino-terminal globular domain of the enzyme. By contrast, binding of the (S)-phosphonate renders the enzyme structure unchanged. Thus, EPSPS may facilitate the tight binding of structurally diverse ligands through conformational flexibility. Molecular docking calculations did not explain why the (R)-phosphonate is the better inhibitor. Therefore, we propose that the structural events during the open-closed transition of EPSPS are altered as a result of inhibitor action.

3-Phosphoshikimate 1-Carboxyvinyltransferase↗

PubMed2006

Novel effect of C75 on carnitine palmitoyltransferase I activity and palmitate oxidation.

C75 is a potential drug for the treatment of obesity. It was first identified as a competitive, irreversible inhibitor of fatty acid synthase (FAS). It has also been described as a malonyl-CoA analogue that antagonizes the allosteric inhibitory effect of malonyl-CoA on carnitine palmitoyltransferase I (CPT I), the main regulatory enzyme involved in fatty acid oxidation. On the basis of MALDI-TOF analysis, we now provide evidence that C75 can be transformed to its C75-CoA derivative. Unlike the activation produced by C75, the CoA derivative is a potent competitive inhibitor that binds tightly but reversibly to CPT I. IC50 values for yeast-overexpressed L- or M-CPT I isoforms, as well as for purified mitochondria from rat liver and muscle, were within the same range as those observed for etomoxiryl-CoA, a potent inhibitor of CPT I. When a pancreatic INS(823/13), muscle L6E9, or kidney HEK293 cell line was incubated directly with C75, fatty acid oxidation was inhibited. This suggests that C75 could be transformed in the cell to its C75-CoA derivative, inhibiting CPT I activity and consequently fatty acid oxidation. In vivo, a single intraperitoneal injection of C75 in mice produced short-term inhibition of CPT I activity in mitochondria from the liver, soleus, and pancreas, indicating that C75 could be transformed to its C75-CoA derivative in these tissues. Finally, in silico molecular docking studies showed that C75-CoA occupies the same pocket in CPT I as palmitoyl-CoA, suggesting an inhibiting mechanism based on mutual exclusion. Overall, our results describe a novel role for C75 in CPT I activity, highlighting the inhibitory effect of its C75-CoA derivative.

4-Butyrolactone↗

PubMed2006

New insights into the inhibition of human neutrophil elastase by heparin.

In the normal feedback mechanism of injury and repair in the lung, fragmented heparan sulfate proteoglycans (HSPGs) from damaged extracellular matrix and cells are believed to interact with elastases to limit their activity. An imbalance in the HSPG-elastase response may play an important role in situations where uncontrolled lung injury leads to diseases such as emphysema. To gain insight into this complex process of heparin and heparan sulfate regulation of elastases, an experimental study was undertaken to resolve the mechanism and structural requirements of heparin inhibition of human neutrophil elastase (HNE). Kinetic analyses were completed using in vitro assays with synthetic and insoluble elastin substrates in the presence of HNE and various heparin preparations (14-15 kDa; 17-19 kDa), heparin-derived oligosaccharides (4-22 saccharides), and chemically modified heparins (2-O-, 6-O-, O-, and N-desulfated). Results showed that heparin inhibits HNE by a tight-binding, hyperbolic, competitive mechanism, contrary to previous reports in the literature. A minimum length of at least 12-14 saccharides is required for inhibition, after which inhibitory activity increases with chain length (or molecular mass). Although all N- and O-sulfate groups contribute to inhibition, 2-O-sulfate groups are less critical than either N- or 6-O-sulfate groups, indicating that inhibitory activity is dependent upon the heparin fine structure. Molecular-docking simulations support the kinetic results and provide a plausible model for the size requirement, whereby positively charged, clamp-like regions at the ends of the interdomain crevice (elastase fold) are used by heparin to bridge the active site and inhibit activity.

Animals↗

PubMedDate not supplied

Comparing the quality and predictiveness between 3D QSAR models obtained from manual and automated alignment.

A set of 113 flexible cyclic urea inhibitors of human immunodeficiency virus protease (HIV-1 PR) was used to compare the quality and predictive power of CoMFA and CoMSIA models for manually or automatically aligned inhibitor set. Inhibitors that were aligned automatically with molecular docking were in agreement with information obtained from existing X-ray structures. Both alignment methods produced statistically significant CoMFA and CoMSIA models, with the best q(2) value being 0.649 and the best predictive r(2) being 0.754. The manual alignment gave statistically higher values, whereas the automated alignment gave more robust models for predicting the activities of an external inhibitor set. Both models utilized similar amino acids in the HIV-1 PR active site, supporting the idea that hydrogen bonds form between an inhibitor and the backbone carbonyl oxygens of Gly48 and Gly48' and also the backbone NH group of Asp30, Gly48, Asp29', and Gly48' of the enzyme. These results suggest that an automated inhibitor alignment can yield predictive 3D QSAR models that are well comparable to manual methods. Thus, an automated alignment method in creating 3D QSAR models is encouragable when a well-characterized structure of the target protein is available.

Binding Sites↗

PubMedDate not supplied

New Born radii deriving method for Generalized Born model.

Here we report a method to calculate Born radii, an important parameter used in a Generalized Born model. Traditional methods to derive Born radii are mostly based on a complicated formula, while our method is easier and more direct. Atoms are classified according to their atom type, and the Born radii of each type are obtained by fitting to experimental solvation free energy. The SMARTS language is used for the exact definition of atoms types, and Ullmann's subgraph isomorphism algorithm is used to deduce the environment. A generic algorithm is used for the parameter fitting because of its efficiency in searching a huge phase space, and its results are then optimized by using the conjugate gradient method. The final parameter set is fitting from a training set containing 357 molecules and is tested using a test set of 44 small organic molecules, and the average error is 0.58 kcal/mol for 36 neutral molecules and is 1.67 kcal/mol for 8 ions. The model is further tested under organic molecules, biopolymers, and a protein-inhibitor complex and yields reliable results in all these cases. This method can be used to accelerate molecular docking calculations.

Journal Article↗

PubMedDate not supplied

On updating torsion angles of molecular conformations.

A conformation of a molecule is defined by the relative positions of atoms and by the chirality of asymmetric atoms in the molecule. The three main representations for conformations of molecules are Cartesian coordinates, a distance geometry descriptor (which consists of a distance matrix and the signs of the volumes of quadruples of atoms), and internal coordinates. In biochemistry, conformational changes of a molecule are usually described in terms of internal coordinates. However, for many applications, such as molecular docking, the Cartesian coordinates of atoms are needed for computation. Although, for each conformational change, the Cartesian coordinates of atoms can be updated in linear time (which is optimal asymptotically), the constant factor becomes significant if a large number of updates are needed. Zhang and Kavraki (J. Chem. Inf. Comput. Sci. 2002, 42, 64-70) examined three methods: the simple rotations, the Denavit-Hartenberg local frames, and the atom-group local frames. On the basis of their implementations, they showed that the atom-group local frames are more efficient than the other two. In this paper, by expressing the torsion-angle change as a composition of translations and rotations, we observe that the simple rotations can be implemented in an efficient way by taking advantage of consecutive operations. Both quantitative and experimental comparisons show that the improved simple rotations, in which rotations are expressed in unit quaternions, are as efficient as the atom-group local frames and, thus, have the advantage of avoiding the need of precomputations of a set of local frames and transformations between them.

Journal Article↗

PubMed2004

Potent and selective farnesyl transferase inhibitors.

We recently described a novel series of CA(1)A(2)X peptidomimetics as farnesyl transferase inhibitors (FTIs). These compounds possess an N-(4-piperidinyl)benzamide scaffold mimicking A(1)A(2) residue. Extensive exploration of structure--activity relationships revealed that replacement of cysteine by substituted benzylimidazoles provided nanomolar FTIs with in vitro activities (18e, IC(50) = 4.60 nM on isolated enzyme, EC(50) = 20.0 nM for growth inhibition on a tumor cell line). The molecular docking of 18e and 19e in the active site of the enzyme provided details of key interactions with the protein and showed that the methionine or phenylalanine residue fits into the aryl binding site.

Alkyl and Aryl Transferases↗

PubMed2005

Impact of induced fit on ligand binding to the androgen receptor: a multidimensional QSAR study to predict endocrine-disrupting effects of environmental chemicals.

We investigated the influence of induced fit of the androgen receptor binding pocket on free energies of ligand binding. On the basis of a novel alignment procedure using flexible docking, molecular dynamics simulations, and linear-interaction energy analysis, we simulated the binding of 119 molecules representing six compound classes. The superposition of the ligand molecules emerging from the combined protocol served as input for Raptor, a receptor-modeling tool based on multidimensional QSAR allowing for ligand-dependent induced fit. Throughout our study, protein flexibility was explicitly accounted for. The model converged at a cross-validated r(2) = 0.858 (88 training compounds) and yielded a predictive r(2) = 0.792 (26 test compounds), thereby predicting the binding affinity of all compounds close to their experimental value. We then challenged the model by testing five molecules not belonging to compound classes used to train the model: the IC(50) values were predicted within a factor of 4.5 compared to the experimental data. The demonstrated predictivity of the model suggests that our approach may well be beneficial for both drug discovery and the screening of environmental chemicals for endocrine-disrupting effects, a problem that has recently become a cause for concern among scientists, environmental advocates, and politicians alike.

Benzhydryl Compounds↗

PubMed2006

Taspine: bioactivity-guided isolation and molecular ligand-target insight of a potent acetylcholinesterase inhibitor from Magnolia x soulangiana.

A bioactivity-guided approach was taken to identify the acetylcholinesterase (AChE, EC 3.1.1.7) inhibitory agent in a Magnolia x soulangiana extract using a microplate enzyme assay with Ellman's reagent. This permitted the isolation of the alkaloids taspine (1) and (-)-asimilobine (2), which were detected for the first time in this species. Compound 1 showed a significantly higher effect on AChE than the positive control galanthamine and selectively inhibited the enzyme in a long-lasting and concentration-dependent fashion with an IC(50) value of 0.33 +/- 0.07 muM. Extensive molecular docking studies were performed with human and Torpedo californica-AChE employing Gold software to rationalize the binding interaction. The results suggested ligand 1 to bind in an alternative binding orientation when compared to galanthamine. While this is located in close vicinity to the catalytic amino acid triad, the 1-AChE complex was found to be stabilized by (i) sandwich-like pi-stacking interactions between the planar aromatic ligand (1) and the Trp84 and Phe330 of the enzyme, (ii) an esteratic site anchoring with the amino side chain, and (iii) a hydrogen-bonding network.

Alkaloids↗

PubMed2024

Deciphering the ghost proteome in ovarian cancer cells by deep proteogenomic characterization.

Proteogenomics is becoming a powerful tool in personalized medicine by linking genomics, transcriptomics and mass spectrometry (MS)-based proteomics. Due to increasing evidence of alternative open reading frame-encoded proteins (AltProts), proteogenomics has a high potential to unravel the characteristics, variants, expression levels of the alternative proteome, in addition to already annotated proteins (RefProts). To obtain a broader view of the proteome of ovarian cancer cells compared to ovarian epithelial cells, cell-specific total RNA-sequencing profiles and customized protein databases were generated. In total, 128 RefProts and 30 AltProts were identified exclusively in SKOV-3 and PEO-4 cells. Among them, an AltProt variant of IP_715944, translated from DHX8, was found mutated (p.Leu44Pro). We show high variation in protein expression levels of RefProts and AltProts in different subcellular compartments. The presence of 117 RefProt and two AltProt variants was described, along with their possible implications in the different physiological/pathological characteristics. To identify the possible involvement of AltProts in cellular processes, cross-linking-MS (XL-MS) was performed in each cell line to identify AltProt-RefProt interactions. This approach revealed an interaction between POLD3 and the AltProt IP_183088, which after molecular docking, was placed between POLD3-POLD2 binding sites, highlighting its possibility of the involvement in DNA replication and repair.

Humans↗

PubMed2026

GFPT1 as a cross-ancestry validated target for degenerative spinal disease: genetic association in a Chinese cohort and functional characterization in zebrafish.

Degenerative spinal disease (DSD), including spinal stenosis and spondylosis, lacks effective pharmacological treatment. To identify druggable targets and assess cross-ancestry applicability, we integrate multi-omics analyses using Summary-data-based Mendelian Randomization (SMR), colocalization, and two-sample Mendelian randomization with European whole-blood, peripheral-blood, and CSF eQTL/pQTL datasets, followed by whole-genome sequencing (WGS) validation in a Chinese cohort. We identify 7 genes/proteins associated with spinal stenosis and 5 with spondylosis, with GFPT1, GPX1, and SERPINA1 shared by both. Two-sample MR further supports the causal associations of these targets with DSD. Phenome-wide MR prioritization selects GFPT1 and GPX1 as favorable candidates with no predicted adverse effects and potential beneficial effects on hypertension. In the Chinese cohort (67 lumbar spinal stenosis patients and 100 controls), WGS identifies 4 GFPT1 cis-eQTL loci (rs13016371, rs35392088, rs12997521, and rs13019789) associated with lumbar spinal stenosis risk; all risk alleles are linked to increased GFPT1 expression, and all 24 variant carriers show L4/L5 stenosis on imaging. Druggability analysis identifies IOX1 as the sole preclinical-stage compound targeting GFPT1, and molecular docking supports robust binding to GFPT1 (- 6.39 kcal/mol). Functional assays show that IOX1 directly inhibits GFPT1 enzymatic activity and induces fructose-6-phosphate accumulation. In zebrafish, IOX1 significantly rescues GFPT1-induced degenerative phenotypes. These findings establish GFPT1 as a cross-ancestry validated therapeutic target for DSD and nominate IOX1 as a promising disease-modifying candidate.

Animals↗