Guess j pro.m

A reasonable attempt at a function that guesses and assigns all J-couplings in a protein.

Syntax

    jmatrix=guess_j_pro(resnums,resnames,labels,coords)

Description

This function assigns J-couplings in proteins from atomic coordinates using semiempirical estimates. It is a graph-theoretical estimator with the following stages:

1. The molecular bonding graph is partitioned into connected subgraphs of size two, and one-bond J-couplings are assigned from a complete database of atom pairs. Our experienceindicates that there are fewer than 100 unique connected atom pairs in regular proteins, and that most one-bond J-couplings within those pairs can be either found in the literature, or measured in individual amino acids, or estimated with sufficient accuracy using electronic structure theory software.
2. The molecular bonding graph is partitioned into connected subgraphs of size three, and two-bond J-couplings assigned from a complete database of connected atom triples. The number of unique connected atom triples in proteins is also reasonable – fewer than 150 in regular proteins, a small enough number for an exhaustive list to be compiled from experiments, literature and electronic structure theory estimates.
3. The molecular bonding graph is partitioned into sequentially connected subgraphs of size four and dihedral angles are computed from atomic coordinates, allowing three-bond J-couplings to be assigned from a complete database of Karplus curves. Karplus curves are a well-researched topic, with specific data available for the backbone and less accurate generic curves available for the rest of the structure. The number of unique sequentially connected atom quartets found in proteins (fewer than 300, many belonging to similar structural types) is sufficiently small for a complete database of Karplus curves to be compiled from literature data, experiments, and electronic structure theory estimates.

J-couplings across more than three bonds are ignored.

Arguments

      aa_num - nspins x 1 vector giving the number of the amino acid to which each spin belongs

      aa_typ - nspins x 1 cell array of strings giving the PDB identifier of the amino acid to 
               which each spin belongs (e.g. 'TYR')

      pdb_id - nspins x 1 cell array of strings giving the PDB identifier of the protein atom 
               type to which each spin belongs (e.g. 'HE2')

      coords - nspins x 1 cell array of 3-vectors giving cartesian coordinates of each spin in Angstrom

Returns

    jmatrix - nspins x nspins sparse matrix of J-couplings in Hz

Examples

The function accepts, for example, the output of read_pdb_pro.m:

    % Parse the PDB file
    [pdb_aa_num,pdb_aa_typ,pdb_atom_id,pdb_coords]=read_pdb_pro('1D3Z.pdb',1);

    % Guess the J-couplings
    jmatrix=guess_j_pro(aa_num,aa_typ,pdb_id,coords)