Rapid simulation of glycoprotein structures by grafting and steric exclusion of glycan conformer libraries.

Most membrane proteins are modified by covalent addition of complex sugars through N- and O-glycosylation. Unlike proteins, glycans do not typically adopt specific secondary structures and remain very mobile, shielding potentially large fractions of protein surface. High glycan conformational freedom hinders complete structural elucidation of glycoproteins. Computer simulations may be used to model glycosylated proteins but require hundreds of thousands of computing hours on supercomputers, thus limiting routine use. Here, we describe GlycoSHIELD, a reductionist method that can be implemented on personal computers to graft realistic ensembles of glycan conformers onto static protein structures in minutes. Using molecular dynamics simulation, small-angle X-ray scattering, cryoelectron microscopy, and mass spectrometry, we show that this open-access toolkit provides enhanced models of glycoprotein structures. Focusing on N-cadherin, human coronavirus spike proteins, and gamma-aminobutyric acid receptors, we show that GlycoSHIELD can shed light on the impact of glycans on the conformation and activity of complex glycoproteins.

Structural mechanisms of α7 nicotinic receptor allosteric modulation and activation.

The α7 nicotinic acetylcholine receptor is a pentameric ligand-gated ion channel that plays an important role in cholinergic signaling throughout the nervous system. Its unique physiological characteristics and implications in neurological disorders and inflammation make it a promising but challenging therapeutic target. Positive allosteric modulators overcome limitations of traditional α7 agonists, but their potentiation mechanisms remain unclear. Here, we present high-resolution structures of α7-modulator complexes, revealing partially overlapping binding sites but varying conformational states. Structure-guided functional and computational tests suggest that differences in modulator activity arise from the stable rotation of a channel gating residue out of the pore. We extend the study using a time-resolved cryoelectron microscopy (cryo-EM) approach to reveal asymmetric state transitions for this homomeric channel and also find that a modulator with allosteric agonist activity exploits a distinct channel-gating mechanism. These results define mechanisms of α7 allosteric modulation and activation with implications across the pentameric receptor superfamily.

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.

Angiotensin Analogs with Divergent Bias Stabilize Distinct Receptor Conformations.

"Biased" G protein-coupled receptor (GPCR) agonists preferentially activate pathways mediated by G proteins or ß-arrestins. Here, we use double electron-electron resonance spectroscopy to probe the changes that ligands induce in the conformational distribution of the angiotensin II type I receptor. Monitoring distances between 10 pairs of nitroxide labels distributed across the intracellular regions enabled mapping of four underlying sets of conformations. Ligands from different functional classes have distinct, characteristic effects on the conformational heterogeneity of the receptor. Compared to angiotensin II, the endogenous agonist, agonists with enhanced Gq coupling more strongly stabilize an "open" conformation with an accessible transducer-binding site. ß-arrestin-biased agonists deficient in Gq coupling do not stabilize this open conformation but instead favor two more occluded conformations. These data suggest a structural mechanism for biased ligand action at the angiotensin receptor that can be exploited to rationally design GPCR-targeting drugs with greater specificity of action.

Isolation of a Structural Mechanism for Uncoupling T Cell Receptor Signaling from Peptide-MHC Binding.

TCR-signaling strength generally correlates with peptide-MHC binding affinity; however, exceptions exist. We find high-affinity, yet non-stimulatory, interactions occur with high frequency in the human T cell repertoire. Here, we studied human TCRs that are refractory to activation by pMHC ligands despite robust binding. Analysis of 3D affinity, 2D dwell time, and crystal structures of stimulatory versus non-stimulatory TCR-pMHC interactions failed to account for their different signaling outcomes. Using yeast pMHC display, we identified peptide agonists of a formerly non-responsive TCR. Single-molecule force measurements demonstrated the emergence of catch bonds in the activating TCR-pMHC interactions, correlating with exclusion of CD45 from the TCR-APC contact site. Molecular dynamics simulations of TCR-pMHC disengagement distinguished agonist from non-agonist ligands based on the acquisition of catch bonds within the TCR-pMHC interface. The isolation of catch bonds as a parameter mediating the coupling of TCR binding and signaling has important implications for TCR and antigen engineering for immunotherapy.

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.

