library(bio3d)
load(file="xrayPCA_EQ.RData")

pdb.closed <- read.pdb("1XCK_chainA_noWAT_noH.pdb")
##inds <- atom.select(pdb, "///500:515///CA/")
ca.inds <- atom.select(pdb.closed, "calpha")

pdb.e434k <- read.pdb("E434K_apo_chainA_noWAT_noH.pdb")
pdb.e461k <- read.pdb("E461K_apo_chainA_noWAT_noH.pdb")
ca.inds.e434k <- atom.select(pdb.e434k, "calpha")
ca.inds.e461k <- atom.select(pdb.e461k, "calpha")


e1 <- gaps.pos$f.inds[atom2xyz(1:134)] ## from Ala3
e2 <- gaps.pos$f.inds[atom2xyz(408:522)] ## to Leu524
eq <- c(e1,e2)

e1 <- atom2xyz(2:135)
e2 <- atom2xyz(409:523)
eq.md <- c(e1,e2)




trj.inds <- seq(1, 1500, by=1)

sim191 <- NULL
prefix191 <- "/net/lutefisk/slars/groel_md/1XCK/191_1XCK_E434K/results/traj_monomer/"
prefix <- "20-50ns_noWAT_1500frames_noATP_noH_chain_"




for ( i in 1:7 ) {
  tmptrj <- read.ncdf(paste(prefix191, prefix, i , ".nc", sep=""))
  sim191$xyzfit[[i]] <- fit.xyz(fixed = pdbs$xyz[1, ],
                                mobile = tmptrj[trj.inds, ca.inds.e434k$xyz],
                                fixed.inds = core$c1A.xyz,
                                mobile.inds = core$c1A.xyz,
                                full.pdbs = FALSE)
}


for ( i in 8:14 ) {
  tmptrj <- read.ncdf(paste(prefix191, prefix, i , ".nc", sep=""))
  sim191$xyzfit[[i]] <- fit.xyz(fixed = pdbs$xyz[1,gaps.pos$f.inds],
                                mobile = tmptrj[trj.inds, ca.inds.e434k$xyz],
                                fixed.inds = core$c1A.xyz,
                                mobile.inds = core$c1A.xyz,
                                full.pdbs = FALSE)
}





tmptrj <- NULL
for ( i in 1:7 ) {
  tmptrj=rbind(tmptrj, sim191$xyzfit[[i]])
}
sim191$cis$proj <- pca.project(tmptrj[,eq.md], pc.xray)


tmptrj <- NULL
for ( i in 8:14 ) {
  tmptrj=rbind(tmptrj, sim191$xyzfit[[i]])
}
sim191$trans$proj <- pca.project(tmptrj[,eq.md], pc.xray)

sim191$xyzfit <- NULL

save(sim191, file="projected_1XCK-E434K.RData")