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.

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.

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.

Molecular dynamics-based identification of binding pathways and two distinct high-affinity sites for succinate in succinate receptor 1/GPR91.

SUCNR1 is an auto- and paracrine sensor of the metabolic stress signal succinate. Using unsupervised molecular dynamics (MD) simulations (170.400 ns) and mutagenesis across human, mouse, and rat SUCNR1, we characterize how a five-arginine motif around the extracellular pole of TM-VI determines the initial capture of succinate in the extracellular vestibule (ECV) to either stay or move down to the orthosteric site. Metadynamics demonstrate low-energy succinate binding in both sites, with an energy barrier corresponding to an intermediate stage during which succinate, with an associated water cluster, unlocks the hydrogen-bond-stabilized conformationally constrained extracellular loop (ECL)-2b. Importantly, simultaneous binding of two succinate molecules through either a "sequential" or "bypassing" mode is a frequent endpoint. The mono-carboxylate NF-56-EJ40 antagonist enters SUCNR1 between TM-I and -II and does not unlock ECL-2b. It is proposed that occupancy of both high-affinity sites is required for selective activation of SUCNR1 by high local succinate concentrations.

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.

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.

The structure of a C. neoformans polysaccharide motif recognized by protective antibodies: A combined NMR and MD study.

Cryptococcus neoformans is a fungal pathogen responsible for cryptococcosis and cryptococcal meningitis. The C. neoformans' capsular polysaccharide and its shed exopolysaccharide function both as key virulence factors and to protect the fungal cell from phagocytosis. Currently, a glycoconjugate of these polysaccharides is being explored as a vaccine to protect against C. neoformans infection. In this study, NOE and J-coupling values from NMR experiments were consistent with a converged structure of the synthetic decasaccharide, GXM10-Ac3, calculated from MD simulations. GXM10-Ac3 was designed as an extension of glucuronoxylomannan (GXM) polysaccharide motif (M2) which is common in the clinically predominant serotype A strains and is recognized by protective forms of GXM-specific monoclonal antibodies. The M2 motif is a hexasaccharide with a three-residue α-mannan backbone, modified by ß-(1→2)-xyloses (Xyl) on the first two mannoses (Man) and a ß-(1→2)-glucuronic acid (GlcA) on the third Man. Combined NMR and MD analyses reveal that GXM10-Ac3 adopts an extended structure, with Xyl/GlcA branches alternating sides along the α-mannan backbone. O-acetyl esters also alternate sides and are grouped in pairs. MD analysis of a twelve M2-repeating unit polymer supports the notion that the GXM10-Ac3 structure is uniformly represented throughout the polysaccharide. This derived GXM model displays high flexibility while maintaining a structural identity, yielding insights to further explore intermolecular interactions between polysaccharides, interactions with anti-GXM mAbs, and the cryptococcal polysaccharide architecture.

Conformational dynamics underlying atypical chemokine receptor 3 activation.

Atypical Chemokine Receptor 3 (ACKR3) belongs to the G protein-coupled receptor family but it does not signal through G proteins. The structural properties that govern the functional selectivity and the conformational dynamics of ACKR3 activation are poorly understood. Here, we combined hydrogen/deuterium exchange mass spectrometry, site-directed mutagenesis, and molecular dynamics simulations to examine the binding mode and mechanism of action of ACKR3 ligands of different efficacies. Our results show that activation or inhibition of ACKR3 is governed by intracellular conformational changes of helix 6, intracellular loop 2, and helix 7, while the DRY motif becomes protected during both processes. Moreover, we identified the binding sites and the allosteric modulation of ACKR3 upon ß-arrestin 1 binding. In summary, this study highlights the structure-function relationship of small ligands, the binding mode of ß-arrestin 1, the activation dynamics, and the atypical dynamic features in ACKR3 that may contribute to its inability to activate G proteins.

Enhancing protein dynamics analysis with hydrophilic polyethylene glycol cross-linkers.

