Structure and dynamics of the lipid modifications of a transmembrane α-helical peptide determined by ²H solid-state NMR spectroscopy

The fusion of biological membranes is mediated by integral membrane proteins with α-helical transmembrane segments. Additionally, those proteins are often modified by the covalent attachment of hydrocarbon chains. Previously, a series of de novo designed α-helical peptides with mixed Leu/Val sequences was presented, mimicking fusiogenically active transmembrane segments in model membranes (Hofmann et al., Proc. Natl. Acad. Sci. USA 101 (2004) 14776-14781). From this series, we have investigated the peptide LV16 (KKKW LVLV LVLV LVLV LVLV KKK), which was synthesized featuring either a free N-terminus or a saturated N-acylation of 2, 8, 12, or 16 carbons. We used 2H and 31P NMR spectroscopy to investigate the structure and dynamics of those peptide lipid modifications in POPC and DLPC bilayers and compared them to the hydrocarbon chains of the surrounding membrane. Except for the C2 chain, all peptide acyl chains were found to insert well into the membrane. This can be explained by the high local lipid concentrations the N-terminal lipid chains experience. Further, the insertion of these peptides did not influence the membrane structure and dynamics as seen from the 2H and 31P NMR data. In spite of the fact that the longer acyl chains insert into the membrane, they do not adapt their lengths to the thickness of the bilayer. Even the C16 lipid chain on the peptide, which could match the length of the POPC palmitoyl chain, exhibited lower order parameters in the upper chain, which get closer and finally reach similar values in the lower chain region. 2H NMR square law plots reveal motions of slightly larger amplitudes for the peptide lipid chains compared to the surrounding phospholipids. In spite of the significantly different chain lengths of the acylations, the fraction of gauche defects in the inserted chains is constant.     

²H solid-state nuclear magnetic resonance investigation of whole Escherichia coli interacting with antimicrobial peptide MSI-78

A key aspect of the activity of antimicrobial peptides (AMPs) is their interaction with membranes. Efforts to elucidate their detailed mechanisms have focused on applying biophysical methods, including nuclear magnetic resonance (NMR), to AMPs in model lipid systems. However, these highly simplified systems fail to capture many of the features of the much more complex cell envelopes with which AMPs interact in vivo. To address this issue, we have designed a procedure to incorporate high levels of (2)H NMR labels specifically into the cell membrane of Escherichia coli and used this approach to study the interactions between the AMP MSI-78 and the membranes of intact bacteria. The (2)H NMR spectra of these membrane-deuterated bacteria can be reproduced in the absence and presence of MSI-78. Because the (2)H NMR data provide a quantitative measure of lipid disorder, they directly report on the lipid bilayer disruption central to the function of AMPs, in the context of intact bacteria. Addition of MSI-78 to the bacteria leads to decreases in the order of the lipid acyl chains. The molar peptide:lipid ratios required to observe the effects of MSI-78 on acyl chain order are approximately 30 times greater than the ratios needed to observe effects in model lipid systems and approximately 100 times less than the ratios required to observe inhibition of cell growth in biological assays. The observations thus suggest that MSI-78 disrupts the bilayer even at sublethal AMP levels and that a large fraction of the peptide does not actually reach the inner membrane.     

Peptide-lipid interactions of the β-hairpin antimicrobial peptide tachyplesin and its linear derivatives from solid-state NMR

The peptide-lipid interaction of a β-hairpin antimicrobial peptide tachyplesin-1 (TP-1) and its linear derivatives are investigated to gain insight into the mechanism of antimicrobial activity. ³¹P and ²H NMR spectra of uniaxially aligned lipid bilayers of varying compositions and peptide concentrations are measured to determine the peptide-induced orientational disorder and the selectivity of membrane disruption by tachyplesin. The disulfide-linked TP-1 does not cause any disorder to the neutral POPC and POPC/cholesterol membranes but induces both micellization and random orientation distribution to the anionic POPE/POPG membranes above a peptide concentration of 2%. In comparison, the anionic POPC/POPG bilayer is completely unaffected by TP-1 binding, suggesting that TP-1 induces negative curvature strain to the membrane as a mechanism of its action. Removal of the disulfide bonds by substitution of Cys residues with Tyr and Ala abolishes the micellization of POPE/POPG bilayers but retains the orientation randomization of both POPC/POPG and POPE/POPG bilayers. Thus, linear tachyplesin derivatives have membrane disruptive abilities but use different mechanisms from the wild-type peptide. The different lipid-peptide interactions between TP-1 and other β-hairpin antimicrobial peptides are discussed in terms of their molecular structure.     

