Lawrence Berkeley National Laboratory Recent Work Title Fully Mechanically Controlled Automated Electron Microscopic Tomography
Jinxin Liu,Hongchang Li,Lei Zhang,Matthew J. Rames,Meng Zhang,Yadong Yu,Bo,Peng,C. Celis,April Xu,Qin Zou,Xu Yang,Xuefeng Chen,G. Ren
2017-01-01
Abstract:Knowledge of three-dimensional (3D) structures for asymmetric and flexible proteins is essential in understanding those proteins’ functions; but their structures are difficult to determine. Electron tomography (ET) provides a tool for imaging a single and unique biological object from a series of tilted angles, but it is challenging to image a single protein for three-dimensional (3D) reconstruction due to the imperfect mechanical control capability of the specimen goniometer under both medium to high magnification (approximately 50,000-160,000×) and under optimized beam coherence conditions. Here, we report a fully mechanical control method for automating ET data acquisition without using beam tilt/shift processes. This method could reduce the accumulation of beam tilt/shift that used to compensate the error from the mechanical control, but downgraded the beam coherence. Our method was developed by minimizing the error of the target object center during the tilting process through a closed-loop proportional-integral (PI) control algorithm. The validations by both negative staining (NS) and cryo-electron microscopy (cryo-EM) suggest that this method has a comparable capability to other ET methods in tracking target proteins while maintaining optimized beam coherence conditions for imaging. Introduction The structural and dynamic characteristics of proteins are essential for understanding their functional activity. The dynamic character of proteins hinders structural determination by conventional approaches, particularly for highly dynamic proteins, such as antibodies, lipoproteins and DNA-protein complexes . Conventional approaches, such as X-ray and electron microscopy (EM) single-particle reconstruction, require thousands to millions of different molecules to average . Averaging these proteins without prior knowledge of the protein dynamics and fluctuations could potentially fail to detect the dynamic characteristics and blur or eliminate any flexible domains. Therefore, a method to reveal the structure from each single and unique molecule is necessary. Electron tomography (ET) is a powerful tool to obtain a snapshot of a single-instance biological object from a series of tilted viewing angles. After computerized image alignment and three-dimensional (3D) reconstruction algorithms, a 3D structure can be revealed from a single and individual object, such as a section of a cell , an individual bacterium , large protein complexes 5 or even a single protein . The 3D reconstruction capability of this technique requires a set of high-resolution and high-quality images. However, imaging a target object from a series of tilted angles under high magnification is challenging, especially for imaging proteins. The imperfect mechanical design and control capability often causes ET data acquisition failure due to a significant shift from the targeted imaging area during the tilting process. For example, if an object is 1 μm away from the Eucentric height of the goniometer, the center of this object can shift away by approximately 0.6 μm at the high tilt angle of 60°, which is often larger than the imaging area under a magnification of 100,000×, resulting in a failure to track and image the object. In the past two decades, several automation-based ET software programs have been developed to control the TEM and allow precise tracking and imaging of a target object and reduce image acquisition time . The early automated method for ET acquisition utilized a pre-illuminated image to calculate the shift from the previous image by cross-correlation and then acquire the real image after compensating for this shift. To reduce the overall illumination dose to the target area induced by the pre-illumination step, later methods were developed by introducing a pre-determined tilting trajectory model of the target area 10,16 to predict and correct the shift before image acquisition. Pre-determination of the tilting trajectory of the object is challenging due to the imperfect mechanical design and control of the goniometer, the unevenness of the specimen and environmental vibrations during tilting, which could cause variance in the determined tilting trajectory. An approach to mitigate the influence of those problems has been developed by predicting specimen movements by using nearby tilt angles , which enables dynamic position tracking of the imaging area. However, the imperfect mechanical control capability of the specimen goniometer still requires compensation by electron beam tilting/shifting, particularly under medium to high magnifications (50,000-160,000×). The accumulation of beam tilt/shift processes could lead to a significant residual beam shift, which could degrade the beam coherence and lower the image quality. Because beam coherence is important for high-resolution imaging, in this study, we propose a method to maintain the optimized beam coherence by only controlling the mechanical stage for tracking and imaging the proteins under medium to high magnification. Mechanical control problems The imperfect alignment of the specimen to the Eucentric height, along with imperfections in the manufactured goniometer design, often causes the target area center to shift away from the imaging area during titling (Fig. 1). The major aspects of the imperfect design and mechanical control of the specimen stage by the goniometer can be categorized by the following three phenomena: uneven moving distances (referred to as moving control error), goniometer backlash (backlash error) and an unrepeatable tilt trajectory (trajectory error). i) Moving control error can be demonstrated by ordering the specimen stage to move a series of identical moving steps, with the actual moving distances measured by crosscorrelated images. The results showed that the actual moving distances were uneven and unrepeatable. For example, a continuous movement with a step length of 400 nm showed an actual moving distance in the range of approximately 320-420 nm (Fig. 2A). This moving control error was even worse when a smaller step was used. For instance, a continuous movement with a step length of 10 nm showed that the actual moving distance could range from 2 to 30 nm (Fig. 2A). ii) Backlash error is introduced by the backlash of mechanical components and occurs when the goniometer changes its movement direction (Fig. 2B left panel). The backlash error can be as large as approximately 1,500 and 400 nm for X and Y, respectively, on our Zeiss Libra 120 Plus TEM (Fig. 2B right panel). iii) Trajectory error may induced by the mechanically imperfect design of the goniometer, environmental vibrations, and misalignments of the tilt axis with the optical axis. For example, by repeating the goniometer tilting three times from −60° to +60° in steps of 1.5°, the trajectories of the same target were neither overlapping nor repeatable (Fig. 2C). Backlash elimination Among the above three major errors, the backlash error is the largest error that is related to the moving direction. Mechanical clearance or lost motion caused by gaps between the gears within the goniometer often generates backlash error (Fig. 2B). The backlash error can typically be significantly reduced by resetting the gear movement direction to the same direction as the previous moving direction. In our program, we move the specimen backward to the targeted moving direction by 5 μm before moving it to the targeted position. This process can reduce the backlash error within a standard deviation for X and Y motions to approximately 27 and 20 nm, respectively (Fig. 2D), which is significantly smaller than the original backlash errors of 1,500 and 400 nm (Fig. 2B). Target Position Tracking Control Tracking of the target object during tilting is still affected by the trajectory error, moving control error, and residual backlash error. The correction of each error is challenging because the errors are convoluted. In our strategy, instead of reducing each error separately, we treated all errors together as an “environmental” disturbance that interrupts the targeting position center during the tilting process. By introducing a closed-loop control system with a proportional-integral (PI) control strategy 18 (Fig. 3), we significantly suppressed those errors and successfully tracked the target position under 160,000 magnification. A closed-loop PI control compares the built-up historical data (integral) weighted against the instantaneous error of each step away from a reference point (proportional) to maintain the system (Fig. 3). This system applied to our XY positional tracking is briefly described as follows. For an example of the n tilting step, the image shift ś [n] is measured by the cross-correlation between the last and current tilting images. The accumulation of all historical image shifts is defined by Ś [n]. The goal of the control system is to maintain the accumulated shift Ś [n] as close to the target ŕ [n] as possible at each tilting step (for XY positional tracking, ŕ [n] is set as zero because a target position of zero corresponds to no positional shift between tilted images, resulting in optimal tracking). However, the actual accumulated shift Ś [n] is away from the target with a residual error of é [n], with the accumulation of such residual error defined by É [n]. Thus, we apply a goniometer motion ú [n], which can reduce both the latest residual error é [n] and its accumulated error É [n] (after the ú [n] motion is complete, we generally wait 10-15 seconds for the stage to stabilize, although this time can be modified). The PI control algorithm works to balance both residual errors via a weight k, resulting in the suggested goniometer motion ú [n]. In practice, k can be adjusted based on different microscopes and specimen holders to a value between 0 and 1. In every tilt, the achieved image shift ś [n] is interrupted by an environmental disturbance d́ [n], which will be controlled in the next loop.