Slow normal modes of a protein using a C-alpha modelΒΆ

 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
29
30
31
32
33
34
35
36
37
38
39
40
41
# This example shows how to calculate approximate low-frequency
# modes for big proteins. For a description of the techniques,
# see
#
#    K. Hinsen
#    Analysis of domain motions by approximate normal mode calculations
#    Proteins 33, 417 (1998)
#
# and
#
#    K. Hinsen, A.J. Petrescu, S. Dellerue, M.C. Bellissent-Funel, G.R. Kneller
#    Harmonicity in slow protein dynamics
#    Chem. Phys. 261, 25 (2000)
#

from MMTK import *
from MMTK.Proteins import Protein
from MMTK.ForceFields import CalphaForceField
from MMTK.FourierBasis import FourierBasis, estimateCutoff
from MMTK.NormalModes import EnergeticModes
from MMTK.Visualization import view

# Construct system
universe = InfiniteUniverse(CalphaForceField(2.5))
universe.protein = Protein('insulin.pdb', model='calpha')

# Find a reasonable basis set size and cutoff
nbasis = max(10, universe.numberOfAtoms()/5)
cutoff, nbasis = estimateCutoff(universe, nbasis)
print "Calculating %d low-frequency modes." % nbasis

if cutoff is None:
    # Do full normal mode calculation
    modes = EnergeticModes(universe, 300.*Units.K)
else:
    # Do subspace mode calculation with Fourier basis
    subspace = FourierBasis(universe, cutoff)
    modes = EnergeticModes(universe, 300.*Units.K, subspace)

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

This Page