in silico
Phase 1: In Silico Validation
- Objective
- Determine computationally whether BHB can bind and inhibit HDAC1/2/3 catalytic zinc at physiologically achievable interstitial concentrations (0.5-2.0 mM), whether the predicted Ki is consistent with ≥30% inhibition, and whether existing transcriptomic/proteomic datasets support the downstream causal chain (HDAC inhibition → H3K9ac/H3K27ac at Bdnf/Arc/SynGAP → spine turnover genes; HDAC3 inhibition → M2 microglial genes).
- Estimated cost
- €0-2000 (compute credits + software licenses)
- Estimated duration
- 4-8 weeks
- Success criteria
- Predicted BHB Ki for HDAC1/2/3 · Ki ≤ 2.0 mM for at least HDAC2 and HDAC3 (consistent with ≥30% inhibition at 0.5-2.0 mM) · (FEP+ free-energy calculation with 3 independent replicates, SEM < 0.5 kcal/mol)
- Overlap between KD-induced and HDAC-inhibition-induced gene sets · ≥30% overlap at FDR < 0.05, with Bdnf and Arc in the intersection · (Fisher exact test on DESeq2/ChIPseeker outputs)
- ODE model predicts spine turnover increase · Predicted increase ≥25% at 1.0 mM interstitial BHB · (COPASI simulation with 1000 Monte Carlo parameter draws)
- Statistical power for Phase 2 · n ≤ 12/group achieves power ≥0.80 for HDAC activity reduction ≥30% · (simr simulation with 1000 iterations)
- Go if
- Predicted Ki ≤ 2.0 mM for HDAC2/3 AND ≥30% gene-set overlap AND ODE predicts ≥25% spine turnover increase AND n ≤ 12/group is sufficient
- No-go if
- Predicted Ki > 5.0 mM for all class I HDACs OR <10% gene-set overlap OR ODE predicts <10% spine turnover increase at 2.0 mM BHB
- Pivot if
- Ki between 2.0-5.0 mM OR 10-30% overlap: pivot to testing acetoacetate as alternative mediator, or to lower BHB threshold hypothesis (HDAC inhibition only at ≥2.0 mM), or to combination with other KD-derived metabolites (e.g., β-hydroxybutyrylation-independent mechanisms)
- Risks
- Docking/FEP cannot accurately model zinc-coordinating inhibitor binding (force-field limitations)Probability: mediumMitigation: Cross-validate with three scoring functions (Vina, Glide, MM-GBSA) and compare to experimental Ki of known HDAC inhibitors; if inconsistent, rely more on meta-analysis and ODE
- Public KD datasets use different brain regions, ages, or KD compositions, confounding meta-analysisProbability: highMitigation: Stratify by region (cortex vs hippocampus) and age; use only adult (P60+) cortex datasets for primary analysis; run sensitivity analysis excluding heterogeneous datasets
- ODE model parameters (HDAC turnover, transcription rates) poorly constrainedProbability: highMitigation: Use Bayesian priors from literature; report credible intervals; identify which parameters need Phase 2 measurement
- Compute resources insufficient for FEP+ on 3 HDAC isoformsProbability: lowMitigation: Use cloud GPU (AWS p4d or Azure NDv4) on demand; or reduce to HDAC2 and HDAC3 only