我试着在微笑中添加原子数：
from rdkit import Chem
mol=Chem.MolFromSmiles('c1ccccc(C(N)=O)1')
for i, atom in enumerate(mol.GetAtoms()):
atom.SetProp('molAtomMapNumber',str(i))
smi=Chem.MolToSmiles(mol)
print(smi)输出为：
[cH:0]1[cH:1][cH:2][cH:3][cH:4][c:5]1C:6=[O:8]然后我想把微笑分解成原子：
from rdkit import Chem
mol=Chem.MolFromSmiles('c1ccccc(C(N)=O)1')
for i, atom in enumerate(mol.GetAtoms()):
atom.SetProp('molAtomMapNumber',str(i))
print(i,atom.GetSymbol())输出为：
0 C
1 C
2 C
3 C
4 C
5 C
6 C
7 N
8 O**但我想要的是这样的东西：
0 cH
1 CH
...
7 NH2
8 O**谁能帮我弄清楚如何从上面的微笑中得到每个原子H？
1条答案
按热度按时间blmhpbnm1#
您可以使用您获得的SMILES
smi变量获得原子沿着氢。