<< Chapter < Page Chapter >> Page >

The only preprocessing step that is necessary to start working with this method, however, is to perform an initial pass on the protein to extract the initial values of the dihedral angles and the constant bond lengths and angles, from the Cartesian coordinates available from PDB files, if the intention is to start from the protein's native state. This is easily done; bond lengths can be obtained by computing the distance between the bonded atoms, and bond angles by computing the angle between the vectors formed by two consecutive bonds (recall that the dot product of two vectors yields the product of their lengths times the cosine of the angle between them). Next, we present the transformations required for the Denavit-Hartenberg method.

Consider three consecutive bonds as in the figure below. Suppose that a local coordinate frame is attached at the beginning of each bond. For example, local coordinate system x i 1 , y i 1 , z i 1 is centered at atom A i 1 . Therefore, imagine that the position of each atom in three-dimensional space is specified in terms ofa frame that is centered at the previous atom. Given the frames at atom A i 2 , and atom A i 1 , one can determine how the frames at atoms A i and atom A i 1 will change in space as a consequence of a rotation around the bond that connects atoms A i 1 and A i with the dihedral angle . The correct transformation can be computed in terms of three primitive operations: two rotations and one translation. The two rotations are a rotation around the dihedral bond by the dihedral angle and a rotation around an axis perpendicular to the bond angle, by the bond angle. The translation refers to the fact that the origins of the frames are on the respective centers of the atoms connected by the bond, thus separated by bond lengths.

The denavit-hartenberg convention

To describe the position of atom i in terms of the coordinate frame centered at atom i-1, two rotations and a translation are composed.
The order in which to compose these 3 transformations, to obtain the total transformation that expresses the position of atom i in terms of frame i-1, is the following:where the rotation axes are the usual x (1,0,0) and z (0,0,1), not to be confused with the DH Local Frames. The resulting homogenous transformation is shown below.

Transformation

Homogeneous transformation to express the coordinates of atom i in terms of the frame centered at i-1
Note that θ i is the dihedral angle on bond b i and α i 1 is the bond angle between bonds b i 1 and b i . d i is the length of bond b i . For a more detailed derivation of this transformation, please read the included material in required readings.The position of any atom in the molecule can be determined by chaining matrices of the form given above. For example, suppose that b i , b i 1 , ..., b 1 , represents the sequence of bonds on the path from a particular atom a to the anchor atom a anch . Then, for atom a , its Cartesian coordinates with respect to the frame attached to the anchor atom is given by:

Equation 1

The coordinates of atom a with respect to the local frame attached to it are 0, 0, 0.
To complete the description, one can allow for rotations or translations of the local frameattached to the anchor atom with respect to some global frame. Rotations of the anchor atom with respect to a global framecause a rigid rotation of the entire polypeptide chain. To do so, one can define the rotation frame as the Euler matrix defined by the Eulerangles of the local frame of the anchor atom to the global frame. As discussed before, there are many conventions to define the Euler matrix. One ofthem, the X-Y-Z convention, defines the Euler matrix as the product of three rotation matrices: rotation around the z axis by angle α ; rotation around the y axis by angle β ; rotation around the x axis by the angle γ . The order of performing these three rotations in the X-Y-Z conventions is:rotation around x axis first, then around y axis second, and around z axis last. The resulting Euler matrix according to this convention is givenbelow:

Euler matrix

α , β , γ are the so-called Euler angles - the angles with respect to each of the Cartesian axes. The convention used here is the XYZ convention. and denote cos( α ) and sin( α ) respectively.
The Euler matrix can be applied last to the accumulating dihedral rotations in order to allow the anchor atom to move with respect to a global frame. For a more detailed explanation, please read the included material in required readings.

Required reading

  • Zhang-Kavraki 2002 [PDF] Zhang, M. and L. E. Kavraki, "A New Method for Fast and Accurate Derivation of Molecular Conformations". Journal of Chemical Information and Computer Sciences, 42:64-70, 2002.
  • Stamati-Shehu-Kavraki. Computing Forward Kinematics for Protein-like linear systems using Denavit-Hartenberg Local Frames [PDF]

    Resources

  • VMD Visual Molecular Dynamics, is an excellent tool for visualization and scripted manipulation of protein structures that uses Tcl scripting - Humphrey, W., Dalke, A. and Schulten, K., "VMD - Visual Molecular Dynamics", J. Molec. Graphics, 1996, vol. 14, pp. 33-38.
  • RasMol is mostly a viewer, but has some built-in tools. - Roger Sayle and E. James Milner-White. "RasMol: Biomolecular graphics for all", Trends in Biochemical Sciences (TIBS), September 1995, Vol. 20, No. 9, p. 374.
  • Chimera is a very powerful visualizer that handles huge structures easily. - Pettersen, E.F., Goddard, T.D., Huang, C.C., Couch, G.S., Greenblatt, D.M., Meng, E.C., and Ferrin, T.E. "UCSF Chimera - A Visualization System for Exploratory Research and Analysis." J. Comput. Chem. 25(13):1605-1612 (2004).
  • InsightII , Cerius2 and Catalyst are products for simulation, discovery and analysis that recently became commercial, and can be found here .
  • CHARMM is a simulation package based on the CHARMM force field. Brooks BR, Bruccoleri RE, Olafson BD, States DJ, Swaminathan S, Karplus M (1983). "CHARMM: A program for macromolecular energy, minimization, and dynamics calculations". J Comp Chem 4: 187217.
  • NAMD is another popular simulation package and can be obtained here . James C. Phillips, Rosemary Braun, Wei Wang, James Gumbart, Emad Tajkhorshid, Elizabeth Villa, Christophe Chipot, Robert D. Skeel, Laxmikant Kale, and Klaus Schulten. Scalable molecular dynamics with NAMD. Journal of Computational Chemistry, 26:1781-1802, 2005.
  • Amber is one of the most widely used molecular dynamics simulators due to its speed. Duan et al. A point-charge force field for molecular mechanics simulations of proteins based on condensed-phase quantum mechanical calculations Journal of Computational Chemistry Vol. 24, Issue 16. Pages 1999-2012 (2003).

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Geometric methods in structural computational biology. OpenStax CNX. Jun 11, 2007 Download for free at http://cnx.org/content/col10344/1.6
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Geometric methods in structural computational biology' conversation and receive update notifications?

Ask