rdkit&python | 通过反应获得断键位点与类型

对反应事先进行mapping

‘CH2:15[Mg+:19].[CH2:20]1[O:21][CH2:22][CH2:23][CH2:24]1.[Cl-:14].[OH:1][c:2]1[n:3][cH:4]c:5[cH:12][cH:13]1>>[OH:1][c:2]1[n:3][cH:4]c:5[cH:12][cH:13]1’

逻辑是根据生成物与产物之间的mapping number的不同，得到断键位置与
键的类型GetBondTypeAsDouble()

def get_changed_bonds(rxn_smi):
    reactants = Chem.MolFromSmiles(rxn_smi.split('>')