Investigating the membrane orientation and transversal distribution of 17β-estradiol in lipid membranes by solid-state NMR

17β-Estradiol (E₂) is a potent estrogen, which modulates many important cellular functions by binding to specific estrogen receptors located in the cell nucleus and also on the plasma membrane. We have studied the membrane interaction of E₂ using a combination of solid-state NMR methods. ²H NMR results indicate that E₂ does not cause a condensation effect of the surrounding phospholipids, which is contrary to the effects of cholesterol, and only very modest E₂ induced alterations of the membrane structure were detected. ¹H magic-angle spinning NMR showed well resolved signals from E₂ as well as of POPC in the membrane-lipid layer. Two-dimensional NOESY spectra revealed intense cross-peaks between E₂ and the membrane lipids indicating that E₂ is stably inserted into the membrane. The determination of intermolecular cross-relaxation rates revealed that E₂ is broadly distributed in the membrane with a maximum of the E₂ distribution function in the upper chain region of the membrane. We conclude that E₂ is highly dynamic in lipid membranes and may undergo rotations as it exhibits two polar hydroxyl groups on either side of the molecule.     

SC²ATmd: a tool for integration of the figure of merit with cluster analysis for gene expression data

Standard and Consensus Clustering Analysis Tool for Microarray Data (SC2ATmd) is a MATLAB-implemented application specifically designed for the exploration of microarray gene expression data via clustering. Implementation of two versions of the clustering validation method figure of merit allows for performance comparisons between different clustering algorithms, and tailors the cluster analysis process to the varying characteristics of each dataset. Along with standard clustering algorithms this application also offers a consensus clustering method that can generate reproducible clusters across replicate experiments or different clustering algorithms. This application was designed specifically for the analysis of gene expression data, but may be used with any numerical data as long as it is in the right format. AVAILABILITY: SC2ATmd may be freely downloaded from http://www.compbiosci.wfu.edu/tools.htm.     

Fluorescence and molecular dynamics studies of the acetylcholine receptor γM4 transmembrane peptide in reconstituted systems

A combination of fluorescence spectroscopy and molecular dynamics (MD) is applied to assess the conformational dynamics of a peptide making up the outermost ring of the nicotinic acetylcholine receptor (AChR) transmembrane region and the effect of membrane thickness and cholesterol on the hydrophobic matching of this peptide. The fluorescence studies exploit the intrinsic fluorescence of the only tryptophan residue in a synthetic peptide corresponding to the fourth transmembrane domain of the AChR γ subunit (γM4-Trp⁶) reconstituted in lipid bilayers of varying thickness, and combine this information with quenching studies using depth-sensitive phosphatidylcholine spin-labeled probes and acrylamide, polarization of fluorescence, and generalized polarization of Laurdan. A direct correlation was found between bilayer width and the depth of insertion of Trp⁶. We further extend our recent MD study of the conformational dynamics of the AChR channel to focus on the crosstalk between M4 and the lipid-belt region. The isolated γM4 peptide is shown to possess considerable orientational flexibility while maintaining a linear α-helical structure, and to vary its tilt depending on bilayer width and cholesterol (Chol) content. MD studies also show that γM4 also establishes contacts with the other TM peptides on its inner face, stabilizing a shorter TM length that is still highly sensitive to the lipid environment. In the native membrane the topology of the M4 ring is likely to exhibit a similar behavior, dynamically modifying its tilt to match the hydrophobic thickness of the bilayer.     

¹H, ¹³C, and ¹⁵N NMR resonance assignments of reduced full length and shortened forms of the Grx domain of Mus musculus TGR

Two forms of the glutaredoxin (Grx) domain (full length Grx domain and short Grx lacking the N-terminal region) of Mus musculus thioredoxin glutathione reductase (TGR) were isotopically labelled with ¹⁵N and ¹³C isotopes, expressed and purified to homogeneity. We report here the ¹H, ¹³C and ¹⁵N NMR assignment for both Grx forms of this mouse TGR. This investigation represents the first NMR analysis of a mammalian TGR.     

High-yield expression and purification of the transmembrane region of ion channel-forming amyloid-β protein for NMR structural studies

