24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 421  |  回复: 0

MolAICal

兑换贵宾

优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!优秀!!有木有!!!

[交流] 使用MolAICal计算药物化学过滤器MCFs及MCE-18描述符

MolAICal官方地址: https://molaical.github.io  或  https://molaical.gitlab.io

1. 引言
药物化学过滤器(MCFs)的作用是剔除含有所谓 “结构警示” 的化合物。含有此类结构单元的分子要么具有不稳定或反应性基团,要么会发生生物转化,进而产生有毒代谢物或中间体。PAINS(泛测定干扰化合物)过滤器由一系列子结构过滤器组成,旨在减少筛选库中的假阳性结果、测定伪影以及非特异性生物活性分子。其设计依据是:结构中的某些片段可能会引发不良性质 —— 如反应性、螯合作用、胶体聚集形成以及染料相关效应等,这些都会干扰测定结果。由于 MCFs 和 PAINS 均涉及反应性,因此在 MolAICal 中,MCFs 的功能涵盖了 MCFs 和 PAINS 所涉及的结构。

“MCE-18”全称为 “药物化学演进,2018”,是一种原创分子描述符,用于有效评估具有药理学相关性分子的新颖性和先导化合物潜力。MCE-18 计算公式:
使用MolAICal计算药物化学过滤器MCFs及MCE-18描述符

参数说明
        AR:芳香环或杂芳香环的存在情况(0 = 不存在,1 = 存在)。
        NAR:脂环或杂脂环的存在情况(0 = 不存在,1 = 存在)。
        CHIRAL:手性中心的存在情况(0 = 不存在,1 = 存在)。
        SPIRO:螺原子的存在情况(0 = 不存在,1 = 存在)。
        sp³:sp³ 杂化碳原子的比例(范围为 0 至 1)。
        Cyc:sp³ 杂化的环碳原子比例(范围为 0 至 1)。
        Acyc:sp³ 杂化的非环碳原子比例(范围为 0 至 1)。
        NCSPTR:复合项,计算公式为(sp³ + Cyc - Acyc)。
        Q¹:归一化二次指数,用于表征结构骨架的支化程度。
如需了解更详细的 MolAICal 相关信息,请访问:https://molaical.github.io


2. 材料
2.1. 软件要求
        MolAICal:https://molaical.github.iohttps://molaical.gitlab.io
注意:确保 MolAICal 安装正确!
2.2. 示例文件
        所有必要的教程文件可从以下网址下载:
https://gitee.com/molaical/tutorials/tree/master/024-mcf_mce18



第一部分:使用 MolAICal 计算化学过滤器(MCFs)
CODE:
1.  为单个 SMILES 字符串计算 MCFs
#> molaical.exe -tool mcf -s "CC(=O)OC1=CC=CC=C1C(=O)O"

2.  为单个 SMILES 字符串计算 MCFs(不输出属性)
#> molaical.exe -tool mcf -s "CC(=O)OC1=CC=CC=C1C(=O)O" -p false

3.  计算 MCFs 并指定输出文件及格式(如 “text”“json” 或 “csv”)
#> molaical.exe -tool mcf -s "CC(=O)OC1=CC=CC=C1C(=O)O" -o result.csv -f csv

4.  计算 MCFs 并指定输出文件(不输出属性)
#> molaical.exe -tool mcf -s "CC(=O)OC1=CC=CC=C1C(=O)O" -o result.csv -f csv -p false

默认情况下,还会输出一些属性用于分子的进一步评估:
        SMILES:输入分子的 SMILES 序列。
        Valid:判断分子是否有效,值为 True(有效)或 False(无效)。
        Passes filters:判断分子是否通过 MCFs 评估,值为 True(通过)或 False(未通过)。若值为 True,表明该分子通过了 MCFs 和 PAINS 过滤器。
        LogP:辛醇 - 水分配系数,指物质在辛醇相中的浓度与在水相中的浓度之比。
        SA Score:即合成可及性分数,用于粗略评估特定分子的合成难度(评分 10 表示难,1 表示易)。该分数由分子片段的综合贡献得出。
        NP Score:即类天然产物分数,是一种贝叶斯指标,用于评估分子与天然产物所涵盖的结构空间的相似程度。
        QED:即药物相似性定量估计,是一个 [0,1] 区间内的值,用于评估分子成为可行药物候选物的可能性。与 SA 分数类似,QED 中的描述符阈值在过去十年中有所变化,当前阈值可能未涵盖最新药物。
        Molecular weight(MW):化合物中原子量的总和。

从上述网址下载 “024-mcf_mce18/molecules.smi”;“molecules.smi” 可包含多个 SMILES 字符串,因此 MolAICal 可一次性为一个或多个分子计算 MCFs。请尝试以下示例:

        从文件中计算 MCFs(文件可包含多行 SMILES,每行代表一个分子的 SMILES)
CODE:
#> molaical.exe -tool mcf -i molecules.smi -o result.csv -f csv

