Missing Wedge Completion via Unsupervised Learning with Coordinate Networks.

Cryogenic electron tomography (cryoET) is a powerful tool in structural biology, enabling detailed 3D imaging of biological specimens at a resolution of nanometers. Despite its potential, cryoET faces challenges such as the missing wedge problem, which limits reconstruction quality due to incomplete data collection angles. Recently, supervised deep learning methods leveraging convolutional neural networks (CNNs) have considerably addressed this issue; however, their pretraining requirements render them susceptible to inaccuracies and artifacts, particularly when representative training data is scarce. To overcome these limitations, we introduce a proof-of-concept unsupervised learning approach using coordinate networks (CNs) that optimizes network weights directly against input projections. This eliminates the need for pretraining, reducing reconstruction runtime by 3-20× compared to supervised methods. Our in silico results show improved shape completion and reduction of missing wedge artifacts, assessed through several voxel-based image quality metrics in real space and a novel directional Fourier Shell Correlation (FSC) metric. Our study illuminates benefits and considerations of both supervised and unsupervised approaches, guiding the development of improved reconstruction strategies.

Addressing preferred specimen orientation in single-particle cryo-EM through tilting.

We present a strategy for tackling preferred specimen orientation in single-particle cryogenic electron microscopy by employing tilts during data collection. We also describe a tool to quantify the resulting directional resolution using 3D Fourier shell correlation volumes. We applied these methods to determine the structures at near-atomic resolution of the influenza hemagglutinin trimer, which adopts a highly preferred specimen orientation, and of ribosomal biogenesis intermediates, which adopt moderately preferred orientations.

Belief about nicotine selectively modulates value and reward prediction error signals in smokers.

Little is known about how prior beliefs impact biophysically described processes in the presence of neuroactive drugs, which presents a profound challenge to the understanding of the mechanisms and treatments of addiction. We engineered smokers' prior beliefs about the presence of nicotine in a cigarette smoked before a functional magnetic resonance imaging session where subjects carried out a sequential choice task. Using a model-based approach, we show that smokers' beliefs about nicotine specifically modulated learning signals (value and reward prediction error) defined by a computational model of mesolimbic dopamine systems. Belief of "no nicotine in cigarette" (compared with "nicotine in cigarette") strongly diminished neural responses in the striatum to value and reward prediction errors and reduced the impact of both on smokers' choices. These effects of belief could not be explained by global changes in visual attention and were specific to value and reward prediction errors. Thus, by modulating the expression of computationally explicit signals important for valuation and choice, beliefs can override the physical presence of a potent neuroactive compound like nicotine. These selective effects of belief demonstrate that belief can modulate model-based parameters important for learning. The implications of these findings may be far ranging because belief-dependent effects on learning signals could impact a host of other behaviors in addiction as well as in other mental health problems.

Genome-wide imaging association study implicates functional activity and glial homeostasis of the caudate in smoking addiction.

BACKGROUND: Nearly 6 million deaths and over a half trillion dollars in healthcare costs worldwide are attributed to tobacco smoking each year. Extensive research efforts have been pursued to elucidate the molecular underpinnings of smoking addiction and facilitate cessation. In this study, we genotyped and obtained both resting state and task-based functional magnetic resonance imaging from 64 non-smokers and 42 smokers. Smokers were imaged after having smoked normally ("sated") and after having not smoked for at least 12 h ("abstinent"). RESULTS: While abstinent smokers did not differ from non-smokers with respect to pairwise resting state functional connectivities (RSFCs) between 12 brain regions of interest, RSFCs involving the caudate and putamen of sated smokers significantly differed from those of non-smokers (P < 0.01). Further analyses of caudate and putamen activity during elicited experiences of reward and disappointment show that caudate activity during reward (CR) correlated with smoking status (P = 0.015). Moreover, abstinent smokers with lower CR experienced greater withdrawal symptoms (P = 0.024), which suggests CR may be related to smoking urges. Associations between genetic variants and CR, adjusted for smoking status, were identified by genome-wide association study (GWAS). Genes containing or exhibiting caudate-specific expression regulation by these variants were enriched within Gene Ontology terms that describe cytoskeleton functions, synaptic organization, and injury response (P < 0.001, FDR < 0.05). CONCLUSIONS: By integrating genomic and imaging data, novel insights into potential mechanisms of caudate activation and homeostasis are revealed that may guide new directions of research toward improving our understanding of addiction pathology.