Numerous studies of APP in the literature to date have focused on its processing at the plasma membrane, its membrane-bound oligomeric state, and the calcium-permeable ion channel formation of non-fibrillar Aβ in the cell membrane. Despite these studies, little is known structurally about the transmembrane region of APP beyond a theoretical model. This is due to challenges in the expression, purification, and sample preparation of eukaryotic membrane proteins. To determine the three-dimensional structures of the intact transmembrane domain from human APP (hAPP-TM) and to elucidate the structure and formation mechanisms of the hAPP channel, the isotopically labeled protein must be produced and purified in sufficient quantity. Here, we describe a procedure whereby the hAPP-TM peptide, comprising residues 692–723 of hAPP, was successfully expressed and purified sufficiently to perform NMR analysis. To increase expression levels of the target protein, we designed a construct containing two tandem repeats of the target gene. The fusion protein was expressed in the form of inclusion bodies, purified on immobilized nickel affinity chromatography, and chemically cleaved by cyanogen bromide. Final purification of hAPP-TM was achieved by preparative reversed-phase high performance liquid chromatography (HPLC). The final yields of purified hAPP-TM were around 5 mg/l of M9 minimal media.     

Deep-UV resonance Raman analysis of the Rhodobacter capsulatus cytochrome bc₁complex reveals a potential marker for the transmembrane peptide backbone

Classical strategies for structure analysis of proteins interacting with a lipid phase typically correlate ensemble secondary structure content measurements with changes in the spectroscopic responses of localized aromatic residues or reporter molecules to map regional solvent environments. Deep-UV resonance Raman (DUVRR) spectroscopy probes the vibrational modes of the peptide backbone itself, is very sensitive to the ensemble secondary structures of a protein, and has been shown to be sensitive to the extent of solvent interaction with the peptide backbone [ Wang , Y. , Purrello , R. , Georgiou , S. , and Spiro , T. G. ( 1991 ) J. Am. Chem. Soc. 113 , 6368 - 6377 ]. Here we show that a large detergent solubilized membrane protein, the Rhodobacter capsulatus cytochrome bc(1) complex, has a distinct DUVRR spectrum versus that of an aqueous soluble protein with similar overall secondary structure content. Cross-section calculations of the amide vibrational modes indicate that the peptide backbone carbonyl stretching modes differ dramatically between these two proteins. Deuterium exchange experiments probing solvent accessibility confirm that the contribution of the backbone vibrational mode differences are derived from the lipid solubilized or transmembrane α-helical portion of the protein complex. These findings indicate that DUVRR is sensitive to both the hydration status of a protein's peptide backbone, regardless of primary sequence, and its secondary structure content. Therefore, DUVRR may be capable of simultaneously measuring protein dynamics and relative water/lipid solvation of the protein.     

IRAP — A system for the retrieval and analysis of orthopaedic implant data

An interactive computer system for the storage, retrieval and analysis of standardized clinical and material characterization data associated with orthopaedic implants is described. The system consists basically of four independent modules. The essence of the system centers on the cross-referencing capabilities which are virtually unlimited. The system has been designed for use by non-computer trained personnel.     

Solid-state ²H NMR shows equivalence of dehydration and osmotic pressures in lipid membrane deformation

Lipid bilayers represent a fascinating class of biomaterials whose properties are altered by changes in pressure or temperature. Functions of cellular membranes can be affected by nonspecific lipid-protein interactions that depend on bilayer material properties. Here we address the changes in lipid bilayer structure induced by external pressure. Solid-state 2H NMR spectroscopy of phospholipid bilayers under osmotic stress allows structural fluctuations and deformation of membranes to be investigated. We highlight the results from NMR experiments utilizing pressure-based force techniques that control membrane structure and tension. Our 2H NMR results using both dehydration pressure (low water activity) and osmotic pressure (poly(ethylene glycol) as osmolyte) show that the segmental order parameters (S(CD)) of DMPC approach very large values of ≈ 0.35 in the liquid-crystalline state. The two stresses are thermodynamically equivalent, because the change in chemical potential when transferring water from the interlamellar space to the bulk water phase corresponds to the induced pressure. This theoretical equivalence is experimentally revealed by considering the solid-state 2H NMR spectrometer as a virtual osmometer. Moreover, we extend this approach to include the correspondence between osmotic pressure and hydrostatic pressure. Our results establish the magnitude of the pressures that lead to significant bilayer deformation including changes in area per lipid and volumetric bilayer thickness. We find that appreciable bilayer structural changes occur with osmotic pressures in the range of 10-100 atm or lower. This research demonstrates the applicability of solid-state 2H NMR spectroscopy together with bilayer stress techniques for investigating the mechanism of pressure sensitivity of membrane proteins.     

HYPER: A hierarchical algorithm for automatic determination of protein dihedral-angle constraints and stereospecific CβH2 resonance assignments from NMR data

