Compute native contacts ratio from a dcd trajectory file with MD-Analysis

Usage:

native_contacts.py file.dcd file.pdb file.psf [any valid selection string]

The selection string is optional and can be any valid selection string for selectAtoms() that produces identical selections in mobile and reference. e.g. "resid 6-91"

I don’t know why but the number of native contacts for the last frame is not computed, for the version of MDAnalyis I use (0.7.3).

If you want to ask me a question or leave me a message add @bougui505 in your comment.