Fast protein dynamics can be missed with techniques that have relatively slow observation times. Using 2D IR spectroscopy and isotope labeling, we have probed the rapid, picosecond dynamics of a membrane protein in its native environment. By measuring the homogeneous and inhomogeneous IR linewidths of 11 amide I modes (backbone carbonyl stretch), we have captured the structural distributions and dynamics of the CD3ζ protein along its transmembrane segment that are lost with slower time-scale techniques. We find that the homogeneous lifetimes and population relaxation times are the same for almost all of the residues. In contrast, the inhomogeneous linewidths vary significantly with the largest inhomogeneous distribution occurring for residues near the N terminus and the narrowest near the center. This behavior is highly consistent with a recently reported experimental model of the protein and water accessibility as observed by molecular dynamics simulations. The data support the proposed CD3ζ peptide structure, and the simulations point to the structural disorder of water and lipid head-groups as the main source of inhomogeneous broadening. Taken together, this rigorous analysis of the vibrational dynamics of a membrane peptide provides experimental insight into a time regime of motions that has so far been largely unexplored.

A molecular-level understanding of the function of a protein requires knowledge of both its structural and dynamic properties. NMR spectroscopy allows the measurement of generalized order parameters that provide an atomistic description of picosecond and nanosecond fluctuations in protein structure. Molecular dynamics (MD) simulation provides a complementary approach to the study of protein dynamics on similar time scales. Comparisons between NMR spectroscopy and MD simulations can be used to interpret experimental results and to improve the quality of simulation-related force fields and integration methods. However, apparent systematic discrepancies between order parameters extracted from simulations and experiments are common, particularly for elements of noncanonical secondary structure. In this paper, results from a 1.2 μs explicit solvent MD simulation of the protein ubiquitin are compared with previously determined backbone order parameters derived from NMR relaxation experiments [Tjandra, N.; Feller, S. E.; Pastor, R. W.; Bax, A. J. Am. Chem. Soc.1995, 117, 12562−12566]. The simulation reveals fluctuations in three loop regions that occur on time scales comparable to or longer than that of the overall rotational diffusion of ubiquitin and whose effects would not be apparent in experimentally derived order parameters. A coupled analysis of internal and overall motion yields simulated order parameters substantially closer to the experimentally determined values than is the case for a conventional analysis of internal motion alone. Improved agreement between simulation and experiment also is encouraging from the viewpoint of assessing the accuracy of long MD simulations.

Statistical analyses were undertaken for putative transmembrane alpha-helices obtained from a database representing the subset of membrane proteins available in Swiss-Prot. The average length of a transmembrane alpha-helix was found to be 22-21 amino acids with a large variation around the mean. The transfer free energy from water to oil of a transmembrane alpha-helix in bitopic proteins, -48 kcal/mol, is higher than that in polytopic proteins, -39 kcal/mol, and is nearly identical to that obtained by assuming a random distribution of solely hydrophobic amino acids in the alpha-helix. The amino acid composition of hydrophobic residues is similar in bitopic and polytopic proteins. In contrast, the more polar the amino acids are, the less likely they are to be found in bitopic proteins compared to polytopic ones. This most likely reflects the ability of alpha-helical bundles to shield the polarity of residues from the hydrophobic bilayer. One half of all amino acids were distributed nonrandomly in both bitopic and polytopic proteins. A preference was found for tyrosine and tryptophan residues to be at the ends of transmembrane alpha-helices. Correlated distribution analysis of amino acid pairs indicated that most amino acids are independently distributed in each helix. Exceptions are cysteine, tyrosine, and tryptophan which appear to cluster closely to one another and glycines which are preferentially found on the same side of alpha-helices.

Glycosylphosphatidylinositol (GPI) anchoring of proteins provides a potential mechanism for targeting to the plant plasma membrane and cell wall. However, relatively few such proteins have been identified. Here, we develop a procedure for database analysis to identify GPI-anchored proteins (GAP) based on their possession of common features. In a comprehensive search of the annotated Arabidopsis genome, we identified 167 novel putative GAP in addition to the 43 previously described candidates. Many of these 210 proteins show similarity to characterized cell surface proteins. The predicted GAP include homologs of β-1,3-glucanases (16), metallo- and aspartyl proteases (13), glycerophosphodiesterases (6), phytocyanins (25), multi-copper oxidases (2), extensins (6), plasma membrane receptors (19), and lipid-transfer-proteins (18). Classical arabinogalactan (AG) proteins (13), AG peptides (9), fasciclin-like proteins (20), COBRA and 10 homologs, and novel potential signaling peptides that we name GAPEPs (8) were also identified. A further 34 proteins of unknown function were predicted to be GPI anchored. A surprising finding was that over 40% of the proteins identified here have probable AG glycosylation modules, suggesting that AG glycosylation of cell surface proteins is widespread. This analysis shows that GPI anchoring is likely to be a major modification in plants that is used to target a specific subset of proteins to the cell surface for extracellular matrix remodeling and signaling.

Na ⁺ /H ⁺ antiporters are central to cellular salt and pH homeostasis. The structure of Escherichia coli NhaA was recently determined, but its mechanisms of transport and pH regulation remain elusive. We performed molecular dynamics simulations of NhaA that, with existing experimental data, enabled us to propose an atomically detailed model of antiporter function. Three conserved aspartates are key to our proposed mechanism: Asp ¹⁶⁴ (D164) is the Na ⁺ -binding site, D163 controls the alternating accessibility of this binding site to the cytoplasm or periplasm, and D133 is crucial for pH regulation. Consistent with experimental stoichiometry, two protons are required to transport a single Na ⁺ ion: D163 protonates to reveal the Na ⁺ -binding site to the periplasm, and subsequent protonation of D164 releases Na ⁺ . Additional mutagenesis experiments further validated the model.