SwarmPS: Rapid, semi-automated single particle selection software

Single particle analysis (SPA) coupled with high-resolution electron cryo-microscopy is emerging as a powerful technique for the structure determination of membrane protein complexes and soluble macromolecular assemblies. Current estimates suggest that ∼104–105 particle projections are required to attain a 3 Å resolution 3D reconstruction (symmetry dependent). Selecting this number of molecular projections differing in size, shape and symmetry is a rate-limiting step for the automation of 3D image reconstruction. Here, we present SwarmPS, a feature rich GUI based software package to manage large scale, semi-automated particle picking projects. The software provides cross-correlation and edge-detection algorithms. Algorithm-specific parameters are transparently and automatically determined through user interaction with the image, rather than by trial and error. Other features include multiple image handling (∼102), local and global particle selection options, interactive image freezing, automatic particle centering, and full manual override to correct false positives and negatives. SwarmPS is user friendly, flexible, extensible, fast, and capable of exporting boxed out projection images, or particle coordinates, compatible with downstream image processing suites.

Position, rotation, and scale invariant recognition of images using higher-order spectra

A new approach to recognition of images using invariant features based on higher-order spectra is presented. Higher-order spectra are translation invariant because translation produces linear phase shifts which cancel. Scale and amplification invariance are satisfied by the phase of the integral of a higher-order spectrum along a radial line in higher-order frequency space because the contour of integration maps onto itself and both the real and imaginary parts are affected equally by the transformation. Rotation invariance is introduced by deriving invariants from the Radon transform of the image and using the cyclic-shift invariance property of the discrete Fourier transform magnitude. Results on synthetic and actual images show isolated, compact clusters in feature space and high classification accuracies

An investigation into factors affecting electron density calibration for a megavoltage cone-beam CT system

There is a growing interest in the use of megavoltage cone-beam computed tomography (MV CBCT) data for radiotherapy treatment planning. To calculate accurate dose distributions, knowledge of the electron density (ED) of the tissues being irradiated is required. In the case of MV CBCT, it is necessary to determine a calibration-relating CT number to ED, utilizing the photon beam produced for MV CBCT. A number of different parameters can affect this calibration. This study was undertaken on the Siemens MV CBCT system, MVision, to evaluate the effect of the following parameters on the reconstructed CT pixel value to ED calibration: the number of monitor units (MUs) used (5, 8, 15 and 60 MUs), the image reconstruction filter (head and neck, and pelvis), reconstruction matrix size (256 by 256 and 512 by 512), and the addition of extra solid water surrounding the ED phantom. A Gammex electron density CT phantom containing EDs from 0.292 to 1.707 was imaged under each of these conditions. The linear relationship between MV CBCT pixel value and ED was demonstrated for all MU settings and over the range of EDs. Changes in MU number did not dramatically alter the MV CBCT ED calibration. The use of different reconstruction filters was found to affect the MV CBCT ED calibration, as was the addition of solid water surrounding the phantom. Dose distributions from treatment plans calculated with simulated image data from a 15 MU head and neck reconstruction filter MV CBCT image and a MV CBCT ED calibration curve from the image data parameters and a 15 MU pelvis reconstruction filter showed small and clinically insignificant differences. Thus, the use of a single MV CBCT ED calibration curve is unlikely to result in any clinical differences. However, to ensure minimal uncertainties in dose reporting, MV CBCT ED calibration measurements could be carried out using parameter-specific calibration measurements.

Ultrasound attenuation computed tomography assessment of PAGAT gel dose

Ultrasound has been previously investigated as an alternative readout method for irradiated polymer gel dosimeters, with authors reporting varying dose responses. We extend previous work utilizing a new computed tomography ultrasound scanner comprising of two identical 5 MHz, 128-element linear-array ultrasound transducers, co-axially aligned and submerged in water as a coupling agent, with rotational of the gel dosimeter between the transducers facilitated by a robotic arm. We have investigated the dose-dependence of both ultrasound bulk attenuation and broadband ultrasound attenuation (BUA) for the PAGAT gel dosimeter. The ultrasound bulk attenuation dose sensitivity was found to be 1.46  ±  0.04 dB m −1 Gy −1, being in agreement with previously published results for PAG and MAGIC gels. BUA was also found to be dose dependent and was measured to be 0.024  ±  0.003 dB MHz −1 Gy −1; the advantage of BUA being its insensitivity to frequency-independent attenuation mechanisms including reflection and refraction, thereby minimizing image reconstruction artefacts.