Cross-linkers play a critical role in capturing protein dynamics in chemical cross-linking mass spectrometry techniques. Various types of cross-linkers with different backbone features are widely used in the study of proteins. However, it is still not clear how the cross-linkers' backbone affect their own structure and their interactions with proteins. In this study, we systematically characterized and compared methylene backbone and polyethylene glycol (PEG) backbone cross-linkers in terms of capturing protein structure and dynamics. The results indicate the cross-linker with PEG backbone have a better ability to capture the inter-domain dynamics of calmodulin, adenylate kinase, maltodextrin binding protein and dual-specificity protein phosphatase. We further conducted quantum chemical calculations and all-atom molecular dynamics simulations to analyze thermodynamic and kinetic properties of PEG backbone and methylene backbone cross-linkers. Solution nuclear magnetic resonance was employed to validate the interaction interface between proteins and cross-linkers. Our findings suggest that the polarity distribution of PEG backbone enhances the accessibility of the cross-linker to the protein surface, facilitating the capture of sites located in dynamic regions. By comprehensively benchmarking with disuccinimidyl suberate (DSS)/bis-sulfosuccinimidyl-suberate(BS3), bis-succinimidyl-(PEG)2 revealed superior advantages in protein dynamic conformation analysis in vitro and in vivo, enabling the capture of a greater number of cross-linking sites and better modeling of protein dynamics. Furthermore, our study provides valuable guidance for the development and application of PEG backbone cross-linkers.

Reconstruction of LPS Transfer Cascade Reveals Structural Determinants within LBP, CD14, and TLR4-MD2 for Efficient LPS Recognition and Transfer.

Lipopolysaccharide (LPS), the major component of the outer membrane of Gram-negative bacteria, binds Toll-like receptor 4 (TLR4)-MD2 complex and activates innate immune responses. LPS transfer to TLR4-MD2 is catalyzed by both LPS binding protein (LBP) and CD14. To define the sequential molecular interactions underlying this transfer, we reconstituted in vitro the entire LPS transfer process from LPS micelles to TLR4-MD2. Using electron microscopy and single-molecule approaches, we characterized the dynamic intermediate complexes for LPS transfer: LBP-LPS micelles, CD14-LBP-LPS micelle, and CD14-LPS-TLR4-MD2 complex. A single LBP molecule bound longitudinally to LPS micelles catalyzed multi-rounds of LPS transfer to CD14s that rapidly dissociated from LPB-LPS complex upon LPS transfer via electrostatic interactions. Subsequently, the single LPS molecule bound to CD14 was transferred to TLR4-MD2 in a TLR4-dependent manner. The definition of the structural determinants of the LPS transfer cascade to TLR4 may enable the development of targeted therapeutics for intervention in LPS-induced sepsis.

The birth and evolution of solvated electrons in the water.

The photo-induced radiolysis of water is an elementary reaction in biology and chemistry, forming solvated electrons, OH radicals, and hydronium cations on fast time scales. Here, we use an optical-pump terahertz-probe spectroscopy setup to trigger the photoionization of water molecules with optical laser pulses at ~400 nm and then time-resolve the transient solvent response with broadband terahertz (THz) fields with a ~90 fs time resolution. We observe three distinct spectral responses. The first is a positive broadband mode that can be attributed to an initial diffuse, delocalized electron with a radius of (22 ± 1) Å, which is short lived (<200 fs) because the absorption is blue-shifting outside of the THz range. The second emerging spectroscopic signature with a lifetime of about 150 ps is attributed to an intermolecular mode associated with a mass rearrangement of solvent molecules due to charge separation of radicals and hydronium cations. After 0.2 ps, we observe a long-lasting THz signature with depleted intensity at 110 cm-1 that is well reproduced by ab initio molecular dynamics. We interpret this negative band at 110 cm-1 as the solvent cage characterized by a weakening of the hydrogen bond network in the first and second hydration shells of the cavity occupied by the localized electron.

Pressure pushes tRNALys3 into excited conformational states.

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

Disulfiram blocks inflammatory TLR4 signaling by targeting MD-2.

Toll-like receptor 4 (TLR4) sensing of lipopolysaccharide (LPS), the most potent pathogen-associated molecular pattern of gram-negative bacteria, activates NF-κB and Irf3, which induces inflammatory cytokines and interferons that trigger an intense inflammatory response, which is critical for host defense but can also cause serious inflammatory pathology, including sepsis. Although TLR4 inhibition is an attractive therapeutic approach for suppressing overexuberant inflammatory signaling, previously identified TLR4 antagonists have not shown any clinical benefit. Here, we identify disulfiram (DSF), an FDA-approved drug for alcoholism, as a specific inhibitor of TLR4-mediated inflammatory signaling. TLR4 cell surface expression, LPS sensing, dimerization and signaling depend on TLR4 binding to MD-2. DSF and other cysteine-reactive drugs, previously shown to block LPS-triggered inflammatory cell death (pyroptosis), inhibit TLR4 signaling by covalently modifying Cys133 of MD-2, a key conserved residue that mediates TLR4 sensing and signaling. DSF blocks LPS-triggered inflammatory cytokine, chemokine, and interferon production by macrophages in vitro. In the aggressive N-methyl-4-phenyl-1,2,3,6-tetrahydropyridine mouse model of Parkinson's disease (PD) in which TLR4 plays an important role, DSF markedly suppresses neuroinflammation and dopaminergic neuron loss, and restores motor function. Our findings identify a role for DSF in curbing TLR4-mediated inflammation and suggest that DSF and other drugs that target MD-2 might be useful for treating PD and other diseases in which inflammation contributes importantly to pathogenesis.

