Shape Estimation of 3-D DNA Molecules from Stereo Cryo-Electron Micro-Graphs

We introduce a 3-D parametric active contour algorithm for the shape estimation of DNA molecules from stereo cryo-electron micrographs. We consider a 3-D filament (consisting of a B-spline skeleton and a specified radial profile) and match its projections with the micrographs using an optimization algorithm. To accelerate the evaluation of the projections, we approximate the global model locally by an elongated blob-like template that is designed to be projection-steerable. This means that the 2-D projections of the template at any 3-D orientation can be expressed as a linear combination of a few basis functions. Thus, the matching of the template projections is reduced to evaluating a weighted sum of the inner-products between the basis functions and the micrographs. We choose an internal energy term that penalizes the total curvature magnitude of the curve. We also use a constraint energy term that forces the curve to have a specified length. The sum of these terms along with the image energy obtained from the matching process is minimized using a conjugate-gradient algorithm. We validate the algorithm using real as well as simulated data.