Electron cryo-microscopy studies of bacteriophage phi8 and archaeal virus SH1

Symmetry is a key principle in viral structures, especially the protein capsid shells. However, symmetry mismatches are very common, and often correlate with dynamic functionality of biological significance. The three-dimensional structures of two isometric viruses, bacteriophage phi8 and the archaeal virus SH1 were reconstructed using electron cryo-microscopy. Two image reconstruction methods were used: the classical icosahedral method yielded high resolution models for the symmetrical parts of the structures, and a novel asymmetric in-situ reconstruction method allowed us to resolve the symmetry mismatches at the vertices of the viruses. Evidence was found that the hexameric packaging enzyme at the vertices of phi8 does not rotate relative to the capsid. The large two-fold symmetric spikes of SH1 were found not to be responsible for infectivity. Both virus structures provided insight into the evolution of viruses. Comparison of the phi8 polymerase complex capsid with those of phi6 and other dsRNA viruses suggests that the quaternary structure in dsRNA bacteriophages differs from other dsRNA viruses. SH1 is unusual because there are two major types of capsomers building up the capsid, both of which seem to be composed mainly of single beta-barrels perpendicular to the capsid surface. This indicates that the beta-barrel may be ancestral to the double beta-barrel fold.

Lipid-Containing Icosahedral dsDNA Bacteriophages: Entry, Exit and Structure

In this thesis three icosahedral lipid-containing double-stranded (ds) deoxyribonucleic acid (DNA) bacteriophages have been studied: PRD1, Bam35 and P23-77. The work focuses on the entry, exit and structure of the viruses. PRD1 is the type member of the Tectiviridae family, infecting a variety of Gram-negative bacteria. The PRD1 receptor binding complex, consisting of the penton protein P31, the spike protein P5 and the receptor binding protein P2 recognizes a specific receptor on the host surface. In this study we found that the transmembrane protein P16 has an important stabilization function as the fourth member of the receptor binding complex and protein P16 may have a role in the formation of a tubular membrane structure, which is needed in the ejection of the genome into the cell. Phage Bam35 (Tectiviridae), which infects Gram-positive hosts, has been earlier found to resemble PRD1 in morphology and genome organization The uncharacterized early and late events in the Bam35 life cycle were studied by electrochemical methods. Physiological changes in the beginning of the infection were found to be similar in both lysogenic and nonlysogenic cell lines, Bam35 inducing a temporal decrease of membrane voltage and K+ efflux. At the end of the infection cycle physiological changes were observed only in the nonlysogenic cell line. The strong K+ efflux 40 min after infection and the induced premature cell lysis propose that Bam35 has a similar holin-endolysin lysis system to that of PRD1. Thermophilic icosahedral dsDNA Thermus phages P23-65H, P23-72 and P23-77 have been proposed to belong to the Tectiviridae family. In this study these phages were compared to each other. Analysis of structural protein patterns and stability revealed these phages to be very similar but not identical. The most stable of the studied viruses, P23-77, was further analyzed in more detail. Cryo-electron microscopy and three-dimensional image reconstruction was used to determine the structure of virus to 14 Å resolution. Results of thin layer chromatography for neutral lipids together with analysis of the three dimensional reconstruction of P23-77 virus particle revealed the presence of an internal lipid membrane. The overall capsid architecture of P23-77 is similar to PRD1 and Bam35, but most closely it resembles the structure of the capsid of archaeal virus SH1. This complicates the classification of dsDNA, internal lipid-containing icosahedral viruses.

Singular value decomposition based computationally efficient algorithm for rapid dynamic near-infrared diffuse optical tomography

