if len(conditions) > 1: element = "(" + element + ")" filtering = element if asa == True: import asa_config import asa_module atoms = asa_module.generate_atom_list(mol) areas_dict = asa_module.calculate_asa(atoms, asa_config.distance) scale_dict = asa_module.normalize_areas(areas_dict) elif asa == False: scale_dict = None (mol, z_1, z_2) = apply_main(mol, filtering, f_structure_a_n, my_normalized_hydrophobicies, mem_window, scale_dict, solvation, max_it) mem1 = dummy_leaflet(pmin, pmax, z_1) mem2 = dummy_leaflet(pmin, pmax, z_2) allmol = mol.copy() allmol.append(mem1) allmol.append(mem2) allmol.center() if opm_molecule_path != None: opm = Molecule(opm_molecule_path) opm.center() if superpose != True: opm.moveBy([50, 0, 0]) allmol.append(opm) allmol.reps.add(sel='protein', style='NewCartoon') allmol.reps.add(sel='resname DUM', style='Licorice')
filtering = filtering + " or (" + element + ")" except NameError: if len(conditions) > 1: element = "(" + element + ")" filtering = element if asa == True: import asa_config import asa_module atoms = asa_module.generate_atom_list(mol) areas_dict = asa_module.calculate_asa(atoms, asa_config.distance) scale_dict = asa_module.normalize_areas(areas_dict) elif asa == False: scale_dict = None (mol, z_1, z_2) = apply_main(mol, filtering, f_structure_a_n, my_normalized_hydrophobicies, mem_window, scale_dict, solvation, max_it, 0) mem1 = dummy_leaflet(pmin,pmax,z_1) mem2 = dummy_leaflet(pmin,pmax,z_2) allmol=mol.copy() allmol.append(mem1) allmol.append(mem2) allmol.center() if opm_molecule_path != None: opm=Molecule(opm_molecule_path) opm.center() if superpose != True: opm.moveBy([50,0,0]) allmol.append(opm) allmol.reps.add(sel='protein', style='NewCartoon') allmol.reps.add(sel='resname DUM', style='Licorice')