A new computer program, HYPER, has been developed for automated analysis of protein dihedral angle values and CH₂ stereospecific assignments from NMR data. HYPER uses a hierarchical grid-search algorithm to determine allowed values of , , and ₁ dihedral angles and CH₂ stereospecific assignments based on a set of NMR-derived distance and/or scalar-coupling constraints. Dihedral-angle constraints are valuable for restricting conformational space and improving convergence in three-dimensional structure calculations. HYPER computes the set of , , and ₁dihedral angles and CH₂ stereospecific assignments that are consistent with up to nine intraresidue and sequential distance bounds, two pairs of relative distance bounds, thirteen homo- and heteronuclear scalar coupling bounds, and two pairs of relative scalar coupling constant bounds. The program is designed to be very flexible, and provides for simple user modification of Karplus equations and standard polypeptide geometries, allowing it to accommodate recent and future improved calibrations of Karplus curves. The C code has been optimized to execute rapidly (0.3–1.5 CPU-sec residue^–1 using a 5° grid) on Silicon Graphics R8000, R10000 and Intel Pentium CPUs, making it useful for interactive evaluation of inconsistent experimental constraints. The HYPER program has been tested for internal consistency and reliability using both simulated and real protein NMR data sets.     

Molecular dynamics simulation approach for the prediction of transmembrane helix–helix heterodimers assembly

Computational methods are useful to identify favorable structures of transmembrane (TM) helix oligomers when experimental data are not available or when they cannot help to interpret helix-helix association. We report here a global search method using molecular dynamics (MD) simulations to predict the structures of transmembrane homo and heterodimers. The present approach is based only on sequence information without any experimental data and is first applied to glycophorin A to validate the protocol and to the HER2-HER3 heterodimer receptor. The method successfully reproduces the experimental structures of the TM domain of glycophorin A (GpA(TM)) with a root mean square deviation of 1.5 A. The search protocol identifies three energetically stable models of the TM domain of HER2-HER3 receptor with favorable helix-helix arrangement, including right-handed and left-handed coiled-coils. The predicted TM structures exhibit the GxxxG-like motif at the dimer interface which is presumed to drive receptor oligomerization. We demonstrate that native structures of TM domain can be predicted without quantitative experimental data. This search protocol could help to predict structures of the TM domain of HER heterodimer family.     

Anesthetic effects on the structure and dynamics of the second transmembrane domains of nAChR α4β2

Channel functions of the neuronal α4β2 nicotinic acetylcholine receptor (nAChR), one of the most widely expressed subtypes in the brain, can be inhibited by volatile anesthetics. Our Na⁺ flux experiments confirmed that the second transmembrane domains (TM2) of α4 and β2 in 2:3 stoichiometry, (α4)₂(β2)₃, could form pentameric channels, whereas the α4 TM2 alone could not. The structure, topology, and dynamics of the α4 TM2 and (α4)₂(β2)₃ TM2 in magnetically aligned phospholipid bicelles were investigated using solid-state NMR spectroscopy in the absence and presence of halothane and isoflurane, two clinically used volatile anesthetics. ²H NMR demonstrated that anesthetics increased lipid conformational heterogeneity. Such anesthetic effects on lipids became more profound in the presence of transmembrane proteins. PISEMA experiments on the selectively ¹⁵N-labeled α4 TM2 showed that the TM2 formed transmembrane helices with tilt angles of 12° ± 1° and 16° ± 1° relative to the bicelle normal for the α4 and (α4)₂(β2)₃ samples, respectively. Anesthetics changed the tilt angle of the α4 TM2 from 12° ± 1° to 14° ± 1°, but had only a subtle effect on the tilt angle of the (α4)₂(β2)₃ TM2. A small degree of wobbling motion of the helix axis occurred in the (α4)₂(β2)₃ TM2. In addition, a subset of the (α4)₂(β2)₃ TM2 exhibited counterclockwise rotational motion around the helix axis on a time scale slower than 10^- 4 s in the presence of anesthetics. Both helical tilting and rotational motions have been identified computationally as critical elements for ion channel functions. This study suggested that anesthetics could alter these motions to modulate channel functions.     

Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain

With improvements in computer speed and algorithm efficiency, MD simulations are sampling larger amounts of molecular and biomolecular conformations. Being able to qualitatively and quantitatively sift these conformations into meaningful groups is a difficult and important task, especially when considering the structure-activity paradigm. Here we present a study that combines two popular techniques, principal component (PC) analysis and clustering, for revealing major conformational changes that occur in molecular dynamics (MD) simulations. Specifically, we explored how clustering different PC subspaces effects the resulting clusters versus clustering the complete trajectory data. As a case example, we used the trajectory data from an explicitly solvated simulation of a bacteria’s L11·23S ribosomal subdomain, which is a target of thiopeptide antibiotics. Clustering was performed, using K-means and average-linkage algorithms, on data involving the first two to the first five PC subspace dimensions. For the average-linkage algorithm we found that data-point membership, cluster shape, and cluster size depended on the selected PC subspace data. In contrast, K-means provided very consistent results regardless of the selected subspace. Since we present results on a single model system, generalization concerning the clustering of different PC subspaces of other molecular systems is currently premature. However, our hope is that this study illustrates a) the complexities in selecting the appropriate clustering algorithm, b) the complexities in interpreting and validating their results, and c) by combining PC analysis with subsequent clustering valuable dynamic and conformational information can be obtained.     

NMR structures of the human α7 nAChR transmembrane domain and associated anesthetic binding sites

The α7 nicotinic acetylcholine receptor (nAChR), assembled as homomeric pentameric ligand-gated ion channels, is one of the most abundant nAChR subtypes in the brain. Despite its importance in memory, learning and cognition, no structure has been determined for the α7 nAChR TM domain, a target for allosteric modulators. Using solution state NMR, we determined the structure of the human α7 nAChR TM domain (PDB ID: 2MAW) and demonstrated that the α7 TM domain formed functional channels in Xenopus oocytes. We identified the associated binding sites for the anesthetics halothane and ketamine; the former cannot sensitively inhibit α7 function, but the latter can. The α7 TM domain folds into the expected four-helical bundle motif, but the intra-subunit cavity at the extracellular end of the α7 TM domain is smaller than the equivalent cavity in the α4β2 nAChRs (PDB IDs: 2LLY; 2LM2). Neither drug binds to the extracellular end of the α7 TM domain, but two halothane molecules or one ketamine molecule binds to the intracellular end of the α7 TM domain. Halothane and ketamine binding sites are partially overlapped. Ketamine, but not halothane, perturbed the α7 channel-gate residue L9′. Furthermore, halothane did not induce profound dynamics changes in the α7 channel as observed in α4β2. The study offers a novel high-resolution structure for the human α7 nAChR TM domain that is invaluable for developing α7-specific therapeutics. It also provides evidence to support the hypothesis: only when anesthetic binding perturbs the channel pore or alters the channel motion, can binding generate functional consequences.     

On the properties of Se ⋯N interaction: the analysis of substituent effects by energy decomposition and orbital interaction

The nature and strength of intermolecular Se ?N interaction between selenium-containing compounds HSeX (X = CH₃, NH₂, CF₃, OCH₃, CN, OH, NO₂, Cl, F), and NH₃ have been investigated at the MP2/aug-cc-pVDZ level. The Se ?N interaction is found to be dependent on the substituent groups, which greatly affect the positive electrostatic potential of Se atoms and the accepting electron ability of X-Se σ ^? antibonding orbital. Energy decomposition of the Se ?N interaction reveals that electrostatic and induction forces are comparable in the weak-bonded complexes and induction becomes more significant in the complexes with strong electron-withdrawing substituents. Natural bond orbital (NBO) analysis indicates that the primary source of the induction is the electron transfer from the N lone pair to the X-Se σ ^? antibonding orbital. The geometry of the complex and the interaction directionality of NH₃ to X-Se bond can be regarded as a consequence of the exchange-repulsion. The topological analysis on the electron density reveals the nature of closed-shell interaction in these X-Se ?N contacts. The Se ?N interaction in the complexes with the strong electron-withdrawing substituent has a partly covalent character.     

Optimization of expression,purification, and NMR measurement for structural studies of transmembrane region of amyloid β protein

The amyloid precursor protein (APP) is an integral transmembrane protein which has been suggested to play a central role in the pathogenesis of Alzheimer’s disease. Despite the enormous amount of research conducted on amyloid protein, the precise mechanism of its toxic effect is not yet fully understood. To better understand the mechanism and function of amyloid protein, it is critical to elucidate the three-dimensional structure of the single transmembrane spanning region of human APP (hAPP-TM). Unfortunately, it is difficult to prepare the peptide sample because hAPP-TM is a membrane-bound protein that transverses the lipid bilayer of the cell membrane. Generally, the preparation of a transmembrane peptide is very difficult and time-consuming. In fact, high yield production of transmembrane peptides has been limited by experimental difficulties related to insufficient yields and the low solubility of such peptides. In this study, we describe experimental processes developed to optimize the expression, purification, and NMR measurement conditions for hAPP-TM transmembrane peptide.