This is a brief overview of the simulations and the analysis. Note that there is some redundency in the files.
setup.py
-Script to run MELD simulation for calmodulin (Trial3). All necesary restraint files are included.
-See github.com/maccallumlab/meld for more information on running MELD.
job.sh
Submits MELD simulation to cluster. Requires simultaneous access to 48 GPUs with CUDA.
See github.com/maccallumlab/meld for more information on running MELD.
Contains restraint files for PRE-derived restraints (split by location and distance)
mdrmsd.py
Clustering protocol. HDBSCAN with min_cluster_size=200
Input | Description |
---|---|
cam_fill.pdb | PDB of native |
CaM_Trial{}.dcd | DCD of trajectory generated using meld "extract_trajectory extract_traj_dcd Cam_Trial{}.dcd --start=0 --end=50000" |
trajectory_1.pdb | PDB of first frame of trajectory.pdb generated manually (functions as topology file) |
Output | Description |
---|---|
k-dist.png | Visualized k-dist to help determine min_cluster_size |
clusters_0{}.dat | RMSDs of clusters |
cluster_0{}.pdb | PDBs of clusters |
plot_rmsd.py
Generates histograms of clusters.
Input | Description |
---|---|
clusters_0{}.dat | RMSDs of clusters (from mdrmsd.py) |
Output | Description |
---|---|
clusters_200h.png | Histogram image |
get_rmsds.py
Determines rmsds to native for plotting KDEs (partially redundant, easier to use when only plotting KDEs than mdrmsd.py).
Excludes flexible tails, includes peptide.
Input | Description |
---|---|
cam_fill.pdb | PDB of native |
CaM_Trial{}.dcd | DCD of trajectory generated using meld "extract_trajectory extract_traj_dcd Cam_Trial{}.dcd --start=0 --end=50000" |
trajectory_1.pdb | PDB of first frame of trajectory.pdb generated manually (used as topology file) |
Output | Description |
---|---|
rmsds_Trial{}.dat | RMSDs of native to each frame in trajectory (1 microsecond) |
plot_kde.py
Plots KDEs using .dat files from get_rmsds.py. Bandwidth of 0.1 for smoothness.
Input | Description |
---|---|
rmsds_Trial{}.dat | RMSDs of native to each frame in trajectory (1 microsecond) (from get_rmsds.py) |
Output | Description |
---|---|
kde_01_{}.png | KDE plot |