Skip to content

Instantly share code, notes, and snippets.

@chenyuxiang0425
Created March 9, 2022 10:27
Show Gist options
  • Save chenyuxiang0425/5d0135d7a018806b5fde67e0b0ce9312 to your computer and use it in GitHub Desktop.
Save chenyuxiang0425/5d0135d7a018806b5fde67e0b0ce9312 to your computer and use it in GitHub Desktop.
covert smiles to .mol files
import pandas as pd
from rdkit.Chem import AllChem as Chem
import argparse
def smiles2molfile(smiles):
""" smiles to mol file
@parma molecule_name
@return None
@obtain a file
"""
m = Chem.MolFromSmiles(smiles)
m2 = Chem.AddHs(m)
print(Chem.MolToMolBlock(m2),file=open("data/"+str(smiles)+".mol",'w+'))
def main():
parser = argparse.ArgumentParser(description='filename')
parser.add_argument('filename')
args = parser.parse_args()
data = pd.read_csv(args.filename)
smiles_lst=data["smiles"]
for i in range(len(smiles_lst)):
smiles = smiles_lst[i]
smiles2molfile(smiles)
if __name__ == "__main__":
main()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment