A graph-based approach identifies dynamic H-bond communication networks in spike protein S of SARS-CoV-2

We apply graph-based approaches to identify H-bond clusters in protein complexes. Three conformations of spike protein S have distinct H-bond clusters at key sites. Hydrogen-bond clusters could govern structural plasticity of spike protein S. Protein S binds to ACE2 receptor via H-bond clusters extending deep across interface.Corona virus spike protein S is a large homo-trimeric protein anchored in the membrane of the virion particle. Protein S binds to angiotensin-converting-enzyme 2, ACE2, of the host cell, followed by proteolysis of the spike protein, drastic protein conformational change with exposure of the fusion peptide of the virus, and entry of the virion into the host cell. The structural elements that govern conformational plasticity of the spike protein are largely unknown. Here, we present a methodology that relies upon graph and centrality analyses, augmented by bioinformatics, to identify and characterize large H-bond clusters in protein structures. We apply this methodology to protein S ectodomain and find that, in the closed conformation, the three protomers of protein S bring the same contribution to an extensive central network of H-bonds, and contribute symmetrically to a relatively large H-bond cluster at the receptor binding domain, and to a cluster near a protease cleavage site. Markedly different H-bonding at these three clusters in open and pre-fusion conformations suggest dynamic H-bond clusters could facilitate structural plasticity and selection of a protein S protomer for binding to the host receptor, and proteolytic cleavage. From analyses of spike protein sequences we identify patches of histidine and carboxylate groups that could be involved in transient proton binding.PSI COVID19 Emergency Science FundSpanish Ministry of Science, Innovation and Universities RTI2018-098983-B-I00Excellence Initiative of the German Federal and State Governments via the Freie Universitat BerlinGerman Research Foundation (DFG) SFB 107

A set of common movements within GPCR-G-protein complexes from variability analysis of cryo-EM datasets

G-protein coupled receptors (GPCRs) are among the most versatile signal transducers in the cell. Once activated, GPCRs sample a large conformational space and couple to G-proteins to initiate distinct signaling pathways. The dynamical behavior of GPCR-G-protein complexes is difficult characterize structurally, and it might hinder obtaining routine high-resolution density maps in single-particle reconstructions. Here, we used variability analysis on the rhodopsin-Gi-Fab16 complex cryo-EM dataset, and the results provide insights into the dynamic nature of the receptor-complex interaction. We compare the outcome of this analysis with recent results obtained on the cannabinoid-Gi- and secretin-Gs-receptor complexes. Despite differences related to the biochemical compositions of the three samples, a set of consensus movements emerges. We anticipate that systematic variability analysis on GPCR-G-protein complexes may provide useful information not only at the biological level, but also for improving the preparation of more stable samples for cryo-EM single-particle analysis.ISSN:1047-8477ISSN:1095-865

The role of water molecules in phototransduction of retinal proteins and G protein-coupled receptors

ISSN:1359-6640ISSN:1364-549

Membrane protein structural biology using X-ray free electron lasers

Membrane protein structural biology has benefitted tremendously from access to micro-focus crystallography at synchrotron radiation sources. X-ray free electron lasers (XFELs) are linear accelerator driven X-ray sources that deliver a jump in peak X-ray brilliance of nine orders of magnitude and represent a disruptive technology with potential to dramatically change the field. Membrane proteins were amongst the first macromolecules to be studied with XFEL radiation and include proof-of-principle demonstrations of serial femtosecond crystallography (SFX), the observation that XFEL data can deliver damage free crystallographic structures, initial experiments towards recording structural information from 2D arrays of membrane proteins, and time-resolved SFX, time- resolved wide angle X-ray scattering and time-resolved X-ray emission spectroscopy studies. Conversely, serial crystallography methods are now being applied using synchrotron radiation. We believe that a context dependent choice of synchrotron or XFEL radiation will accelerate progress towards novel insights in understanding membrane protein structure and dynamics.ISSN:0959-440XISSN:1879-033

Structural basis of the partially open central gate in the human CNGA1/CNGB1 channel explained by additional density for calmodulin in cryo-EM map

The recently reported structure of the human CNGA1/CNGB1 CNG channel in the open state (Xue et al., 2021a) shows that one CNGA1 and one CNGB1 subunit do not open the central hydrophobic gate completely upon cGMP binding. This is different from what has been reported for CNGA homomeric channels (Xue et al., 2021b; Zheng et al., 2020). In seeking to understand how this difference is due to the presence of the CNGB1 subunit, we find that the deposited density map (Xue et al., 2021a) (EMDB 24465) contains an additional density not reported in the images of the original publication. This additional density fits well the structure of calmodulin (CaM), and it unambiguously connects the newly identified D-helix of CNGB1 to one of the CNGA1 helices (A1R) participating in the coiled-coil region. Interestingly, the CNGA1 subunit that engages in the interaction with this additional density is the one that, together with CNGB1, does not open completely the central gate. The sequence of the D-helix of CNGB1 contains a known CaM-binding site of exquisitely high affinity – named CaM2 (Weitz et al., 1998) –, and thus the presence of CaM in that region is not surprising. The mechanism through which CaM reduces currents across the membrane by acting on the native channel (Bauer, 1996; Hsu and Molday, 1993; Weitz et al., 1998) remains unclear. We suggest that the presence of CaM may explain the partially open central gate reported by Xue et al. (2021a). The structure of the open and closed states of the CNGA1/CNGB1 channel may be different with and without CaM present.ISSN:1047-8477ISSN:1095-865

Molecular architecture of a retinal cGMP-gated channel: the arrangement of the cytoplasmic domains

