NMR assignments of a 48 kDa tetramer of the T1 domain of the mammalian voltage gated potassium channel Kv1.4

The N-terminal cytosolic T1 domain of the mammalian voltage gated potassium channel Kv1.4 is strongly involved in the tetramerization of the Kv1.4 subunit that is required for forming a functional ion channel. The T1 domain forms a stable tetramer of 48 kDa in solution that cannot be dissociated into monomers. In spite of the high molecular mass it was possible to completely assign the backbone and part of the side chain resonances by multidimensional NMR spectroscopy on uniformly ²H, ¹³C, ¹⁵N enriched protein. The secondary structure analysis derived from the chemical shifts is in line with the expectations from X-ray structures of related proteins.     

Self-assembly of a simple membrane protein: coarse-grained molecular dynamics simulations of the influenza M2 channel

The transmembrane (TM) domain of the M2 channel protein from influenza A is a homotetrameric bundle of α-helices and provides a model system for computational approaches to self-assembly of membrane proteins. Coarse-grained molecular dynamics (CG-MD) simulations have been used to explore partitioning into a membrane of M2 TM helices during bilayer self-assembly from lipids. CG-MD is also used to explore tetramerization of preinserted M2 TM helices. The M2 helix monomer adopts a membrane spanning orientation in a lipid (DPPC) bilayer. Multiple extended CG-MD simulations (5 × 5 μs) were used to study the tetramerization of inserted M2 helices. The resultant tetramers were evaluated in terms of the most populated conformations and the dynamics of their interconversion. This analysis reveals that the M2 tetramer has 2× rotationally symmetrical packing of the helices. The helices form a left-handed bundle, with a helix tilt angle of ∼16°. The M2 helix bundle generated by CG-MD was converted to an atomistic model. Simulations of this model reveal that the bundle's stability depends on the assumed protonation state of the H37 side chains. These simulations alongside comparison with recent x-ray (3BKD) and NMR (2RLF) structures of the M2 bundle suggest that the model yielded by CG-MD may correspond to a closed state of the channel.     

Protonation of histidine and histidine-tryptophan interaction in the activation of the M2 ion channel from influenza a virus

The M2 protein of influenza A virus forms a homotetramer ion channel in the lipid membrane. The channel is specific for proton conductance and is activated by low pH with a transition midpoint at pH 5.7. We have studied the structure of the transmembrane domain of the M2 ion channel by using UV resonance Raman spectroscopy, with special attention to the side chains of histidine (His37) and tryptophan (Trp41) residues. The Raman spectra provide direct evidence that the imidazole ring of His37 is protonated upon channel activation at low pH. Concomitantly, the UV resonance Raman scattering from Trp41 shows an unusual intensity change, which is ascribed to a cation-pi interaction between the protonated (cationic) imidazole ring of His37 and the indole ring of Trp41. The protonation of His37 and the Raman intensity change of Trp41 do not occur in the presence of amantadine that blocks the M2 ion channel. These observations clearly show that the protonation of His37 and concomitant cation-pi interaction with Trp41 is a key step in the activation of the M2 ion channel. The His37-Trp41 interaction associated with the channel activation is explained by assuming a conformational transition of His37 induced by electrostatic repulsion among the protonated imidazole rings of four His37 residues in the tetramer channel. Trp41 may play a role in stabilizing the channel open state through cation-pi interaction with His37. A molecular model for the activation of M2 ion channel is proposed on the basis of the gating mechanism.     

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Comparison of solution and crystal structures of maize nonspecific lipid transfer protein: A model for a potential in vivo lipid carrier protein

The three-dimensional solution structure of maize nonspecific lipid transfer protein (nsLTP) obtained by nuclear magnetic resonance (NMR) is compared to the X-ray structure. Although both structures are very similar, some local structural differences are observed in the first and the fourth helices and in several side-chain conformations. These discrepancies arise partly from intermolecular contacts in the crystal lattice. The main characteristic of nsLTP structures is the presence of an internal hydrophobic cavity whose volume was found to vary from 237 to 513 ų without major variations in the 15 solution structures. Comparison of crystal and NMR structures shows the existence of another small hollow at the periphery of the protein containing a water molecule in the X-ray structure, which could play an important structural role. A model of the complexed form of maize nsLTP by α-lysopalmitoylphosphatidylcholine was built by docking the lipid inside the protein cavity of the NMR structure. The main structural feature is a hydrogen bond found also in the X-ray structure of the complex maize nsLTP/palmitate between the hydroxyl of Tyr81 and the carbonyl of the lipid. Comparison of 12 primary sequences of nsLTPs emphasizes that all residues delineating the cavities calculated on solution and X-ray structures are conserved, which suggests that this large cavity is a common feature of all compared plant nsLTPs. Furthermore several conserved basic residues seem to be involved in the stabilization of the protein architecture. Proteins 31:160–171, 1998. © 1998 Wiley-Liss, Inc.     

