Structure-based discovery of CFTR potentiators and inhibitors.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a crucial ion channel whose loss of function leads to cystic fibrosis, whereas its hyperactivation leads to secretory diarrhea. Small molecules that improve CFTR folding (correctors) or function (potentiators) are clinically available. However, the only potentiator, ivacaftor, has suboptimal pharmacokinetics and inhibitors have yet to be clinically developed. Here, we combine molecular docking, electrophysiology, cryo-EM, and medicinal chemistry to identify CFTR modulators. We docked ∼155 million molecules into the potentiator site on CFTR, synthesized 53 test ligands, and used structure-based optimization to identify candidate modulators. This approach uncovered mid-nanomolar potentiators, as well as inhibitors, that bind to the same allosteric site. These molecules represent potential leads for the development of more effective drugs for cystic fibrosis and secretory diarrhea, demonstrating the feasibility of large-scale docking for ion channel drug discovery.

Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

High-precision coding in visual cortex.

Individual neurons in visual cortex provide the brain with unreliable estimates of visual features. It is not known whether the single-neuron variability is correlated across large neural populations, thus impairing the global encoding of stimuli. We recorded simultaneously from up to 50,000 neurons in mouse primary visual cortex (V1) and in higher order visual areas and measured stimulus discrimination thresholds of 0.35° and 0.37°, respectively, in an orientation decoding task. These neural thresholds were almost 100 times smaller than the behavioral discrimination thresholds reported in mice. This discrepancy could not be explained by stimulus properties or arousal states. Furthermore, behavioral variability during a sensory discrimination task could not be explained by neural variability in V1. Instead, behavior-related neural activity arose dynamically across a network of non-sensory brain areas. These results imply that perceptual discrimination in mice is limited by downstream decoders, not by neural noise in sensory representations.

Changes of Cell Biochemical States Are Revealed in Protein Homomeric Complex Dynamics.

We report here a simple and global strategy to map out gene functions and target pathways of drugs, toxins, or other small molecules based on "homomer dynamics" protein-fragment complementation assays (hdPCA). hdPCA measures changes in self-association (homomerization) of over 3,500 yeast proteins in yeast grown under different conditions. hdPCA complements genetic interaction measurements while eliminating the confounding effects of gene ablation. We demonstrate that hdPCA accurately predicts the effects of two longevity and health span-affecting drugs, the immunosuppressant rapamycin and the type 2 diabetes drug metformin, on cellular pathways. We also discovered an unsuspected global cellular response to metformin that resembles iron deficiency and includes a change in protein-bound iron levels. This discovery opens a new avenue to investigate molecular mechanisms for the prevention or treatment of diabetes, cancers, and other chronic diseases of aging.

Target-Specific Precision of CRISPR-Mediated Genome Editing.

The CRISPR-Cas9 system has successfully been adapted to edit the genome of various organisms. However, our ability to predict the editing outcome at specific sites is limited. Here, we examined indel profiles at over 1,000 genomic sites in human cells and uncovered general principles guiding CRISPR-mediated DNA editing. We find that precision of DNA editing (i.e., recurrence of a specific indel) varies considerably among sites, with some targets showing one highly preferred indel and others displaying numerous infrequent indels. Editing precision correlates with editing efficiency and a preference for single-nucleotide homologous insertions. Precise targets and editing outcome can be predicted based on simple rules that mainly depend on the fourth nucleotide upstream of the protospacer adjacent motif (PAM). Indel profiles are robust, but they can be influenced by chromatin features. Our findings have important implications for clinical applications of CRISPR technology and reveal general patterns of broken end joining that can provide insights into DNA repair mechanisms.

Space wandering in the rodent default mode network.

