next up previous contents index
Next: Labels Submenu Up: The Popup Menu Previous: Display Modes Submenu

Fit Submenu

   

When one has two similar structures, one often wants to compare them. What's the difference between two X-ray structures? How much did the structure change during a simulation? To answer these questions, you must first figure out how to compare two structures, which usually means that you must find the root mean square deviation (RMSD).

Formally, given N atom positions from structure x and the corresponding N atoms from structure y with a weighting factor , the RMSD is defined as:

Using this equation by itself probably won't give you the answer you are looking for. Imagine two identical structures offset by some distance. The RMSD should be 0, but the offset prevents that from happening. What you really want is the minimum RMSD between two given structures; the best fit. There are many ways to do this, but for VMD we have implemented the method of Kabsch (Acta Cryst. (1978) A34, 827-828 or see file Measure.C in the VMD source code). This algorithm computes the transformation, needed to move one structure onto another in order to minimize the RMSD.

With the mathematical prerequisites behind us, we still need to be able to specify how to choose the atoms to compare. If you want to compare all the atoms in both structures, and they both have the same number of atoms, then the problem is easy - N is everything. This occurs most often in MD simulations when the only thing different between two structures are the coordinates.

But what about homologous sequences? In this case, the number of atoms differ because while the number of residues is the same, the sidechains have different numbers of atoms. The usual solution is to determine the RMSD based solely on the backbone atoms or, in some X-ray structures where only the atoms have been determined, based on the atoms. In addition, VMD allows two other methods for fitting. Fitting by heavy atoms omits the hydrogens, since their positions are often not well determined. Fitting by ``picked atoms'' performs only a translation to bring one atom directly on top of another molecule.

Hopefully the previous discussion revealed the importance of the options available in the fit submenu. Before examining each of them in turn, you should be aware of some VMD definitions. In VMD, a ``molecule'' refers to all the atoms from a structure file. The file may contain multiple molecules under standard chemical usage, but VMD still thinks of them as one molecule. Instead, those individual parts are called ``fragments,'' for lack of a better term. With this said, the fit options in the submenu include:


next up previous contents index
Next: Labels Submenu Up: The Popup Menu Previous: Display Modes Submenu

Justin Gullingsrud
Tue Apr 6 09:22:39 CDT 1999