a <- read.table("distance_1.dat")
b <- read.table("distance_2.dat")
c <- read.table("distance_3.dat")
d <- read.table("distance_4.dat")
e <- read.table("distance_5.dat")
f <- read.table("distance_6.dat")
g <- read.table("distance_7.dat")


cex=.7

pdf("interring_distance.pdf", w=4,h=3)
par(mgp=c(2,.75,0), mar=c(3,3,1,1))
plot(a[,"V2"], type='l', col="grey40", cex=cex,
     xlab="Time [ns]", ylab="Distance [Å]",
     xlim=c(0,900), ylim=c(28.6, 35), axes=F)
lines(b[,"V2"], col="grey80", lty=3)
lines(c[,"V2"], col="grey80", lty=1)
#lines(d[,"V2"], col="grey30", lty=2)
lines(e[,"V2"], col="black", lty=1)
lines(f[,"V2"], col="red", lty=3)
lines(g[,"V2"], col="red", lty=4)

abline(v=500, col="grey60", lty=2)
legend("topright", c("E461K", "E434K", "WT", "WT-ATP"), cex=cex,
       col=c("grey40","grey80","black", "red"), lty=c(1,1,1,3), ncol=2)
mtext(1, line=-1.5, at=c(50,850), text=c("300K", "350K"), cex=cex)
axis(2)
axis(1, at=c(1,200,400,600,800), labels=c(1,20,40,60,80))
box()

dev.off()