Structural basis for proton conduction and inhibition by the influenza M2 protein

The influenza M2 protein forms an acid‐activated and drug‐sensitive proton channel in the virus envelope that is important for the virus lifecycle. The functional properties and high‐resolution structures of this proton channel have been extensively studied to understand the mechanisms of proton conduction and drug inhibition. We review biochemical and electrophysiological studies of M2 and discuss how high‐resolution structures have transformed our understanding of this proton channel. Comparison of structures obtained in different membrane‐mimetic solvents and under different pH using X‐ray crystallography, solution NMR, and solid‐state NMR spectroscopy revealed how the M2 structure depends on the environment and showed that the pharmacologically relevant drug‐binding site lies in the transmembrane (TM) pore. Competing models of proton conduction have been evaluated using biochemical experiments, high‐resolution structural methods, and computational modeling. These results are converging to a model in which a histidine residue in the TM domain mediates proton relay with water, aided by microsecond conformational dynamics of the imidazole ring. These mechanistic insights are guiding the design of new inhibitors that target drug‐resistant M2 variants and may be relevant for other proton channels.     

Identifying inter-residue resonances in crowded 2D 13C–13C chemical shift correlation spectra of membrane proteins by solid-state MAS NMR difference spectroscopy

The feasibility of using difference spectroscopy, i.e. subtraction of two correlation spectra at different mixing times, for substantially enhanced resolution in crowded two-dimensional ¹³C–¹³C chemical shift correlation spectra is presented. With the analyses of ¹³C–¹³C spin diffusion in simple spin systems, difference spectroscopy is proposed to partially separate the spin diffusion resonances of relatively short intra-residue distances from the longer inter-residue distances, leading to a better identification of the inter-residue resonances. Here solid-state magic-angle-spinning NMR spectra of the full length M2 protein embedded in synthetic lipid bilayers have been used to illustrate the resolution enhancement in the difference spectra. The integral membrane M2 protein of Influenza A virus assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH and is essential for the viral lifecycle. Based on known amino acid resonance assignments from amino acid specific labeled samples of truncated M2 sequences or from time-consuming 3D experiments of uniformly labeled samples, some inter-residue resonances of the full length M2 protein can be identified in the difference spectra of uniformly ¹³C labeled protein that are consistent with the high resolution structure of the M2 (22–62) protein (Sharma et al., Science 330(6003):509–512, 2010).     

Binding Hot Spots and Amantadine Orientation in the Influenza A Virus M2 Proton Channel

Structures of truncated versions of the influenza A virus M2 proton channel have been determined recently by x-ray crystallography in the open conformation of the channel, and by NMR in the closed state. The structures differ in the position of the bound inhibitors. The x-ray structure shows a single amantadine molecule in the middle of the channel, whereas in the NMR structure four drug molecules bind at the channel's outer surface. To study this controversy we applied computational solvent mapping, a technique developed for the identification of the most druggable binding hot spots of proteins. The method moves molecular probes—small organic molecules containing various functional groups—around the protein surface, finds favorable positions using empirical free energy functions, clusters the conformations, and ranks the clusters on the basis of the average free energy. The results of the mapping show that in both structures the primary hot spot is an internal cavity overlapping the amantadine binding site seen in the x-ray structure. However, both structures also have weaker hot spots at the exterior locations that bind rimantadine in the NMR structure, although these sites are partially due to the favorable interactions with the interfacial region of the lipid bilayer. As confirmed by docking calculations, the open channel binds amantadine at the more favorable internal site, in good agreement with the x-ray structure. In contrast, the NMR structure is based on a peptide/micelle construct that is able to accommodate the small molecular probes used for the mapping, but has a too narrow pore for the rimantadine to access the internal hot spot, and hence the drug can bind only at the exterior sites.     

Biological membrane structure by solid-state NMR