How acidic amino acid residues facilitate DNA target site selection.

Despite the negative charge of the DNA backbone, acidic residues (Asp/Glu) commonly participate in the base readout, with a strong preference for cytosine. In fact, in the solved DNA/protein structures, cytosine is recognized almost exclusively by Asp/Glu through a direct hydrogen bond, while at the same time, adenine, regardless of its amino group, shows no propensity for Asp/Glu. Here, we analyzed the contribution of Asp/Glu to sequence-specific DNA binding using classical and ab initio simulations of selected transcription factors and found that it is governed by a fine balance between the repulsion from backbone phosphates and attractive interactions with cytosine. Specifically, Asp/Glu lower the affinity for noncytosine sites and thus act as negative selectors preventing off-target binding. At cytosine-containing sites, the favorable contribution does not merely rely on the formation of a single H-bond but usually requires the presence of positive potential generated by multiple cytosines, consistently with the observed excess of cytosine in the target sites. Finally, we show that the preference of Asp/Glu for cytosine over adenine is a result of the repulsion from the adenine imidazole ring and a tendency of purine-purine dinucleotides to adopt the BII conformation.

Road-blocker HSP disease mutation disrupts pre-organization for ATP hydrolysis in kinesin through a second sphere control.

Kinesin motor proteins perform several essential cellular functions powered by the adenosine triphosphate (ATP) hydrolysis reaction. Several single-point mutations in the kinesin motor protein KIF5A have been implicated to hereditary spastic paraplegia disease (HSP), a lethal neurodegenerative disease in humans. In earlier studies, we have shown that a series of HSP-related mutations can impair the kinesin's long-distance displacement or processivity by modulating the order-disorder transition of the linker connecting the heads to the coiled coil. On the other hand, the reduction of kinesin's ATP hydrolysis reaction rate by a distal asparagine-to-serine mutation is also known to cause HSP disease. However, the molecular mechanism of the ATP hydrolysis reaction in kinesin by this distal mutation is still not fully understood. Using classical molecular dynamics simulations combined with quantum mechanics/molecular mechanics calculations, the pre-organization geometry required for optimal hydrolysis in kinesin motor bound to α/ß-tubulin is determined. This optimal geometry has only a single salt-bridge (of the possible two) between Arg203-Glu236, putting a reactive water molecule at a perfect position for hydrolysis. Such geometry is also needed to create the appropriate configuration for proton translocation during ATP hydrolysis. The distal asparagine-to-serine mutation is found to disrupt this optimal geometry. Therefore, the current study along with our previous one demonstrates how two different effects on kinesin dynamics (processivity and ATP hydrolysis), caused by a different set of genotypes, can give rise to the same phenotype leading to HSP disease.

Specific residues and conformational plasticity define the substrate specificity of short-chain dehydrogenases/reductases.

Short-chain dehydrogenases/reductases (SDRs) are one of the most prevalent enzyme families distributed among the sequenced microorganisms. Despite the presence of a conserved catalytic tetrad and high structural similarity, these enzymes exhibit different substrate specificities. The insufficient knowledge regarding the amino acids underlying substrate specificity hinders the understanding of the SDRs' roles in diverse and significant biological processes. Here, we performed bioinformatic analysis, molecular modeling, and mutagenesis studies to identify the key residues that regulate the substrate specificities of two homologous microbial SDRs (i.e., DesE and KduD). Further, we investigated the impact of altering the physicochemical properties of these amino acids on enzyme activity. Interestingly, molecular dynamics simulations also suggest a critical role of enzyme conformational flexibility in substrate recognition and catalysis. Overall, our findings improve the understanding of microbial SDR substrate specificity and shed light on future rational design of more efficient and effective biocatalysts.

Helical domain of hGBP3 cannot stimulate the second phosphate cleavage of GTP.

