#NEWLINE# def salt_bridge_plus_minus( mol1 , mol2 , cutoff = 4 ) : #NEWLINE# #INDENT# 'Returns pairs of plus-mins atoms, which meet salt bridge criteria\n\n    Parameters\n    ----------\n        mol1, mol2 : oddt.toolkit.Molecule object\n            Molecules to compute plus and minus pairs\n\n        cutoff : float, (default=4)\n            Distance cutoff for A-H pairs\n\n    Returns\n    -------\n        plus, minus : atom_dict-type numpy array\n            Aligned arrays of atoms forming salt bridge, firstly plus, secondly minus\n\n    ' #NEWLINE# ( m1_plus , m2_minus ) = close_contacts ( mol1 . atom_dict [ mol1 . atom_dict [ 'isplus' ] ] , mol2 . atom_dict [ mol2 . atom_dict [ 'isminus' ] ] , m2_minus ) #NEWLINE# return ( m1_plus , m2_minus )