Nuclear magnetic resonance (NMR) spectroscopy, and particularly solid-state NMR spectroscopy, is a method of choice to study the structure and dynamics of both the lipid and the protein components of model and biological membranes. Different approaches have been developed to study these systems in which the restricted molecular motions result in broad NMR spectra. This contribution will first present an overview of the different techniques used to study lipid bilayers, namely 31p, 2H and 13C solid-state NMR spectroscopy. On the other hand, the study of the structure of membrane peptides and proteins is a rapidly growing field and several methods developed in the last two decades will be presented. These methods allow the investigation of protein systems for which structural information is often difficult to obtain by techniques such as X-ray diffraction and multidimensional solution NMR.     

pH-driven helix rotations in the influenza M2 H+ channel: a potential gating mechanism

The pH activated M2 H⁺ channel from influenza A has been a subject of numerous studies due to following: (1) It serves as a target for the aminoadamantane drugs that block its channel activity. (2) M2’s small size makes it amenable to biophysical scrutiny. (3) A single histidine residue is thought to control the pH gating of the channel. Recent FTIR analysis proposed that the helices of the channel rotate about their directors during pH activation. Herein, we report on molecular dynamics simulations of the X-ray structure of the protein with three charged histidine residues, representing the open form of the protein and two rotated forms with neutral histidines, representing its closed form. We compare the channel stability, convergence, interaction with water and hydration of the histidine residues that have been implicated in channel gating. Taken together, we show that both forms of the protein are stable during the course of the MD simulation and that indeed a rotation of the helices leads to channel closure. Finally, we propose a mechanism for channel gating that involves protonation of the histidine residues that necessities their increased solvation.     

Thermodynamic and structural behaviour of equimolar POPC/CnE4 (n = 8, 12, 16) mixtures by sorption gravimetry, 2H NMR spectroscopy and X-ray diffraction

The hydration behaviour of equimolar mixtures of phospholipids and nonionic surfactants with different chain length was investigated by gravimetric sorption, NMR spectroscopy and X-ray diffraction. At the most hydration degrees investigated, the incorporation of nonionic surfactants in a phospholipid bilayer leads to an increase of the hydrophilicity, which can be shown by the presence of excess hydration. The increased hydrophilicity could be explained by the excavation of additional water binding sites due to the “dilution” of the dipole field of the phospholipid bilayer. Another related contribution arises from the increase of the accessible surface area due to the increase of gauche conformers that result from the steric mismatch when surfactants are incorporated into the phospholipid matrix. ²H NMR spectroscopy was used to determine the quadrupolar splitting representing a measure of the order state of water. The swelling behaviour could be assessed by small-angle X-ray diffraction. ³¹P NMR spectroscopy was applied for the assignment of phase structures to the respective hydration range.     

Structure of Amantadine-Bound M2 Transmembrane Peptide of Influenza A in Lipid Bilayers from Magic-Angle-Spinning Solid-State NMR: The Role of Ser31 in Amantadine Binding

The M2 proton channel of influenza A is the target of the antiviral drugs amantadine and rimantadine, whose effectiveness has been abolished by a single-site mutation of Ser31 to Asn in the transmembrane domain of the protein. Recent high-resolution structures of the M2 transmembrane domain obtained from detergent-solubilized protein in solution and crystal environments gave conflicting drug binding sites. We present magic-angle-spinning solid-state NMR results of Ser31 and a number of other residues in the M2 transmembrane peptide (M2TMP) bound to lipid bilayers. Comparison of the spectra of the membrane-bound apo and complexed M2TMP indicates that Ser31 is the site of the largest chemical shift perturbation by amantadine. The chemical shift constraints lead to a monomer structure with a small kink of the helical axis at Gly34. A tetramer model is then constructed using the helix tilt angle and several interhelical distances previously measured on unoriented bilayer samples. This tetramer model differs from the solution and crystal structures in terms of the openness of the N-terminus of the channel, the constriction at Ser31, and the side-chain conformations of Trp41, a residue important for channel gating. Moreover, the tetramer model suggests that Ser31 may interact with amantadine amine via hydrogen bonding. While the apo and drug-bound M2TMP have similar average structures, the complexed peptide has much narrower linewidths at physiological temperature, indicating drug-induced changes of the protein dynamics in the membrane. Further, at low temperature, several residues show narrower lines in the complexed peptide than the apo peptide, indicating that amantadine binding reduces the conformational heterogeneity of specific residues. The differences of the current solid-state NMR structure of the bilayer-bound M2TMP from the detergent-based M2 structures suggest that the M2 conformation is sensitive to the environment, and care must be taken when interpreting structural findings from non-bilayer samples.     

Differential Lipid Dependence of the Function of Bacterial Sodium Channels

