Source code for rxn.chemutils.miscellaneous

import logging
import re
import typing
from collections import Counter
from functools import partial
from typing import Callable, List

from rdkit.Chem import AddHs, Atom, Mol
from rxn.utilities.files import (
    PathLike,
    dump_list_to_file,
    iterate_lines_from_file,
    raise_if_paths_are_identical,
)

from .conversion import canonicalize_smiles, smiles_to_mol
from .exceptions import InvalidSmiles
from .multicomponent_smiles import (
    apply_to_multicomponent_smiles,
    list_to_multicomponent_smiles,
    multicomponent_smiles_to_list,
    sort_multicomponent_smiles,
)
from .reaction_equation import (
    ReactionEquation,
    apply_to_compound_groups,
    apply_to_compounds,
    sort_compounds,
)
from .reaction_smiles import (
    determine_format,
    parse_any_reaction_smiles,
    parse_reaction_smiles,
    to_reaction_smiles,
)

logger = logging.getLogger(__name__)
logger.addHandler(logging.NullHandler())

CHIRAL_CENTER_PATTERN = re.compile(
    r"\[([^],@]+)@+([^]]*)]"
)  # Matches stereo centres, and groups what comes before and after "@"


[docs]def is_valid_smiles(smiles: str, check_valence: bool = True) -> bool: """ Whether a given string corresponds to a valid SMILES string. Args: smiles: string to check. check_valence: whether to check the valence. Returns: True if the given SMILES is valid, else False. """ try: canonicalize_smiles(smiles, check_valence=check_valence) return True except InvalidSmiles: return False
[docs]def equivalent_smiles(*smiles: str, check_valence: bool = False) -> bool: """ Returns true if all the given SMILES strings are equivalent. Will catch the exceptions for invalid SMILES and return false in that case. Args: smiles: multiple SMILES to check for equivalence. check_valence: if True, molecules with invalid valence will be invalidated. """ try: canonical_smiles = [ canonicalize_smiles(s, check_valence=check_valence) for s in smiles ] return len(set(canonical_smiles)) == 1 except InvalidSmiles: return False
[docs]def atom_type_counter(smiles: str) -> typing.Counter[str]: """ Return a counter of atom types (as symbols). """ mol: Mol = AddHs(smiles_to_mol(smiles, sanitize=False)) atoms: List[Atom] = mol.GetAtoms() # type: ignore[call-arg,no-untyped-call] return Counter(atom.GetSymbol() for atom in atoms)
[docs]def remove_chiral_centers(smiles: str) -> str: """ Return SMILES where all the chiral centres are removed. Args: smiles: non-atom-mapped SMILES string. Returns: SMILES with no chiral information. It is not canonical. """ return re.sub(CHIRAL_CENTER_PATTERN, r"[\g<1>\g<2>]", smiles)
[docs]def remove_double_bond_stereochemistry(smiles: str) -> str: """ Return SMILES where all the E/Z information on double bonds is removed. Args: smiles: SMILES string. Returns: SMILES with no sterochemical information for double bonds. The SMILES is not guaranteed to be canonical. """ return smiles.replace("/", "").replace("\\", "")
[docs]def apply_to_any_smiles( any_smiles: str, fn: Callable[[str], str], force_multicomponent: bool = False ) -> str: """ Apply a given function to individual compound SMILES strings given in any kind of SMILES string (molecule SMILES, multicomponent SMILES, reaction SMILES). In the case of reaction SMILES, the format is kept. Args: any_smiles: any kind of SMILES string. fn: callback to apply to every compound SMILES. force_multicomponent: by default, when a SMILES string contains no ">" or "~", it is assumed to just be a normal single-component SMILES string. Providing force_multicomponent=True leads to an interpretation as a multismiles string, i.e. splitting at all the dots. Raises: Exception: different kinds of exception may be raised during parsing, or during execution of the callback. Returns: the new (molecule, multicomponent, or reaction) SMILES string after application of the callback to all the component SMILES. """ if ">" in any_smiles: # we have a reaction SMILES reaction_format = determine_format(any_smiles) reaction = parse_reaction_smiles(any_smiles, reaction_format) reaction = apply_to_compounds(reaction, fn) return to_reaction_smiles(reaction, reaction_format) elif "~" in any_smiles or force_multicomponent: # we have a multicomponent SMILES return apply_to_multicomponent_smiles(any_smiles, fn=fn, fragment_bond="~") else: # we have a single-component SMILES return fn(any_smiles)
[docs]def apply_to_smiles_groups( any_smiles: str, fn: Callable[[List[str]], List[str]] ) -> str: """ Apply a given function to groups of SMILES strings given in any multicomponent SMILES or reaction SMILES. This function can typically be used for sorting or shuffling precursors, products, etc. In the case of reaction SMILES, the format is kept. Args: any_smiles: any kind of SMILES string. fn: callback to apply to every compound SMILES. Raises: Exception: different kinds of exception may be raised during parsing, or during execution of the callback. Returns: the new (multicomponent, or reaction) SMILES string after application of the callback to all the groups of SMILES. """ if ">" in any_smiles: # we have a reaction SMILES reaction_format = determine_format(any_smiles) reaction = parse_reaction_smiles(any_smiles, reaction_format) reaction = apply_to_compound_groups(reaction, fn) return to_reaction_smiles(reaction, reaction_format) else: # we have a multicomponent SMILES compounds = multicomponent_smiles_to_list(any_smiles, fragment_bond="~") return list_to_multicomponent_smiles(fn(compounds), fragment_bond="~")
[docs]def canonicalize_any( any_smiles: str, check_valence: bool = True, sort_molecules: bool = False, fallback_value: typing.Optional[str] = None, ) -> str: """ Canonicalize any SMILES string (molecule SMILES, multicomponent SMILES, reaction SMILES). In the case of reaction SMILES, the format is kept. Args: any_smiles: any kind of SMILES string. check_valence: if False, will not do any valence check. sort_molecules: whether to sort the compounds alphabetically at the same time. fallback_value: what value to returns when the canonicalization is unsuccessful. Default: no fallback, will propagate the exception. Raises: Exception: different kinds of exception may be raised during parsing. InvalidSmiles: for canonicalization errors. Returns: the canonical (molecule, multicomponent, or reaction) SMILES string. """ try: fn = partial(canonicalize_smiles, check_valence=check_valence) canonical_smiles = apply_to_any_smiles(any_smiles, fn) if sort_molecules: canonical_smiles = sort_any(canonical_smiles) return canonical_smiles except Exception as e: if fallback_value is not None: logger.debug(f'Error when canonicalizing "{any_smiles}": {e}') return fallback_value raise
[docs]def canonicalize_file( input_file: PathLike, output_file: PathLike, check_valence: bool = True, fallback_value: str = "", sort_molecules: bool = False, ) -> None: raise_if_paths_are_identical(input_file, output_file) logger.info(f'Canonicalizing file "{input_file}" -> "{output_file}".') # We formulate it as a generator, so that the file below is written directly canonical = ( canonicalize_any( line, check_valence=check_valence, fallback_value=fallback_value, sort_molecules=sort_molecules, ) for line in iterate_lines_from_file(input_file) ) dump_list_to_file(canonical, output_file)
[docs]def sort_any(any_smiles: str) -> str: """ Sort any SMILES string (molecule SMILES, multicomponent SMILES, reaction SMILES). For single-component SMILES, the fragments will be reorderd. In the case of reaction SMILES, the format is kept. Args: any_smiles: any kind of SMILES string. Raises: Exception: different kinds of exception may be raised during parsing. Returns: the sorted SMILES string. """ if ">" in any_smiles: # we have a reaction SMILES reaction_format = determine_format(any_smiles) reaction = parse_reaction_smiles(any_smiles, reaction_format) reaction = sort_compounds(reaction) return to_reaction_smiles(reaction, reaction_format) else: # we call the same function for single- and multi-component SMILES return sort_multicomponent_smiles(any_smiles)
[docs]def get_individual_compounds(any_smiles: str) -> List[str]: """ Get the individual compound SMILES strings starting from any SMILES string (multicomponent SMILES, reaction SMILES). Single-component SMILES with dots are interpreted as multicomponent SMILES strings. Args: any_smiles: any kind of SMILES string. Raises: Exception: different kinds of exception may be raised during parsing. Returns: List of individual compound SMILES. """ if ">" in any_smiles: # We have a reaction SMILES reaction = parse_any_reaction_smiles(any_smiles) return list(reaction.iter_all_smiles()) else: # We interpret it as a multicomponent SMILES. # We use "~" as a fragment bond even if it is not actually needed. return multicomponent_smiles_to_list(any_smiles, fragment_bond="~")
[docs]def merge_reactions(*reactions: ReactionEquation) -> ReactionEquation: """Merge several reactions into one. Useful when ReactionEquation is used to store partial equations.""" reactants = [] agents = [] products = [] for reaction in reactions: reactants.extend(reaction.reactants) agents.extend(reaction.agents) products.extend(reaction.products) return ReactionEquation(reactants=reactants, agents=agents, products=products)
[docs]def mol_has_atom_mapping(mol: Mol) -> bool: """ Whether at least one atom of an RDKit Mol contains an atom map number. Args: mol: RDKit Mol. """ atom: Atom for atom in mol.GetAtoms(): # type: ignore[call-arg,no-untyped-call] if atom.GetAtomMapNum() != 0: return True return False
[docs]def smiles_has_atom_mapping(smiles: str) -> bool: """ Whether at least one atom of a compound SMILES contains an atom map number. Args: smiles: compound SMILES. """ mol = smiles_to_mol(smiles, sanitize=False) return mol_has_atom_mapping(mol)