Lipid-Sorting Specificity Encoded in K-Ras Membrane Anchor Regulates Signal Output.

K-Ras is targeted to the plasma membrane by a C-terminal membrane anchor that comprises a farnesyl-cysteine-methyl-ester and a polybasic domain. We used quantitative spatial imaging and atomistic molecular dynamics simulations to examine molecular details of K-Ras plasma membrane binding. We found that the K-Ras anchor binds selected plasma membrane anionic lipids with defined head groups and lipid side chains. The precise amino acid sequence and prenyl group define a combinatorial code for lipid binding that extends beyond simple electrostatics; within this code lysine and arginine residues are non-equivalent and prenyl chain length modifies nascent polybasic domain lipid preferences. The code is realized by distinct dynamic tertiary structures of the anchor on the plasma membrane that govern amino acid side-chain-lipid interactions. An important consequence of this specificity is the ability of such anchors when aggregated to sort subsets of phospholipids into nanoclusters with defined lipid compositions that determine K-Ras signaling output.

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.

Structural and Functional Analysis of a ß2-Adrenergic Receptor Complex with GRK5.

The phosphorylation of agonist-occupied G-protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) functions to turn off G-protein signaling and turn on arrestin-mediated signaling. While a structural understanding of GPCR/G-protein and GPCR/arrestin complexes has emerged in recent years, the molecular architecture of a GPCR/GRK complex remains poorly defined. We used a comprehensive integrated approach of cross-linking, hydrogen-deuterium exchange mass spectrometry (MS), electron microscopy, mutagenesis, molecular dynamics simulations, and computational docking to analyze GRK5 interaction with the ß2-adrenergic receptor (ß2AR). These studies revealed a dynamic mechanism of complex formation that involves large conformational changes in the GRK5 RH/catalytic domain interface upon receptor binding. These changes facilitate contacts between intracellular loops 2 and 3 and the C terminus of the ß2AR with the GRK5 RH bundle subdomain, membrane-binding surface, and kinase catalytic cleft, respectively. These studies significantly contribute to our understanding of the mechanism by which GRKs regulate the function of activated GPCRs. PAPERCLIP.

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.

The solute carrier SPNS2 recruits PI(4,5)P2 to synergistically regulate transport of sphingosine-1-phosphate.

Solute carrier spinster homolog 2 (SPNS2), one of only four known major facilitator superfamily (MFS) lysolipid transporters in humans, exports sphingosine-1-phosphate (S1P) across cell membranes. Here, we explore the synergistic effects of lipid binding and conformational dynamics on SPNS2's transport mechanism. Using mass spectrometry, we discovered that SPNS2 interacts preferentially with PI(4,5)P2. Together with functional studies and molecular dynamics (MD) simulations, we identified potential PI(4,5)P2 binding sites. Mutagenesis of proposed lipid binding sites and inhibition of PI(4,5)P2 synthesis reduce S1P transport, whereas the absence of the N terminus renders the transporter essentially inactive. Probing the conformational dynamics of SPNS2, we show how synergistic binding of PI(4,5)P2 and S1P facilitates transport, increases dynamics of the extracellular gate, and stabilizes the intracellular gate. Given that SPNS2 transports a key signaling lipid, our results have implications for therapeutic targeting and also illustrate a regulatory mechanism for MFS transporters.

Regulation of GTPase function by autophosphorylation.

A unifying feature of the RAS superfamily is a conserved GTPase cycle by which these proteins transition between active and inactive states. We demonstrate that autophosphorylation of some GTPases is an intrinsic regulatory mechanism that reduces nucleotide hydrolysis and enhances nucleotide exchange, altering the on/off switch that forms the basis for their signaling functions. Using X-ray crystallography, nuclear magnetic resonance spectroscopy, binding assays, and molecular dynamics on autophosphorylated mutants of H-RAS and K-RAS, we show that phosphoryl transfer from GTP requires dynamic movement of the switch II region and that autophosphorylation promotes nucleotide exchange by opening the active site and extracting the stabilizing Mg2+. Finally, we demonstrate that autophosphorylated K-RAS exhibits altered effector interactions, including a reduced affinity for RAF proteins in mammalian cells. Thus, autophosphorylation leads to altered active site dynamics and effector interaction properties, creating a pool of GTPases that are functionally distinct from their non-phosphorylated counterparts.