Interferon-gamma-inducible large GTPases, hGBPs, possess antipathogenic and antitumor activities in human cells. Like hGBP1, its closest homolog, hGBP3 has two domains; an N-terminal catalytic domain and a C-terminal helical domain, connected by an intermediate region. The biochemical function of this protein and the role of its domains in substrate hydrolysis have not yet been investigated. Here, we report that while hGBP3 can produce both GDP and GMP, GMP is the minor product, 30% (unlike 85% in hGBP1), indicating that hGBP3 is unable to produce enhanced GMP. To understand which domain(s) are responsible for this deficiency, we created hGBP3 truncated variants. Surprisingly, GMP production was similar upon deletion of the helical domain, suggesting that in contrast to hGBP1, the helical domain of hGBP3 cannot stimulate the second phosphate cleavage of GTP. We conducted computational and solution studies to understand the underlying basis. We found that the regulatory residue W79, present in the catalytic domain, forms an H-bond with the backbone carbonyl of K76 (located in the catalytic loop) of the substrate-bound hGBP3. However, after gamma-phosphate cleavage of GTP, the W79-containing region does not undergo a conformational change, failing to redirect the catalytic loop toward the beta-phosphate. This is necessary for efficient GMP formation because hGBP homologs utilize the same catalytic residue for both phosphate cleavages. We suggest that the lack of specific interdomain contacts mediated by the helical domain prevents the catalytic loop movement, resulting in reduced GMP formation. These findings may provide insight into how hGBP3 contributes to immunity.

Interaction networks within disease-associated GαS variants characterized by an integrative biophysical approach.

Activation of G proteins through nucleotide exchange initiates intracellular signaling cascades essential for life processes. Under normal conditions, nucleotide exchange is regulated by the formation of G protein-G protein-coupled receptor complexes. Single point mutations in the Gα subunit of G proteins bypass this interaction, leading to loss of function or constitutive gain of function, which is closely linked with the onset of multiple diseases. Despite the recognized significance of Gα mutations in disease pathology, structural information for most variants is lacking, potentially due to inherent protein dynamics that pose challenges for crystallography. To address this, we leveraged an integrative spectroscopic and computational approach to structurally characterize seven of the most frequently observed and clinically relevant mutations in the stimulatory Gα subunit, GαS. A previously proposed allosteric model of Gα activation linked structural changes in the nucleotide-binding pocket with functionally important changes in interactions between switch regions. We investigated this allosteric connection in GαS by integrating data from variable temperature CD spectroscopy, which measured changes in global protein structure and stability, and molecular dynamics simulations, which observed changes in interaction networks between GαS switch regions. Additionally, saturation-transfer difference NMR spectroscopy was applied to observe changes in nucleotide interactions with residues within the nucleotide binding site. These data have enabled testing of predictions regarding how mutations in GαS result in loss or gain of function and evaluation of proposed structural mechanisms. The integration of experimental and computational data allowed us to propose a more nuanced classification of mechanisms underlying GαS gain-of-function and loss-of-function mutations.

Recovery from desensitization in GluA2 AMPA receptors is affected by a single mutation in the N-terminal domain interface.

AMPA-type ionotropic glutamate receptors (AMPARs) are central to various neurological processes, including memory and learning. They assemble as homo- or heterotetramers of GluA1, GluA2, GluA3, and GluA4 subunits, each consisting of an N-terminal domain (NTD), a ligand-binding domain, a transmembrane domain, and a C-terminal domain. While AMPAR gating is primarily controlled by reconfiguration in the ligand-binding domain layer, our study focuses on the NTDs, which also influence gating, yet the underlying mechanism remains enigmatic. In this investigation, we employ molecular dynamics simulations to evaluate the NTD interface strength in GluA1, GluA2, and NTD mutants GluA2-H229N and GluA1-N222H. Our findings reveal that GluA1 has a significantly weaker NTD interface than GluA2. The NTD interface of GluA2 can be weakened by a single point mutation in the NTD dimer-of-dimer interface, namely H229N, which renders GluA2 more GluA1-like. Electrophysiology recordings demonstrate that this mutation also leads to slower recovery from desensitization. Moreover, we observe that lowering the pH induces more splayed NTD states and enhances desensitization in GluA2. We hypothesized that H229 was responsible for this pH sensitivity; however, GluA2-H229N was also affected by pH, meaning that H229 is not solely responsible and that protons exert their effect across multiple domains of the AMPAR. In summary, our work unveils an allosteric connection between the NTD interface strength and AMPAR desensitization.

Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling.

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work, we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi, and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.