[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [ccp4bb]: How to make a movie
*** For details on how to be removed from this list visit the ***
*** CCP4 home page http://www.ccp4.ac.uk ***
Jeff Taylor wrote:
>
> Edward Berry wrote:
> >
> > That sounds almost too simple. Are bond lengths/angles preserved during
> > the morph? if a Histidine rotated 180* about CHI2, wouldn't the NE1 and
> > CE2 atoms pass through each other on the way to their new positions?
> > Do any morphing programs preserve bond lengths and angles during the
> > morph? Perhaps one could express the structure in terms of dihedral
> > angles and interpolate linearly between the starting and final values
> > for the dihedral angles. Still in the case of large-scale movements
> > you might get one domain passing through another on the way to its
> > final position, where in nature it makes a two-stage movement to go
> > around.
> >
> > ed
>
> Yes it is too simple. I just make a linear interpolation between the
> old and new coordinates. This doesn't preserve any physical aspects of
> the molecule and the intermediate structures can be severely mangled. A
> 180 degree rotation of a Histidine (or tyrosine, phenylalanine, etc.)
> will look really funny. To make a physically reasonable morph one would
> need to employ some kind of directed molecular dynamics, or at least
> incorporate some kind of energy minimization for the intermediate
> steps. What a mess. On the bright side, the linear interpolations
> often work quite well, provided that the motions aren't too large or
> complex, AND they aren't taken too seriously.
>
yes, I can see that for rotations of 10-20 degrees the linear interpolation
would be a good approximation - the distortion would not be discernible.
and this would be the case for the majority of conformational changes
observed in different crystal structures. Sorry for nit-picking.
Still I'm curious if there are slightly more sophisticated approaches
that don't take as much CPU-time or know-how as directed molecular dynamics
a la Izrailev et al. (Biophys J 1999 77:1753-68 Steered molecular dynamics
simulation of the Rieske subunit motion in the cytochrome bc(1) complex.)
http://www.ncbi.nlm.nih.gov:80/entrez/query.fcgi?cmd=Retrieve&db=PubMed&list_uids=10512801
Ed