## Loading from the file system
fn <- system.file("extdata", "1AJC_chainA.pdb", package="Rknots")
protein <- loadProtein(fn)
protein <- newKnot(protein$A)
protein.cp <- closeAndProject(protein)
# Plot the results
par(mfrow = c(1,2))
plot(protein, lwd = 2)
plot(protein.cp, lwd = 2)
Run the code above in your browser using DataLab