The default mode network (DMN) is a large-scale brain network known to be suppressed during a wide range of cognitive tasks. However, our comprehension of its role in naturalistic and unconstrained behaviors has remained elusive because most research on the DMN has been conducted within the restrictive confines of MRI scanners. Here, we use multisite GCaMP (a genetically encoded calcium indicator) fiber photometry with simultaneous videography to probe DMN function in awake, freely exploring rats. We examined neural dynamics in three core DMN nodes-the retrosplenial cortex, cingulate cortex, and prelimbic cortex-as well as the anterior insula node of the salience network, and their association with the rats' spatial exploration behaviors. We found that DMN nodes displayed a hierarchical functional organization during spatial exploration, characterized by stronger coupling with each other than with the anterior insula. Crucially, these DMN nodes encoded the kinematics of spatial exploration, including linear and angular velocity. Additionally, we identified latent brain states that encoded distinct patterns of time-varying exploration behaviors and found that higher linear velocity was associated with enhanced DMN activity, heightened synchronization among DMN nodes, and increased anticorrelation between the DMN and anterior insula. Our findings highlight the involvement of the DMN in collectively and dynamically encoding spatial exploration in a real-world setting. Our findings challenge the notion that the DMN is primarily a "task-negative" network disengaged from the external world. By illuminating the DMN's role in naturalistic behaviors, our study underscores the importance of investigating brain network function in ecologically valid contexts.

Multiscale Flow for robust and optimal cosmological analysis.

We propose Multiscale Flow, a generative Normalizing Flow that creates samples and models the field-level likelihood of two-dimensional cosmological data such as weak lensing. Multiscale Flow uses hierarchical decomposition of cosmological fields via a wavelet basis and then models different wavelet components separately as Normalizing Flows. The log-likelihood of the original cosmological field can be recovered by summing over the log-likelihood of each wavelet term. This decomposition allows us to separate the information from different scales and identify distribution shifts in the data such as unknown scale-dependent systematics. The resulting likelihood analysis can not only identify these types of systematics, but can also be made optimal, in the sense that the Multiscale Flow can learn the full likelihood at the field without any dimensionality reduction. We apply Multiscale Flow to weak lensing mock datasets for cosmological inference and show that it significantly outperforms traditional summary statistics such as power spectrum and peak counts, as well as machine learning-based summary statistics such as scattering transform and convolutional neural networks. We further show that Multiscale Flow is able to identify distribution shifts not in the training data such as baryonic effects. Finally, we demonstrate that Multiscale Flow can be used to generate realistic samples of weak lensing data.

Multifractality in spin glasses.

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics.

A Guide to Emerging Technologies for Large-Scale and Whole-Brain Optical Imaging of Neuronal Activity.

The mammalian brain is a densely interconnected network that consists of millions to billions of neurons. Decoding how information is represented and processed by this neural circuitry requires the ability to capture and manipulate the dynamics of large populations at high speed and high resolution over a large area of the brain. Although the use of optical approaches by the neuroscience community has rapidly increased over the past two decades, most microscopy approaches are unable to record the activity of all neurons comprising a functional network across the mammalian brain at relevant temporal and spatial resolutions. In this review, we survey the recent development in optical technologies for Ca2+ imaging in this regard and provide an overview of the strengths and limitations of each modality and its potential for scalability. We provide guidance from the perspective of a biological user driven by the typical biological applications and sample conditions. We also discuss the potential for future advances and synergies that could be obtained through hybrid approaches or other modalities.

Biolinguistic graph fusion model for circRNA-miRNA association prediction.

Emerging clinical evidence suggests that sophisticated associations with circular ribonucleic acids (RNAs) (circRNAs) and microRNAs (miRNAs) are a critical regulatory factor of various pathological processes and play a critical role in most intricate human diseases. Nonetheless, the above correlations via wet experiments are error-prone and labor-intensive, and the underlying novel circRNA-miRNA association (CMA) has been validated by numerous existing computational methods that rely only on single correlation data. Considering the inadequacy of existing machine learning models, we propose a new model named BGF-CMAP, which combines the gradient boosting decision tree with natural language processing and graph embedding methods to infer associations between circRNAs and miRNAs. Specifically, BGF-CMAP extracts sequence attribute features and interaction behavior features by Word2vec and two homogeneous graph embedding algorithms, large-scale information network embedding and graph factorization, respectively. Multitudinous comprehensive experimental analysis revealed that BGF-CMAP successfully predicted the complex relationship between circRNAs and miRNAs with an accuracy of 82.90% and an area under receiver operating characteristic of 0.9075. Furthermore, 23 of the top 30 miRNA-associated circRNAs of the studies on data were confirmed in relevant experiences, showing that the BGF-CMAP model is superior to others. BGF-CMAP can serve as a helpful model to provide a scientific theoretical basis for the study of CMA prediction.

