sapiopycommons 2024.12.4a375__py3-none-any.whl → 2024.12.9a377__py3-none-any.whl

sapiopycommons/chem/IndigoMolecules.py

@@ -10,7 +10,50 @@ indigo.setOption("render-stereo-style", "ext")
 indigo.setOption("aromaticity-model", "generic")
 indigo.setOption("render-coloring", True)
 indigo.setOption("molfile-saving-mode", "3000")
-indigo_inchi = IndigoInchi(indigo);
+indigo_inchi = IndigoInchi(indigo)
+
+
+# Function to process dative bonds in a molecule
+# Returns True if at least one dative bond (_BOND_COORDINATION) was removed
+def remove_dative_bonds_in_mol(molecule: IndigoObject) -> bool:
+    """
+    Remove all dative bonds in a molecule or a query molecule.
+    :param molecule: The molecule to remove.
+    :return: Whether there are any dative bonds in the molecule that were removed.
+    """
+    dative_bond_removed = False  # Flag to track if any dative bond was removed
+
+    bonds_idx_to_remove = []
+    for bond in molecule.iterateBonds():
+        # Check if the bond is of a dative type (_BOND_COORDINATION = 9)
+        if bond.bondOrder() == 9:  # _BOND_COORDINATION
+            atom1 = bond.source()
+            atom2 = bond.destination()
+
+            # Print bond details for debugging
+            print(f"Processing dative bond between atoms {atom1.index()} and {atom2.index()}")
+
+            # Cache bond information
+            # bond.setBondOrder(1)
+            bonds_idx_to_remove.append(bond.index())
+            dative_bond_removed = True  # Set flag to True
+
+    if not dative_bond_removed:
+        return False
+
+    molecule.removeBonds(bonds_idx_to_remove)
+    return True  # Return whether any dative bond was removed
+
+
+def remove_dative_in_reaction(reaction: IndigoObject) -> bool:
+    """
+    Remove all dative bonds in a reaction or a query reaction, from all reactants and products.
+    :param reaction: The reaction to remove dative bonds.
+    :return: Whether there are any dative bonds in the reaction that were removed.
+    """
+    reactant_dative_removed: bool = any(remove_dative_bonds_in_mol(reactant) for reactant in reaction.iterateReactants())
+    product_dative_removed: bool = any(remove_dative_bonds_in_mol(product) for product in reaction.iterateProducts())
+    return reactant_dative_removed or product_dative_removed
 
 
 def highlight_mol_substructure(query: IndigoObject, sub_match: IndigoObject):

sapiopycommons/chem/Molecules.py

@@ -9,7 +9,8 @@ from rdkit.Chem import rdMolDescriptors
 from rdkit.Chem.EnumerateStereoisomers import StereoEnumerationOptions, EnumerateStereoisomers
 from rdkit.Chem.MolStandardize import rdMolStandardize
 from rdkit.Chem.SaltRemover import SaltRemover
-from rdkit.Chem.rdchem import Mol
+from rdkit.Chem.rdChemReactions import ChemicalReaction
+from rdkit.Chem.rdchem import Mol, RWMol, Bond
 
 from sapiopycommons.chem.IndigoMolecules import indigo, renderer, indigo_inchi
 
@@ -22,6 +23,24 @@ tautomer_params.tautomerRemoveIsotopicHs = True
 enumerator = rdMolStandardize.TautomerEnumerator(tautomer_params)
 
 
