RDKit 给3D信息缺失的sdf生成三维结构

ops/2025/2/13 16:52:37/

要生成包含三维结构的 SDF 文件,可以使用 RDKit 等化学信息学工具。以下是一个 Python 脚本示例,使用 RDKit 读取 SDF 文件、生成三维结构并保存。

### 安装 RDKit
如果尚未安装 RDKit,可以通过以下命令安装:

```bash
conda install -c conda-forge rdkit
```

### Python 脚本

```python
import os
from rdkit import Chem
from rdkit.Chem import AllChem

def generate_3d_structure(input_sdf, output_sdf):
    # 读取 SDF 文件
    suppl = Chem.SDMolSupplier(input_sdf)
    writer = Chem.SDWriter(output_sdf)

    for mol in suppl:
        if mol is not None:
            # 生成三维结构
            AllChem.EmbedMolecule(mol)
            AllChem.UFFOptimizeMolecule(mol)
            # 写入新的 SDF 文件
            writer.write(mol)

    writer.close()

def process_directory(directory):
    for filename in os.listdir(directory):
        if filename.endswith(".sdf"):
            input_sdf = os.path.join(directory, filename)
            output_sdf = os.path.join(directory, f"3d_{filename}")
            generate_3d_structure(input_sdf, output_sdf)
            print(f"Processed {filename} and saved to {output_sdf}")

if __name__ == "__main__":
    # 指定目录
    directory = "path/to/your/directory"
    process_directory(directory)
```

### 脚本说明

1. **generate_3d_structure**: 读取 SDF 文件,生成三维结构并保存。
2. **process_directory**: 遍历指定目录,处理所有 `.sdf` 文件。
3. **EmbedMolecule**: 生成三维坐标。
4. **UFFOptimizeMolecule**: 使用 UFF 力场优化结构。

### 使用步骤

1. 将脚本保存为 `generate_3d_structures.py`。
2. 修改 `directory` 变量为目标目录路径。
3. 运行脚本:

```bash
python generate_3d_structures.py
```

脚本会为每个 `.sdf` 文件生成一个包含三维结构的新文件,文件名以 `3d_` 开头。


http://www.ppmy.cn/ops/158087.html

相关文章

JUC并发—1.Java集合包底层源码剖析

大纲 1.为什么要对JDK源码剖析 2.ArrayList源码一:基本原理以及优缺点 3.ArrayList源码二:核心方法的原理 4.ArrayList源码三:数组扩容以及元素拷贝 5.LinkedList源码一:优缺点和使用场景 6.LinkedList源码二:双…

数字游牧时代:IT人力外包的范式革命与文明重构

当英国工业革命时期的企业主们将生产环节外包给家庭作坊时,他们不会想到这种生产组织方式会演变为21世纪最复杂的商业形态。IT人力外包行业在经历三十年爆炸式增长后,正在经历一场静默的范式革命。这场革命不仅重构着全球IT产业链的拓扑结构,…

DatePicker 实现:日期范围截止时间为23:59:59

文章目录 需求描述实现逻辑 需求描述 在使用 Element Plus 的 el-date-picker 组件进行日期范围选择时,如果你希望选择的日期范围截止时间为所选时间的23:59:59,你可以通过设置 type 属性为 daterange,并结合使用 value-format 属性来控制时间…

Python基础语法精要

文章目录 一、Python的起源二、Python的用途三、Python的优缺点优点缺点 四、基础语法(1)常量和表达式(2)变量变量的语法(i)定义变量(ii)变量命名的规则 (3)变…

HTML5--网页前端编程(下)

HTML5–网页前端编程(下) 9.常用标签下 (1)表格标签 用来展示数据,显示数据,规整条理,可读性好 基本语法 <table><tr> <td>单元格内的文字</td> <td>单元格内的文字</td>… </tr> <tr> <td>单元格内的文字&l…

verilog练习:i2c slave 模块设计

文章目录 前言1.结构2.代码2.1 iic_slave.v2.2 sync.v2.3 wr_fsm.v2.3.1 状态机状态解释 2.4 ram.v 3. 波形展示4. 建议5. 资料总结 前言 首先就不啰嗦iic协议了&#xff0c;网上有不少资料都是叙述此协议的。 下面将是我本次设计的一些局部设计汇总&#xff0c;如果对读者有…

c/c++蓝桥杯经典编程题100道(18)括号匹配

括号匹配 ->返回c/c蓝桥杯经典编程题100道-目录 目录 括号匹配 一、题型解释 二、例题问题描述 三、C语言实现 解法1&#xff1a;栈匹配法&#xff08;难度★&#xff09; 解法2&#xff1a;计数器法&#xff08;仅限单一括号类型&#xff0c;难度★☆&#xff09; …

Java面试题(11) 整理Java面试题及参考答案

下面是最近翻阅各类博客网站收集整理的一些Java面试题&#xff0c;您值得拥有&#xff1a; 史上最全Java面试题&#xff08;带全部答案&#xff09; 2018JAVA面试题附答案(长期更新) 这是我见过最有用的java面试题&#xff0c;面试了无数公司总结的 2017 最新java面试题&…