Farmland practices are driving bird population decline across Europe.

Declines in European bird populations are reported for decades but the direct effect of major anthropogenic pressures on such declines remains unquantified. Causal relationships between pressures and bird population responses are difficult to identify as pressures interact at different spatial scales and responses vary among species. Here, we uncover direct relationships between population time-series of 170 common bird species, monitored at more than 20,000 sites in 28 European countries, over 37 y, and four widespread anthropogenic pressures: agricultural intensification, change in forest cover, urbanisation and temperature change over the last decades. We quantify the influence of each pressure on population time-series and its importance relative to other pressures, and we identify traits of most affected species. We find that agricultural intensification, in particular pesticides and fertiliser use, is the main pressure for most bird population declines, especially for invertebrate feeders. Responses to changes in forest cover, urbanisation and temperature are more species-specific. Specifically, forest cover is associated with a positive effect and growing urbanisation with a negative effect on population dynamics, while temperature change has an effect on the dynamics of a large number of bird populations, the magnitude and direction of which depend on species' thermal preferences. Our results not only confirm the pervasive and strong effects of anthropogenic pressures on common breeding birds, but quantify the relative strength of these effects stressing the urgent need for transformative changes in the way of inhabiting the world in European countries, if bird populations shall have a chance of recovering.

Large-scale land acquisitions exacerbate local farmland inequalities in Tanzania.

Land inequality stalls economic development, entrenches poverty, and is associated with environmental degradation. Yet, rigorous assessments of land-use interventions attend to inequality only rarely. A land inequality lens is especially important to understand how recent large-scale land acquisitions (LSLAs) affect smallholder and indigenous communities across as much as 100 million hectares around the world. This paper studies inequalities in land assets, specifically landholdings and farm size, to derive insights into the distributional outcomes of LSLAs. Using a household survey covering four pairs of land acquisition and control sites in Tanzania, we use a quasi-experimental design to characterize changes in land inequality and subsequent impacts on well-being. We find convincing evidence that LSLAs in Tanzania lead to both reduced landholdings and greater farmland inequality among smallholders. Households in proximity to LSLAs are associated with 21.1% (P = 0.02) smaller landholdings while evidence, although insignificant, is suggestive that farm sizes are also declining. Aggregate estimates, however, hide that households in the bottom quartiles of farm size suffer the brunt of landlessness and land loss induced by LSLAs that combine to generate greater farmland inequality. Additional analyses find that land inequality is not offset by improvements in other livelihood dimensions, rather farm size decreases among households near LSLAs are associated with no income improvements, lower wealth, increased poverty, and higher food insecurity. The results demonstrate that without explicit consideration of distributional outcomes, land-use policies can systematically reinforce existing inequalities.

Data-driven large-scale genomic analysis reveals an intricate phylogenetic and functional landscape in J-domain proteins.

The 70-kD heat shock protein (Hsp70) chaperone system is a central hub of the proteostasis network that helps maintain protein homeostasis in all organisms. The recruitment of Hsp70 to perform different and specific cellular functions is regulated by the J-domain protein (JDP) co-chaperone family carrying the small namesake J-domain, required to interact and drive the ATPase cycle of Hsp70s. Besides the J-domain, prokaryotic and eukaryotic JDPs display a staggering diversity in domain architecture, function, and cellular localization. Very little is known about the overall JDP family, despite their essential role in cellular proteostasis, development, and its link to a broad range of human diseases. In this work, we leverage the exponentially increasing number of JDP gene sequences identified across all kingdoms owing to the advancements in sequencing technology and provide a broad overview of the JDP repertoire. Using an automated classification scheme based on artificial neural networks (ANNs), we demonstrate that the sequences of J-domains carry sufficient discriminatory information to reliably recover the phylogeny, localization, and domain composition of the corresponding full-length JDP. By harnessing the interpretability of the ANNs, we find that many of the discriminatory sequence positions match residues that form the interaction interface between the J-domain and Hsp70. This reveals that key residues within the J-domains have coevolved with their obligatory Hsp70 partners to build chaperone circuits for specific functions in cells.

