Structural mobility tunes signalling of the GluA1 AMPA glutamate receptor.

AMPA glutamate receptors (AMPARs), the primary mediators of excitatory neurotransmission in the brain, are either GluA2 subunit-containing and thus Ca2+-impermeable, or GluA2-lacking and Ca2+-permeable1. Despite their prominent expression throughout interneurons and glia, their role in long-term potentiation and their involvement in a range of neuropathologies2, structural information for GluA2-lacking receptors is currently absent. Here we determine and characterize cryo-electron microscopy structures of the GluA1 homotetramer, fully occupied with TARPγ3 auxiliary subunits (GluA1/γ3). The gating core of both resting and open-state GluA1/γ3 closely resembles GluA2-containing receptors. However, the sequence-diverse N-terminal domains (NTDs) give rise to a highly mobile assembly, enabling domain swapping and subunit re-alignments in the ligand-binding domain tier that are pronounced in desensitized states. These transitions underlie the unique kinetic properties of GluA1. A GluA2 mutant (F231A) increasing NTD dynamics phenocopies this behaviour, and exhibits reduced synaptic responses, reflecting the anchoring function of the AMPAR NTD at the synapse. Together, this work underscores how the subunit-diverse NTDs determine subunit arrangement, gating properties and ultimately synaptic signalling efficiency among AMPAR subtypes.

Scipion-EM-ProDy: A Graphical Interface for the ProDy Python Package within the Scipion Workflow Engine Enabling Integration of Databases, Simulations and Cryo-Electron Microscopy Image Processing.

Macromolecular assemblies, such as protein complexes, undergo continuous structural dynamics, including global reconfigurations critical for their function. Two fast analytical methods are widely used to study these global dynamics, namely elastic network model normal mode analysis and principal component analysis of ensembles of structures. These approaches have found wide use in various computational studies, driving the development of complex pipelines in several software packages. One common theme has been conformational sampling through hybrid simulations incorporating all-atom molecular dynamics and global modes of motion. However, wide functionality is only available for experienced programmers with limited capabilities for other users. We have, therefore, integrated one popular and extensively developed software for such analyses, the ProDy Python application programming interface, into the Scipion workflow engine. This enables a wider range of users to access a complete range of macromolecular dynamics pipelines beyond the core functionalities available in its command-line applications and the normal mode wizard in VMD. The new protocols and pipelines can be further expanded and integrated into larger workflows, together with other software packages for cryo-electron microscopy image analysis and molecular simulations. We present the resulting plugin, Scipion-EM-ProDy, in detail, highlighting the rich functionality made available by its development.

GPCR kinases generate an APH1A phosphorylation barcode to regulate amyloid-ß generation.

Emerging evidence suggests that G protein-coupled receptor (GPCR) kinases (GRKs) are associated with the pathophysiology of Alzheimer's disease (AD). However, GRKs have not been directly implicated in regulation of the amyloid-ß (Aß) pathogenic cascade in AD. Here, we determine that GRKs phosphorylate a non-canonical substrate, anterior pharynx-defective 1A (APH1A), an integral component of the γ-secretase complex. Significantly, we show that GRKs generate distinct phosphorylation barcodes in intracellular loop 2 (ICL2) and the C terminus of APH1A, which differentially regulate recruitment of the scaffolding protein ß-arrestin 2 (ßarr2) to APH1A and γ-secretase-mediated Aß generation. Further molecular dynamics simulation studies reveal an interaction between the ßarr2 finger loop domain and ICL2 and ICL3 of APH1A, similar to a GPCR-ß-arrestin complex, which regulates γ-secretase activity. Collectively, these studies provide insight into the molecular and structural determinants of the APH1A-ßarr2 interaction that critically regulate Aß generation.

Mechanisms underlying TARP modulation of the GluA1/2-γ8 AMPA receptor.

AMPA-type glutamate receptors (AMPARs) mediate rapid signal transmission at excitatory synapses in the brain. Glutamate binding to the receptor's ligand-binding domains (LBDs) leads to ion channel activation and desensitization. Gating kinetics shape synaptic transmission and are strongly modulated by transmembrane AMPAR regulatory proteins (TARPs) through currently incompletely resolved mechanisms. Here, electron cryo-microscopy structures of the GluA1/2 TARP-γ8 complex, in both open and desensitized states (at 3.5 Å), reveal state-selective engagement of the LBDs by the large TARP-γ8 loop ('ß1'), elucidating how this TARP stabilizes specific gating states. We further show how TARPs alter channel rectification, by interacting with the pore helix of the selectivity filter. Lastly, we reveal that the Q/R-editing site couples the channel constriction at the filter entrance to the gate, and forms the major cation binding site in the conduction path. Our results provide a mechanistic framework of how TARPs modulate AMPAR gating and conductance.

Impact of new variants on SARS-CoV-2 infectivity and neutralization: A molecular assessment of the alterations in the spike-host protein interactions.