Purpose: A computationally efficient algorithm (linear iterative type) based on singular value decomposition (SVD) of the Jacobian has been developed that can be used in rapid dynamic near-infrared (NIR) diffuse optical tomography. Methods: Numerical and experimental studies have been conducted to prove the computational efficacy of this SVD-based algorithm over conventional optical image reconstruction algorithms. Results: These studies indicate that the performance of linear iterative algorithms in terms of contrast recovery (quantitation of optical images) is better compared to nonlinear iterative (conventional) algorithms, provided the initial guess is close to the actual solution. The nonlinear algorithms can provide better quality images compared to the linear iterative type algorithms. Moreover, the analytical and numerical equivalence of the SVD-based algorithm to linear iterative algorithms was also established as a part of this work. It is also demonstrated that the SVD-based image reconstruction typically requires O(NN2) operations per iteration, as contrasted with linear and nonlinear iterative methods that, respectively, requir O(NN3) and O(NN6) operations, with ``NN'' being the number of unknown parameters in the optical image reconstruction procedure. Conclusions: This SVD-based computationally efficient algorithm can make the integration of image reconstruction procedure with the data acquisition feasible, in turn making the rapid dynamic NIR tomography viable in the clinic to continuously monitor hemodynamic changes in the tissue pathophysiology.

A FEM-Based Forward Solver for Studying the Forward Problem of Electrical Impedance Tomography (EIT) with A Practical Biological Phantom

A Finite Element Method based forward solver is developed for solving the forward problem of a 2D-Electrical Impedance Tomography. The Method of Weighted Residual technique with a Galerkin approach is used for the FEM formulation of EIT forward problem. The algorithm is written in MatLAB7.0 and the forward problem is studied with a practical biological phantom developed. EIT governing equation is numerically solved to calculate the surface potentials at the phantom boundary for a uniform conductivity. An EIT-phantom is developed with an array of 16 electrodes placed on the inner surface of the phantom tank filled with KCl solution. A sinusoidal current is injected through the current electrodes and the differential potentials across the voltage electrodes are measured. Measured data is compared with the differential potential calculated for known current and solution conductivity. Comparing measured voltage with the calculated data it is attempted to find the sources of errors to improve data quality for better image reconstruction.

A Simple Instrumentation Calibration Technique for Electrical Impedance Tomography (EIT) Using A 16-Electrode Phantom

A simple analog instrumentation for Electrical Impedance Tomography is developed and calibrated using the practical phantoms. A constant current injector consisting of a modified Howland voltage controlled current source fed by a voltage controlled oscillator is developed to inject a constant current to the phantom boundary. An instrumentation amplifier, 50 Hz notch filter and a narrow band pass filter are developed and used for signal conditioning. Practical biological phantoms are developed and the forward problem is studied to calibrate the EIT-instrumentation. An array of sixteen stainless steel electrodes is developed and placed inside the phantom tank filled with KCl solution. 1 mA, 50 kHz sinusoidal current is injected at the phantom boundary using adjacent current injection protocol. The differential potentials developed at the voltage electrodes are measured for sixteen current injections. Differential voltage signal is passed through an instrumentation amplifier and a filtering block and measured by a digital multimeter. A forward solver is developed using Finite Element Method in MATLAB7.0 for solving the EIT governing equation. Differential potentials are numerically calculated using the forward solver with a simulated current and bathing solution conductivity. Measured potential data is compared with the differential potentials calculated for calibrating the instrumentation to acquire the voltage data suitable for better image reconstruction.

Structural studies on members of the picornavirus superfamily

The pathogenic members of the picornavirus superfamily have adverse effects on humans, their crops and their livestock. As structure is related to function, detailed structural studies on these viruses are important not only for fundamental understanding of the viral life cycle, but also for the rational design of vaccines and inhibitors for disease control. These viruses have positive sense, single-stranded RNA genomes enclosed in a protein capsid. X-ray crystallography and cryo-electron microscopy studies have revealed that the isometric members of this group have icosahedrally-symmetric capsids made up of 60 copies of each of the structural proteins. The members that infect animal cells often employ one or more cellular receptors to facilitate cell entry which in some cases is known to initiate the uncoating sequence of the genome. The nature of the interactions between individual viruses and alternative cellular receptors has rarely been probed. The capsid assembly of the members of the picornavirus superfamily is considered to be cooperative and the interactions of RNA and capsid proteins are thought to play an important role in orchestrating virus assembly. The major aims of this thesis were to solve the structures of blackcurrant reversion virus (BRV), human parechovirus 1 (HPEV1) and coxsackievirus A7 (CAV7), as well as the structure of HPEV1 complexed with two of its cellular receptors using cryo-electron microscopy, three-dimensional image reconstruction and homology modeling. Each of the selected viruses represents a taxonomic group where little or no structural data was previously available. The results enabled the detailed comparison of the new structures to those of known picornaviruses, the identification of surface-exposed epitopes potentially important for host interaction, the mapping of RNA-capsid protein interactions and the elucidation of the basis for the specificity of two different receptor molecules for the same capsid. This work will form the basis for further studies on the influence of RNA on parechovirus assembly as a potential target for drug design.