A Novel Approach to Identifying a Neuroimaging Biomarker for Patients With Serious Mental Illness.

Serious mental illness (SMI) is disabling, and current interventions are ineffective for many. This exploratory study sought to demonstrate the feasibility of applying topological data analysis (TDA) to resting-state functional connectivity data obtained from a heterogeneous sample of 235 adult inpatients to identify a biomarker of treatment response. TDA identified two groups based on connectivity between the prefrontal cortex and striatal regions: patients admitted with greater functional connectivity between these regions evidenced less improvement from admission to discharge than patients with lesser connectivity between them. TDA identified a potential biomarker of an attenuated treatment response among inpatients with SMI. Insofar as the observed pattern of resting-state functional connectivity collected early during treatment is replicable, this potential biomarker may indicate the need to modify standard of care for a small, albeit meaningful, percentage of patients.

High resolution single particle refinement in EMAN2.1.

EMAN2.1 is a complete image processing suite for quantitative analysis of grayscale images, with a primary focus on transmission electron microscopy, with complete workflows for performing high resolution single particle reconstruction, 2-D and 3-D heterogeneity analysis, random conical tilt reconstruction and subtomogram averaging, among other tasks. In this manuscript we provide the first detailed description of the high resolution single particle analysis pipeline and the philosophy behind its approach to the reconstruction problem. High resolution refinement is a fully automated process, and involves an advanced set of heuristics to select optimal algorithms for each specific refinement task. A gold standard FSC is produced automatically as part of refinement, providing a robust resolution estimate for the final map, and this is used to optimally filter the final CTF phase and amplitude corrected structure. Additional methods are in-place to reduce model bias during refinement, and to permit cross-validation using other computational methods.

Increased habenular connectivity in opioid users is associated with an α5 subunit nicotinic receptor genetic variant.

BACKGROUND AND OBJECTIVES: Opioid use disorder (OUD) is a chronic disorder with relapse based on both desire for reinforcement (craving) and avoidance of withdrawal. The aversive aspect of dependence and relapse has been associated with a small brain structure called the habenula, which expresses large numbers of both opioid and nicotinic receptors. Additionally, opioid withdrawal symptoms can be induced in opioid-treated rodents by blocking not only opioid, but also nicotinic receptors. This receptor co-localization and cross-induction of withdrawal therefore might lead to genetic variation in the nicotinic receptor influencing development of human opioid dependence through its impact on the aversive components of opioid dependence. METHODS: We studied habenular resting state functional connectivity with related brain structures, specifically the striatum. We compared abstinent psychiatric patients who use opioids (N = 51) to psychiatric patients who do not (N = 254) to identify an endophenotype of opioid use that focused on withdrawal avoidance and aversion rather than the more commonly examined craving aspects of relapse. RESULTS: We found that habenula-striatal connectivity was stronger in opioid-using patients. Increased habenula-striatum connectivity was observed in opioid-using patients with the low risk rs16969968 GG genotype, but not in patients carrying the high risk AG or AA genotypes. CONCLUSIONS: We propose that increased habenula-striatum functional connectivity may be modulated by the nicotinic receptor variant rs16969968 and may lead to increased opioid use. SCIENTIFIC SIGNIFICANCE: Our data uncovered a promising brain target for development of novel anti-addiction therapies and may help the development of personalized therapies against opioid abuse. (Am J Addict 2017;26:751-759).

De Novo modeling in cryo-EM density maps with Pathwalking.