The emergence of SARS-CoV-2 variants necessitates rational assessment of their impact on the recognition and neutralization of the virus by the host cell. We present a comparative analysis of the interactions of Alpha, Beta, Gamma, and Delta variants with cognate molecules (ACE2 and/or furin), neutralizing nanobodies (Nbs), and monoclonal antibodies (mAbs) using in silico methods, in addition to Nb-binding assays. Our study elucidates the molecular origin of the ability of Beta and Delta variants to evade selected antibodies, such as REGN10933, LY-CoV555, B38, C105, or H11-H4, while being insensitive to others including REGN10987. Experiments confirm that nanobody Nb20 retains neutralizing activity against the Delta variant. The substitutions T478K and L452R in the Delta variant enhance associations with ACE2, whereas P681R promotes recognition by proteases, thus facilitating viral entry. The Ab-specific responses of variants highlight how full-atomic structure and dynamics analyses are required for assessing the response to newly emerging variants.

Sampling of Protein Conformational Space Using Hybrid Simulations: A Critical Assessment of Recent Methods.

Recent years have seen several hybrid simulation methods for exploring the conformational space of proteins and their complexes or assemblies. These methods often combine fast analytical approaches with computationally expensive full atomic molecular dynamics (MD) simulations with the goal of rapidly sampling large and cooperative conformational changes at full atomic resolution. We present here a systematic comparison of the utility and limits of four such hybrid methods that have been introduced in recent years: MD with excited normal modes (MDeNM), collective modes-driven MD (CoMD), and elastic network model (ENM)-based generation, clustering, and relaxation of conformations (ClustENM) as well as its updated version integrated with MD simulations (ClustENMD). We analyzed the predicted conformational spaces using each of these four hybrid methods, applied to four well-studied proteins, triosephosphate isomerase (TIM), 3-phosphoglycerate kinase (PGK), HIV-1 protease (PR) and HIV-1 reverse transcriptase (RT), which provide extensive ensembles of experimental structures for benchmarking and comparing the methods. We show that a rigorous multi-faceted comparison and multiple metrics are necessary to properly assess the differences between conformational ensembles and provide an optimal protocol for achieving good agreement with experimental data. While all four hybrid methods perform well in general, being especially useful as computationally efficient methods that retain atomic resolution, the systematic analysis of the same systems by these four hybrid methods highlights the strengths and limitations of the methods and provides guidance for parameters and protocols to be adopted in future studies.

A monoclonal antibody against staphylococcal enterotoxin B superantigen inhibits SARS-CoV-2 entry in vitro.

We recently discovered a superantigen-like motif sequentially and structurally similar to a staphylococcal enterotoxin B (SEB) segment, near the S1/S2 cleavage site of the SARS-CoV-2 spike protein, which might explain the multisystem inflammatory syndrome (MIS-C) observed in children and the cytokine storm in severe COVID-19 patients. We show here that an anti-SEB monoclonal antibody (mAb), 6D3, can bind this viral motif at its polybasic (PRRA) insert to inhibit infection in live virus assays. The overlap between the superantigenic site of the spike and its proteolytic cleavage site suggests that the mAb prevents viral entry by interfering with the proteolytic activity of cell proteases (furin and TMPRSS2). The high affinity of 6D3 for this site originates from a polyacidic segment at its heavy chain CDR2. The study points to the potential utility of 6D3 for possibly treating COVID-19, MIS-C, or common colds caused by human coronaviruses that also possess a furin-like cleavage site.

ProDy 2.0: increased scale and scope after 10 years of protein dynamics modelling with Python.

SUMMARY: ProDy, an integrated application programming interface developed for modelling and analysing protein dynamics, has significantly evolved in recent years in response to the growing data and needs of the computational biology community. We present major developments that led to ProDy 2.0: (i) improved interfacing with databases and parsing new file formats, (ii) SignDy for signature dynamics of protein families, (iii) CryoDy for collective dynamics of supramolecular systems using cryo-EM density maps and (iv) essential site scanning analysis for identifying sites essential to modulating global dynamics. AVAILABILITY AND IMPLEMENTATION: ProDy is open-source and freely available under MIT License from https://github.com/prody/ProDy. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A monoclonal antibody against staphylococcal enterotoxin B superantigen inhibits SARS-CoV-2 entry in vitro.

We recently discovered a superantigen-like motif, similar to Staphylococcal enterotoxin B (SEB), near the S1/S2 cleavage site of SARS-CoV-2 Spike protein, which might explain the multisystem-inflammatory syndrome (MIS-C) observed in children and cytokine storm in severe COVID-19 patients. We show here that an anti-SEB monoclonal antibody (mAb), 6D3, can bind this viral motif, and in particular its PRRA insert, to inhibit infection by blocking the access of host cell proteases, TMPRSS2 or furin, to the cleavage site. The high affinity of 6D3 for the furin-cleavage site originates from a poly-acidic segment at its heavy chain CDR2, a feature shared with SARS-CoV-2-neutralizing mAb 4A8. The affinity of 6D3 and 4A8 for this site points to their potential utility as therapeutics for treating COVID-19, MIS-C, or common cold caused by human coronaviruses (HCoVs) that possess a furin-like cleavage site.