Elimination of cross-talk in diffuse optical tomographic reconstructions through a weighted update procedure: results of simulation study - art. no. 61640S

Reconstructions in optical tomography involve obtaining the images of absorption and reduced scattering coefficients. The integrated intensity data has greater sensitivity to absorption coefficient variations than scattering coefficient. However, the sensitivity of intensity data to scattering coefficient is not zero. We considered an object with two inhomogeneities (one in absorption and the other in scattering coefficient). The standard iterative reconstruction techniques produced results, which were plagued by cross talk, i.e., the absorption coefficient reconstruction has a false positive corresponding to the location of scattering inhomogeneity, and vice-versa. We present a method to remove cross talk in the reconstruction, by generating a weight matrix and weighting the update vector during the iteration. The weight matrix is created by the following method: we first perform a simple backprojection of the difference between the experimental and corresponding homogeneous intensity data. The built up image has greater weightage towards absorption inhomogeneity than the scattering inhomogeneity and its appropriate inverse is weighted towards the scattering inhomogeneity. These two weight matrices are used as multiplication factors in the update vectors, normalized backprojected image of difference intensity for absorption inhomogeneity and the inverse of the above for the scattering inhomogeneity, during the image reconstruction procedure. We demonstrate through numerical simulations, that cross-talk is fully eliminated through this modified reconstruction procedure.

Implementation of a phase array diffuse optical tomographic imager

Diffuse optical tomography (DOT) using near-infrared (NIR) light is a promising tool for noninvasive imaging of deep tissue. This technique is capable of quantitative reconstructions of absorption coefficient inhomogeneities of tissue. The motivation for reconstructing the optical property variation is that it, and, in particular, the absorption coefficient variation, can be used to diagnose different metabolic and disease states of tissue. In DOT, like any other medical imaging modality, the aim is to produce a reconstruction with good spatial resolution and accuracy from noisy measurements. We study the performance of a phase array system for detection of optical inhomogeneities in tissue. The light transport through a tissue is diffusive in nature and can be modeled using diffusion equation if the optical parameters of the inhomogeneity are close to the optical properties of the background. The amplitude cancellation method that uses dual out-of-phase sources (phase array) can detect and locate small objects in turbid medium. The inverse problem is solved using model based iterative image reconstruction. Diffusion equation is solved using finite element method for providing the forward model for photon transport. The solution of the forward problem is used for computing the Jacobian and the simultaneous equation is solved using conjugate gradient search. The simulation studies have been carried out and the results show that a phase array system can resolve inhomogeneities with sizes of 5 mm when the absorption coefficient of the inhomogeneity is twice that of the background tissue. To validate this result, a prototype model for performing a dual-source system has been developed. Experiments are carried out by inserting an inhomogeneity of high optical absorption coefficient in an otherwise homogeneous phantom while keeping the scattering coefficient same. The high frequency (100 MHz) modulated dual out-of-phase laser source light is propagated through the phantom. The interference of these sources creates an amplitude null and a phase shift of 180° along a plane between the two sources with a homogeneous object. A solid resin phantom with inhomogeneities simulating the tumor is used in our experiment. The amplitude and phase changes are found to be disturbed by the presence of the inhomogeneity in the object. The experimental data (amplitude and the phase measured at the detector) are used for reconstruction. The results show that the method is able to detect multiple inhomogeneities with sizes of 4 mm. The localization error for a 5 mm inhomogeneity is found to be approximately 1 mm.

Studying the resistivity imaging of chicken tissue phantoms with different current patterns in Electrical Impedance Tomography (EIT)