从文件中计算 MCFs(不输出属性,文件可包含多行 SMILES,每行代表一个分子的 SMILES)
CODE:
#> molaical.exe -tool mcf -i molecules.smi -o result.csv -f csv -p false

注意:用户可选择上述任一方式计算 MCFs。如需更详细的帮助,请参考 MolAICal 手册。


第二部分:使用 MolAICal 计算 MCE-18
CODE:
1.  为单个 SMILES 计算 MCE-18
#> molaical.exe -tool mce18 -s "CC(=O)OC1=CC=CC=C1C(=O)O"

2.  仅输出 SMILES 和 MCE-18 值
#> molaical.exe -tool mce18 -s "CC(=O)OC1=CC=CC=C1C(=O)O" -sm true

3.  输入 SMILES 并输出 CSV 格式文件
#> molaical.exe -tool mce18 -s "CC(=O)OC1=CC=CC=C1C(=O)O" -o output.csv

4.  输入 SMILES 并输出 CSV 格式文件(不输出详细内容)
#> molaical.exe -tool mce18 -s "CC(=O)OC1=CC=CC=C1C(=O)O" -o output.csv -q true

5.  输入 SMILES 并输出 CSV 格式文件(仅包含 SMILES 和 MCE-18 值)
#> molaical.exe -tool mce18 -s "CC(=O)OC1=CC=CC=C1C(=O)O" -o output.csv -sm true

6.  输入包含 SMILES 字符串的文件(每行一个)并输出 CSV 格式文件
#> molaical.exe -tool mce18 -i molecules.smi -o results.csv

7.  输入包含 SMILES 字符串的文件(每行一个)并输出 CSV 格式文件(仅包含 SMILES 和 MCE-18 值)
#> molaical.exe -tool mce18 -i molecules.smi -o results.csv -sm true

对MCE-18的 值进行分析:
从实际应用场景来看,45-78 是一个相对理想的范围
对供应商化合物库(如 ChemDiv、Enamine)的分析显示:
MCE-18 <45:分子结构简单、新颖性低,多为 “旧骨架”,吸引力有限;
MCE-18: 45-78:新颖性足够,符合当前药物化学趋势(45-63),与制药公司近期专利化合物的结构相似性高(63−78);
MCE-18 > 78:需要手动评估靶点适应性和药物相似性,因为过度复杂可能导致开发困难。


注意:用户可选择上述任一方式计算 MCE-18。如需更详细的帮助,请参考 MolAICal 手册。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 MolAICal 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 319求调剂 +3 太容易1018 2026-04-01 3/150 2026-04-01 16:40 by JourneyLucky
[考研] 找调剂 +5 楚乔乔 2026-04-01 5/250 2026-04-01 16:12 by w虫虫123
[考研] 0703求调剂 +4 zizimo 2026-03-31 4/200 2026-04-01 16:04 by yanflower7133
[考研] 08工科275求调剂,可跨考。 +5 AaAa7420 2026-03-31 5/250 2026-04-01 15:21 by 159357hjz
[考研] 349求调剂 +6 zwjjjjjj 2026-03-31 6/300 2026-04-01 09:16 by JourneyLucky
[考研] 初试301,代码085701环境工程,本硕一致,四六级已过,有二区一作,共发表5篇论文 +3 axibli 2026-04-01 3/150 2026-04-01 08:43 by i_cooler
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +12 Koxui 2026-03-28 12/600 2026-03-31 23:17 by wwytracy
[考研] 一志愿:西北大学,英一数一408-284分求调剂 +7 12.27 2026-03-27 7/350 2026-03-31 21:59 by lbsjt
[考研] 化学工程085602 305分求调剂 +28 RichLi_ 2026-03-25 36/1800 2026-03-31 14:56 by JourneyLucky
[考研] 一志愿哈尔滨工业大学材料与化工方向336分 +13 辰沐5211314 2026-03-26 13/650 2026-03-31 14:37 by 记事本2026
[考研] 274求调剂 +6 xiao爱同学 2026-03-30 6/300 2026-03-31 10:04 by cal0306
[考研] 一志愿南开大学0710生物学359求调剂 +5 兔兔兔111223314 2026-03-29 7/350 2026-03-30 18:29 by 兔兔兔111223314
[考研] 332求92调剂 +8 蕉蕉123 2026-03-28 8/400 2026-03-29 10:46 by 周梓丹
[考研] 305求调剂 +8 RuiFairyrui 2026-03-28 8/400 2026-03-29 08:22 by fmesaito
[考研] 071000生物学求调剂,初试成绩343 +7 小小甜面团 2026-03-25 7/350 2026-03-28 20:25 by 唐沐儿
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 315分求调剂 +7 26考研上岸版26 2026-03-26 7/350 2026-03-28 04:05 by fmesaito
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 一志愿上海理工能源动力(085800)310分求调剂 +3 zhangmingc 2026-03-27 4/200 2026-03-27 19:01 by 给你你注意休息
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
信息提示
请填处理意见