The lipid bilayer is important for maintaining the integrity of cellular compartments and plays a vital role in providing the hydrophobic and charged interactions necessary for membrane protein structure, conformational flexibility and function. To directly assess the lipid dependence of activity for voltage-gated sodium channels, we compared the activity of three bacterial sodium channel homologues (NaChBac, NavMs, and NavSp) by cumulative ²²Na⁺ uptake into proteoliposomes containing a 3∶1 ratio of 1-palmitoyl 2-oleoyl phosphatidylethanolamine and different “guest” glycerophospholipids. We observed a unique lipid profile for each channel tested. NavMs and NavSp showed strong preference for different negatively-charged lipids (phosphatidylinositol and phosphatidylglycerol, respectively), whilst NaChBac exhibited a more modest variation with lipid type. To investigate the molecular bases of these differences we used synchrotron radiation circular dichroism spectroscopy to compare structures in liposomes of different composition, and molecular modeling and electrostatics calculations to rationalize the functional differences seen. We then examined pore-only constructs (with voltage sensor subdomains removed) and found that in these channels the lipid specificity was drastically reduced, suggesting that the specific lipid influences on voltage-gated sodium channels arise primarily from their abilities to interact with the voltage-sensing subdomains.     

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.     

Protein-protein binding is often associated with changes in protonation state

pK(a) values of ionizable residues have been calculated using the PROPKA method and structures of 75 protein-protein complexes and their corresponding free forms. These pK(a) values were used to compute changes in protonation state of individual residues, net changes in protonation state of the complex relative to the uncomplexed proteins, and the correction to a binding energy calculated assuming standard protonation states at pH 7. For each complex, two different structures for the uncomplexed form of the proteins were used: the X-ray structures determined for the proteins in the absence of the other protein and the individual protein structures taken from the structure of the complex (referred to as unbound and bound structures, respectively). In 28 and 77% of the cases considered here, protein-protein binding is accompanied by a complete (>95%) or significant (>50%) change in protonation state of at least one residue using unbound structures. Furthermore, in 36 and 61% of the cases, protein-protein binding is accompanied by a complete or significant net change in protonation state of the complex relative to the separated monomers. Using bound structures, the corresponding values are 12, 51, 20, and 48%. Comparison to experimental data suggest that using unbound and bound structures lead to over- and underestimation of binding-induced protonation state changes, respectively. Thus, we conclude that protein-protein binding is often associated with changes in protonation state of amino acid residues and with changes in the net protonation state of the proteins. The pH-dependent correction to the binding energy contributes at least one order of magnitude to the binding constant in 45 and 23%, using unbound and bound structures, respectively.     

Solution structure and dynamics of GCN4 cognate DNA: NMR investigations

A 12 bp long GCN4-binding, self-complementary duplex DNA d(CATGACGTCATG)₂ has been investigated by NMR spectroscopy to study the structure and dynamics of the molecule in aqueous solution. The NMR structure of the DNA obtained using simulated annealing and iterative relaxation matrix calculations compares quite closely with the X-ray structure of ATF/CREB DNA in complex with GCN4 protein (DNA-binding domain). The DNA is also seen to be curved in the free state and this has a significant bearing on recognition by the protein. The dynamic characteristics of the molecule have been studied by ¹³C relaxation measurements at natural abundance. A correlation has been observed between sequence-dependent dynamics and recognition by GCN4 protein.     

Auxiliary Ca2+ binding sites can influence the structure of CIB1

Recent X-ray crystal structures and solution NMR spectroscopy data for calcium- and integrin-binding protein 1 (CIB1) have all revealed a common EF-hand domain structure for the protein. However, the orientation of the two protein domains, the oligomerization state, and the conformations of the N- and C-terminal extensions differ among the structures. In this study, we examine whether the binding of glutathione or auxiliary Ca²⁺ ions as observed in the crystal structures, occur in solution, and whether these interactions can influence the structure or dimerization of CIB1. In addition, we test the potential phosphatase activity of CIB1, which was hypothesized based on the glutathione binding site geometry observed in one of the crystal structures of the protein. Biophysical and biochemical experiments failed to detect glutathione binding, protein dimerization, or phosphatase activity for CIB1 under several solution conditions. However, our data identify low affinity (K_d, 10⁻²M) Ca²⁺ binding events that influence the structures of the N- and C-terminal extensions of CIB1 under high (300 mM) Ca²⁺ crystallization conditions. In addition to providing a rationale for differences amongst the various solution and crystal structures of CIB1, our results show that the impact of low affinity Ca²⁺ binding events should be considered when analyzing and interpreting protein crystallographic structures determined in the presence of very high Ca²⁺ concentrations.     