As electron cryo-microscopy (cryo-EM) can now frequently achieve near atomic resolution, accurate interpretation of these density maps in terms of atomistic detail has become paramount in deciphering macromolecular structure and function. However, there are few software tools for modeling protein structure from cryo-EM density maps in this resolution range. Here, we present an extension of our original Pathwalking protocol, which can automatically trace a protein backbone directly from a near-atomic resolution (3-6Å) density map. The original Pathwalking approach utilized a Traveling Salesman Problem solver for backbone tracing, but manual adjustment was still required during modeling. In the new version, human intervention is minimized and we provide a more robust approach for backbone modeling. This includes iterative secondary structure identification, termini detection and the ability to model multiple subunits without prior segmentation. Overall, the new Pathwalking procedure provides a more complete and robust tool for annotating protein structure function in near-atomic resolution density maps.

Alignment algorithms and per-particle CTF correction for single particle cryo-electron tomography.

Single particle cryo-electron tomography (cryoSPT) extracts features from cryo-electron tomograms, followed by 3D classification, alignment and averaging to generate improved 3D density maps of such features. Robust methods to correct for the contrast transfer function (CTF) of the electron microscope are necessary for cryoSPT to reach its resolution potential. Many factors can make CTF correction for cryoSPT challenging, such as lack of eucentricity of the specimen stage, inherent low dose per image, specimen charging, beam-induced specimen motions, and defocus gradients resulting both from specimen tilting and from unpredictable ice thickness variations. Current CTF correction methods for cryoET make at least one of the following assumptions: that the defocus at the center of the image is the same across the images of a tiltseries, that the particles all lie at the same Z-height in the embedding ice, and/or that the specimen, the cryo-electron microscopy (cryoEM) grid and/or the carbon support are flat. These experimental conditions are not always met. We have developed a CTF correction algorithm for cryoSPT without making any of the aforementioned assumptions. We also introduce speed and accuracy improvements and a higher degree of automation to the subtomogram averaging algorithms available in EMAN2. Using motion-corrected images of isolated virus particles as a benchmark specimen, recorded with a DE20 direct detection camera, we show that our CTF correction and subtomogram alignment routines can yield subtomogram averages close to 4/5 Nyquist frequency of the detector under our experimental conditions.

Anterior cingulum white matter is altered in tobacco smokers.

BACKGROUND AND OBJECTIVES: The anterior cingulate cortex (ACC) is hypothesized to be involved in decision making and emotion regulation. Previous observations of drug dependent individuals indicate that substance dependence may be associated with cingulum white matter abnormalities. The present study evaluated cingulum white matter in cigarette smokers. METHODS: Diffusion tensor imaging (DTI) in adult tobacco smokers and healthy non-smoker controls (total N = 70) was performed in a 3T Siemens Trio MRI scanner. RESULTS: Analyses of DTI tractography of the cingulum in tobacco-smoking individuals and controls indicated that tobacco abusers have significantly reduced fractional anisotropy (FA) in the right cingulum. In addition, FA in the left cingulum white matter was negatively associated with the number of cigarettes smoked per day and the Fagerstrom test for nicotine dependence, a self-report measure of tobacco dependence severity. CONCLUSIONS: The white matter of the cingulum is altered in a non-symmetrical way in tobacco smokers. An inverse relationship between FA and reported number of cigarettes per day was observed. Previous studies have also noted altered neural connectivity in cigarette smokers using similar methods. Similar white matter differences in the cingulum have been observed in methamphetamine dependent individuals and patients with dementia, which suggests that the cingulum may be altered by mechanisms not specific to tobacco exposure. SCIENTIFIC SIGNIFICANCE: By better understanding the effects of tobacco abuse on the brain, we hope to gain insight into how drug dependence influences the neurological foundations of behavior.

Synaptic Plasticity in the Injured Brain Depends on the Temporal Pattern of Stimulation.

