Structural basis for ion selectivity in potassium-selective channelrhodopsins.

KCR channelrhodopsins (K+-selective light-gated ion channels) have received attention as potential inhibitory optogenetic tools but more broadly pose a fundamental mystery regarding how their K+ selectivity is achieved. Here, we present 2.5-2.7 Å cryo-electron microscopy structures of HcKCR1 and HcKCR2 and of a structure-guided mutant with enhanced K+ selectivity. Structural, electrophysiological, computational, spectroscopic, and biochemical analyses reveal a distinctive mechanism for K+ selectivity; rather than forming the symmetrical filter of canonical K+ channels achieving both selectivity and dehydration, instead, three extracellular-vestibule residues within each monomer form a flexible asymmetric selectivity gate, while a distinct dehydration pathway extends intracellularly. Structural comparisons reveal a retinal-binding pocket that induces retinal rotation (accounting for HcKCR1/HcKCR2 spectral differences), and design of corresponding KCR variants with increased K+ selectivity (KALI-1/KALI-2) provides key advantages for optogenetic inhibition in vitro and in vivo. Thus, discovery of a mechanism for ion-channel K+ selectivity also provides a framework for next-generation optogenetics.

Impaired neural replay of inferred relationships in schizophrenia.

An ability to build structured mental maps of the world underpins our capacity to imagine relationships between objects that extend beyond experience. In rodents, such representations are supported by sequential place cell reactivations during rest, known as replay. Schizophrenia is proposed to reflect a compromise in structured mental representations, with animal models reporting abnormalities in hippocampal replay and associated ripple activity during rest. Here, utilizing magnetoencephalography (MEG), we tasked patients with schizophrenia and control participants to infer unobserved relationships between objects by reorganizing visual experiences containing these objects. During a post-task rest session, controls exhibited fast spontaneous neural reactivation of presented objects that replayed inferred relationships. Replay was coincident with increased ripple power in hippocampus. Patients showed both reduced replay and augmented ripple power relative to controls, convergent with findings in animal models. These abnormalities are linked to impairments in behavioral acquisition and subsequent neural representation of task structure.

Antidepressant drugs act by directly binding to TRKB neurotrophin receptors.

It is unclear how binding of antidepressant drugs to their targets gives rise to the clinical antidepressant effect. We discovered that the transmembrane domain of tyrosine kinase receptor 2 (TRKB), the brain-derived neurotrophic factor (BDNF) receptor that promotes neuronal plasticity and antidepressant responses, has a cholesterol-sensing function that mediates synaptic effects of cholesterol. We then found that both typical and fast-acting antidepressants directly bind to TRKB, thereby facilitating synaptic localization of TRKB and its activation by BDNF. Extensive computational approaches including atomistic molecular dynamics simulations revealed a binding site at the transmembrane region of TRKB dimers. Mutation of the TRKB antidepressant-binding motif impaired cellular, behavioral, and plasticity-promoting responses to antidepressants in vitro and in vivo. We suggest that binding to TRKB and allosteric facilitation of BDNF signaling is the common mechanism for antidepressant action, which may explain why typical antidepressants act slowly and how molecular effects of antidepressants are translated into clinical mood recovery.

How GPCR Phosphorylation Patterns Orchestrate Arrestin-Mediated Signaling.

Binding of arrestin to phosphorylated G-protein-coupled receptors (GPCRs) controls many aspects of cell signaling. The number and arrangement of phosphates may vary substantially for a given GPCR, and different phosphorylation patterns trigger different arrestin-mediated effects. Here, we determine how GPCR phosphorylation influences arrestin behavior by using atomic-level simulations and site-directed spectroscopy to reveal the effects of phosphorylation patterns on arrestin binding and conformation. We find that patterns favoring binding differ from those favoring activation-associated conformational change. Both binding and conformation depend more on arrangement of phosphates than on their total number, with phosphorylation at different positions sometimes exerting opposite effects. Phosphorylation patterns selectively favor a wide variety of arrestin conformations, differently affecting arrestin sites implicated in scaffolding distinct signaling proteins. We also reveal molecular mechanisms of these phenomena. Our work reveals the structural basis for the long-standing "barcode" hypothesis and has important implications for design of functionally selective GPCR-targeted drugs.

Mechanobiology of T Cell Activation: To Catch a Bond.

T cell activation is a critical event in the adaptive immune response, indispensable for cell-mediated and humoral immunity as well as for immune regulation. Recent years have witnessed an emerging trend emphasizing the essential role that physical force and mechanical properties play at the T cell interface. In this review, we integrate current knowledge of T cell antigen recognition and the different models of T cell activation from the perspective of mechanobiology, focusing on the interaction between the T cell receptor (TCR) and the peptide-major histocompatibility complex (pMHC) antigen. We address the shortcomings of TCR affinity alone in explaining T cell functional outcomes and the rising status of force-regulated TCR bond lifetimes, most notably the TCR catch bond. Ultimately, T cell activation and the ensuing physiological responses result from mechanical interaction between TCRs and the pMHC.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

Determinants of Polar versus Nematic Organization in Networks of Dynamic Microtubules and Mitotic Motors.

During cell division, mitotic motors organize microtubules in the bipolar spindle into either polar arrays at the spindle poles or a "nematic" network of aligned microtubules at the spindle center. The reasons for the distinct self-organizing capacities of dynamic microtubules and different motors are not understood. Using in vitro reconstitution experiments and computer simulations, we show that the human mitotic motors kinesin-5 KIF11 and kinesin-14 HSET, despite opposite directionalities, can both organize dynamic microtubules into either polar or nematic networks. We show that in addition to the motor properties the natural asymmetry between microtubule plus- and minus-end growth critically contributes to the organizational potential of the motors. We identify two control parameters that capture system composition and kinetic properties and predict the outcome of microtubule network organization. These results elucidate a fundamental design principle of spindle bipolarity and establish general rules for active filament network organization.

Structural Basis of Membrane Protein Chaperoning through the Mitochondrial Intermembrane Space.

The exchange of metabolites between the mitochondrial matrix and the cytosol depends on ß-barrel channels in the outer membrane and α-helical carrier proteins in the inner membrane. The essential translocase of the inner membrane (TIM) chaperones escort these proteins through the intermembrane space, but the structural and mechanistic details remain elusive. We have used an integrated structural biology approach to reveal the functional principle of TIM chaperones. Multiple clamp-like binding sites hold the mitochondrial membrane proteins in a translocation-competent elongated form, thus mimicking characteristics of co-translational membrane insertion. The bound preprotein undergoes conformational dynamics within the chaperone binding clefts, pointing to a multitude of dynamic local binding events. Mutations in these binding sites cause cell death or growth defects associated with impairment of carrier and ß-barrel protein biogenesis. Our work reveals how a single mitochondrial "transfer-chaperone" system is able to guide α-helical and ß-barrel membrane proteins in a "nascent chain-like" conformation through a ribosome-free compartment.

Whole-Cell Models and Simulations in Molecular Detail.

Comprehensive data about the composition and structure of cellular components have enabled the construction of quantitative whole-cell models. While kinetic network-type models have been established, it is also becoming possible to build physical, molecular-level models of cellular environments. This review outlines challenges in constructing and simulating such models and discusses near- and long-term opportunities for developing physical whole-cell models that can connect molecular structure with biological function.

Mechanism of Substrate Translocation in an Alternating Access Transporter.

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

A universal allosteric mechanism for G protein activation.

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gßγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gßγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

Cryo-EM reveals mechanistic insights into lipid-facilitated polyamine export by human ATP13A2.

The cytoplasmic polyamine maintains cellular homeostasis by chelating toxic metal cations, regulating transcriptional activity, and protecting DNA. ATP13A2 was identified as a lysosomal polyamine exporter responsible for polyamine release into the cytosol, and its dysfunction is associated with Alzheimer's disease and other neural degradation diseases. ATP13A2 belongs to the P5 subfamily of the P-type ATPase family, but its mechanisms remain unknown. Here, we report the cryoelectron microscopy (cryo-EM) structures of human ATP13A2 under four different conditions, revealing the structural coupling between the polyamine binding and the dephosphorylation. Polyamine is bound at the luminal tunnel and recognized through numerous electrostatic and π-cation interactions, explaining its broad specificity. The unique N-terminal domain is anchored to the lipid membrane to stabilize the E2P conformation, thereby accelerating the E1P-to-E2P transition. These findings reveal the distinct mechanism of P5B ATPases, thereby paving the way for neuroprotective therapy by activating ATP13A2.

Physical models of plant development.

The definition of shape in multicellular organisms is a major issue of developmental biology. It is well established that morphogenesis relies on genetic regulation. However, cells, tissues, and organism behaviors are also bound by the laws of physics, which limit the range of possible deformations organisms can undergo but also define what organisms must do to achieve specific shapes. Besides experiments, theoretical models and numerical simulations of growing tissues are powerful tools to investigate the link between genetic regulation and mechanics. Here, we provide an overview of the main mechanical models of plant morphogenesis developed so far, from subcellular scales to whole tissues. The common concepts and discrepancies between the various models are discussed.

A Zn-dependent structural transition of SOD1 modulates its ability to undergo phase separation.

The misfolding and mutation of Cu/Zn superoxide dismutase (SOD1) is commonly associated with amyotrophic lateral sclerosis (ALS). SOD1 can accumulate within stress granules (SGs), a type of membraneless organelle, which is believed to form via liquid-liquid phase separation (LLPS). Using wild-type, metal-deficient, and different ALS disease mutants of SOD1 and computer simulations, we report here that the absence of Zn leads to structural disorder within two loop regions of SOD1, triggering SOD1 LLPS and amyloid formation. The addition of exogenous Zn to either metal-free SOD1 or to the severe ALS mutation I113T leads to the stabilization of the loops and impairs SOD1 LLPS and aggregation. Moreover, partial Zn-mediated inhibition of LLPS was observed for another severe ALS mutant, G85R, which shows perturbed Zn-binding. By contrast, the ALS mutant G37R, which shows reduced Cu-binding, does not undergo LLPS. In addition, SOD1 condensates induced by Zn-depletion exhibit greater cellular toxicity than aggregates formed by prolonged incubation under aggregating conditions. Overall, our work establishes a role for Zn-dependent modulation of SOD1 conformation and LLPS properties that may contribute to amyloid formation.

Adhesion and shrinkage transform the rounded pupal horn into an angular adult horn in Japanese rhinoceros beetle.

Clarifying the mechanisms underlying shape alterations during insect metamorphosis is important for understanding exoskeletal morphogenesis. The large horn of the Japanese rhinoceros beetle Trypoxylus dichotomus is the result of drastic metamorphosis, wherein it appears as a rounded shape during pupation and then undergoes remodeling into an angular adult shape. However, the mechanical mechanisms underlying this remodeling process remain unknown. In this study, we investigated the remodeling mechanisms of the Japanese rhinoceros beetle horn by developing a physical simulation. We identified three factors contributing to remodeling by biological experiments - ventral adhesion, uneven shrinkage, and volume reduction - which were demonstrated to be crucial for transformation using a physical simulation. Furthermore, we corroborated our findings by applying the simulation to the mandibular remodeling of stag beetles. These results indicated that physical simulation applies to pupal remodeling in other beetles, and the morphogenic mechanism could explain various exoskeletal shapes.

Accreditation of analogue quantum simulators.

We present an accreditation protocol for analogue, i.e., continuous-time, quantum simulators. For a given simulation task, it provides an upper bound on the variation distance between the probability distributions at the output of an erroneous and error-free analogue quantum simulator. As its overheads are independent of the size and nature of the simulation, the protocol is ready for immediate usage and practical for the long term. It builds on the recent theoretical advances of strongly universal Hamiltonians and quantum accreditation as well as experimental progress toward the realization of programmable hybrid analogue-digital quantum simulators.

