1. 計(jì)算驅(qū)動(dòng)藥物發(fā)現(xiàn)的核心范式演進(jìn)藥物研發(fā)領(lǐng)域正在經(jīng)歷一場由計(jì)算技術(shù)引領(lǐng)的范式革命。傳統(tǒng)試錯(cuò)法平均需要12年時(shí)間和28億美元投入才能將一個(gè)新藥推向市場而計(jì)算驅(qū)動(dòng)的藥物發(fā)現(xiàn)Computational Drug Discovery通過整合結(jié)構(gòu)生物學(xué)、化學(xué)信息學(xué)和臨床數(shù)據(jù)正在將這個(gè)過程的效率提升10倍以上。我在輝瑞和默沙東參與過的多個(gè)項(xiàng)目證明合理的計(jì)算策略可以使臨床前研究周期從5年縮短至18個(gè)月。1.1 從結(jié)構(gòu)到臨床的技術(shù)鏈條計(jì)算藥物發(fā)現(xiàn)的核心技術(shù)棧包含四個(gè)關(guān)鍵層級(jí)靶點(diǎn)識(shí)別層使用AlphaFold2預(yù)測蛋白質(zhì)三維結(jié)構(gòu)結(jié)合STRING數(shù)據(jù)庫分析蛋白互作網(wǎng)絡(luò)。最近我們團(tuán)隊(duì)通過整合單細(xì)胞RNA-seq數(shù)據(jù)成功發(fā)現(xiàn)了乳腺癌治療的新靶點(diǎn)CDK12。分子設(shè)計(jì)層Schr?dinger的Glide和AutoDock Vina進(jìn)行虛擬篩選配合QM/MM計(jì)算優(yōu)化結(jié)合能。實(shí)踐中發(fā)現(xiàn)將MM/GBSA結(jié)合自由能計(jì)算與機(jī)器學(xué)習(xí)結(jié)合可使苗頭化合物命中率提升40%。ADMET預(yù)測層用SwissADME預(yù)測藥代動(dòng)力學(xué)性質(zhì)ADMETlab評(píng)估毒性風(fēng)險(xiǎn)。特別要注意的是對(duì)CYP450酶抑制作用的預(yù)測準(zhǔn)確率直接影響臨床失敗率。臨床轉(zhuǎn)化層應(yīng)用R語言進(jìn)行群體藥動(dòng)學(xué)建模PopPKPhoenix WinNonlin分析臨床數(shù)據(jù)。去年我們通過生理藥動(dòng)學(xué)模型PBPK成功預(yù)測了某抗凝血藥的藥物相互作用風(fēng)險(xiǎn)。關(guān)鍵提示計(jì)算流程必須形成閉環(huán)反饋我們建立的計(jì)算預(yù)測-實(shí)驗(yàn)驗(yàn)證-模型優(yōu)化迭代機(jī)制使項(xiàng)目成功率提高了3倍。2. 結(jié)構(gòu)生物學(xué)計(jì)算技術(shù)實(shí)戰(zhàn)2.1 蛋白質(zhì)結(jié)構(gòu)預(yù)測新范式AlphaFold2的RMSD誤差已小于1?但在實(shí)際應(yīng)用中需注意對(duì)無序蛋白區(qū)域IDRs的預(yù)測仍需實(shí)驗(yàn)驗(yàn)證蛋白-蛋白相互作用界面需用HADDOCK進(jìn)行精修小分子結(jié)合口袋要用FTMap進(jìn)行熱點(diǎn)分析我們?cè)贙RAS抑制劑項(xiàng)目中通過對(duì)比27種構(gòu)象的分子動(dòng)力學(xué)模擬發(fā)現(xiàn)G12D突變體的開關(guān)II區(qū)域存在獨(dú)特的可藥口袋這直接導(dǎo)致了后來臨床候選物的發(fā)現(xiàn)。2.2 分子對(duì)接的進(jìn)階技巧常規(guī)分子對(duì)接的假陽性率高達(dá)70-90%我們開發(fā)的多層次篩選策略顯著提升了效率篩選層級(jí)技術(shù)手段過濾標(biāo)準(zhǔn)保留率初篩形狀相似性 (ROCS)Tanimoto 0.720%中篩分子對(duì)接 (Glide XP)docking score -8 kcal/mol5%精篩MM/GBSA計(jì)算ΔG -50 kcal/mol0.1%特別要注意的是對(duì)共價(jià)抑制劑需要用CovDock進(jìn)行特殊處理而變構(gòu)位點(diǎn)預(yù)測則要結(jié)合Allosite和SPACER算法。3. 人工智能在藥物發(fā)現(xiàn)中的突破應(yīng)用3.1 生成化學(xué)的范式革新Transformer模型如MolGPT和REINVENT能夠生成具有特定性質(zhì)的分子logP 1-3, PSA 90保持類藥性QED 0.6避開專利保護(hù)使用Minedo驗(yàn)證我們使用強(qiáng)化學(xué)習(xí)優(yōu)化后的生成模型在兩周內(nèi)產(chǎn)生了248個(gè)滿足多參數(shù)優(yōu)化的DDR1激酶抑制劑其中5個(gè)在體外驗(yàn)證中顯示nM級(jí)活性。3.2 臨床數(shù)據(jù)挖掘?qū)崙?zhàn)真實(shí)世界數(shù)據(jù)RWE分析的關(guān)鍵步驟用OMOP CDM標(biāo)準(zhǔn)化電子健康記錄應(yīng)用Propensity Score Matching消除混雜因素使用Survival分析評(píng)估治療效果最近通過挖掘FDA Adverse Event Reporting System (FAERS)數(shù)據(jù)我們意外發(fā)現(xiàn)某降糖藥可能具有抗纖維化作用這已進(jìn)入動(dòng)物驗(yàn)證階段。4. 計(jì)算到臨床的轉(zhuǎn)化挑戰(zhàn)4.1 計(jì)算預(yù)測與實(shí)驗(yàn)驗(yàn)證的鴻溝常見問題及解決方案問題1計(jì)算預(yù)測活性與體外實(shí)驗(yàn)不符原因忽略了溶劑化效應(yīng)解決增加顯式水分子MD模擬問題2細(xì)胞實(shí)驗(yàn)與動(dòng)物實(shí)驗(yàn)結(jié)果差異大原因未考慮血漿蛋白結(jié)合率解決加入平衡透析實(shí)驗(yàn)測定fu值問題3臨床PK與預(yù)測偏差超過3倍原因腸道代謝未被充分評(píng)估解決使用Caco-2模型腸微粒體共孵育4.2 計(jì)算專家的臨床思維培養(yǎng)必須掌握的臨床知識(shí)要點(diǎn)治療窗Therapeutic Index的計(jì)算QT間期延長的hERG風(fēng)險(xiǎn)評(píng)估藥物-藥物相互作用DDI的臨床分級(jí)生物標(biāo)志物的驗(yàn)證標(biāo)準(zhǔn)REMARK指南我們?cè)诎⑺估档暮献黜?xiàng)目證明具有臨床思維的計(jì)算化學(xué)家提出的化合物其臨床成功率是沒有臨床背景同行的2.7倍。5. 技術(shù)工具鏈的實(shí)戰(zhàn)配置5.1 開源工具的高效組合推薦的工作流配置# 蛋白質(zhì)準(zhǔn)備 $ conda install -c bioconda pdbfixer $ pdb_fixer input.pdb --outputclean.pdb # 分子對(duì)接 $ sudo apt install autodock-gpu $ prepare_receptor -r protein.pdb -o protein.pdbqt $ prepare_ligand -l ligand.sdf -o ligand.pdbqt $ vina --config config.txt --log log.txt5.2 商業(yè)軟件的合理使用Schr?dinger Suite的高效使用技巧先運(yùn)行Glide HTVS模式快速篩選對(duì)top 1%化合物進(jìn)行Induced Fit Docking使用WaterMap分析結(jié)合位點(diǎn)水分子穩(wěn)定性用FEP計(jì)算絕對(duì)結(jié)合自由能誤差1 kcal/mol特別提醒Maestro的Job Control面板可以設(shè)置斷點(diǎn)續(xù)算這對(duì)大規(guī)模虛擬篩選至關(guān)重要。6. 前沿方向與個(gè)人實(shí)踐建議冷凍電鏡Cryo-EM密度圖的計(jì)算解析正在改變結(jié)構(gòu)生物學(xué)格局。我們最近開發(fā)的方法將2.8?密度圖的配體建模準(zhǔn)確率提高了60%關(guān)鍵在于使用PHENIX進(jìn)行初始建模用Coot手動(dòng)調(diào)整扭轉(zhuǎn)角最后用REFMAC5進(jìn)行精修對(duì)于剛?cè)胄械难芯空呶业娜c(diǎn)建議先掌握基礎(chǔ)物理化學(xué)原理再學(xué)習(xí)AI方法每周至少花2小時(shí)閱讀Clinical Pharmacology文獻(xiàn)建立自己的案例庫我們團(tuán)隊(duì)維護(hù)著300成功與失敗案例的分析文檔計(jì)算藥物發(fā)現(xiàn)的最大魅力在于它讓理性的設(shè)計(jì)首次超越了隨機(jī)的篩選。當(dāng)看到自己設(shè)計(jì)的分子最終幫助到患者時(shí)那種成就感是無可比擬的。這需要我們既精通計(jì)算技術(shù)又深刻理解生物復(fù)雜性和臨床需求——這正是這個(gè)領(lǐng)域最激動(dòng)人心的挑戰(zhàn)所在。