Towards gaining sight of multiscale events: utilizing network models and normal modes in hybrid methods.

With the explosion of normal mode analyses (NMAs) based on elastic network models (ENMs) in the last decade, and the proven precision of MD simulations for visualizing interactions at atomic scale, many hybrid methods have been proposed in recent years. These aim at exploiting the best of both worlds: the atomic precision of MD that often fall short of exploring time and length scales of biological interest, and the capability of ENM-NMA to predict the cooperative and often functional rearrangements of large structures and assemblies, albeit at low resolution. We present an overview of recent progress in the field with examples of successful applications highlighting the utility of such hybrid methods and pointing to emerging future directions guided by advances in experimental characterization of biomolecular systems structure and dynamics.

Pharmmaker: Pharmacophore modeling and hit identification based on druggability simulations.

Recent years have seen progress in druggability simulations, that is, molecular dynamics simulations of target proteins in solutions containing drug-like probe molecules to characterize their drug-binding abilities, if any. An important consecutive step is to analyze the trajectories to construct pharmacophore models (PMs) to use for virtual screening of libraries of small molecules. While considerable success has been observed in this type of computer-aided drug discovery, a systematic tool encompassing multiple steps from druggability simulations to pharmacophore modeling, to identifying hits by virtual screening of libraries of compounds, has been lacking. We address this need here by developing a new tool, Pharmmaker, building on the DruGUI module of our ProDy application programming interface. Pharmmaker is composed of a suite of steps: (Step 1) identification of high affinity residues for each probe molecule type; (Step 2) selecting high affinity residues and hot spots in the vicinity of sites identified by DruGUI; (Step 3) ranking of the interactions between high affinity residues and specific probes; (Step 4) obtaining probe binding poses and corresponding protein conformations by collecting top-ranked snapshots; and (Step 5) using those snapshots for constructing PMs. The PMs are then used as filters for identifying hits in structure-based virtual screening. Pharmmaker, accessible online at http://prody.csb.pitt.edu/pharmmaker, can be used in conjunction with other tools available in ProDy.

Shared Signature Dynamics Tempered by Local Fluctuations Enables Fold Adaptability and Specificity.

Recent studies have drawn attention to the evolution of protein dynamics, in addition to sequence and structure, based on the premise structure-encodes-dynamics-encodes-function. Of interest is to understand how functional differentiation is accomplished while maintaining the fold, or how intrinsic dynamics plays out in the evolution of structural variations and functional specificity. We performed a systematic computational analysis of 26,899 proteins belonging to 116 CATH superfamilies. Characterizing cooperative mechanisms and convergent/divergent features that underlie the shared/differentiated dynamics of family members required a methodology that lends itself to efficient analyses of large ensembles of proteins. We therefore introduced, SignDy, an integrated pipeline for evaluating the signature dynamics of families based on elastic network models. Our analysis confirmed that family members share conserved, highly cooperative (global) modes of motion. Importantly, our analysis discloses a subset of motions that sharply distinguishes subfamilies, which lie in a low-to-intermediate frequency regime of the mode spectrum. This regime has maximal impact on functional differentiation of families into subfamilies, while being evolutionarily conserved among subfamily members. Notably, the high-frequency end of the spectrum also reveals evolutionary conserved features across and within subfamilies; but in sharp contrast to global motions, high-frequency modes are minimally collective. Modulation of robust/conserved global dynamics by low-to-intermediate frequency fluctuations thus emerges as a versatile mechanism ensuring the adaptability of selected folds and the specificity of their subfamilies. SignDy further allows for dynamics-based categorization as a new layer of information relevant to distinctive mechanisms of action of subfamilies, beyond sequence or structural classifications.

Conformational recognition of an intrinsically disordered protein.

There is a growing interest in understanding the properties of intrinsically disordered proteins (IDPs); however, the characterization of these states remains an open challenge. IDPs appear to have functional roles that diverge from those of folded proteins and revolve around their ability to act as hubs for protein-protein interactions. To gain a better understanding of the modes of binding of IDPs, we combined statistical mechanics, calorimetry, and NMR spectroscopy to investigate the recognition and binding of a fragment from the disordered protein Gab2 by the growth factor receptor-bound protein 2 (Grb2), a key interaction for normal cell signaling and cancer development. Structural ensemble refinement by NMR chemical shifts, thermodynamics measurements, and analysis of point mutations indicated that the population of preexisting bound conformations in the free-state ensemble of Gab2 is an essential determinant for recognition and binding by Grb2. A key role was found for transient polyproline II (PPII) structures and extended conformations. Our findings are likely to have very general implications for the biological behavior of IDPs in light of the evidence that a large fraction of these proteins possess a specific propensity to form PPII and to adopt conformations that are more extended than the typical random-coil states.