Cyclic nucleotide-gated (CNG) channels play a central role in the conversion of sensory information, such as light and scent, into primary electrical signals. We have purified the CNG channel from bovine retina and have studied it using electron microscopy and image processing. We present the structure of the channel to 35Å resolution. This three-dimensional reconstruction provides insight into the architecture of the protein, suggesting that the cyclic nucleotide-binding domains, which initiate the response to ligand, ‘hang’ below the pore-forming part of the channel, attached by narrow linkers. The structure also suggests that the four cyclic nucleotide-binding domains present in each channel form two distinct domains, lending structural weight to the suggestion that the four subunits of the CNG channels are arranged as a pair of dimers

Dynamics retrieval from stochastically weighted incomplete data by low-pass spectral analysis

ABSTRACT Time-resolved serial femtosecond crystallography (TR-SFX) provides access to protein dynamics on sub-picosecond timescales, and with atomic resolution. Due to the nature of the experiment, these datasets are often highly incomplete and the measured diffracted intensities are affected by partiality. To tackle these issues, one established procedure is that of splitting the data into time bins, and averaging the multiple measurements of equivalent reflections within each bin. This binning and averaging often involve a loss of information. Here, we propose an alternative approach, which we call low-pass spectral analysis (LPSA). In this method, the data are projected onto the subspace defined by a set of trigonometric functions, with frequencies up to a certain cutoff. This approach attenuates undesirable high-frequency features and facilitates retrieving the underlying dynamics. A time-lagged embedding step can be included prior to subspace projection to improve the stability of the results with respect to the parameters involved. Subsequent modal decomposition allows to produce a low-rank description of the system's evolution. Using a synthetic time-evolving model with incomplete and partial observations, we analyze the LPSA results in terms of quality of the retrieved signal, as a function of the parameters involved. We compare the performance of LPSA to that of a range of other sophisticated data analysis techniques. We show that LPSA allows to achieve excellent dynamics reconstruction at modest computational cost. Finally, we demonstrate the superiority of dynamics retrieval by LPSA compared to time binning and merging, which is, to date, the most commonly used method to extract dynamical information from TR-SFX data.ISSN:2329-777

Correction of rhodopsin serial crystallography diffraction intensities for a lattice-translocation defect

Rhodopsin is a G-protein-coupled receptor that detects light and initiates the intracellular signalling cascades that underpin vertebrate vision. Light sensitivity is achieved by covalent linkage to 11-cis retinal, which isomerizes upon photoabsorption. Serial femtosecond crystallography data collected from rhodopsin microcrystals grown in the lipidic cubic phase were used to solve the roomtemperature structure of the receptor. Although the diffraction data showed high completeness and good consistency to 1.8 A degrees resolution, prominent electron-density features remained unaccounted for throughout the unit cell after model building and refinement. A deeper analysis of the diffraction intensities uncovered the presence of a lattice-translocation defect (LTD) within the crystals. The procedure followed to correct the diffraction intensities for this pathology enabled the building of an improved resting-state model. The correction was essential to both confidently model the structure of the unilluminated state and interpret the light-activated data collected after photoexcitation of the crystals. It is expected that similar cases of LTD will be observed in other serial crystallography experiments and that correction will be required in a variety of systems.ISSN:2059-798

The counterion–retinylidene Schiff base interaction of an invertebrate rhodopsin rearranges upon light activation

Animals sense light using photosensitive proteins—rhodopsins—containing a chromophore—retinal—that intrinsically absorbs in the ultraviolet. Visible light-sensitivity depends primarily on protonation of the retinylidene Schiff base (SB), which requires a negatively-charged amino acid residue—counterion—for stabilization. Little is known about how the most common counterion among varied rhodopsins, Glu181, functions. Here, we demonstrate that in a spider visual rhodopsin, orthologue of mammal melanopsins relevant to circadian rhythms, the Glu181 counterion functions likely by forming a hydrogen-bonding network, where Ser186 is a key mediator of the Glu181–SB interaction. We also suggest that upon light activation, the Glu181–SB interaction rearranges while Ser186 changes its contribution. This is in contrast to how the counterion of vertebrate visual rhodopsins, Glu113, functions, which forms a salt bridge with the SB. Our results shed light on the molecular mechanisms of visible light-sensitivity relevant to invertebrate vision and vertebrate non-visual photoreception.ISSN:2399-364

GPCR Activation States Induced by Nanobodies and Mini-G Proteins Compared by NMR Spectroscopy

In this work, we examine methyl nuclear magnetic resonance (NMR) spectra of the methionine ε-[13CH3] labelled thermostabilized β1 adrenergic receptor from turkey in association with a variety of different effectors, including mini-Gs and nanobody 60 (Nb60), which have not been previously studied in complex with β1 adrenergic receptor (β1AR) by NMR. Complexes with pindolol and Nb60 induce highly similar inactive states of the receptor, closely resembling the resting state conformational ensemble. We show that, upon binding of mini-Gs or nanobody 80 (Nb80), large allosteric changes throughout the receptor take place. The conformation of tβ1AR stabilized by the native-like mini-Gs protein is highly similar to the conformation induced by the currently used surrogate Nb80. Interestingly, in both cases residual dynamics are present, which were not observed in the resting states. Finally, we reproduce a pharmaceutically relevant situation, where an antagonist abolishes the interaction of the receptor with the mini-G protein in a competitive manner, validating the functional integrity of our preparation. The presented system is therefore well suited for reproducing the individual steps of the activation cycle of a G protein-coupled receptor (GPCR) in vitro and serves as a basis for functional and pharmacological characterizations of more native-like systems in the future.ISSN:1420-304