SwissParam命令行完全指南：从小分子参数化到结果获取

本文的主体翻译自：https://www.swissparam.ch/command-line.php

本文信息

• 工具名称: SwissParam Command Line Interface
• 官方网站: https://www.swissparam.ch

什么是SwissParam？

SwissParam是一个基于网络的自动参数化工具，专门为小分子生成CHARMM力场（MATCH）和MMFF力场参数。它通过命令行接口提供了灵活的参数化方式，支持非共价和共价小分子的处理，是目前分子模拟中常用的参数化工具之一。

基础使用流程

1. 检查服务器状态

在开始使用之前，首先确认SwissParam服务器是否正常运行：

curl "https://www.swissparam.ch:8443/"

如果服务器正常运行，你将收到”Hello World!”消息。如果没有响应，请联系SwissParam团队。

2. 启动参数化任务

a. 非共价小分子参数化

对于普通的非共价小分子，可以使用以下命令启动参数化：

curl -F "myMol2=@molecule.mol2" "https://www.swissparam.ch:8443/startparam?approach=both"

其中：

• molecule.mol2 是小分子的mol2文件，可以是任意文件名
• approach 是参数化方法的选择

可用的参数化方法包括：

• both (默认方法)
• mmff-based
• match

注意：使用mmff-based方法时，可以通过添加&c22或&c27来使用CHARMM22/27替代CHARMM36生成参数。

如果mol2文件不包含氢原子，可以添加&addH来在pH 7.4条件下质子化分子：

curl -F "myMol2=@molecule.mol2" "https://www.swissparam.ch:8443/startparam?approach=both&addH"

如果想要使用SMILES字符串替代mol2文件：

curl -g "https://www.swissparam.ch:8443/startparam?mySMILES=NC(=N)NC1=CC=CC=C1&approach=both"

如果没有问题，计算将被提交到服务器队列。用户将获得一个随机分配的会话编号(Session Number)，这个编号允许用户检查计算状态，并在计算成功后检索结果。

示例：使用GF1.mol2文件运行参数化，命令为：

curl -F "myMol2=@GF1.mol2" "https://www.swissparam.ch:8443/startparam?approach=both"

这里，65720367是提交的参数化任务的会话编号。

b. 共价小分子参数化

要参数化共价小分子，需要使用以下命令并指定一些参数：

curl -F "myMol2=@molecule.mol2" "https://www.swissparam.ch:8443/startparam?ligsite=l&reaction=r&protres=p&topology=t"

其中：

• molecule.mol2 是小分子的mol2文件，可以是任意文件名
• ligsite 是共价连接的配体位点（原子名称）
• reaction 是反应命名空间
• protres 是进行共价连接的蛋白质残基，可以是CYS、SER、LYS、ASP、GLU、THR、TYR
• topology 是配体的拓扑结构（反应后或反应前）

可用的反应类型包括：

示例：使用92V.mol2文件运行参数化，其中配体位点是S24，蛋白质残基是CYS，反应是disulf_form，拓扑是反应后，命令为：

curl -F "myMol2=@92V.mol2" "https://www.swissparam.ch:8443/startparam?ligsite=S24&reaction=disulf_form&protres=CYS&topology=post"

使用的参数化方法将自动选择为MMFF-based。

注意：同样可以通过添加&c22或&c27来使用CHARMM22/27替代CHARMM36。

重要提示：使用反应后拓扑时，可以指定必须删除哪些原子以获得反应前拓扑。如果这些原子没有”官方PDB名称”，请通过添加&delete=atom1,atom2来指定它们。

例如，使用CB0000002.mol2文件：

curl -F "myMol2=@CB0000002.mol2" "https://www.swissparam.ch:8443/startparam?delete=SG,H49&reaction=carbonyl_add&topology=post-cap&protres=CYS&ligsite=C32"

3. 检查参数化状态

你可以使用提交时收到的会话编号来检查作业状态。如果计算正在队列中等待轮到它，你将收到相关信息，并会被告知在它之前队列中等待的作业数量。如果作业正在运行，你将收到运行信息，并会报告运行时间。如果参数化已完成，你将被告知作业已完成。

curl "https://www.swissparam.ch:8443/checksession?sessionNumber=65720367"

4. 取消参数化任务

你可以取消当前正在运行或在队列中等待的参数化任务。以下命令将从服务器队列中移除计算：

curl "https://www.swissparam.ch:8443/cancelsession?sessionNumber=1742524"

5. 获取参数化结果

确认提交的作业已完成（见上文）后，你可以获取结果：

curl "https://www.swissparam.ch:8443/retrievesession?sessionNumber=65720367"

直接运行给定命令来获取你的结果：

curl "https://www.swissparam.ch:8443/retrievesession?sessionNumber=65720367" -o results.tar.gz

你将在你的机器上下载gzip压缩的结果文件。

实用技巧与最佳实践

📋 完整工作流程示例

# 1. 检查服务器状态
curl "https://www.swissparam.ch:8443/"

# 2. 提交参数化任务（普通小分子）
curl -F "myMol2=@ligand.mol2" "https://www.swissparam.ch:8443/startparam?approach=both&addH"

# 3. 定期检查状态（假设会话编号为12345678）
curl "https://www.swissparam.ch:8443/checksession?sessionNumber=12345678"

# 4. 下载结果
curl "https://www.swissparam.ch:8443/retrievesession?sessionNumber=12345678" -o results.tar.gz

# 5. 解压结果
tar -xzf results.tar.gz

⚡ 批量处理建议

对于多个分子的批量参数化，建议：

1. 编写脚本：使用shell脚本或Python脚本自动化处理流程
2. 会话管理：保存所有会话编号，便于后续状态检查
3. 错误处理：添加适当的错误处理机制
4. 结果整理：建立清晰的结果文件命名和组织系统

🔄 参数化方法选择指南

常见问题解答

Q1: 如何知道我的参数化任务是否成功？

A1: 使用checksession命令检查状态。如果显示作业完成，且下载的结果文件中包含了参数文件（.rtf, .par, .str），则表示参数化成功。

Q2: 参数化失败的原因有哪些？

A2: 常见失败原因包括：

• mol2文件格式错误
• 分子结构过于复杂或特殊
• 服务器负载过高
• 网络连接问题

Q3: 共价小分子参数化时如何选择正确的反应类型？

A3: 根据你的分子和目标蛋白质之间形成的共价键类型来选择。例如，如果形成的是二硫键，选择disulf_form；如果是Michael加成，选择michael_add。

Q4: 可以自定义力场参数吗？

A4: SwissParam主要提供基于CHARMM力场的标准参数。如果需要高度自定义的参数，建议使用其他专门的力场开发工具。

Q5: 结果文件的格式有哪些？

A5: 主要结果文件包括：

• .rtf - 残基拓扑文件
• .par - 参数文件
• .str - 结构文件
• .log - 日志文件

总结

SwissParam命令行工具为分子模拟研究者提供了一个强大而灵活的小分子参数化解决方案。通过其直观的命令行接口，用户可以轻松地完成从普通小分子到复杂共价分子的参数化工作。掌握这些命令行操作将大大提高分子动力学模拟前处理的效率和准确性。

无论是学术研究还是药物开发，SwissParam都是一个值得信赖的参数化工具，它让力场参数生成变得简单而可靠。