Add contacts#
import molsysmt as msm
from molsysmt import pyunitwizard as puw
import nglview as nv
import numpy as np
molsys = msm.convert(msm.systems['chicken villin HP35']['chicken_villin_HP35.h5msm'])
view = msm.view(molsys, standard=True)
msm.thirds.nglview.add_contacts(view, selection='atom_name=="CA"', threshold='9 angstroms')
view
contacts = msm.structure.get_contacts(molsys, selection='atom_name=="CA"', threshold='9 angstroms',
output_type='pairs', output_indices='atom')
view = msm.view(molsys, standard=True)
msm.thirds.nglview.add_contacts(view, contact_pairs=contacts[0])
view