Lilien Lab
Department of Computer Science
Centre for Cellular and Biomolecular Research
University of Toronto

CryoBayes
Navdeep Jaitly, Marcus Brubaker, John Rubinstein, Ryan Lilien

CryoBayes is a Bayesian method for inferring 3D structure of macromolecules imaged by Electron Cryo-Microscopy experiments. We encourage experimental Structural Biologists to download and try our code. We hope that our technique will be useful as a method of automatically generating intial models given only class average images.

Our method performs ab initio inference of the three-dimensional structures of macromolecules from single particle electron cryo-microscopy experiments using class average images. The main aspects of our approach are as follows:
  • We model electron density in real space using tri-linear interpolation over a grid of points.
  • Our model is projected on to a set of rotations and the model projections are compared to class averages
  • The set of rotations are chosen based on our prior belief about how the macromolecule under study orients itself in the imaging plane.
  • A Bayesian framework is used to assess the probability of observing the class averages from a given model
  • We optimize our model to find the model that has the highest probability of generating the class averages

Sample Reconstructions

Figure: An illustration of the algorithm's progress towards its 'final' structural model across multiple cycles and iterations for three test systems: (a) A phantom synthetic data set, (b) experimental class averages for ATP synthase, and (c) experimental class averages for GroEL. Cx,y indicates cycle x, iteration y. The columns 'Final Model' show two views of the model returned by our algorithm. This model should serve as a starting point for traditional refinement.

Author Summary: Electron Cryo-Microscopy (Cryo-EM) is an experimental technique for determining the high resolution structures of proteins and protein complexes. The experimental data in a Cryo-EM experiment is a set of extremely noisy projections of the complex along unknown orientations. State of the art methods for structural inference typically consist of two phases, the creation of an initial structural model and then the refinement of this model into a final structure. Here, we present a statistical method to create an initial structural model using a small number of ensemble images. Each ensemble image (or class average) is the average of multiple projections each hypothesized to be taken from the same orientation. Our technique uses a Bayesian or probabilistic approach to infer the molecular structure which has the highest probability of generating the observed images. In other words, our technique identifies the molecular structure that is most consistent with the observed data. Our method is a novel automated solution to a critical step in the Cryo-EM structure determination pipeline. By automatically generating an initial structural model, our software can expedite the structure determination process. To facilitate its use by experimentalists, we have made our software is freely available to the structural biology community.

 

 

News

Spring 2010
CryoBayes beta has been released