Membrane protein assembly patterns reflect selection for non-proliferative structures

Membrane proteins that regulate solute movement are often built from multiple copies of an identical polypeptide chain. These complexes represent striking examples of self-assembling systems that recruit monomers only until a prescribed level for function is reached. Here we report that three modes of assembly - distinguished by sequence and stoichiometry - describe all helical membrane protein complexes currently solved to high resolution. Using the 13 presently available non-redundant homo-oligomeric structures, we show that two of these types segregate with protein function: one produces energy-dependent transporters, while the other builds channels for passive diffusion. Given such limited routes to functional complexes, membrane proteins that self-assemble exist on the edge of aggregation, susceptible to mutations that may underlie human diseases.     

Biophysical approaches to membrane protein structure determination.

Recently, there have been several technical advances in the use of solution and solid-state NMR spectroscopy to determine the structures of membrane proteins. The structures of several isolated transmembrane (TM) helices and pairs of TM helices have been solved by solution NMR methods. Similarly, the complete folds of two TM beta-barrel proteins with molecular weights of 16 and 19 kDa have been determined by solution NMR in detergent micelles. Solution NMR has also provided a first glimpse at the dynamics of an integral membrane protein. Structures of individual TM helices have also been determined by solid-state NMR. A combination of NMR with site-directed spin-label electron paramagnetic resonance or Fourier transform IR spectroscopy allows one to assemble quite detailed protein structures in the membrane.     

Deriving correlated motions in proteins from X-ray structure refinement by using TLS parameters

Dynamic information in proteins may provide valuable information for understanding allosteric regulation of protein complexes or long-range effects of the mutations on enzyme activity. Experimental data such as X-ray B-factors or NMR order parameters provide a convenient estimate of atomic fluctuations (or atomic auto-correlated motions) in proteins. However, it is not as straightforward to obtain atomic cross-correlated motions in proteins — one usually resorts to more sophisticated computational methods such as Molecular Dynamics, normal mode analysis or atomic network models. In this report, we show that atomic cross-correlations can be reliably obtained directly from protein structure using X-ray refinement data. We have derived an analytic form of atomic correlated motions in terms of the original TLS parameters used to refine the B-factors of X-ray structures. The correlated maps computed using this equation are well correlated with those of the method based on a mechanical model (the correlation coefficient is 0.75) for a non-homologous dataset comprising 100 structures. We have developed an approach to compute atomic cross-correlations directly from X-ray protein structure. Being in analytic form, it is fast and provides a feasible way to compute correlated motions in proteins in a high throughput way. In addition, avoiding sophisticated computational operations; it provides a quick, reliable way, especially for non-computational biologists, to obtain dynamics information directly from protein structure relevant to its function.     

Solid-state NMR structures of integral membrane proteins

Abstract

Solid-state NMR is unique for its ability to obtain three-dimensional structures and to measure atomic-resolution structural and dynamic information for membrane proteins in native lipid bilayers. An increasing number and complexity of integral membrane protein structures have been determined by solid-state NMR using two main methods. Oriented sample solid-state NMR uses macroscopically aligned lipid bilayers to obtain orientational restraints that define secondary structure and global fold of embedded peptides and proteins and their orientation and topology in lipid bilayers. Magic angle spinning (MAS) solid-state NMR uses unoriented rapidly spinning samples to obtain distance and torsion angle restraints that define tertiary structure and helix packing arrangements. Details of all current protein structures are described, highlighting developments in experimental strategy and other technological advancements. Some structures originate from combining solid- and solution-state NMR information and some have used solid-state NMR to refine X-ray crystal structures. Solid-state NMR has also validated the structures of proteins determined in different membrane mimetics by solution-state NMR and X-ray crystallography and is therefore complementary to other structural biology techniques. By continuing efforts in identifying membrane protein targets and developing expression, isotope labelling and sample preparation strategies, probe technology, NMR experiments, calculation and modelling methods and combination with other techniques, it should be feasible to determine the structures of many more membrane proteins of biological and biomedical importance using solid-state NMR. This will provide three-dimensional structures and atomic-resolution structural information for characterising ligand and drug interactions, dynamics and molecular mechanisms of membrane proteins under physiological lipid bilayer conditions.     

Structural Interactions between Lipids, Water and S1-S4 Voltage-Sensing Domains

Membrane proteins serve crucial signaling and transport functions, yet relatively little is known about their structures in membrane environments or how lipids interact with these proteins. For voltage-activated ion channels, X-ray structures suggest that the mobile voltage-sensing S4 helix would be exposed to the membrane, and functional studies reveal that lipid modification can profoundly alter channel activity. Here, we use solid-state NMR to investigate structural interactions of lipids and water with S1-S4 voltage-sensing domains and to explore whether lipids influence the structure of the protein. Our results demonstrate that S1-S4 domains exhibit extensive interactions with lipids and that these domains are heavily hydrated when embedded in a membrane. We also find evidence for preferential interactions of anionic lipids with S1-S4 domains and that these interactions have lifetimes on the timescale of ≤ 10^− 3 s. Arg residues within S1-S4 domains are well hydrated and are positioned in close proximity to lipids, exhibiting local interactions with both lipid headgroups and acyl chains. Comparative studies with a positively charged lipid lacking a phosphodiester group reveal that this lipid modification has only modest effects on the structure and hydration of S1-S4 domains. Taken together, our results demonstrate that Arg residues in S1-S4 voltage-sensing domains reside in close proximity to the hydrophobic interior of the membrane yet are well hydrated, a requirement for carrying charge and driving protein motions in response to changes in membrane voltage.     

Structure and dynamics of the influenza A M2 Channel: a comparison of three structures

The M2 protein is an essential component of the Influenza virus’ infectivity cycle. It is a homo-tetrameric bundle forming a pH-gated H⁺ channel. The structure of M2 was solved by three different groups, using different techniques, protein sequences and pH environment. For example, solid-state NMR spectroscopy was used on a protein in lipid bilayers, while X-ray crystallography and solution NMR spectroscopy were applied on a protein in detergent micelles. The resulting structures from the above efforts are rather distinct. Herein, we examine the different structures under uniform conditions such as a lipid bilayer and specified protonation state. We employ extensive molecular dynamics simulations, in several protonation states, representing both closed and open forms of the channel. Exploring the properties of each of these structures has shown that the X-ray structure is more stable than the other structures according to various criteria, although its water conductance and water-wire formation do not correlate to the protonation state of the channel.     

Statistical analysis of interface similarity in crystals of homologous proteins

Many proteins function as homo-oligomers and are regulated via their oligomeric state. For some proteins, the stoichiometry of homo-oligomeric states under various conditions has been studied using gel filtration or analytical ultracentrifugation experiments. The interfaces involved in these assemblies may be identified using cross-linking and mass spectrometry, solution-state NMR, and other experiments. However, for most proteins, the actual interfaces that are involved in oligomerization are inferred from X-ray crystallographic structures using assumptions about interface surface areas and physical properties. Examination of interfaces across different Protein Data Bank (PDB) entries in a protein family reveals several important features. First, similarities in space group, asymmetric unit size, and cell dimensions and angles (within 1%) do not guarantee that two crystals are actually the same crystal form, containing similar relative orientations and interactions within the crystal. Conversely, two crystals in different space groups may be quite similar in terms of all the interfaces within each crystal. Second, NMR structures and an existing benchmark of PDB crystallographic entries consisting of 126 dimers as well as larger structures and 132 monomers were used to determine whether the existence or lack of common interfaces across multiple crystal forms can be used to predict whether a protein is an oligomer or not. Monomeric proteins tend to have common interfaces across only a minority of crystal forms, whereas higher-order structures exhibit common interfaces across a majority of available crystal forms. The data can be used to estimate the probability that an interface is biological if two or more crystal forms are available. Finally, the Protein Interfaces, Surfaces, and Assemblies (PISA) database available from the European Bioinformatics Institute is more consistent in identifying interfaces observed in many crystal forms compared with the PDB and the European Bioinformatics Institute's Protein Quaternary Server (PQS). The PDB, in particular, is missing highly likely biological interfaces in its biological unit files for about 10% of PDB entries.     

NMR studies on fully hydrated membrane proteins, with emphasis on bacteriorhodopsin as a typical and prototype membrane protein

The 3D structures or dynamic feature of fully hydrated membrane proteins are very important at ambient temperature, in relation to understanding their biological activities, although their data, especially from the flexible portions such as surface regions, are unavailable from X-ray diffraction or cryoelectron microscope at low temperature. In contrast, high-resolution solid-state NMR spectroscopy has proved to be a very convenient alternative means to be able to reveal their dynamic structures. To clarify this problem, we describe here how we are able to reveal such structures and dynamic features, based on intrinsic probes from high-resolution solid-state NMR studies on bacteriorhodopsin (bR) as a typical membrane protein in 2D crystal, regenerated preparation in lipid bilayer and detergents. It turned out that their dynamic features are substantially altered upon their environments where bR is present. We further review NMR applications to study structure and dynamics of a variety of membrane proteins, including sensory rhodopsin, rhodopsin, photoreaction centers, diacylglycerol kinases, etc.     

Solid-State NMR-Restrained Ensemble Dynamics of a Membrane Protein in Explicit Membranes

Solid-state NMR has been used to determine the structures of membrane proteins in native-like lipid bilayer environments. Most structure calculations based on solid-state NMR observables are performed using simulated annealing with restrained molecular dynamics and an energy function, where all nonbonded interactions are represented by a single, purely repulsive term with no contributions from van der Waals attractive, electrostatic, or solvation energy. To our knowledge, this is the first application of an ensemble dynamics technique performed in explicit membranes that uses experimental solid-state NMR observables to obtain the refined structure of a membrane protein together with information about its dynamics and its interactions with lipids. Using the membrane-bound form of the fd coat protein as a model membrane protein and its experimental solid-state NMR data, we performed restrained ensemble dynamics simulations with different ensemble sizes in explicit membranes. For comparison, a molecular dynamics simulation of fd coat protein was also performed without any restraints. The average orientation of each protein helix is similar to a structure determined by traditional single-conformer approaches. However, their variations are limited in the resulting ensemble of structures with one or two replicas, as they are under the strong influence of solid-state NMR restraints. Although highly consistent with all solid-state NMR observables, the ensembles of more than two replicas show larger orientational variations similar to those observed in the molecular dynamics simulation without restraints. In particular, in these explicit membrane simulations, Lys⁴⁰, residing at the C-terminal side of the transmembrane helix, is observed to cause local membrane curvature. Therefore, compared to traditional single-conformer approaches in implicit environments, solid-state NMR restrained ensemble simulations in explicit membranes readily characterize not only protein dynamics but also protein-lipid interactions in detail.     

Cover Image,Volume 88, Issue 9

There have been several studies suggesting that protein structures solved by NMR spectroscopy and X-ray crystallography show significant differences. To understand the origin of these differences, we assembled a database of high-quality protein structures solved by both methods. We also find significant differences between NMR and crystal structures—in the root-mean-square deviations of the C _α atomic positions, identities of core amino acids, backbone, and side-chain dihedral angles, and packing fraction of core residues. In contrast to prior studies, we identify the physical basis for these differences by modeling protein cores as jammed packings of amino acid-shaped particles. We find that we can tune the jammed packing fraction by varying the degree of thermalization used to generate the packings. For an athermal protocol, we find that the average jammed packing fraction is identical to that observed in the cores of protein structures solved by X-ray crystallography. In contrast, highly thermalized packing-generation protocols yield jammed packing fractions that are even higher than those observed in NMR structures. These results indicate that thermalized systems can pack more densely than athermal systems, which suggests a physical basis for the structural differences between protein structures solved by NMR and X-ray crystallography.     

Improved technologies now routinely provide protein NMR structures useful for molecular replacement

Molecular replacement (MR) is widely used for addressing the phase problem in X-ray crystallography. Historically, crystallographers have had limited success using NMR structures as MR search models. Here, we report a comprehensive investigation of the utility of protein NMR ensembles as MR search models, using data for 25 pairs of X-ray and NMR structures solved and refined using modern NMR methods. Starting from NMR ensembles prepared by an improved protocol, FindCore, correct MR solutions were obtained for 22 targets. Based on these solutions, automatic model rebuilding could be done successfully. Rosetta refinement of NMR structures provided MR solutions for another two proteins. We also demonstrate that such properly prepared NMR ensembles and X-ray crystal structures have similar performance when used as MR search models for homologous structures, particularly for targets with sequence identity >40%.     

Evaluating protein structures determined by structural genomics consortia

Structural genomics projects are providing large quantities of new 3D structural data for proteins. To monitor the quality of these data, we have developed the protein structure validation software suite (PSVS), for assessment of protein structures generated by NMR or X-ray crystallographic methods. PSVS is broadly applicable for structure quality assessment in structural biology projects. The software integrates under a single interface analyses from several widely-used structure quality evaluation tools, including PROCHECK (Laskowski et al., J Appl Crystallog 1993;26:283-291), MolProbity (Lovell et al., Proteins 2003;50:437-450), Verify3D (Luthy et al., Nature 1992;356:83-85), ProsaII (Sippl, Proteins 1993;17: 355-362), the PDB validation software, and various structure-validation tools developed in our own laboratory. PSVS provides standard constraint analyses, statistics on goodness-of-fit between structures and experimental data, and knowledge-based structure quality scores in standardized format suitable for database integration. The analysis provides both global and site-specific measures of protein structure quality. Global quality measures are reported as Z scores, based on calibration with a set of high-resolution X-ray crystal structures. PSVS is particularly useful in assessing protein structures determined by NMR methods, but is also valuable for assessing X-ray crystal structures or homology models. Using these tools, we assessed protein structures generated by the Northeast Structural Genomics Consortium and other international structural genomics projects, over a 5-year period. Protein structures produced from structural genomics projects exhibit quality score distributions similar to those of structures produced in traditional structural biology projects during the same time period. However, while some NMR structures have structure quality scores similar to those seen in higher-resolution X-ray crystal structures, the majority of NMR structures have lower scores. Potential reasons for this "structure quality score gap" between NMR and X-ray crystal structures are discussed.     

Comparison of multiple crystal structures with NMR data for engrailed homeodomain

Two methods are currently available to solve high resolution protein structures—X-ray crystallography and nuclear magnetic resonance (NMR). Both methods usually produce highly similar structures, but small differences between both solutions are always observed. Here the raw NMR data as well as the solved NMR structure were compared to the multiple crystal structures solved for the WT 60 residue three helix bundle engrailed homeodomain (EnHD) and single point mutants. There was excellent agreement between TALOS-predicted and crystal structure-observed dihedral angles and a good agreement for the ³ J(H ^N H ^α ) couplings for the multiple crystal structures. Around 1% of NOEs were violated for any crystal structure, but no NOE was inconsistent with all of the crystal structures. Violations usually occurred for surface residues or for residues for which multiple discreet conformations were observed between the crystal structures. Comparison of the disorder shown in the multiple crystal structures shows little correlation with dynamics under native conditions for this protein.     

Mapping of protein-protein interaction sites by the 'absence of interference' approach

Protein-protein interactions are critical to most biological processes, and locating protein-protein interfaces on protein structures is an important task in molecular biology. We developed a new experimental strategy called the ‘absence of interference’ approach to determine surface residues involved in protein-protein interaction of established yeast two-hybrid pairs of interacting proteins. One of the proteins is subjected to high-level randomization by error-prone PCR. The resulting library is selected by yeast two-hybrid system for interacting clones that are isolated and sequenced. The interaction region can be identified by an absence or depletion of mutations. For data analysis and presentation, we developed a Web interface that analyzes the mutational spectrum and displays the mutational frequency on the surface of the structure (or a structural model) of the randomized protein†. Additionally, this interface might be of use for the display of mutational distributions determined by other types of random mutagenesis experiments. We applied the approach to map the interface of the catalytic domain of the DNA methyltransferase Dnmt3a with its regulatory factor Dnmt3L. Dnmt3a was randomized with high mutational load. A total of 76 interacting clones were isolated and sequenced, and 648 mutations were identified. The mutational pattern allowed to identify a unique interaction region on the surface of Dnmt3a, which comprises about 500-600 Å². The results were confirmed by site-directed mutagenesis and structural analysis. The absence-of-interference approach will allow high-throughput mapping of protein interaction sites suitable for functional studies and protein docking.     

NMR Resonance Assignment Methodology: Characterizing Large Sparsely Labeled Glycoproteins

Characterization of proteins using NMR methods begins with assignment of resonances to specific residues. This is usually accomplished using sequential connectivities between nuclear pairs in proteins uniformly labeled with NMR active isotopes. This becomes impractical for larger proteins, and especially for proteins that are best expressed in mammalian cells, including glycoproteins. Here an alternate protocol for the assignment of NMR resonances of sparsely labeled proteins, namely, the ones labeled with a single amino acid type, or a limited subset of types, isotopically enriched with ¹⁵N or ¹³C, is described. The protocol is based on comparison of data collected using extensions of simple two-dimensional NMR experiments (correlated chemical shifts, nuclear Overhauser effects, residual dipolar couplings) to predictions from molecular dynamics trajectories that begin with known protein structures. Optimal pairing of predicted and experimental values is facilitated by a software package that employs a genetic algorithm, ASSIGN_SLP_MD. The approach is applied to the 36-kDa luminal domain of the sialyltransferase, rST6Gal1, in which all phenylalanines are labeled with ¹⁵N, and the results are validated by elimination of resonances via single-point mutations of selected phenylalanines to tyrosines. Assignment allows the use of previously published paramagnetic relaxation enhancements to evaluate placement of a substrate analog in the active site of this protein. The protocol will open the way to structural characterization of the many glycosylated and other proteins that are best expressed in mammalian cells.     

Crystallographic and NMR spectroscopic protein structures: Interresidue contacts

Interresidue pair contacts were analyzed in detail for four pairs of protein structures solved using X-ray analysis (X-ray) and nuclear magnetic resonance (NMR). In the four NMR structures, at distances of ≤4.0 Å, the total number of pair contacts was 4–9% lower and, in general, the pair contacts were 0.02–0.16 Å shorter compared to the X-ray structures. Each of the four structural pairs contained 83–94% common pair contacts (CPCs), which were formed by identical residues in both structures; the other 6–17% were longer intrinsic pair contacts (IPCs) formed by different residues in NMR and X-ray structures, while the latter contained more IPC. Every NMR structure contained three types of CPC that were shorter, longer, or equal to the identical contact pairs in the X-ray structure of this protein. Methodologically different short CPCs prevailed at a known distance dependence of the interresidue contact density in 60–61 pairs of NMR/X-ray structures. Among the analyzed four structural pairs, contact shortening appeared upon the energy minimization of the crambin NMR structure and upon solving the ubiquitin, hen lysozyme, and monomeric hemoglobin NMR structures using X-PLOR software with decreased van der Waals atomic radii. The degree of contact shortening in the NMR structures diminished with an increase in the NMR data used to solve these structures. Among the 60 pairs of NMR/X-ray structures, the major difference between α-helical and β-structural proteins in the dependences on interresidue distances of average contact density appeared due to strong α/β differences in the backbone local geometry.     

Membranome 3.0: Database of single‐pass membrane proteins with AlphaFold models

The Membranome database provides comprehensive structural information on single‐pass (i.e., bitopic) membrane proteins from six evolutionarily distant organisms, including protein–protein interactions, complexes, mutations, experimental structures, and models of transmembrane α‐helical dimers. We present a new version of this database, Membranome 3.0, which was significantly updated by revising the set of 5,758 bitopic proteins and incorporating models generated by AlphaFold 2 in the database. The AlphaFold models were parsed into structural domains located at the different membrane sides, modified to exclude low‐confidence unstructured terminal regions and signal sequences, validated through comparison with available experimental structures, and positioned with respect to membrane boundaries. Membranome 3.0 was re‐developed to facilitate visualization and comparative analysis of multiple 3D structures of proteins that belong to a specified family, complex, biological pathway, or membrane type. New tools for advanced search and analysis of proteins, their interactions, complexes, and mutations were included. The database is freely accessible at https://membranome.org.     

Tricks of the trade used to accelerate high-resolution structure determination of membrane proteins

The rate at which X-ray structures of membrane proteins are solved is on a par with that of soluble proteins in the late 1970s. There are still many obstacles facing the membrane protein structural community. Recently, there have been several technical achievements in the field that have started to dramatically accelerate structural studies. Here, we summarize these so-called ‘tricks-of-the-trade’ and include case studies of several mammalian transporters.