ArXiv最新文献

Recent advances in techniques for monitoring and perturbing neural populations have greatly enhanced our ability to study circuits in the brain. In particular, two-photon holographic optogenetics now enables precise photostimulation of experimenter-specified groups of individual neurons, while simultaneous two-photon calcium imaging enables the measurement of ongoing and induced activity across the neural population. Despite the enormous space of potential photostimulation patterns and the time-consuming nature of photostimulation experiments, very little algorithmic work has been done to determine the most effective photostimulation patterns for identifying the neural population dynamics. Here, we develop methods to efficiently select which neurons to stimulate such that the resulting neural responses will best inform a dynamical model of the neural population activity. Using neural population responses to photostimulation in mouse motor cortex, we demonstrate the efficacy of a low-rank linear dynamical systems model, and develop an active learning procedure which takes advantage of low-rank structure to determine informative photostimulation patterns. We demonstrate our approach on both real and synthetic data, obtaining in some cases as much as a two-fold reduction in the amount of data required to reach a given predictive power. Our active stimulation design method is based on a novel active learning procedure for low-rank regression, which may be of independent interest.

It is widely recognized that reciprocal interactions between cells and their microenvironment, via mechanical forces and biochemical signaling pathways, regulate cell behaviors during normal development, homeostasis and disease progression such as cancer. However, it is still not well understood how complex patterns of tissue growth emerge. Here, we propose a framework for the chemomechanical regulation of growth based on thermodynamics of continua and growth-elasticity to predict growth patterns. Combining the elastic and chemical energies, we use an energy variational approach to derive a novel formulation that incorporates an energy-dissipating stress relaxation and biochemomechanical regulation of the volumetric growth rate. We validate the model using experimental data from growth of tumor spheroids in confined environments. We also investigate the influence of model parameters, including tissue rearrangement rate, tissue compressibility, strength of mechanical feedback and external mechanical stimuli, on the growth patterns of tumor spheroids.

F"orster resonance energy transfer (FRET) is a quantum mechanical phenomenon involving the non-radiative transfer of energy between coupled electric dipoles. Due to the strong dependence of FRET on the distance between the dipoles, it is frequently used as a ``molecular ruler" in biology, chemistry, and physics. This is done by placing dipolar molecules called dyes on molecules of interest. In time-resolved confocal single-molecule FRET (smFRET) experiments, the joint distribution of the FRET efficiency and the donor fluorescence lifetime can reveal underlying molecular conformational dynamics via deviation from their theoretical F"orster relationship. This deviation is referred to as a dynamic shift. Quantifying the dynamic shift caused by the motion of the fluorescent dyes is essential to decoupling the dynamics of the studied molecules and the dyes. We develop novel Langevin models for the dye linker dynamics, including rotational dynamics, based on first principle physics and proper dye linker chemistry to match accessible volumes predicted by molecular dynamics simulations. By simulating the dyes' stochastic translational and rotational dynamics, we show that the observed dynamic shift can largely be attributed to the mutual orientational dynamics of the electric dipole moments associated with the dyes, not their accessible volume. Our models provide the most up-to-date and accurate estimation of FRET.

Human braingraphs or connectomes are widely studied in the last decade to understand the structural and functional properties of our brain. In the last several years our research group has computed and deposited thousands of human braingraphs to the braingraph.org site, by applying public structural (diffusion) MRI data from young and healthy subjects. Here we describe a recent addition to the {tt braingraph.org} site, which contains connectomes from healthy and demented subjects between 42 and 95 years of age, based on the public release of the OASIS-3 dataset. The diffusion MRI data was processed with the Connectome Mapper Toolkit v.3.1. We believe that the new addition to the braingraph.org site will become a useful resource for enlightening the aging circuitry of the human brain in healthy and diseased subjects, including those with Alzheimer's disease in several stages.

In certain three-dimensional (3D) applications of photoacoustic computed tomography (PACT), including textit{in vivo} breast imaging, hemispherical measurement apertures that enclose the object within their convex hull are employed for data acquisition. Data acquired with such measurement geometries are referred to as textit{half-scan} data, as only half of a complete spherical measurement aperture is employed. Although previous studies have demonstrated that half-scan data can uniquely and stably reconstruct the sought-after object, no closed-form reconstruction formula for use with half-scan data has been reported. To address this, a semi-analytic reconstruction method in the form of filtered backprojection (FBP), referred to as the half-scan FBP method, is developed in this work. Because the explicit form of the filtering operation in the half-scan FBP method is not currently known, a learning-based method is proposed to approximate it. The proposed method is systematically investigated by use of virtual imaging studies of 3D breast PACT that employ ensembles of numerical breast phantoms and a physics-based model of the data acquisition process. The method is subsequently applied to experimental data acquired in an textit{in vivo} breast PACT study. The results confirm that the half-scan FBP method can accurately reconstruct 3D images from half-scan data. Importantly, because the sought-after inverse mapping is well-posed, the reconstruction method remains accurate even when applied to data that differ considerably from those employed to learn the filtering operation.

