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
-
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
-
Abstract The Thai Meteorological Department (TMD) seismic network began development in 2008. There are a total of 71 seismic stations consisting of 26 borehole stations and 45 surface stations currently installed. The three-component data from the TMD seismic network have been widely used in previous seismological studies. In a recent analysis, we have found that sensor orientation as reported in the site metadata is sometimes significantly incorrect, especially for borehole stations. In this study, we analyze P-wave polarization data from regional and teleseismic earthquakes recorded in the network to estimate the true instrument orientation relative to geographic north and compare that to station metadata. Of the 45 surface stations, we found that at present, ~ 82% are well oriented (i.e., aligned within 0–15° of true north). However, 8 sites have sensors misoriented by more than 15°, and some stations had a temporal change in sensor orientation during an upgrade to the seismic system with replacement of the sensor. We also evaluated sensor orientations for 26 TMD borehole seismic stations, from 2018 to the 2022. For many of the borehole stations, the actual sensor orientation differs significantly from the TMD metadata, especially at short-period stations. Many of those stations have sensor misorientations approaching 180°, due to errors in the ambient noise analysis calibration techniques used during installation. We have also investigated how this sensor misorientation affects previous seismic studies, such as regional moment tensor inversion of earthquakes sources and receiver function stacking. We have found that the large deviations in sensor orientation can result in erroneous results and/or large measurement errors. A cause of the orientation error for borehole sites could be a combination of strong background surface ambient seismic noise coupled with an incorrect reference instrument response.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Conference Paper:
- https://doi.org/10.1117/12.2543813
