Dihedrals and Ramachandran plots

Dihderal angles can be computed with the dihedral function, and an application of this function is the computation of Ramachandran plots.

MolSimToolkitShared.dihedral — Function

dihedral(at1::Atom, at2::Atom, at3::Atom, at4::Atom)

Computes the dihderal angle given four atoms of type PDBTools.Atom.

Example

julia> using PDBTools

julia> pdb = read_pdb(PDBTools.TESTPDB);

julia> C1 = pdb[11]; N2 = pdb[13]; CA2 = pdb[15]; C2 = pdb[22];

julia> phi = dihedral(C1, N2, CA2, C2) 
-36.70359f0

source

PDBTools.Ramachandran — Type

Ramachandran(prot::AbstractVector{<:PDBTools.Atom})
Ramachandran # type

The Ramachandran function receives a vector of atoms of a protein and and returns a Ramachandran object, with two fields phi and psi, containing the lists of corresponding angles, that is:

phi: C(-1) - N - CA - C
psi: N - CA - C - N(+1)

If any of the above atoms is missing, the function errors. The residues are expected to belong to a single chain and consecutive.

The resulting Ramachandran object can be plotted with the Plots.scatter function.

Example

julia> using PDBTools

julia> prot = read_pdb(PDBTools.TESTPDB, "protein");

julia> ram = Ramachandran(prot)
Ramachandran data: phi, psi vectors with 102 angles.

source

Plots.scatter — Method

scatter(ram::Ramachandran; kargs...)

Creates a Ramachandran plot given a Ramachandran object.

Arguments

ram::Ramachandran: the Ramachandran object, containing ϕ and ψ angles, resulting from the the Ramachandan function.

All other arguments are default keywords of Plots.scatter and can be adjusted to customize the plot.

Example

julia> using PDBTools, Plots

julia> prot = read_pdb(PDBTools.TESTPDB, "protein");

julia> ram = Ramachandran(prot)
Ramachandran data: phi, psi vectors with 102 angles.

julia> # plt = scatter(map) # uncomment to plot

source

Dihedral angles

The dihderal function computes the dihderal angle given four atoms:

julia> using PDBTools

julia> prot = read_pdb(PDBTools.TESTPDB, "protein");

julia> dihedral(prot[1], prot[5], prot[11], prot[13])
64.07296f0

The Ramachandran function and object are used to compute and plot Ramachandran plots for a protein structure. The call to Ramachandran(vec) where vec is a vector of Atoms returns a Ramachandran object, with fields phi and psi, containing the list of dihedral angles:

using PDBTools
prot = read_pdb(PDBTools.TESTPDB, "protein");
ram = Ramachandran(prot)

Ramachandran data: phi, psi vectors with 102 angles.

Given the ram::Ramachandran object, the scatter function from Plots can be used to produce the Ramachandran plot:

using Plots
scatter(ram)

All scatter parametes can be customized using the Plots keyword syntax.