Normal modes of a protein using a simplified force fieldΒΆ

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
# Normal mode calculation using a simplified harmonic force field.
# For a description of the force field, see
#
#     K. Hinsen & G.R. Kneller
#     A simplified force field for describing vibrational protein dynamics
#     over the whole frequency range
#     J. Chem. Phys. 111, 10766 (1999)
#

from MMTK import *
from MMTK.Proteins import Protein
from MMTK.ForceFields import HarmonicForceField
from MMTK.NormalModes import VibrationalModes
from MMTK.Visualization import view

# Construct system
universe = InfiniteUniverse(HarmonicForceField())
universe.protein = Protein('bala1')

# Calculate normal modes
modes = VibrationalModes(universe)

# Print frequencies
for mode in modes:
    print mode

# Show animation of the first non-trivial mode 
view(modes[6])

This Page