Dynamics and timekeeping in biological systems.

This article introduces three reviews on the theme of circadian rhythms.

Molecular insights into the biased signaling mechanism of the µ-opioid receptor.

GPCR functional selectivity opens new opportunities for the design of safer drugs. Ligands orchestrate GPCR signaling cascades by modulating the receptor conformational landscape. Our study provides insights into the dynamic mechanism enabling opioid ligands to preferentially activate the G protein over the ß-arrestin pathways through the µ-opioid receptor (µOR). We combine functional assays in living cells, solution NMR spectroscopy, and enhanced-sampling molecular dynamic simulations to identify the specific µOR conformations induced by G protein-biased agonists. In particular, we describe the dynamic and allosteric communications between the ligand-binding pocket and the receptor intracellular domains, through conserved motifs in class A GPCRs. Most strikingly, the biased agonists trigger µOR conformational changes in the intracellular loop 1 and helix 8 domains, which may impair ß-arrestin binding or signaling. The findings may apply to other GPCR families and provide key molecular information that could facilitate the design of biased ligands.

Gating movements and ion permeation in HCN4 pacemaker channels.

The HCN1-4 channel family is responsible for the hyperpolarization-activated cation current If/Ih that controls automaticity in cardiac and neuronal pacemaker cells. We present cryoelectron microscopy (cryo-EM) structures of HCN4 in the presence or absence of bound cAMP, displaying the pore domain in closed and open conformations. Analysis of cAMP-bound and -unbound structures sheds light on how ligand-induced transitions in the channel cytosolic portion mediate the effect of cAMP on channel gating and highlights the regulatory role of a Mg2+ coordination site formed between the C-linker and the S4-S5 linker. Comparison of open/closed pore states shows that the cytosolic gate opens through concerted movements of the S5 and S6 transmembrane helices. Furthermore, in combination with molecular dynamics analyses, the open pore structures provide insights into the mechanisms of K+/Na+ permeation. Our results contribute mechanistic understanding on HCN channel gating, cyclic nucleotide-dependent modulation, and ion permeation.

Structure of a GRK5-Calmodulin Complex Reveals Molecular Mechanism of GRK Activation and Substrate Targeting.

The phosphorylation of G protein-coupled receptors (GPCRs) by GPCR kinases (GRKs) facilitates arrestin binding and receptor desensitization. Although this process can be regulated by Ca2+-binding proteins such as calmodulin (CaM) and recoverin, the molecular mechanisms are poorly understood. Here, we report structural, computational, and biochemical analysis of a CaM complex with GRK5, revealing how CaM shapes GRK5 response to calcium. The CaM N and C domains bind independently to two helical regions at the GRK5 N and C termini to inhibit GPCR phosphorylation, though only the C domain interaction disrupts GRK5 membrane association, thereby facilitating cytoplasmic translocation. The CaM N domain strongly activates GRK5 via ordering of the amphipathic αN-helix of GRK5 and allosteric disruption of kinase-RH domain interaction for phosphorylation of cytoplasmic GRK5 substrates. These results provide a framework for understanding how two functional effects, GRK5 activation and localization, can cooperate under control of CaM for selective substrate targeting by GRK5.

Computationally reconstructing cotranscriptional RNA folding from experimental data reveals rearrangement of non-native folding intermediates.

The series of RNA folding events that occur during transcription can critically influence cellular RNA function. Here, we present reconstructing RNA dynamics from data (R2D2), a method to uncover details of cotranscriptional RNA folding. We model the folding of the Escherichia coli signal recognition particle (SRP) RNA and show that it requires specific local structural fluctuations within a key hairpin to engender efficient cotranscriptional conformational rearrangement into the functional structure. All-atom molecular dynamics simulations suggest that this rearrangement proceeds through an internal toehold-mediated strand-displacement mechanism, which can be disrupted with a point mutation that limits local structural fluctuations and rescued with compensating mutations that restore these fluctuations. Moreover, a cotranscriptional folding intermediate could be cleaved in vitro by recombinant E. coli RNase P, suggesting potential cotranscriptional processing. These results from experiment-guided multi-scale modeling demonstrate that even an RNA with a simple functional structure can undergo complex folding and processing during synthesis.