Neurostimulation protocols are increasingly used as therapeutic interventions, including for brain injury. In addition to the direct activation of neurons, these stimulation protocols are also likely to have downstream effects on those neurons' synaptic outputs. It is well known that alterations in the strength of synaptic connections (long-term potentiation, LTP; long-term depression, LTD) are sensitive to the frequency of stimulation used for induction; however, little is known about the contribution of the temporal pattern of stimulation to the downstream synaptic plasticity that may be induced by neurostimulation in the injured brain. We explored interactions of the temporal pattern and frequency of neurostimulation in the normal cerebral cortex and after mild traumatic brain injury (mTBI), to inform therapies to strengthen or weaken neural circuits in injured brains, as well as to better understand the role of these factors in normal brain plasticity. Whole-cell (WC) patch-clamp recordings of evoked postsynaptic potentials in individual neurons, as well as field potential (FP) recordings, were made from layer 2/3 of visual cortex in response to stimulation of layer 4, in acute slices from control (naive), sham operated, and mTBI rats. We compared synaptic plasticity induced by different stimulation protocols, each consisting of a specific frequency (1 Hz, 10 Hz, or 100 Hz), continuity (continuous or discontinuous), and temporal pattern (perfectly regular, slightly irregular, or highly irregular). At the individual neuron level, dramatic differences in plasticity outcome occurred when the highly irregular stimulation protocol was used at 1 Hz or 10 Hz, producing an overall LTD in controls and shams, but a robust overall LTP after mTBI. Consistent with the individual neuron results, the plasticity outcomes for simultaneous FP recordings were similar, indicative of our results generalizing to a larger scale synaptic network than can be sampled by individual WC recordings alone. In addition to the differences in plasticity outcome between control (naive or sham) and injured brains, the dynamics of the changes in synaptic responses that developed during stimulation were predictive of the final plasticity outcome. Our results demonstrate that the temporal pattern of stimulation plays a role in the polarity and magnitude of synaptic plasticity induced in the cerebral cortex while highlighting differences between normal and injured brain responses. Moreover, these results may be useful for optimization of neurostimulation therapies to treat mTBI and other brain disorders, in addition to providing new insights into downstream plasticity signaling mechanisms in the normal brain.

Missing Wedge Completion via Unsupervised Learning with Coordinate Networks.

Cryogenic electron tomography (cryoET) is a powerful tool in structural biology, enabling detailed 3D imaging of biological specimens at a resolution of nanometers. Despite its potential, cryoET faces challenges such as the missing wedge problem, which limits reconstruction quality due to incomplete data collection angles. Recently, supervised deep learning methods leveraging convolutional neural networks (CNNs) have considerably addressed this issue; however, their pretraining requirements render them susceptible to inaccuracies and artifacts, particularly when representative training data is scarce. To overcome these limitations, we introduce a proof-of-concept unsupervised learning approach using coordinate networks (CNs) that optimizes network weights directly against input projections. This eliminates the need for pretraining, reducing reconstruction runtime by 3 - 20× compared to supervised methods. Our in silico results show improved shape completion and reduction of missing wedge artifacts, assessed through several voxel-based image quality metrics in real space and a novel directional Fourier Shell Correlation (FSC) metric. Our study illuminates benefits and considerations of both supervised and unsupervised approaches, guiding the development of improved reconstruction strategies.

Overcoming resolution attenuation during tilted cryo-EM data collection.

Structural biology efforts using cryogenic electron microscopy are frequently stifled by specimens adopting "preferred orientations" on grids, leading to anisotropic map resolution and impeding structure determination. Tilting the specimen stage during data collection is a generalizable solution but has historically led to substantial resolution attenuation. Here, we develop updated data collection and image processing workflows and demonstrate, using multiple specimens, that resolution attenuation is negligible or significantly reduced across tilt angles. Reconstructions with and without the stage tilted as high as 60° are virtually indistinguishable. These strategies allowed the reconstruction to 3 Å resolution of a bacterial RNA polymerase with preferred orientation, containing an unnatural nucleotide for studying novel base pair recognition. Furthermore, we present a quantitative framework that allows cryo-EM practitioners to define an optimal tilt angle during data acquisition. These results reinforce the utility of employing stage tilt for data collection and provide quantitative metrics to obtain isotropic maps.

