An official website of the United States government
Simultaneous measurements of single-molecule positions and orientations provide critical insight into a variety of biological and chemical processes. Various engineered point spread functions (PSFs) have been introduced for measuring the orientation and rotational diffusion of dipole-like emitters, but the widely used Cramér-Rao bound (CRB) only evaluates performance for one specific orientation at a time. Here, we report a performance metric, termed variance upper bound (VUB), that yields a global maximum CRB for all possible molecular orientations, thereby enabling the measurement performance of any PSF to be computed efficiently ( faster than calculating average CRB). Our VUB reveals that the simple polarized standard PSF provides robust and precise orientation measurements if emitters are near a refractive index interface. Using this PSF, we measure the orientations and positions of Nile red (NR) molecules transiently bound to amyloid aggregates. Our super-resolved images reveal the main binding mode of NR on amyloid fiber surfaces, as well as structural heterogeneities along amyloid fibrillar networks, that cannot be resolved by single-molecule localization alone.
more »
« less
A computationally-efficient bound for the variance of measuring the orientation of single molecules
Modulating the polarization of excitation light, resolving the polarization of emitted fluorescence, and point spread function (PSF) engineering have been widely leveraged for measuring the orientation of single molecules. Typically, the performance of these techniques is optimized and quantified using the Cramér-Rao bound (CRB), which describes the best possible measurement variance of an unbiased estimator. However, CRB is a local measure and requires exhaustive sampling across the measurement space to fully characterize measurement precision. We develop a global variance upper bound (VUB) for fast quantification and comparison of orientation measurement techniques. Our VUB tightly bounds the diagonal elements of the CRB matrix from above; VUB overestimates the mean CRB by ~34%. However, compared to directly calculating the mean CRB over orientation space, we are able to calculate VUB ~1000 times faster.
more »
« less
- Award ID(s):
- 1653777
- PAR ID:
- 10135046
- Date Published:
- Journal Name:
- Proc. SPIE
- Volume:
- 11246
- Page Range / eLocation ID:
- 1124616
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
PurposeTo improve the performance of neural networks for parameter estimation in quantitative MRI, in particular when the noise propagation varies throughout the space of biophysical parameters. Theory and MethodsA theoretically well‐founded loss function is proposed that normalizes the squared error of each estimate with respective Cramér–Rao bound (CRB)—a theoretical lower bound for the variance of an unbiased estimator. This avoids a dominance of hard‐to‐estimate parameters and areas in parameter space, which are often of little interest. The normalization with corresponding CRB balances the large errors of fundamentally more noisy estimates and the small errors of fundamentally less noisy estimates, allowing the network to better learn to estimate the latter. Further, proposed loss function provides an absolute evaluation metric for performance: A network has an average loss of 1 if it is a maximally efficient unbiased estimator, which can be considered the ideal performance. The performance gain with proposed loss function is demonstrated at the example of an eight‐parameter magnetization transfer model that is fitted to phantom and in vivo data. ResultsNetworks trained with proposed loss function perform close to optimal, that is, their loss converges to approximately 1, and their performance is superior to networks trained with the standard mean‐squared error (MSE). The proposed loss function reduces the bias of the estimates compared to the MSE loss, and improves the match of the noise variance to the CRB. This performance gain translates to in vivo maps that align better with the literature. ConclusionNormalizing the squared error with the CRB during the training of neural networks improves their performance in estimating biophysical parameters.more » « less
-
Collagen quantity and integrity play an important role in understanding diseases such as myelofibrosis (MF). Label-free mid-infrared spectroscopic imaging (MIRSI) has the potential to quantify collagen while minimizing the subjective variance observed with conventional histopathology. Infrared (IR) spectroscopy with polarization sensitivity provides chemical information while also estimating tissue dichroism. This can potentially aid MF grading by revealing the structure and orientation of collagen fibers. Simultaneous measurement of collagen structure and biochemical properties can translate clinically into improved diagnosis and enhance our understanding of disease progression. In this paper, we present the first report of polarization-dependent spectroscopic variations in collagen from human bone marrow samples. We build on prior work with animal models and extend it to human clinical biopsies with a practical method for high-resolution chemical and structural imaging of bone marrow on clinical glass slides. This is done using a new polarization-sensitive photothermal mid-infrared spectroscopic imaging scheme that enables sample and source independent polarization control. This technology provides 0.5 µm spatial resolution, enabling the identification of thin (≈1 µm) collagen fibers that were not separable using Fourier Transform Infrared (FT-IR) imaging in the fingerprint region at diffraction-limited resolution ( ≈ 5 µm). Finally, we propose quantitative metrics to identify fiber orientation from discrete band images (amide I and amide II) measured under three polarizations. Previous studies have used a pair of orthogonal polarization measurements, which is insufficient for clinical samples since human bone biopsies contain collagen fibers with multiple orientations. Here, we address this challenge and demonstrate that three polarization measurements are necessary to resolve orientation ambiguity in clinical bone marrow samples. This is also the first study to demonstrate the ability to spectroscopically identify thin collagen fibers (≈1 µm diameter) and their orientations, which is critical for accurate grading of human bone marrow fibrosis.more » « less
-
Abstract Polarization observations of the Milky Way and many other spiral galaxies have found a close correspondence between the orientation of spiral arms and magnetic field lines on scales of hundreds of parsecs. This paper presents polarization measurements at 214μm toward 10 filamentary candidate “bones” in the Milky Way using the High-resolution Airborne Wide-band Camera on the Stratospheric Observatory for Infrared Astronomy. These data were taken as part of the Filaments Extremely Long and Dark: A Magnetic Polarization Survey and represent the first study to resolve the magnetic field in spiral arms at parsec scales. We describe the complex yet well-defined polarization structure of all 10 candidate bones, and we find a mean difference and standard deviation of −74° ± 32° between their filament axis and the plane-of-sky magnetic field, closer to a field perpendicular to their length rather than parallel. By contrast, the 850μm polarization data from Planck on scales greater than 10 pc show a nearly parallel mean difference of 3° ± 21°. These findings provide further evidence that magnetic fields can change orientation at the scale of dense molecular clouds, even along spiral arms. Finally, we use a power law to fit the dust polarization fraction as a function of total intensity on a cloud-by-cloud basis and find indices between −0.6 and −0.9, with a mean and standard deviation of −0.7 ± 0.1. The polarization, dust temperature, and column density data presented in this work are publicly available online.more » « less
-
In this paper, we study the geometry of quadratic covariance bounds on the estimation error covariance, in a properly defined Hilbert space of random variables. We show that a lower bound on the error covariance may be represented by the Grammian of the error score after projection onto the subspace spanned by the measurement scores. The Grammian is defined with respect to inner products in a Hilbert space of second order random variables. This geometric result holds for a large class of quadratic covariance bounds including the Barankin, Cramer-Rao, and Bhattacharyya bounds, where each bound is characterized by its corresponding measurement scores. When parameters consist of essential parameters and nuisance parameters, the Cram´er-Rao covariance bound is the inverse of the Grammian of essential scores after projection onto the subspace orthogonal to the subspace spanned by the nuisance scores. In two examples, we show that for complex multivariate normal measurements with parameterized mean or covariance, there exist well-known Euclidean space geometries for the general Hilbert space geometry derived in this paper.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1117/12.2543813