Designing novel functional proteins crucially depends on accurately modeling their fitness landscape. Given the limited availability of functional annotations from wet-lab experiments, previous methods have primarily relied on self-supervised models trained on vast, unlabeled protein sequence or structure datasets. While initial protein representation learning studies solely focused on either sequence or structural features, recent hybrid architectures have sought to merge these modalities to harness their respective strengths. However, these sequence-structure models have so far achieved only incremental improvements when compared to the leading sequence-only approaches, highlighting unresolved challenges effectively leveraging these modalities together. Moreover, the function of certain proteins is highly dependent on the granular aspects of their surface topology, which have been overlooked by prior models. To address these limitations, we introduce the Sequence-Structure-Surface Fitness (S3F) model - a novel multimodal representation learning framework that integrates protein features across several scales. Our approach combines sequence representations from a protein language model with Geometric Vector Perceptron networks encoding protein backbone and detailed surface topology. The proposed method achieves state-of-the-art fitness prediction on the ProteinGym benchmark encompassing 217 substitution deep mutational scanning assays, and provides insights into the determinants of protein function. Our code is at https://github.com/DeepGraphLearning/S3F.

Motivation: With the rapid expansion of large-scale biological datasets, DNA and protein sequence alignments have become essential for comparative genomics and proteomics. These alignments facilitate the exploration of sequence similarity patterns, providing valuable insights into sequence conservation, evolutionary relationships and for functional analyses. Typically, sequence alignments are stored in formats such as the Multiple Alignment Format (MAF). Counting k-mer occurrences is a crucial task in many computational biology applications, but currently, there is no algorithm designed for k-mer counting in alignment files.

Results: We have developed MAFcounter, the first k-mer counter dedicated to alignment files. MAFcounter is multithreaded, fast, and memory efficient, enabling k-mer counting in DNA and protein sequence alignment files.

Availability: The MAFcounter package and its Python bindings are released under GPL license as a multi-platform application and are available at: https://github.com/Georgakopoulos-Soares-lab/MAFcounter.

Motivated by elastography that utilizes tissue mechanical properties as biomarkers for liver disease, with the eventual objective of quantitatively linking histopathology and bulk mechanical properties, we develop a micromechanical modeling approach to capture the effects of fat and collagen deposition in the liver. Specifically, we utilize computational homogenization to convert the microstructural changes in hepatic lobule to the effective viscoelastic modulus of the liver tissue, i.e., predict the bulk material properties by analyzing the deformation of repeating unit cell. The lipid and collagen deposition is simulated with the help of ad hoc algorithms informed by histological observations. Collagen deposition is directly included in the computational model, while composite material theory is used to convert fat content to the microscopic mechanical properties, which in turn is included in the computational model. The results illustrate the model's ability to capture the effect of both fat and collagen deposition on the viscoelastic moduli and represents a step towards linking histopathological changes in the liver to its bulk mechanical properties, which can eventually provide insights for accurate diagnosis with elastography.

Magnetic resonance (MR)-guided radiation therapy (RT) is enhancing head and neck cancer (HNC) treatment through superior soft tissue contrast and longitudinal imaging capabilities. However, manual tumor segmentation remains a significant challenge, spurring interest in artificial intelligence (AI)-driven automation. To accelerate innovation in this field, we present the Head and Neck Tumor Segmentation for MR-Guided Applications (HNTS-MRG) 2024 Challenge, a satellite event of the 27th International Conference on Medical Image Computing and Computer Assisted Intervention. This challenge addresses the scarcity of large, publicly available AI-ready adaptive RT datasets in HNC and explores the potential of incorporating multi-timepoint data to enhance RT auto-segmentation performance. Participants tackled two HNC segmentation tasks: automatic delineation of primary gross tumor volume (GTVp) and gross metastatic regional lymph nodes (GTVn) on pre-RT (Task 1) and mid-RT (Task 2) T2-weighted scans. The challenge provided 150 HNC cases for training and 50 for testing, hosted on grand-challenge.org using a Docker submission framework. In total, 19 independent teams from across the world qualified by submitting both their algorithms and corresponding papers, resulting in 18 submissions for Task 1 and 15 submissions for Task 2. Evaluation using the mean aggregated Dice Similarity Coefficient showed top-performing AI methods achieved scores of 0.825 in Task 1 and 0.733 in Task 2. These results surpassed clinician interobserver variability benchmarks, marking significant strides in automated tumor segmentation for MR-guided RT applications in HNC.