RESUMO
An increasingly common viewpoint is that protein dynamics datasets reside in a nonlinear subspace of low conformational energy. Ideal data analysis tools should therefore account for such nonlinear geometry. The Riemannian geometry setting can be suitable for a variety of reasons. First, it comes with a rich mathematical structure to account for a wide range of geometries that can be modeled after an energy landscape. Second, many standard data analysis tools developed for data in Euclidean space can be generalized to Riemannian manifolds. In the context of protein dynamics, a conceptual challenge comes from the lack of guidelines for constructing a smooth Riemannian structure based on an energy landscape. In addition, computational feasibility in computing geodesics and related mappings poses a major challenge. This work considers these challenges. The first part of the paper develops a local approximation technique for computing geodesics and related mappings on Riemannian manifolds in a computationally feasible manner. The second part constructs a smooth manifold and a Riemannian structure that is based on an energy landscape for protein conformations. The resulting Riemannian geometry is tested on several data analysis tasks relevant for protein dynamics data. In particular, the geodesics with given start- and end-points approximately recover corresponding molecular dynamics trajectories for proteins that undergo relatively ordered transitions with medium-sized deformations. The Riemannian protein geometry also gives physically realistic summary statistics and retrieves the underlying dimension even for large-sized deformations within seconds on a laptop.
Assuntos
Conformação Proteica , Proteínas , Proteínas/química , Algoritmos , Simulação de Dinâmica MolecularRESUMO
We introduce a reconstruction framework that can account for shape related a priori information in ill-posed linear inverse problems in imaging. It is a variational scheme that uses a shape functional defined using deformable templates machinery from shape theory. As proof of concept, we apply the proposed shape based reconstruction to 2D tomography with very sparse measurements, and demonstrate strong empirical results.
RESUMO
Accurate modeling of image formation in cryo-electron microscopy is an important requirement for quantitative image interpretation and optimization of the data acquisition strategy. Here we present a forward model that accounts for the specimen's scattering properties, microscope optics, and detector response. The specimen interaction potential is calculated with the isolated atom superposition approximation (IASA) and extended with the influences of solvent's dielectric and ionic properties as well as the molecular electrostatic distribution. We account for an effective charge redistribution via the Poisson-Boltzmann approach and find that the IASA-based potential forms the dominant part of the interaction potential, as the contribution of the redistribution is less than 10%. The electron wave is propagated through the specimen by a multislice approach and the influence of the optics is included via the contrast transfer function. We incorporate the detective quantum efficiency of the camera due to the difference between signal and noise transfer characteristics, instead of using only the modulation transfer function. The full model was validated against experimental images of 20S proteasome, hemoglobin, and GroEL. The simulations adequately predict the effects of phase contrast, changes due to the integrated electron flux, thickness, inelastic scattering, detective quantum efficiency and acceleration voltage. We suggest that beam-induced specimen movements are relevant in the experiments whereas the influence of the solvent amorphousness can be neglected. All simulation parameters are based on physical principles and, when necessary, experimentally determined.
Assuntos
Chaperonina 60/ultraestrutura , Microscopia Crioeletrônica/métodos , Hemoglobinas/ultraestrutura , Modelos Moleculares , Complexo de Endopeptidases do Proteassoma/ultraestrutura , Chaperonina 60/química , Hemoglobinas/química , Processamento de Imagem Assistida por Computador , Distribuição de Poisson , Complexo de Endopeptidases do Proteassoma/química , Software , Eletricidade EstáticaRESUMO
Uncertainty quantification in automated image analysis is highly desired in many applications. Typically, machine learning models in classification or segmentation are only developed to provide binary answers; however, quantifying the uncertainty of the models can play a critical role for example in active learning or machine human interaction. Uncertainty quantification is especially difficult when using deep learning-based models, which are the state-of-the-art in many imaging applications. The current uncertainty quantification approaches do not scale well in high-dimensional real-world problems. Scalable solutions often rely on classical techniques, such as dropout, during inference or training ensembles of identical models with different random seeds to obtain a posterior distribution. In this paper, we present the following contributions. First, we show that the classical approaches fail to approximate the classification probability. Second, we propose a scalable and intuitive framework for uncertainty quantification in medical image segmentation that yields measurements that approximate the classification probability. Third, we suggest the usage of k-fold cross-validation to overcome the need for held out calibration data. Lastly, we motivate the adoption of our method in active learning, creating pseudo-labels to learn from unlabeled images and human-machine collaboration.
Assuntos
Aprendizado Profundo , Humanos , Incerteza , Probabilidade , Calibragem , Processamento de Imagem Assistida por ComputadorRESUMO
Artificial intelligence (AI) methods applied to healthcare problems have shown enormous potential to alleviate the burden of health services worldwide and to improve the accuracy and reproducibility of predictions. In particular, developments in computer vision are creating a paradigm shift in the analysis of radiological images, where AI tools are already capable of automatically detecting and precisely delineating tumours. However, such tools are generally developed in technical departments that continue to be siloed from where the real benefit would be achieved with their usage. Significant effort still needs to be made to make these advancements available, first in academic clinical research and ultimately in the clinical setting. In this paper, we demonstrate a prototype pipeline based entirely on open-source software and free of cost to bridge this gap, simplifying the integration of tools and models developed within the AI community into the clinical research setting, ensuring an accessible platform with visualisation applications that allow end-users such as radiologists to view and interact with the outcome of these AI tools.
RESUMO
PURPOSE: To determine if pelvic/ovarian and omental lesions of ovarian cancer can be reliably segmented on computed tomography (CT) using fully automated deep learning-based methods. METHODS: A deep learning model for the two most common disease sites of high-grade serous ovarian cancer lesions (pelvis/ovaries and omentum) was developed and compared against the well-established "no-new-Net" framework and unrevised trainee radiologist segmentations. A total of 451 CT scans collected from four different institutions were used for training (n = 276), evaluation (n = 104) and testing (n = 71) of the methods. The performance was evaluated using the Dice similarity coefficient (DSC) and compared using a Wilcoxon test. RESULTS: Our model outperformed no-new-Net for the pelvic/ovarian lesions in cross-validation, on the evaluation and test set by a significant margin (p values being 4 × 10-7, 3 × 10-4, 4 × 10-2, respectively), and for the omental lesions on the evaluation set (p = 1 × 10-3). Our model did not perform significantly differently in segmenting pelvic/ovarian lesions (p = 0.371) compared to a trainee radiologist. On an independent test set, the model achieved a DSC performance of 71 ± 20 (mean ± standard deviation) for pelvic/ovarian and 61 ± 24 for omental lesions. CONCLUSION: Automated ovarian cancer segmentation on CT scans using deep neural networks is feasible and achieves performance close to a trainee-level radiologist for pelvic/ovarian lesions. RELEVANCE STATEMENT: Automated segmentation of ovarian cancer may be used by clinicians for CT-based volumetric assessments and researchers for building complex analysis pipelines. KEY POINTS: ⢠The first automated approach for pelvic/ovarian and omental ovarian cancer lesion segmentation on CT images has been presented. ⢠Automated segmentation of ovarian cancer lesions can be comparable with manual segmentation of trainee radiologists. ⢠Careful hyperparameter tuning can provide models significantly outperforming strong state-of-the-art baselines.
Assuntos
Aprendizado Profundo , Cistos Ovarianos , Neoplasias Ovarianas , Humanos , Feminino , Neoplasias Ovarianas/diagnóstico por imagem , Redes Neurais de Computação , Tomografia Computadorizada por Raios XRESUMO
Filtered back-projection and weighted back-projection have long been the methods of choice within the electron microscopy community for reconstructing the structure of macromolecular assemblies from electron tomography data. Here, we describe electron lambda-tomography, a reconstruction method that enjoys the benefits of the above mentioned methods, namely speed and ease of implementation, but also addresses some of their shortcomings. In particular, compared to these standard methods, electron lambda-tomography is less sensitive to artifacts that come from structures outside the region that is being reconstructed, and it can sharpen boundaries.
RESUMO
Three-dimensional reconstruction of the electron-scattering potential of biological macromolecules from electron cryo-microscopy (cryo-EM) projection images is an ill-posed problem. The most popular cryo-EM software solutions to date rely on a regularization approach that is based on the prior assumption that the scattering potential varies smoothly over three-dimensional space. Although this approach has been hugely successful in recent years, the amount of prior knowledge that it exploits compares unfavorably with the knowledge about biological structures that has been accumulated over decades of research in structural biology. Here, a regularization framework for cryo-EM structure determination is presented that exploits prior knowledge about biological structures through a convolutional neural network that is trained on known macromolecular structures. This neural network is inserted into the iterative cryo-EM structure-determination process through an approach that is inspired by regularization by denoising. It is shown that the new regularization approach yields better reconstructions than the current state of the art for simulated data, and options to extend this work for application to experimental cryo-EM data are discussed.
RESUMO
Model-based learned iterative reconstruction methods have recently been shown to outperform classical reconstruction algorithms. Applicability of these methods to large scale inverse problems is however limited by the available memory for training and extensive training times, the latter due to computationally expensive forward models. As a possible solution to these restrictions we propose a multi-scale learned iterative reconstruction scheme that computes iterates on discretisations of increasing resolution. This procedure does not only reduce memory requirements, it also considerably speeds up reconstruction and training times, but most importantly is scalable to large scale inverse problems with non-trivial forward operators, such as those that arise in many 3D tomographic applications. In particular, we propose a hybrid network that combines the multiscale iterative approach with a particularly expressive network architecture which in combination exhibits excellent scalability in 3D. Applicability of the algorithm is demonstrated for 3D cone beam computed tomography from real measurement data of an organic phantom. Additionally, we examine scalability and reconstruction quality in comparison to established learned reconstruction methods in two dimensions for low dose computed tomography on human phantoms.
RESUMO
Semantic edge detection has recently gained a lot of attention as an image-processing task, mainly because of its wide range of real-world applications. This is based on the fact that edges in images contain most of the semantic information. Semantic edge detection involves two tasks, namely pure edge detection and edge classification. Those are in fact fundamentally distinct in terms of the level of abstraction that each task requires. This fact is known as the distracted supervision paradox and limits the possible performance of a supervised model in semantic edge detection. In this work, we will present a novel hybrid method that is based on a combination of the model-based concept of shearlets, which provides probably optimally sparse approximations of a model class of images, and the data-driven method of a suitably designed convolutional neural network. We show that it avoids the distracted supervision paradox and achieves high performance in semantic edge detection. In addition, our approach requires significantly fewer parameters than a pure data-driven approach. Finally, we present several applications such as tomographic reconstruction and show that our approach significantly outperforms former methods, thereby also indicating the value of such hybrid methods for biomedical imaging.
RESUMO
The task of separating an image into distinct components that represent different features plays an important role in many applications. Traditionally, such separation techniques are applied once the image in question has been reconstructed from measured data. We propose an efficient iterative algorithm, where reconstruction is performed jointly with the task of separation. A key assumption is that the image components have different sparse representations. The algorithm is based on a scheme that minimizes a functional composed of a data discrepancy term and the l1 -norm of the coefficients of the different components with respect to their corresponding dictionaries. The performance is demonstrated for joint 2D deconvolution and separation into curve- and point-like components, and tests are performed on synthetic data as well as experimental stimulated emission depletion and confocal microscopy data. Experiments show that such a joint approach outperforms a sequential approach, where one first deconvolves data and then applies image separation.
RESUMO
Protein contents in platelets are frequently changed upon tumor development and metastasis. However, how cancer cells can influence protein-selective redistribution and release within platelets, thereby promoting tumor development, remains largely elusive. With fluorescence-based super-resolution stimulated emission depletion (STED) imaging we reveal how specific proteins, implicated in tumor progression and metastasis, re-distribute within platelets, when subject to soluble activators (thrombin, adenosine diphosphate and thromboxane A2), and when incubated with cancer (MCF-7, MDA-MB-231, EFO21) or non-cancer cells (184A1, MCF10A). Upon cancer cell incubation, the cell-adhesion protein P-selectin was found to re-distribute into circular nano-structures, consistent with accumulation into the membrane of protein-storing alpha-granules within the platelets. These changes were to a significantly lesser extent, if at all, found in platelets incubated with normal cells, or in platelets subject to soluble platelet activators. From these patterns, we developed a classification procedure, whereby platelets exposed to cancer cells, to non-cancer cells, soluble activators, as well as non-activated platelets all could be identified in an automatic, objective manner. We demonstrate that STED imaging, in contrast to electron and confocal microscopy, has the necessary spatial resolution and labelling efficiency to identify protein distribution patterns in platelets and can resolve how they specifically change upon different activations. Combined with image analyses of specific protein distribution patterns within the platelets, STED imaging can thus have a role in future platelet-based cancer diagnostics and therapeutic monitoring. The presented approach can also bring further clarity into fundamental mechanisms for cancer cell-platelet interactions, and into non-contact cell-to-cell interactions in general.
Assuntos
Plaquetas/metabolismo , Microscopia de Fluorescência , Plaquetas/citologia , Plaquetas/ultraestrutura , Linhagem Celular Tumoral , Técnicas de Cocultura , Fibrinogênio/química , Fibrinogênio/metabolismo , Proteínas de Choque Térmico/química , Proteínas de Choque Térmico/metabolismo , Humanos , Nanoestruturas/química , Selectina-P/química , Selectina-P/metabolismo , Fator A de Crescimento do Endotélio Vascular/química , Fator A de Crescimento do Endotélio Vascular/metabolismoRESUMO
We propose the Learned Primal-Dual algorithm for tomographic reconstruction. The algorithm accounts for a (possibly non-linear) forward operator in a deep neural network by unrolling a proximal primal-dual optimization method, but where the proximal operators have been replaced with convolutional neural networks. The algorithm is trained end-to-end, working directly from raw measured data and it does not depend on any initial reconstruction such as filtered back-projection (FBP). We compare performance of the proposed method on low dose computed tomography reconstruction against FBP, total variation (TV), and deep learning based post-processing of FBP. For the Shepp-Logan phantom we obtain >6 dB peak signal to noise ratio improvement against all compared methods. For human phantoms the corresponding improvement is 6.6 dB over TV and 2.2 dB over learned post-processing along with a substantial improvement in the structural similarity index. Finally, our algorithm involves only ten forward-back-projection computations, making the method feasible for time critical clinical applications.
Assuntos
Algoritmos , Aprendizado Profundo , Intensificação de Imagem Radiográfica/métodos , Humanos , Imagens de Fantasmas , Tomografia Computadorizada por Raios XRESUMO
A novel device that can be used as a tunable support-free phase plate for transmission electron microscopy of weakly scattering specimens is described. The device relies on the generation of a controlled phase shift by the magnetic field of a segment of current-carrying wire that is oriented parallel or antiparallel to the electron beam. The validity of the concept is established using both experimental electron holographic measurements and a theoretical model based on Ampere's law. Computer simulations are used to illustrate the resulting contrast enhancement for studies of biological cells and macromolecules.
RESUMO
We introduce a tomographic reconstruction method implemented using a shape-based regularization technique. Spatial models of known features in the structure being reconstructed are integrated into the reconstruction process as regularizers. Our regularization scheme is driven locally through shape information obtained from segmentation and compared with a known spatial model. We demonstrated our method on tomography data from digital phantoms, simulated data, and experimental electron tomography (ET) data of virus complexes. Our reconstruction showed reduced blurring and an improvement in the resolution of the reconstructed volume was also measured. This method also produced improved demarcation of spike boundaries in viral membranes when compared with popular techniques like weighted back projection and the algebraic reconstruction technique. Improved ET reconstructions will provide better structure elucidation and improved feature visualization, which can aid in solving key biological issues. Our method can also be generalized to other tomographic modalities.