name = 'Cyanide' C = Atom('C') N = Atom('N') b = Bond(C, N, 3) mm3_atom_type = { C : 4, N : 10 } db = DanglingBond(C, Dihedral(None, N, C, None, 0.0)) hook1 = Hook(None, [db])