A current injection pattern in Electrical Impedance Tomography (EIT) has its own current distribution profile within the domain under test. Hence, different current patterns have different sensitivity, spatial resolution and distinguishability. Image reconstruction studies with practical phantoms are essential to assess the performance of EIT systems for their validation, calibration and comparison purposes. Impedance imaging of real tissue phantoms with different current injection methods is also essential for better assessment of the biomedical EIT systems. Chicken tissue paste phantoms and chicken tissue block phantoms are developed and the resistivity image reconstruction is studied with different current injection methods. A 16-electrode array is placed inside the phantom tank and the tank is filled with chicken muscle tissue paste or chicken tissue blocks as the background mediums. Chicken fat tissue, chicken bone, air hole and nylon cylinders are used as the inhomogeneity to obtained different phantom configurations. A low magnitude low frequency constant sinusoidal current is injected at the phantom boundary with opposite and neighboring current patterns and the boundary potentials are measured. Resistivity images are reconstructed from the boundary data using EIDORS and the reconstructed images are analyzed with the contrast parameters calculated from their elemental resistivity profiles. Results show that the resistivity profiles of all the phantom domains are successfully reconstructed with a proper background resistivity and high inhomogeneity resistivity for both the current injection methods. Reconstructed images show that, for all the chicken tissue phantoms, the inhomogeneities are suitably reconstructed with both the current injection protocols though the chicken tissue block phantom and opposite method are found more suitable. It is observed that the boundary potentials of the chicken tissue block phantoms are higher than the chicken tissue paste phantom. SNR of the chicken tissue block phantoms are found comparatively more and hence the chicken tissue block phantom is found more suitable for its lower noise performance. The background noise is found less in opposite method for all the phantom configurations which yields the better resistivity images with high PCR and COC and proper IRMean and IRMax neighboring method showed higher noise level for both the chicken tissue paste phantoms and chicken tissue block phantoms with all the inhomogeneities. Opposite method is found more suitable for both the chicken tissue phantoms, and also, chicken tissue block phantoms are found more suitable compared to the chicken tissue paste phantom. (C) 2012 Elsevier Ltd. All rights reserved.

Mesh Simplification Based on Edge Collapsing Could Improve Computational Efficiency in Near Infrared Optical Tomographic Imaging

The diffusion equation-based modeling of near infrared light propagation in tissue is achieved by using finite-element mesh for imaging real-tissue types, such as breast and brain. The finite-element mesh size (number of nodes) dictates the parameter space in the optical tomographic imaging. Most commonly used finite-element meshing algorithms do not provide the flexibility of distinct nodal spacing in different regions of imaging domain to take the sensitivity of the problem into consideration. This study aims to present a computationally efficient mesh simplification method that can be used as a preprocessing step to iterative image reconstruction, where the finite-element mesh is simplified by using an edge collapsing algorithm to reduce the parameter space at regions where the sensitivity of the problem is relatively low. It is shown, using simulations and experimental phantom data for simple meshes/domains, that a significant reduction in parameter space could be achieved without compromising on the reconstructed image quality. The maximum errors observed by using the simplified meshes were less than 0.27% in the forward problem and 5% for inverse problem.

Data-resolution based optimization of the data-collection strategy for near infrared diffuse optical tomography

Purpose: To optimize the data-collection strategy for diffuse optical tomography and to obtain a set of independent measurements among the total measurements using the model based data-resolution matrix characteristics. Methods: The data-resolution matrix is computed based on the sensitivity matrix and the regularization scheme used in the reconstruction procedure by matching the predicted data with the actual one. The diagonal values of data-resolution matrix show the importance of a particular measurement and the magnitude of off-diagonal entries shows the dependence among measurements. Based on the closeness of diagonal value magnitude to off-diagonal entries, the independent measurements choice is made. The reconstruction results obtained using all measurements were compared to the ones obtained using only independent measurements in both numerical and experimental phantom cases. The traditional singular value analysis was also performed to compare the results obtained using the proposed method. Results: The results indicate that choosing only independent measurements based on data-resolution matrix characteristics for the image reconstruction does not compromise the reconstructed image quality significantly, in turn reduces the data-collection time associated with the procedure. When the same number of measurements (equivalent to independent ones) are chosen at random, the reconstruction results were having poor quality with major boundary artifacts. The number of independent measurements obtained using data-resolution matrix analysis is much higher compared to that obtained using the singular value analysis. Conclusions: The data-resolution matrix analysis is able to provide the high level of optimization needed for effective data-collection in diffuse optical imaging. The analysis itself is independent of noise characteristics in the data, resulting in an universal framework to characterize and optimize a given data-collection strategy. (C) 2012 American Association of Physicists in Medicine. http://dx.doi.org/10.1118/1.4736820]