Community recommendations on cryoEM data archiving and validation: Outcomes of a wwPDB/EMDB workshop on cryoEM data management, deposition and validation.

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and consensus recommendations resulting from the workshop. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Community recommendations on cryoEM data archiving and validation.

In January 2020, a workshop was held at EMBL-EBI (Hinxton, UK) to discuss data requirements for the deposition and validation of cryoEM structures, with a focus on single-particle analysis. The meeting was attended by 47 experts in data processing, model building and refinement, validation, and archiving of such structures. This report describes the workshop's motivation and history, the topics discussed, and the resulting consensus recommendations. Some challenges for future methods-development efforts in this area are also highlighted, as is the implementation to date of some of the recommendations.

Transformations between rotational and translational invariants formulated in reciprocal spaces.

Correlation functions play an important role in the theoretical underpinnings of many disparate areas of the physical sciences: in particular, scattering theory. More recently, they have become useful in the classification of objects in areas such as computer vision and our area of cryoEM. Our primary classification scheme in the cryoEM image processing system, EMAN2, is now based on third order invariants formulated in Fourier space. This allows a factor of 8 speed up in the two classification procedures inherent in our software pipeline, because it allows for classification without the need for computationally costly alignment procedures. In this work, we address several formal and practical aspects of such multispectral invariants. We show that we can formulate such invariants in the representation in which the original signal is most compact. We explicitly construct transformations between invariants in different orientations for arbitrary order of correlation functions and dimension. We demonstrate that third order invariants distinguish 2D mirrored patterns (unlike the radial power spectrum), which is a fundamental aspects of its classification efficacy. We show the limitations of 3rd order invariants also, by giving an example of a wide family of patterns with identical (vanishing) set of 3rd order invariants. For sufficiently rich patterns, the third order invariants should distinguish typical images, textures and patterns.

Overcoming Resolution Attenuation During Tilted Cryo-EM Data Collection.

Structural biology efforts using cryogenic electron microscopy are frequently stifled by specimens adopting "preferred orientations" on grids, leading to anisotropic map resolution and impeding structure determination. Tilting the specimen stage during data collection is a generalizable solution but has historically led to substantial resolution attenuation. Here, we develop updated data collection and image processing workflows and demonstrate, using multiple specimens, that resolution attenuation is negligible or significantly reduced across tilt angles. Reconstructions with and without the stage tilted as high as 60° are virtually indistinguishable. These strategies allowed the reconstruction to 3 Å resolution of a bacterial RNA polymerase with preferred orientation. Furthermore, we present a quantitative framework that allows cryo-EM practitioners to define an optimal tilt angle for dataset acquisition. These data reinforce the utility of employing stage tilt for data collection and provide quantitative metrics to obtain isotropic maps.

Tools for visualizing and analyzing Fourier space sampling in Cryo-EM.

A complete understanding of how an orientation distribution contributes to a cryo-EM reconstruction remains lacking. It is necessary to begin critically assessing the set of views to gain an understanding of its effect on experimental reconstructions. Toward that end, we recently suggested that the type of orientation distribution may alter resolution measures in a systematic manner. We introduced the sampling compensation factor (SCF), which incorporates how the collection geometry might change the spectral signal-to-noise ratio (SSNR), irrespective of the other experimental aspects. We show here that knowledge of the sampling restricted to spherical surfaces of sufficiently large radii in Fourier space is equivalent to knowledge of the set of projection views. Moreover, the SCF geometrical factor may be calculated from one such surface. To aid cryo-EM practitioners, we developed a graphical user interface (GUI) tool that evaluates experimental orientation distributions. The GUI returns plots of projection directions, sampling constrained to the surface of a sphere, the SCF value, the fraction of the empty region of Fourier space, and a histogram of the sampling values over the points on a sphere. Finally, a fixed tilt angle may be incorporated to determine how tilting the grid during collection may improve the distribution of views and Fourier space sampling. We advocate this simple conception of sampling and the use of such tools as a complement to the distribution of views to capture the different aspects of the effect of projection directions on cryo-EM reconstructions.