You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Hello,@FreyrS .Could you answer the following two questions for me?Thank you in advance.
① In "01-pdb_extract_and_triangulate.py",what are the functions oflines 90-110?
② What is the meaning of d >= 2.0?Can it be set to other values?
③ What's the meaning of 'iface'.Can I set the value of "masif_opts ['compute_iface ']" to False?
Here are 90-110 lines of code:
if 'compute_iface' in masif_opts and masif_opts['compute_iface']:
# Compute the surface of the entire complex and from that compute the interface.
v3, f3, _, _, _ = computeMSMS(pdb_filename,
protonate=True)
# Regularize the mesh
mesh = pymesh.form_mesh(v3, f3)
# I believe It is not necessary to regularize the full mesh. This can speed up things by a lot.
full_regular_mesh = mesh
# Find the vertices that are in the iface.
v3 = full_regular_mesh.vertices
# Find the distance between every vertex in regular_mesh.vertices and those in the full complex.
kdt = KDTree(v3)
d, r = kdt.query(regular_mesh.vertices)
d = np.square(d) # Square d, because this is how it was in the pyflann version.
assert(len(d) == len(regular_mesh.vertices))
iface_v = np.where(d >= 2.0)[0]
iface[iface_v] = 1.0
# Convert to ply and save.
save_ply(out_filename1+".ply", regular_mesh.vertices,
regular_mesh.faces, normals=vertex_normal, charges=vertex_charges,
normalize_charges=True, hbond=vertex_hbond, hphob=vertex_hphobicity,
iface=iface)
The text was updated successfully, but these errors were encountered:
Hi! 'iface' stands for 'interaction face', for Q2, d is the distance between the two interacting subelements. According to my using experience, for your Q3, if you are running masif on a single chain, it's ok to set masif_opts ['compute_iface '] = False, as for a single chain the output will be all zeros. Hope this is helpful.
Hello,@FreyrS .Could you answer the following two questions for me?Thank you in advance.
① In "01-pdb_extract_and_triangulate.py",what are the functions oflines 90-110?
② What is the meaning of d >= 2.0?Can it be set to other values?
③ What's the meaning of 'iface'.Can I set the value of "masif_opts ['compute_iface ']" to False?
Here are 90-110 lines of code:
if 'compute_iface' in masif_opts and masif_opts['compute_iface']:
# Compute the surface of the entire complex and from that compute the interface.
v3, f3, _, _, _ = computeMSMS(pdb_filename,
protonate=True)
# Regularize the mesh
mesh = pymesh.form_mesh(v3, f3)
# I believe It is not necessary to regularize the full mesh. This can speed up things by a lot.
full_regular_mesh = mesh
# Find the vertices that are in the iface.
v3 = full_regular_mesh.vertices
# Find the distance between every vertex in regular_mesh.vertices and those in the full complex.
kdt = KDTree(v3)
d, r = kdt.query(regular_mesh.vertices)
d = np.square(d) # Square d, because this is how it was in the pyflann version.
assert(len(d) == len(regular_mesh.vertices))
iface_v = np.where(d >= 2.0)[0]
iface[iface_v] = 1.0
# Convert to ply and save.
save_ply(out_filename1+".ply", regular_mesh.vertices,
regular_mesh.faces, normals=vertex_normal, charges=vertex_charges,
normalize_charges=True, hbond=vertex_hbond, hphob=vertex_hphobicity,
iface=iface)
The text was updated successfully, but these errors were encountered: