Search

Abstract Voids—the nothingness—broadly exist within nanomaterials and impact properties ranging from catalysis to mechanical response. However, understanding nanovoids is challenging due to lack of imaging methods with the needed penetration depth and spatial resolution. Here, we integrate electron tomography, morphometry, graph theory and coarse-grained molecular dynamics simulation to study the formation of interconnected nanovoids in polymer films and their impacts on permeance and nanomechanical behaviour. Using polyamide membranes for molecular separation as a representative system, three-dimensional electron tomography at nanometre resolution reveals nanovoid formation from coalescence of oligomers, supported by coarse-grained molecular dynamics simulations. Void analysis provides otherwise inaccessible inputs for accurate fittings of methanol permeance for polyamide membranes. Three-dimensional structural graphs accounting for the tortuous nanovoids within, measure higher apparent moduli with polyamide membranes of higher graph rigidity. Our study elucidates the significance of nanovoids beyond the nothingness, impacting the synthesis‒morphology‒function relationships of complex nanomaterials.

Abstract The defining feature of Parkinson disease (PD) and Lewy body dementia (LBD) is the accumulation of alpha-synuclein (Asyn) fibrils in Lewy bodies and Lewy neurites. Here we develop and validate a method to amplify Asyn fibrils extracted from LBD postmortem tissue samples and use solid state nuclear magnetic resonance (SSNMR) studies to determine atomic resolution structure. Amplified LBD Asyn fibrils comprise a mixture of single protofilament and two protofilament fibrils with very low twist. The protofilament fold is highly similar to the fold determined by a recent cryo-electron microscopy study for a minority population of twisted single protofilament fibrils extracted from LBD tissue. These results expand the structural characterization of LBD Asyn fibrils and approaches for studying disease mechanisms, imaging agents and therapeutics targeting Asyn.

Abstract Metal-organic frameworks (MOFs) exhibit great promise for CO2 capture. However, finding the best performing materials poses computational and experimental grand challenges in view of the vast chemical space of potential building blocks. Here, we introduce GHP-MOFassemble, a generative artificial intelligence (AI), high performance framework for the rational and accelerated design of MOFs with high CO2 adsorption capacity and synthesizable linkers. GHP-MOFassemble generates novel linkers, assembled with one of three pre-selected metal nodes (Cu paddlewheel, Zn paddlewheel, Zn tetramer) into MOFs in a primitive cubic topology. GHP-MOFassemble screens and validates AI-generated MOFs for uniqueness, synthesizability, structural validity, uses molecular dynamics simulations to study their stability and chemical consistency, and crystal graph neural networks and Grand Canonical Monte Carlo simulations to quantify their CO2 adsorption capacities. We present the top six AI-generated MOFs with CO2 capacities greater than 2m mol g−1, i.e., higher than 96.9% of structures in the hypothetical MOF dataset.

Abstract Despite a vaccine, hepatitis B virus (HBV) remains a world-wide source of infections and deaths. We develop a whole-cell computational platform combining spatial and kinetic models describing the infection cycle of HBV in a hepatocyte host. We simulate key parts of the infection cycle with this whole-cell platform for 10 min of biological time, to predict infection progression, map out virus-host and virus-drug interactions. We find that starting from an established infection, decreasing the copy number of the viral envelope proteins shifts the dominant infection pathway from capsid secretion to re-importing the capsids into the nucleus, resulting in more nuclear-localized viral covalently closed circular DNA (cccDNA) and boosting transcription. This scenario can mimic the consequence of drugs designed to manipulate viral gene expression. Mutating capsid proteins facilitates capsid destabilization and disassembly at nuclear pore complexes, resulting in an increase in cccDNA copy number. However, excessive destabilization leads to premature cytoplasmic disassembly and does not increase the cccDNA counts. Finally, our simulations can predict the best drug dosage and its administration timing to reduce the cccDNA counts. Our adaptable computational platform can be parameterized to study other viruses and identify the most central viral pathways that can be targeted by drugs.

Abstract Here we used cryo-electron microscopy (cryo-EM), double electron-electron resonance spectroscopy (DEER), and molecular dynamics (MD) simulations, to capture and characterize ATP- and substrate-bound inward-facing (IF) and occluded (OC) conformational states of the heterodimeric ATP binding cassette (ABC) multidrug exporter BmrCD in lipid nanodiscs. Supported by DEER analysis, the structures reveal that ATP-powered isomerization entails changes in the relative symmetry of the BmrC and BmrD subunits that propagates from the transmembrane domain to the nucleotide binding domain. The structures uncover asymmetric substrate and Mg2+ binding which we hypothesize are required for triggering ATP hydrolysis preferentially in one of the nucleotide-binding sites. MD simulations demonstrate that multiple lipid molecules differentially bind the IF versus the OC conformation thus establishing that lipid interactions modulate BmrCD energy landscape. Our findings are framed in a model that highlights the role of asymmetric conformations in the ATP-coupled transport with general implications to the mechanism of ABC transporters.

Abstract ADP-ribosylation factor 1 (Arf1) interacts with multiple cellular partners and membranes to regulate intracellular traffic, organelle structure and actin dynamics. Defining the dynamic conformational landscape of Arf1 in its active form, when bound to the membrane, is of high functional relevance and key to understanding how Arf1 can alter diverse cellular processes. Through concerted application of nuclear magnetic resonance (NMR), neutron reflectometry (NR) and molecular dynamics (MD) simulations, we show that, while Arf1 is anchored to the membrane through its N-terminal myristoylated amphipathic helix, the G domain explores a large conformational space, existing in a dynamic equilibrium between membrane-associated and membrane-distal conformations. These configurational dynamics expose different interfaces for interaction with effectors. Interaction with the Pleckstrin homology domain of ASAP1, an Arf-GTPase activating protein (ArfGAP), restricts motions of the G domain to lock it in what seems to be a conformation exposing functionally relevant regions.

Aromatic amino acids in proteins support ligand binding and protein stability. To parse the physiocochemical roles of aromatic interactions, here Galles, Infield and co-authors identify pyrrolysine-based aminoacyl-tRNA synthetases that enable the encoding of fluorinated phenylalanine amino acids.

CLC-type channels selectively transport Cl− across biological membranes, but it is unclear how discrimination between anions is maintained. Here, authors use a combination of non-natural amino acid substitutions, electrophysiology, and molecular dynamics simulations to determine Cl− specificity within this family of ion channels.

Prestin, a motor protein, plays a major role in sound amplification. Using molecular dynamics simulations, the authors show that prestin causes membrane deformation patterns thereby achieving a particular lipid-mediated alignment in the membrane.

The Spns lipid transporters use ion gradients to drive substrate transport, including bioactive sphingolipids. Here, Dastvan et al. investigated how binding of protons powers the conformational changes that enable broad transport by a bacterial Spns.

Pant et al. describe the mechanism by which PIP2 might regulate homomeric Kv7.2 channels. They identify sites important in the binding of the PIP2 lipid to Kv7.2 channels and propose that the PIP2 binding to a specific site results in the coupling between the voltage sensor domain (VSD) and pore domain (PD), which stabilizes the open state of the channel.

Haloi et al. investigate the role of the hydrophobic N-terminus of hexokinase II (HKII) in their binding to mitochondrial VDAC1. They use a hybrid approach combining molecular and Brownian dynamics to show that HKII first inserts into the outer membrane of mitochondria (OMM) and then interacts with VDAC1 on the cytosolic surface of OMM to form a binary complex.

Endocannabinoids are ligands of cannabinoid receptors and a promising target for pain management. Here, the authors report a class of lipids formed by the epoxidation of N-arachidonoyl dopamine and N-arachidonoyl serotonin by cytochrome P450 epoxygenases, which reciprocally regulate canabinoid receptors and display anti-inflammatory activity.

Calmodulin (CaM) is a ubiquitous Ca2+ sensing protein that binds to and modulates numerous target proteins and enzymes during cellular signaling processes. A large number of CaM-target complexes have been identified and structurally characterized, revealing a wide diversity of CaM-binding modes. A newly identified target is creatine kinase (CK), a central enzyme in cellular energy homeostasis. This study reports two high-resolution X-ray structures, determined to 1.24 ​Å and 1.43 ​Å resolution, of calmodulin in complex with peptides from human brain and muscle CK, respectively. Both complexes adopt a rare extended binding mode with an observed stoichiometry of 1:2 CaM:peptide, confirmed by isothermal titration calorimetry, suggesting that each CaM domain independently binds one CK peptide in a Ca2+-depended manner. While the overall binding mode is similar between the structures with muscle or brain-type CK peptides, the most significant difference is the opposite binding orientation of the peptides in the N-terminal domain. This may extrapolate into distinct binding modes and regulation of the full-length CK isoforms. The structural insights gained in this study strengthen the link between cellular energy homeostasis and Ca2+-mediated cell signaling and may shed light on ways by which cells can ‘fine tune’ their energy levels to match the spatial and temporal demands.

Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system (CNS). Excitatory amino acid transporters (EAATs) regulate extracellular glutamate by transporting it into cells, mostly glia, to terminate neurotransmission and to avoid neurotoxicity. EAATs are also chloride (Cl–) channels, but the physiological role of Cl– conductance through EAATs is poorly understood. Mutations of human EAAT1 (hEAAT1) have been identified in patients with episodic ataxia type 6 (EA6). One mutation showed increased Cl– channel activity and decreased glutamate transport, but the relative contributions of each function of hEAAT1 to mechanisms underlying the pathology of EA6 remain unclear. Here we investigated the effects of 5 additional EA6-related mutations on hEAAT1 function in Xenopus laevis oocytes, and on CNS function in a Drosophila melanogaster model of locomotor behavior. Our results indicate that mutations resulting in decreased hEAAT1 Cl– channel activity but with functional glutamate transport can also contribute to the pathology of EA6, highlighting the importance of Cl– homeostasis in glial cells for proper CNS function. We also identified what we believe is a novel mechanism involving an ectopic sodium (Na+) leak conductance in glial cells. Together, these results strongly support the idea that EA6 is primarily an ion channelopathy of CNS glia.

XylE is a bacterial xylose transporter and homologue of human glucose transporters GLUTs 1-4. HDX-MS, mutagenesis and MD simulations suggest that protonation of a conserved aspartate triggers conformational transition from outward- to inward facing state only in the presence of substrate xylose. In contrast, inhibitor glucose locks the transporter in the outward facing state.

Formation of amyloid-beta (Aβ) oligomer pores in the membrane of neurons has been proposed to explain neurotoxicity in Alzheimer´s disease. Here authors present the 3D- structure of an Aβ oligomer formed in a membrane mimicking environment and observe that Aβ tetramers and octamers inserted into lipid bilayers as well-defined pores.

We introduce an end-to-end computational framework that allows for hyperparameter optimization using the DeepHyper library, accelerated model training, and interpretable AI inference. The framework is based on state-of-the-art AI models including CGCNN , PhysNet , SchNet , MPNN , MPNN-transformer , and TorchMD-NET . We employ these AI models along with the benchmark QM9 , hMOF , and MD17 datasets to showcase how the models can predict user-specified material properties within modern computing environments. We demonstrate transferable applications in the modeling of small molecules, inorganic crystals and nanoporous metal organic frameworks with a unified, standalone framework. We have deployed and tested this framework in the ThetaGPU supercomputer at the Argonne Leadership Computing Facility, and in the Delta supercomputer at the National Center for Supercomputing Applications to provide researchers with modern tools to conduct accelerated AI-driven discovery in leadership-class computing environments. We release these digital assets as open source scientific software in GitLab, and ready-to-use Jupyter notebooks in Google Colab.

Secretory (S) Immunoglobulin (Ig) A is the predominant mucosal antibody, which binds pathogens and commensal microbes. SIgA is a polymeric antibody, typically containing two copies of IgA that assemble with one joining-chain (JC) to form dimeric (d) IgA that is bound by the polymeric Ig-receptor ectodomain, called secretory component (SC). Here, we report the cryo-electron microscopy structures of murine SIgA and dIgA. Structures reveal two IgAs conjoined through four heavy-chain tailpieces and the JC that together form a β-sandwich-like fold. The two IgAs are bent and tilted with respect to each other, forming distinct concave and convex surfaces. In SIgA, SC is bound to one face, asymmetrically contacting both IgAs and JC. The bent and tilted arrangement of complex components limits the possible positions of both sets of antigen-binding fragments (Fabs) and preserves steric accessibility to receptor-binding sites, likely influencing antigen binding and effector functions.

Citrate, α-ketoglutarate and succinate are TCA cycle intermediates that also play essential roles in metabolic signaling and cellular regulation. These di- and tricarboxylates are imported into the cell by the divalent anion sodium symporter (DASS) family of plasma membrane transporters, which contains both cotransporters and exchangers. While DASS proteins transport substrates via an elevator mechanism, to date structures are only available for a single DASS cotransporter protein in a substrate-bound, inward-facing state. We report multiple cryo-EM and X-ray structures in four different states, including three hitherto unseen states, along with molecular dynamics simulations, of both a cotransporter and an exchanger. Comparison of these outward- and inward-facing structures reveal how the transport domain translates and rotates within the framework of the scaffold domain through the transport cycle. Additionally, we propose that DASS transporters ensure substrate coupling by a charge-compensation mechanism, and by structural changes upon substrate release.

Xenopus oocytes expressing human aquaporin-7 (AQP7) exhibit greater osmotic water permeability and 3H-glycerol uptake vs. those expressing the bacterial glycerol facilitator GlpF. AQP7-expressing oocytes exposed to increasing extracellular [glycerol] under isosmolal conditions exhibit increasing swelling rates, whereas GlpF-expressing oocytes do not swell at all. To provide a structural basis for these observed physiological differences, we performed X-ray crystallographic structure determination of AQP7 and molecular-dynamics simulations on AQP7 and GlpF. The structure reveals AQP7 tetramers containing two monomers with 3 glycerols, and two monomers with 2 glycerols in the pore. In contrast to GlpF, no glycerol is bound at the AQP7 selectivity filter (SF), comprising residues F74, G222, Y223, and R229. The AQP7 SF is resolved in its closed state because F74 blocks the passage of small solutes. Molecular dynamics simulations demonstrate that F74 undergoes large and rapid conformational changes, allowing glycerol molecules to permeate without orientational restriction. The more rigid GlpF imposes orientational constraints on glycerol molecules passing through the SF. Moreover, GlpF-W48 (analogous to AQP7-F74) undergoes rare but long-lasting conformational changes that block the pore to H2O and glycerol.

Among coupled exchangers, CLCs uniquely catalyze the exchange of oppositely charged ions (Cl– for H+). Transport-cycle models to describe and explain this unusual mechanism have been proposed based on known CLC structures. While the proposed models harmonize with many experimental findings, gaps and inconsistencies in our understanding have remained. One limitation has been that global conformational change – which occurs in all conventional transporter mechanisms – has not been observed in any high-resolution structure. Here, we describe the 2.6 Å structure of a CLC mutant designed to mimic the fully H+-loaded transporter. This structure reveals a global conformational change to improve accessibility for the Cl– substrate from the extracellular side and new conformations for two key glutamate residues. Together with DEER measurements, MD simulations, and functional studies, this new structure provides evidence for a unified model of H+/Cl– transport that reconciles existing data on all CLC-type proteins.

The helical transmembrane segments of membrane proteins play central roles in sensory transduction but the mechanistic basis for their function remains unresolved. Here the authors identify regions in the S4 voltage-sensing segment of Shaker potassium channels where local helical structure is reliant upon backbone amide support.

Secondary transporters catalyse substrate translocation across the cell membrane but the role of lipids during the transport cycle remains unclear. Here authors used hydrogen-deuterium exchange mass spectrometry and molecular dynamics simulations to understand how lipids regulate the conformational dynamics of secondary transporters.

Acquired resistance to endocrine therapy remains a significant clinical burden for breast cancer patients. Somatic mutations in the ESR1 (estrogen receptor alpha (ERα)) gene ligand-binding domain (LBD) represent a recognized mechanism of acquired resistance. Antiestrogens with improved efficacy versus tamoxifen might overcome the resistant phenotype in ER +breast cancers. Bazedoxifene (BZA) is a potent antiestrogen that is clinically approved for use in hormone replacement therapies. We found that BZA possesses improved inhibitory potency against the Y537S and D538G ERα mutants compared to tamoxifen and has additional inhibitory activity in combination with the CDK4/6 inhibitor palbociclib. In addition, comprehensive biophysical and structural biology studies show BZA’s selective estrogen receptor degrading (SERD) properties that override the stabilizing effects of the Y537S and D538G ERα mutations.

From bacteria to mammals, different phospholipid species are segregated between the inner and outer leaflets of the plasma membrane by ATP-dependent lipid transporters. Disruption of this asymmetry by ATP-independent phospholipid scrambling is important in cellular signaling, but its mechanism remains incompletely understood. Using MD simulations coupled with experimental assays, we show that the surface hydrophilic transmembrane cavity exposed to the lipid bilayer on the fungal scramblase nhTMEM16 serves as the pathway for both lipid translocation and ion conduction across the membrane. Ca2+ binding stimulates its open conformation by altering the structure of transmembrane helices that line the cavity. We have identified key amino acids necessary for phospholipid scrambling and validated the idea that ions permeate TMEM16 Cl- channels via a structurally homologous pathway by showing that mutation of two residues in the pore region of the TMEM16A Ca2+-activated Cl- channel convert it into a robust scramblase.

Mitochondria are the key source of ATP that fuels cellular functions, and they are also central in cellular signaling, cell division and apoptosis. Dysfunction of mitochondria has been implicated in a wide range of diseases, including neurodegenerative and cardiac diseases, and various types of cancer. One of the key proteins that regulate mitochondrial function is the voltage-dependent anion channel 1 (VDAC1), the most abundant protein on the outer membrane of mitochondria. VDAC1 is the gatekeeper for the passages of metabolites, nucleotides, and ions; it plays a crucial role in regulating apoptosis due to its interaction with apoptotic and anti-apoptotic proteins, namely members of the Bcl-2 family of proteins and hexokinase. Therefore, regulation of VDAC1 is crucial not only for metabolic functions of mitochondria, but also for cell survival. In fact, multiple lines of evidence have confirmed the involvement of VDAC1 in several diseases. Consequently, modulation or dysregulation of VDAC1 function can potentially attenuate or exacerbate pathophysiological conditions. Understanding the role of VDAC1 in health and disease could lead to selective protection of cells in different tissues and diverse diseases. The purpose of this review is to discuss the role of VDAC1 in the pathogenesis of diseases and as a potentially effective target for therapeutic management of various pathologies.

Somatic mutations in the estrogen receptor alpha (ERα) gene (ESR1), especially Y537S and D538G, have been linked to acquired resistance to endocrine therapies. Cell-based studies demonstrated that these mutants confer ERα constitutive activity and antiestrogen resistance and suggest that ligand-binding domain dysfunction leads to endocrine therapy resistance. Here, we integrate biophysical and structural biology data to reveal how these mutations lead to a constitutively active and antiestrogen-resistant ERα. We show that these mutant ERs recruit coactivator in the absence of hormone while their affinities for estrogen agonist (estradiol) and antagonist (4-hydroxytamoxifen) are reduced. Further, they confer antiestrogen resistance by altering the conformational dynamics of the loop connecting Helix 11 and Helix 12 in the ligand-binding domain of ERα, which leads to a stabilized agonist state and an altered antagonist state that resists inhibition.

CLC secondary active transporters exchange Cl- for H+. Crystal structures have suggested that the conformational change from occluded to outward-facing states is unusually simple, involving only the rotation of a conserved glutamate (Gluex) upon its protonation. Using 19F NMR, we show that as [H+] is increased to protonate Gluex and enrich the outward-facing state, a residue ~20 Å away from Gluex, near the subunit interface, moves from buried to solvent-exposed. Consistent with functional relevance of this motion, constriction via inter-subunit cross-linking reduces transport. Molecular dynamics simulations indicate that the cross-link dampens extracellular gate-opening motions. In support of this model, mutations that decrease steric contact between Helix N (part of the extracellular gate) and Helix P (at the subunit interface) remove the inhibitory effect of the cross-link. Together, these results demonstrate the formation of a previously uncharacterized 'outward-facing open' state, and highlight the relevance of global structural changes in CLC function.

Multidrug ATP binding cassette (ABC) exporters are ubiquitous ABC transporters that extrude cytotoxic molecules across cell membranes. Despite recent progress in structure determination of these transporters, the conformational motion that transduces the energy of ATP hydrolysis to the work of substrate translocation remains undefined. Here, we have investigated the conformational cycle of BmrCD, a representative of the heterodimer family of ABC exporters that have an intrinsically impaired nucleotide binding site. We measured distances between pairs of spin labels monitoring the movement of the nucleotide binding (NBD) and transmembrane domains (TMD). The results expose previously unobserved structural intermediates of the NBDs arising from asymmetric configuration of catalytically inequivalent nucleotide binding sites. The two-state transition of the TMD, from an inward- to an outward-facing conformation, is driven exclusively by ATP hydrolysis. These findings provide direct evidence of divergence in the mechanism of ABC exporters.

The leucine transporter (LeuT) has recently commanded exceptional attention due mainly to two distinctions; it provides the only crystal structures available for a protein homologous to the pharmacologically relevant neurotransmitter: sodium symporters (NSS), and, it exhibits a hallmark 5-TM inverted repeat ("LeuT-fold"), a fold recently discovered to also exist in several secondary transporter families, underscoring its general role in transporter function. Constructing the transport cycle of "LeuT-fold" transporters requires detailed structural and dynamic descriptions of the outward-facing (OF) and inward-facing (IF) states, as well as the intermediate states. To this end, we have modeled the structurally unknown IF state of LeuT, based on the known crystal structures of the OF state of LeuT and the IF state of vSGLT, a "LeuT-fold" transporter. The detailed methodology developed for the study combines structure-based alignment, threading, targeted MD and equilibrium MD, and can be applied to other proteins. The resulting IF-state models maintain the secondary structural features of LeuT. Water penetration and solvent accessibility calculations show that TM1, TM3, TM6 and TM8 line the substrate binding/unbinding pathway with TM10 and its pseudosymmetric partner, TM5, participating in the extracellular and intracellular halves of the lumen, respectively. We report conformational hotspots where notable changes in interactions occur between the IF and OF states. We observe Na2 exiting the LeuT-substrate- complex in the IF state, mainly due to TM1 bending. Inducing a transition in only one of the two pseudosymmetric domains, while allowing the second to respond dynamically, is found to be sufficient to induce the formation of the IF state. We also propose that TM2 and TM7 may be facilitators of TM1 and TM6 motion. Thus, this study not only presents a novel modeling methodology applied to obtain the IF state of LeuT, but also describes structural elements involved in a possibly general transport mechanism in transporters adopting the "LeuT-fold".