Census of exposed aggregation-prone regions in proteomes.

Loss of solubility usually leads to the detrimental elimination of protein function. In some cases, the protein aggregation is also required for beneficial functions. Given the duality of this phenomenon, it remains a fundamental question how natural selection controls the aggregation. The exponential growth of genomic sequence data and recent progress with in silico predictors of the aggregation allows approaching this problem by a large-scale bioinformatics analysis. Most of the aggregation-prone regions are hidden within the 3D structure, rendering them inaccessible for the intermolecular interactions responsible for aggregation. Thus, the most realistic census of the aggregation-prone regions requires crossing aggregation prediction with information about the location of the natively unfolded regions. This allows us to detect so-called 'exposed aggregation-prone regions' (EARs). Here, we analyzed the occurrence and distribution of the EARs in 76 reference proteomes from the three kingdoms of life. For this purpose, we used a bioinformatics pipeline, which provides a consensual result based on several predictors of aggregation. Our analysis revealed a number of new statistically significant correlations about the presence of EARs in different organisms, their dependence on protein length, cellular localizations, co-occurrence with short linear motifs and the level of protein expression. We also obtained a list of proteins with the conserved aggregation-prone sequences for further experimental tests. Insights gained from this work led to a deeper understanding of the relationship between protein evolution and aggregation.

CIForm as a Transformer-based model for cell-type annotation of large-scale single-cell RNA-seq data.

Single-cell omics technologies have made it possible to analyze the individual cells within a biological sample, providing a more detailed understanding of biological systems. Accurately determining the cell type of each cell is a crucial goal in single-cell RNA-seq (scRNA-seq) analysis. Apart from overcoming the batch effects arising from various factors, single-cell annotation methods also face the challenge of effectively processing large-scale datasets. With the availability of an increase in the scRNA-seq datasets, integrating multiple datasets and addressing batch effects originating from diverse sources are also challenges in cell-type annotation. In this work, to overcome the challenges, we developed a supervised method called CIForm based on the Transformer for cell-type annotation of large-scale scRNA-seq data. To assess the effectiveness and robustness of CIForm, we have compared it with some leading tools on benchmark datasets. Through the systematic comparisons under various cell-type annotation scenarios, we exhibit that the effectiveness of CIForm is particularly pronounced in cell-type annotation. The source code and data are available at https://github.com/zhanglab-wbgcas/CIForm.

Large-scale parameters framework with large convolutional kernel for encoding visual fMRI activity information.

Visual encoding models often use deep neural networks to describe the brain's visual cortex response to external stimuli. Inspired by biological findings, researchers found that large receptive fields built with large convolutional kernels improve convolutional encoding model performance. Inspired by scaling laws in recent years, this article investigates the performance of large convolutional kernel encoding models on larger parameter scales. This paper proposes a large-scale parameters framework with a sizeable convolutional kernel for encoding visual functional magnetic resonance imaging activity information. The proposed framework consists of three parts: First, the stimulus image feature extraction module is constructed using a large-kernel convolutional network while increasing channel numbers to expand the parameter size of the framework. Second, enlarging the input data during the training stage through the multi-subject fusion module to accommodate the increase in parameters. Third, the voxel mapping module maps from stimulus image features to functional magnetic resonance imaging signals. Compared to sizeable convolutional kernel visual encoding networks with base parameter scale, our visual encoding framework improves by approximately 7% on the Natural Scenes Dataset, the dedicated dataset for the Algonauts 2023 Challenge. We further analyze that our encoding framework made a trade-off between encoding performance and trainability. This paper confirms that expanding parameters in visual coding can bring performance improvements.

Facilely tuning the intrinsic catalytic sites of the spinel oxide for peroxymonosulfate activation: From fundamental investigation to pilot-scale demonstration.

