From 2f16896dc2978f69b4ed7aab672291de53ac44b0 Mon Sep 17 00:00:00 2001
From: hyunp2 <42776897+hyunp2@users.noreply.github.com>
Date: Mon, 18 Sep 2023 10:40:11 -0500
Subject: [PATCH] Update prepare_data_from_sdf.py

---
 utils/prepare_data_from_sdf.py | 75 +++++++++++++++++++++++++++++++++-
 1 file changed, 74 insertions(+), 1 deletion(-)

diff --git a/utils/prepare_data_from_sdf.py b/utils/prepare_data_from_sdf.py
index 31505b4d..4de13df0 100644
--- a/utils/prepare_data_from_sdf.py
+++ b/utils/prepare_data_from_sdf.py
@@ -18,6 +18,79 @@
 
 from docopt import docopt
 
+def prepare_sdf():
+    # Parse args
+    args = docopt(__doc__)
+    sdf_path = args.get('--sdf_path')
+    output_path = args.get('--output_path')
+    no_filters = args.get('--no_filters')
+    verbose = args.get('--verbose')
+    
+    # Load data
+    conformers = Chem.SDMolSupplier(sdf_path)
+    smiles = []
+    errors = 0
+    for num, sdf in enumerate(conformers):
+        try:
+            smiles.append(Chem.MolToSmiles(sdf))
+        except:
+            errors +=1
+    if verbose:
+        print("Original num entries: \t%d" % len(conformers))
+        print("Number parsed by RDKit: %d" % len(smiles))
+        print("Number of errors: \t%d" % errors)
+
+    # Filter SMILES for permitted atom types
+    smiles_filt = []
+    errors = 0
+    for i, smi in enumerate(smiles):
+        if frag_utils.check_smi_atom_types(smi):
+            smiles_filt.append(smi)
+        else:
+            errors +=1
+    
+    if i % 1000 == 0 and verbose:
+        print("\rProcessed smiles: %d" % i, end='')
+        
+    if verbose:
+        print("Original num entries: \t\t\t%d" % len(smiles))
+        print("Number with permitted atom types: \t%d" % len(smiles_filt))
+        print("Number of errors: \t\t\t%d" % errors)
+
+    # Fragment dataset
+    fragmentations = frag_utils.fragment_dataset(smiles_filt, linker_min=3, fragment_min=5, min_path_length=2, linker_leq_frags=True, verbose=False)
+
+    if verbose:
+        print("Processed smiles: \t%d" % len(smiles_filt))
+        print("Num fragmentations: \t%d" % len(fragmentations))
+
+    # Filter fragmentions based on 2D properties
+    if no_filters:
+        fragmentations_filt = fragmentations
+    else:
+        fragmentations_filt = frag_utils.check_2d_filters_dataset(fragmentations, n_cores=3)
+
+    print("Number fragmentations: \t\t%d" % len(fragmentations))
+    print("Number passed 2D filters: \t%d" % len(fragmentations_filt))
+
+    # Calculate structural information
+    fragmentations_new, distances, angles, fails = frag_utils.compute_distance_and_angle_dataset(fragmentations_filt, sdf_path, dataset="CASF", verbose=False)
+
+    if verbose:
+        print("Number of successful fragmentations: \t\t%d" % len(distances))
+        print("Number failed fragmentations: \t\t\t%d" % fails[0]) 
+        print("Number of SDF errors: \t\t\t%d" % fails[1])
+        print("Number of fragmentations without conformers: \t%d" % (len(fragmentations_filt)-len(fragmentations_new)))
+
+
+    # Write data to file
+    # Format: full_mol (SMILES), linker (SMILES), fragments (SMILES), distance (Angstrom), angle (Radians)
+    with open(output_path, 'w') as f:
+        for fragmentation, dist, ang in zip(fragmentations_new, distances, angles):
+            f.write("%s %s %s %s %s\n" % (fragmentation[0], fragmentation[1], fragmentation[2], dist, ang))
+            
+    print("Done")
+
 if __name__ == "__main__":
     # Parse args
     args = docopt(__doc__)
@@ -89,4 +162,4 @@
         for fragmentation, dist, ang in zip(fragmentations_new, distances, angles):
             f.write("%s %s %s %s %s\n" % (fragmentation[0], fragmentation[1], fragmentation[2], dist, ang))
             
-    print("Done")
\ No newline at end of file
+    print("Done")