小米10pro支持存储卡扩展吗,华为matepad10.4支持多大内存卡扩展
文章目录一、引入所需库二、芳香性三、配位键四、闭环五、通过原子序数指定原子六、ChemAxon微笑拓展CXSMILES扩展七、使用坐标和异构笑容从molfile中识别手性中心
rdkit涵盖了日光微笑所有的标准功能以及一些有用的扩展,下面是扩展的部分内容
一、引入所需库#!/usr/傻傻的小蚂蚁/python#编码:utf-8从rdkit导入旧历法从rdkit导入化学.化学进口图纸二、芳香性
和氧同族的Te(碲,拼音:d,原子序数52 , 是银白色的类金属)元素也可能具有芳香性,当其连接2个芳香原子时,它贡献2个圆周率电子。
m=化学莫尔弗洛西莱斯( OC(=O)C1[te]cc C1 )平局.MolToImageFile( m,/Users/zeoy/ST/drug _ development/ST _ rdcit/img/mol 48。jpg ,legend=碲基苯-2-羧酸)#碲吩-2甲酸分子# Te原子的编号是4,下面检查其芳香性芳香原子4=m . GetAtomWithIdx(4).GetIsAromatic()print(atom4是芳香的芳族原子4) #原子4是芳香的真实的三、配位键rdkit通过-和-来支持配位键表示,箭头的方向非常重要,代表了谁提供电子
配位键不会影响起始原子的价态,只会影响指向原子的价态
铜摩尔=化学molfromselles([Cu](Cl)Cl )bipy=Chem .C1(C2=NC=CC=C2)=CC=CC=N1).C1 cccn-2 C1-c1n-3c CCC 1 .[Cu]23(Cl)Cl )mol=[Cu _ mol,bipy,bipycu]img=Draw .MolsToGridImage( mols,molsPerRow=3,subImgSize=(400,400),legends=[CuCl2 , bipy , bipycu ])img。保存(/Users/zeoy/ST/drug _ development/ST _ rdcit/img/mol 49。jpg’)
# 获取铜原子铜原子=铜摩尔GetAtoms()[0]# Cu原子在分子中的价位铜原子价=铜原子GetExplicitValence()print(第一个分子中的铜原子价,铜_化合价)#第一个分子中的铜原子价2n_atom=bipy .get atoms()[2]n _价=n _原子GetExplicitValence()print(第二个分子中的普通原子价,n _价)#第二个分子中的普通原子价3n_atom=bipycu .get atoms()[4]n _价=n _原子GetExplicitValence()print(第三分子中的普通原子价,n _价)#第三分子中的普通原子价3cu_atom=bipycu .get atoms()[12]Cu _ valence=Cu _ atom .GetExplicitValence()print(第三分子中的铜原子价,铜_化合价)#第三分子中的铜原子价四有上述的数据我们看到,在配位后铜的价态由二价变为四价,N的价态不变还是三价。
四、闭环rdkit除了支持[原子]n的方式表示闭环外,还支持%(n)的表示方法…N .的方位是0-99999;n可以不从0开始,可以是任意数字。
mol1=化学来自微笑的摩尔数( C %(1000)OC %(1000))摩尔数2=化学.摩尔来自微笑( c2oc 2 ).molfromselles( c1oc 1 )平局10 . MolToImageFile(mol 1,/Users/zeoy/ST/drug _ development/ST _ rdcit/img/mol 50。jpg ,legend=氧化丙烯)环氧丙烷
五、通过原子序数指定原子rdkit除了直接指定原子标志也支持通过原子序数来指定原子[原子序数].原子序数默认是在聪明中使用的,微笑也支持这种形式。
mol1=化学2005年,molfromselles( c1oc 1 )
mol2=化学molfromselles([# 6]1[# 8][# 6]1 )
六、ChemAxon微笑拓展CXSMILES扩展RDKit支持部分ChemAxon拓展的笑容语法功能。
1 .原子坐标原子坐标
2 .原子价值观原子值
3 .原子标签原子标签
4 .原子属性原子性质
5 .配位键(这些被翻译成双键)
6 .自由基
7 .增强立体声(这些被转换成立体声组)
rdkit .Chem.rdmolfiles.MolToCXSmiles_导出ChemAxon格式的微笑。前6个性质可通过方法rdkit .化学。rdmolfiles。moltocxsmiles()https://www . rd kit . org/docs/source/rd kit。Chem.rdmolfiles.html
m=化学molfromselles( OC )m . GetAtomWithIdx(0).SetProp(p1 , 2)m.GetAtomWithIdx(1).SetProp(p1 , 5)m.GetAtomWithIdx(1).SetProp(p2 , A1)m.GetAtomWithIdx(0).SetProp(atomLabel , O1)m.GetAtomWithIdx(1).SetProp(atomLabel , C2)sm=Chem .MolToCXSmiles(m)print(sm)# CO $ C2;O1$,atomProp:0。p 1.5:0。p2。a1:1。p 1.2 七、使用坐标和异构笑容从molfile中识别手性中心摩尔=化学2D 6 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2.5981-1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 1 6 2 3 1 0 2 4 1 0 4 0 4 5 0 2 0 0 4 6 1 0M END )平局10 . MolToImageFile(mol,/Users/zeoy/ST/drug _ development/ST _ rdcit/img/mol 51。 jpg ,legend= L-丙氨酸)我-丙氨酸
# 根据结构指定原子手性标记化学assignatomchiraltagsfrostructure(mol)#找到分子的手性中心手性中心=化学findmolchiral centers(mol)print(chiral _ center)#[(1, S)]#在笑容中也有该性质sm=化学MolToSmiles(mol)print(sm)# C[C @ H](N)C(=O)Om2=Chem .莫尔弗罗斯(sm)化学公司assignatomchiraltagsfrostructure(平方米)手性中心=化学findmolchiral centers(m2)print(chiral _ center)#[(1, S)]#当以非异构体的形式输出读取时,因为分子不在有构象所以手性信息会丢失m3=化学100毫摩尔(mol,异构miles=False)print(m3)# CC(N)C(=O)Om4=Chem .摩尔弗洛姆塞尔(立方米)化学公司106 . assignatomchiraltagsfrostructure(M4)手性中心=化学findmolchiral centers(M4)print(chiral _ center)#[]
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。