+def remove_dative_bonds_from_mol(mol: Mol) -> RWMol:
+    """
+    Create a new copy of RWMol molecule and remove all dative bonds in the molecule.
+    :param mol: The original molecule
+    :return: The new molecule with dative bonds removed.
+    """
+    ret: RWMol = Chem.RWMol(mol)
+    bonds_to_remove = []
+    bond: Bond
+    for bond in ret.GetBonds():
+        if bond.GetBondType() in [Chem.BondType.DATIVER, Chem.BondType.DATIVE, Chem.BondType.DATIVEL,
+                                  Chem.BondType.DATIVEONE]:
+            bonds_to_remove.append((bond.GetBeginAtomIdx(), bond.GetEndAtomIdx()))
+    for atom1_idx, atom2_idx in bonds_to_remove:
+        ret.RemoveBond(atom1_idx, atom2_idx)
+    return ret
+
+
 def get_enhanced_stereo_reg_hash(mol: Mol, enhanced_stereo: bool) -> str:
     """
     Get the Registration Hash for the molecule by the current registration configuration.
@@ -217,7 +236,9 @@ def mol_to_sapio_substance(mol: Mol, include_stereoisomers=False,
     # This is number of H-Bond Donor
     molecule["numHBonds"] = rdMolDescriptors.CalcNumHBD(mol)
     molecule["molBlock"] = molBlock
-    rdkit_inchi = MolToInchi(mol)
+    # Create a copy of molecule before modifying it for InChI generation.
+    inchi_mol: Mol = remove_dative_bonds_from_mol(mol)
+    rdkit_inchi = MolToInchi(inchi_mol)
     # If INCHI is completely invalid, we fail this molecule.
     if not rdkit_inchi:
         MolToInchi(mol, treatWarningAsError=True)
@@ -231,13 +252,13 @@ def mol_to_sapio_substance(mol: Mol, include_stereoisomers=False,
     indigo_mol.aromatize()
     if enhanced_stereo:
         # Remove enhanced stereo layer when generating InChI as the stereo hash is generated separately for reg.
-        mol_copy: Mol = Chem.Mol(mol)
-        Chem.CanonicalizeEnhancedStereo(mol_copy)
-        molecule["inchi"] = Chem.MolToInchi(mol_copy)
-        molecule["inchiKey"] = Chem.MolToInchiKey(mol_copy)
+        Chem.CanonicalizeEnhancedStereo(inchi_mol)
+        molecule["inchi"] = Chem.MolToInchi(inchi_mol)
+        molecule["inchiKey"] = Chem.MolToInchiKey(inchi_mol)
     else:
         indigo_inchi.resetOptions()
-        indigo_inchi_str = indigo_inchi.getInchi(indigo_mol)
+        indigo_inchi_mol = indigo.loadMolecule(Chem.MolToMolBlock(inchi_mol, forceV3000=True))
+        indigo_inchi_str = indigo_inchi.getInchi(indigo_inchi_mol)
         molecule["inchi"] = indigo_inchi_str
         indigo_inchi_key_str = indigo_inchi.getInchiKey(indigo_inchi_str)
         molecule["inchiKey"] = indigo_inchi_key_str
@@ -245,7 +266,6 @@ def mol_to_sapio_substance(mol: Mol, include_stereoisomers=False,
     molecule["reg_hash"] = get_enhanced_stereo_reg_hash(mol, enhanced_stereo=enhanced_stereo)
     molecule["cxsmiles_hash"] = get_cxs_smiles_hash(mol, enhanced_stereo=enhanced_stereo)
     molecule["has_or_group"] = get_has_or_group(mol, enhanced_stereo=enhanced_stereo)
-
     return molecule
 
 

{sapiopycommons-2024.12.4a375.dist-info → sapiopycommons-2024.12.9a377.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: sapiopycommons
-Version: 2024.12.4a375
+Version: 2024.12.9a377
 Summary: Official Sapio Python API Utilities Package
 Project-URL: Homepage, https://github.com/sapiosciences
 Author-email: Jonathan Steck <jsteck@sapiosciences.com>, Yechen Qiao <yqiao@sapiosciences.com>

{sapiopycommons-2024.12.4a375.dist-info → sapiopycommons-2024.12.9a377.dist-info}/RECORD

@@ -2,8 +2,8 @@ sapiopycommons/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 sapiopycommons/callbacks/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 sapiopycommons/callbacks/callback_util.py,sha256=nb6cXK8yFq96gtG0Z2NiK-qdNaRh88bavUH-ZoBjh18,67953
 sapiopycommons/callbacks/field_builder.py,sha256=p2XacN99MuKk3ite8GAqstUMpixqugul2CsC4gB83-o,38620
-sapiopycommons/chem/IndigoMolecules.py,sha256=3f-aig3AJkKJhRmhlQ0cI-5G8oeaQk_3foJTDZCvoko,2040
-sapiopycommons/chem/Molecules.py,sha256=SQKnqdZnhYj_6HGtEZmE_1DormonRR1-nBAQ__z4gms,11485
+sapiopycommons/chem/IndigoMolecules.py,sha256=slM2y39zZFHc468c366EqR8T-GYJ24UnM9HWAqWFEwQ,3900
+sapiopycommons/chem/Molecules.py,sha256=9B4sbwbvYs50XHRn0TZiu-D1Oa3pKrI9qE5vNW8zv-U,12464
 sapiopycommons/chem/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 sapiopycommons/customreport/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 sapiopycommons/customreport/column_builder.py,sha256=0RO53e9rKPZ07C--KcepN6_tpRw_FxF3O9vdG0ilKG8,3014
@@ -53,7 +53,7 @@ sapiopycommons/webhook/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3
 sapiopycommons/webhook/webhook_context.py,sha256=D793uLsb1691SalaPnBUk3rOSxn_hYLhdvkaIxjNXss,1909
 sapiopycommons/webhook/webhook_handlers.py,sha256=MdsVK4bHffkMNmNWl0_qvu-5Lz8-qGu4Ryi7lZO1BZs,18586
 sapiopycommons/webhook/webservice_handlers.py,sha256=AFM2Va9Fpb2BvlFycIUUOdghtGiEv1Ab5jf44yjHSgU,17218
-sapiopycommons-2024.12.4a375.dist-info/METADATA,sha256=K76hHDwZ7i68OAHXdFl42xXA49M50LxMhGimBZV1PaM,3143
-sapiopycommons-2024.12.4a375.dist-info/WHEEL,sha256=C2FUgwZgiLbznR-k0b_5k3Ai_1aASOXDss3lzCUsUug,87
-sapiopycommons-2024.12.4a375.dist-info/licenses/LICENSE,sha256=HyVuytGSiAUQ6ErWBHTqt1iSGHhLmlC8fO7jTCuR8dU,16725
-sapiopycommons-2024.12.4a375.dist-info/RECORD,,
+sapiopycommons-2024.12.9a377.dist-info/METADATA,sha256=R-mI82Ozt9hd8UYMs_9XcFF6aRhT4HqUKJ7Ge88fN7k,3143
+sapiopycommons-2024.12.9a377.dist-info/WHEEL,sha256=C2FUgwZgiLbznR-k0b_5k3Ai_1aASOXDss3lzCUsUug,87
+sapiopycommons-2024.12.9a377.dist-info/licenses/LICENSE,sha256=HyVuytGSiAUQ6ErWBHTqt1iSGHhLmlC8fO7jTCuR8dU,16725
+sapiopycommons-2024.12.9a377.dist-info/RECORD,,