The chemical and dynamical influence of the anti-viral drug amantadine on the M2 proton channel transmembrane domain

The M(2) proton channel plays a vital role in the life cycle of the influenza A virus. His(37), the key residue in the M(2) transmembrane domain (M(2)-TMD), plays a central role in the proton conductance mechanism. The anti-influenza drug, amantadine, inhibits the channel activity through binding to the pore of the M(2) channel. The nuclear spin relaxation data and polarization inversion spin exchange at the magic angle spectra show that both the polypeptide backbone and His(37) side chain are more constrained in the presence of amantadine. Using (15)N cross polarization magic-angle spinning NMR spectroscopy, the protonation of His(37) of M(2)-TMD in lipid bilayers was monitored in the absence and presence of amantadine as a function of pH. Binding amantadine lowers the His(37) pK(a) values by approximately three orders of magnitude compared with the first pK(a) of histidine in amantadine-free M(2)-TMD. Amantadine's influence on the His(37) chemical properties suggests a novel mechanism by which amantadine may inhibit proton conductance.     

Conformational plasticity of the influenza A M2 transmembrane helix in lipid bilayers under varying pH, drug binding, and membrane thickness

Membrane proteins change their conformations to respond to environmental cues, thus conformational plasticity is important for function. The influenza A M2 protein forms an acid-activated proton channel important for the virus lifecycle. Here we have used solid-state NMR spectroscopy to examine the conformational plasticity of membrane-bound transmembrane domain of M2 (M2TM). (13)C and (15)N chemical shifts indicate coupled conformational changes of several pore-facing residues due to changes in bilayer thickness, drug binding, and pH. The structural changes are attributed to the formation of a well-defined helical kink at G34 in the drug-bound state and in thick lipid bilayers, nonideal backbone conformation of the secondary-gate residue V27 in the presence of drug, and nonideal conformation of the proton-sensing residue H37 at high pH. The chemical shifts constrained the (?, ψ) torsion angles for three "basis" states, the equilibrium among which explains the multiple resonances per site in the NMR spectra under different combinations of bilayer thickness, drug binding, and pH conditions. Thus, conformational plasticity is important for the proton conduction and inhibition of M2TM. The study illustrates the utility of NMR chemical shifts for probing the structural plasticity and folding of membrane proteins.     

The solution structure of desulforedoxin,a simple iron-sulfur protein

Desulforedoxin is a simple dimeric protein isolated from Desulfovibrio gigas containing a distorted rubredoxin-like center with one iron coordinated by four cysteinyl residues (7.9?kDa with a 36-amino-acid monomer). ¹H NMR spectra of the oxidized Dx(Fe³⁺) and reduced Dx(Fe²⁺) forms were analyzed. The spectra show substantial line broadening due to the paramagnetism of iron. However, very low-field-shifted resonances, assigned to Hβ protons, were observed in the reduced state and their temperature dependence analyzed. The active site of Dx was reconstituted with zinc, and its solution structure was determined using 2D NMR methods. This diamagnetic form gave high-resolution NMR data enabling the identification of all the amino acid spin systems. Sequential assignment and the determination of secondary structural elements was attempted using 2D NOESY experiments. However, because of the symmetrical dimer nature of the protein standard, NMR sequential assignment methods could not resolve all cross peaks due to inter- and intra-chain effects. The X-ray structure enabled the spatial relationship between the monomers to be obtained, and resolved the assignment problems. Secondary structural features could be identified from the NMR data; an antiparallel β-sheet running from D5 to V18 with a well-defined β-turn around cysteines C9 and C12. The section G22 to T25 is poorly defined by the NMR data and is followed by a turn around V27-C29. The C-terminus ends up near residues V6 and Y7. Distance geometry (DG) calculations allowed families of structures to be generated from the NMR data. A family of structures with a low target function violation for the Dx monomer and dimer were found to have secondary structural elements identical to those seen in the X-ray structure. The amide protons for G4, D5, G13, L11 NH and Q14 NHε amide protons, H-bonded in the X-ray structure, were not seen by NMR as slowly exchanging, while structural disorder at the N-terminus, for the backbone at E10 and for the section G22–T25, was observed. Comparison between the Fe and Zn forms of Dx suggests that metal substitution does not have an effect on the structure of the protein.