## Fetch stucture
pdb <- read.pdb("1hel")
## Calculate (vibrational) normal modes
modes <- nma(pdb)
## Fluctuations
f <- fluct.nma(modes)
## Fluctuations of first non-trivial mode
f <- fluct.nma(modes, mode.inds=c(7,8))
Run the code above in your browser using DataLab