Heterogeneous peroxymonosulfate (PMS)-based advanced oxidation processes (AOPs) have shown a great potential for pollutant degradation, but their feasibility for large-scale water treatment application has not been demonstrated. Herein, we develop a facile coprecipitation method for the scalable production (∼10 kg) of the Cu-Fe-Mn spinel oxide (CuFeMnO). Such a catalyst has rich oxygen vacancies and symmetry-breaking sites, which endorse it with a superior PMS-catalytic capacity. We find that the working reactive species and their contributions are highly dependent on the properties of target organic pollutants. For the organics with electron-donating group (e.g., -OH), high-valent metal species are mainly responsible for the pollutant degradation, whereas for the organics with electron-withdrawing group (e.g., -COOH and -NO2), hydroxyl radical (•OH) as the secondary oxidant also plays an important role. We demonstrate that the CuFeMnO-PMS system is able to achieve efficient and stable removal of the pollutants in the secondary effluent from a municipal wastewater plant at both bench and pilot scales. Moreover, we explore the application prospect of this PMS-based AOP process for large-scale wastewater treatment. This work describes an opportunity to scalably prepare robust spinel oxide catalysts for water purification and is beneficial to the practical applications of the heterogeneous PMS-AOPs.

A pantropical assessment of deforestation caused by industrial mining.

Growing demand for minerals continues to drive deforestation worldwide. Tropical forests are particularly vulnerable to the environmental impacts of mining and mineral processing. Many local- to regional-scale studies document extensive, long-lasting impacts of mining on biodiversity and ecosystem services. However, the full scope of deforestation induced by industrial mining across the tropics is yet unknown. Here, we present a biome-wide assessment to show where industrial mine expansion has caused the most deforestation from 2000 to 2019. We find that 3,264 km2 of forest was directly lost due to industrial mining, with 80% occurring in only four countries: Indonesia, Brazil, Ghana, and Suriname. Additionally, controlling for other nonmining determinants of deforestation, we find that mining caused indirect forest loss in two-thirds of the investigated countries. Our results illustrate significant yet unevenly distributed and often unmanaged impacts on these biodiverse ecosystems. Impact assessments and mitigation plans of industrial mining activities must address direct and indirect impacts to support conservation of the world's tropical forests.

Hierarchical timescales in the neocortex: Mathematical mechanism and biological insights.

A cardinal feature of the neocortex is the progressive increase of the spatial receptive fields along the cortical hierarchy. Recently, theoretical and experimental findings have shown that the temporal response windows also gradually enlarge, so that early sensory neural circuits operate on short timescales whereas higher-association areas are capable of integrating information over a long period of time. While an increased receptive field is accounted for by spatial summation of inputs from neurons in an upstream area, the emergence of timescale hierarchy cannot be readily explained, especially given the dense interareal cortical connectivity known in the modern connectome. To uncover the required neurobiological properties, we carried out a rigorous analysis of an anatomically based large-scale cortex model of macaque monkeys. Using a perturbation method, we show that the segregation of disparate timescales is defined in terms of the localization of eigenvectors of the connectivity matrix, which depends on three circuit properties: 1) a macroscopic gradient of synaptic excitation, 2) distinct electrophysiological properties between excitatory and inhibitory neuronal populations, and 3) a detailed balance between long-range excitatory inputs and local inhibitory inputs for each area-to-area pathway. Our work thus provides a quantitative understanding of the mechanism underlying the emergence of timescale hierarchy in large-scale primate cortical networks.

Concentration-Driven Interfacial Amorphization toward Highly Stable and High-Rate Zn Metal Batteries.

The interfacial structure holds great promise in suppressing dendrite growth and parasitic reactions of zinc metal in aqueous media. Current advancements prioritize novel component fabrication, yet the local crystal structure significantly impacts the interfacial properties. In addition, there is still a critical need for scalable synthesis methods for expediting the commercialization of aqueous zinc metal batteries (AZMBs). Herein, we propose a scalable concentration-controlled method for realizing crystalline to amorphous transformation of the Zn metal interface with exceptional scalability (>1 m2) and processing consistency (>30 trials). Theoretical and experimental analyses highlight the advantages of amorphous ZnO, which exhibits moderate adsorption energy, strong desolvation ability, and hydrophilicity. Employing the amorphous ZnO-coated zinc metal anode (AZO-Zn) significantly enhances the cycling performance, impressively maintaining 1000 cycles at 100 mA cm-2. The prototype AZO-Zn||MnO2@CNT pouch cell demonstrates a capacity of 15.7 mAh and maintains 91% of its highest capacity over 100 cycles, presenting promising avenues for the future commercialization of AZMBs.