Package schrodinger :: Package trajectory :: Module jcoupling :: Class Jcoupling

Class Jcoupling

               object --+    
                        |    
analysistool.AnalysisTool --+
                            |
                           Jcoupling


A class to calculate NMR J-coupling, using Karplus' relation.  This gives
an empirical relationship between NMR J-coupling constants and the torsion
dihedrals.
    3J(phi) = A cos^2(phi+theta) + B cos(phi+theta) + C
A, B and C are the constants parameterized by:
    Hu & Bax; JACS (1997)
More recent parameters are by:
    Schmidt, et al.; Biomol NMR (199)
theta is an offset. if phi is measured from the backbone (C'-N-CA-C'),
then offset needs to be provided, otherwise if offset is zero.

Instance Methods

[hide private]

__init__(self, atomsel, jtype='JHn_Ha', from_phi=True, *args, **kargs)
Initialize frame slicing information.

get_quartets(self, st, pdb_atom_names)

find_dihedrals(self, st)

doFrameInit(self, frame)
Initialize data that requires a frame instance for reference.

processFrame(self, frame)
Collect data for the provided frame.

convert_phi_to_j(self, params)

getFrameResults(self)
Return results as a numpy array for all values and frames with the first index being across frames.

Inherited from analysistool.AnalysisTool: getDataLabels, getLabel, getResults, plot, plot2D, postProcess, setSlice, shouldProcessFrame

Inherited from object: __delattr__, __format__, __getattribute__, __hash__, __new__, __reduce__, __reduce_ex__, __repr__, __setattr__, __sizeof__, __str__, __subclasshook__

Class Variables

[hide private]

PARAMS = {'JHa_C': {'A': 3.76, 'B': -1.63, 'C': 0.89, 'theta':...

DIHED_ATOMS = {'JHa_C': ['C', 'N', 'CA', 'HA'], 'JHa_N': ['HA'...

backbone_asl = ' and (backbone or (sidechain and atom.ptype CB))'

Properties

[hide private]

Inherited from object: __class__

Method Details

[hide private]

init(self, atomsel, jtype=`'JHn_Ha'`, from_phi=True, *args, **kargs)
(Constructor)

Initialize frame slicing information.

The behavior of 'start', 'stop', 'step' is the same as the python built-in range function. That is, 'start' is included in the values to be processed, but 'stop' is not. Only the largest value satisfying 'start' and 'step' that is less than 'stop' will be processed.

Overrides: object.__init__: (inherited documentation)

doFrameInit(self, frame)

Initialize data that requires a frame instance for reference.

This method will be called once, with the first frame that is processed.

Overrides: analysistool.AnalysisTool.doFrameInit: (inherited documentation)

processFrame(self, frame)

Collect data for the provided frame.

Overrides: analysistool.AnalysisTool.processFrame: (inherited documentation)

getFrameResults(self)

Return results as a numpy array for all values and frames with the first index being across frames.

So, if fr = analysis.getFrameResults(), and the analysis instance calculates three values for each frame, then fr[0] will be the three values for the first frame.

Overrides: analysistool.AnalysisTool.getFrameResults: (inherited documentation)

Class Variable Details

[hide private]

PARAMS

Value:

{'JHa_C': {'A': 3.76, 'B': -1.63, 'C': 0.89, 'theta': 120.0},
 'JHa_N': {'A': -0.88, 'B': -0.61, 'C': -0.27, 'theta': 60.0},
 'JHn_C': {'A': 4.41, 'B': -1.36, 'C': 0.24, 'theta': 180.0},
 'JHn_Cb': {'A': 2.9, 'B': -0.56, 'C': 0.18, 'theta': 60.0},
 'JHn_Ha': {'A': 7.9, 'B': -1.05, 'C': 0.65, 'theta': -60.0},
 'Jc_c': {'A': 1.51, 'B': -1.09, 'C': 0.52, 'theta': 0.0}}

DIHED_ATOMS

Value:

{'JHa_C': ['C', 'N', 'CA', 'HA'],
 'JHa_N': ['HA', 'CA', 'C', 'N'],
 'JHn_C': ['H', 'N', 'CA', 'C'],
 'JHn_Cb': ['H', 'N', 'CA', 'CB'],
 'JHn_Ha': ['H', 'N', 'CA', 'HA'],
 'Jc_c': ['C', 'N', 'CA', 'C']}

Class Jcoupling

__init__(self, atomsel, jtype='JHn_Ha', from_phi=True, *args, **kargs) (Constructor)

doFrameInit(self, frame)

processFrame(self, frame)

getFrameResults(self)

PARAMS

DIHED_ATOMS

init(self, atomsel, jtype=`'JHn_Ha'`, from_phi=True, *args, **kargs)
(Constructor)