Unveiling the effect of Ni on the formation and structure of Earth's inner core.

Ni is the second most abundant element in the Earth's core. Yet, its effects on the inner core's structure and formation process are usually disregarded because of its electronic and size similarity with Fe. Using ab initio molecular dynamics simulations, we find that the bcc phase can spontaneously crystallize in liquid Ni at temperatures above Fe's melting point at inner core pressures. The melting temperature of Ni is shown to be 700 to 800 K higher than that of Fe at 323 to 360 GPa. hcp, bcc, and liquid phase relations differ for Fe and Ni. Ni can be a bcc stabilizer for Fe at high temperatures and inner core pressures. A small amount of Ni can accelerate Fe's crystallization at core pressures. These results suggest that Ni may substantially impact the structure and formation process of the solid inner core.

Atomic view of photosynthetic metabolite permeability pathways and confinement in synthetic carboxysome shells.

Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate ([Formula: see text]), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic [Formula: see text]-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to [Formula: see text] relative to CO2. To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO2-concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO2 assimilation in a crowded model carboxysome. Despite a relatively high CO2 permeability of approximately 10-2 cm/s across the carboxysome shell, the shell proteins reflect enough CO2 back toward rubisco that 2,650 CO2 molecules can be fixed by rubisco for every 1 CO2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.

Entropy drives the ligand recognition in G-protein-coupled receptor subtypes.

Achieving ligand subtype selectivity within highly homologous subtypes of G-protein-coupled receptor (GPCR) is critical yet challenging for GPCR drug discovery, primarily due to the unclear mechanism underlying ligand subtype selectivity, which hampers the rational design of subtype-selective ligands. Herein, we disclose an unusual molecular mechanism of entropy-driven ligand recognition in cannabinoid (CB) receptor subtypes, revealed through atomic-level molecular dynamics simulations, cryoelectron microscopy structure, and mutagenesis experiments. This mechanism is attributed to the distinct conformational dynamics of the receptor's orthosteric pocket, leading to variations in ligand binding entropy and consequently, differential binding affinities, which culminate in specific ligand recognition. We experimentally validated this mechanism and leveraged it to design ligands with enhanced or ablated subtype selectivity. One such ligand demonstrated favorable pharmacokinetic properties and significant efficacy in rodent inflammatory analgesic models. More importantly, it is precisely due to the high subtype selectivity obtained based on this mechanism that this ligand does not show addictive properties in animal models. Our findings elucidate the unconventional role of entropy in CB receptor subtype selectivity and suggest a strategy for rational design of ligands to achieve entropy-driven subtype selectivity for many pharmaceutically important GPCRs.

Impedance of nanocapacitors from molecular simulations to understand the dynamics of confined electrolytes.

Nanoelectrochemical devices have become a promising candidate technology across various applications, including sensing and energy storage, and provide new platforms for studying fundamental properties of electrode/electrolyte interfaces. In this work, we employ constant-potential molecular dynamics simulations to investigate the impedance of gold-aqueous electrolyte nanocapacitors, exploiting a recently introduced fluctuation-dissipation relation. In particular, we relate the frequency-dependent impedance of these nanocapacitors to the complex conductivity of the bulk electrolyte in different regimes, and use this connection to design simple but accurate equivalent circuit models. We show that the electrode/electrolyte interfacial contribution is essentially capacitive and that the electrolyte response is bulk-like even when the interelectrode distance is only a few nanometers, provided that the latter is sufficiently large compared to the Debye screening length. We extensively compare our simulation results with spectroscopy experiments and predictions from analytical theories. In contrast to experiments, direct access in simulations to the ionic and solvent contributions to the polarization allows us to highlight their significant and persistent anticorrelation and to investigate the microscopic origin of the timescales observed in the impedance spectrum. This work opens avenues for the molecular interpretation of impedance measurements, and offers valuable contributions for future developments of accurate coarse-grained representations of confined electrolytes.