Ez, a Depth-dependent Potential for Assessing the Energies of Insertion of Amino Acid Side-chains into Membranes: Derivation and Applications to Determining the Orientation of Transmembrane and Interfacial Helices

Mutational scanning reveals the determinants of protein insertion and association energetics in the plasma membrane

Insertion of helix-forming segments into the membrane and their association determines the structure, function, and expression levels of all plasma membrane proteins. However, systematic and reliable quantification of membrane-protein energetics has been challenging. We developed a deep mutational scanning method to monitor the effects of hundreds of point mutations on helix insertion and self-association within the bacterial inner membrane. The assay quantifies insertion energetics for all natural amino acids at 27 positions across the membrane, revealing that the hydrophobicity of biological membranes is significantly higher than appreciated. We further quantitate the contributions to membrane-protein insertion from positively charged residues at the cytoplasm-membrane interface and reveal large and unanticipated differences among these residues. Finally, we derive comprehensive mutational landscapes in the membrane domains of Glycophorin A and the ErbB2 oncogene, and find that insertion and self-association are strongly coupled in receptor homodimers.

Membrane positioning for high- and low-resolution protein structures through a binary classification approach

The critical importance of algorithms for orienting proteins in the lipid bilayer stems from the extreme difficulty in obtaining experimental data about the membrane boundaries. Here, we present a computational method for positioning protein structures in the membrane, based on the sole alpha carbon coordinates and, therefore, compatible with both high and low structural resolutions. Our algorithm follows a new and simple approach, by treating the membrane assignment problem as a binary classification. Compared with the state-of-the-art algorithms, our method achieves similar accuracy, while being faster. Finally, our open-source software is also capable of processing coarse-grained models of protein structures.

An Integrated Framework Advancing Membrane Protein Modeling and Design

Membrane proteins are critical functional molecules in the human body, constituting more than 30% of open reading frames in the human genome. Unfortunately, a myriad of difficulties in overexpression and reconstitution into membrane mimetics severely limit our ability to determine their structures. Computational tools are therefore instrumental to membrane protein structure prediction, consequently increasing our understanding of membrane protein function and their role in disease. Here, we describe a general framework facilitating membrane protein modeling and design that combines the scientific principles for membrane protein modeling with the flexible software architecture of Rosetta3. This new framework, called RosettaMP, provides a general membrane representation that interfaces with scoring, conformational sampling, and mutation routines that can be easily combined to create new protocols. To demonstrate the capabilities of this implementation, we developed four proof-of-concept applications for (1) prediction of free energy changes upon mutation; (2) high-resolution structural refinement; (3) protein-protein docking; and (4) assembly of symmetric protein complexes, all in the membrane environment. Preliminary data show that these algorithms can produce meaningful scores and structures. The data also suggest needed improvements to both sampling routines and score functions. Importantly, the applications collectively demonstrate the potential of combining the flexible nature of RosettaMP with the power of Rosetta algorithms to facilitate membrane protein modeling and design.

Computational redesign of the lipid-facing surface of the outer membrane protein OmpA

Advances in computational design methods have made possible extensive engineering of soluble proteins, but designed β-barrel membrane proteins await improvements in our understanding of the sequence determinants of folding and stability. A subset of the amino acid residues of membrane proteins interact with the cell membrane, and the design rules that govern this lipid-facing surface are poorly understood. We applied a residue-level depth potential for β-barrel membrane proteins to the complete redesign of the lipid-facing surface of Escherichia coli OmpA. Initial designs failed to fold correctly, but reversion of a small number of mutations indicated by backcross experiments yielded designs with substitutions to up to 60% of the surface that did support folding and membrane insertion.

Amphipathic solvation of indole: implications for the role of tryptophan in membrane proteins

The microscopic structure of the tryptophan side chain, indole, in an amphiphilic environment has been investigated using a combination of neutron diffraction measurements and simulations in solution. The results show that indole is preferentially solvated by hydrogen bonding interactions between water and alcohol -OH groups rather than the interaction being dominated by indole-methyl interactions. This has implications for understanding how tryptophan interacts with the amphipathic membrane environment to anchor proteins into membranes, where the results here suggest that the benzene ring of tryptophan interacts directly with the interfacial water at the membrane surface rather than being buried into the hydrophobic regions of the membrane bilayer.

De novo design of a transmembrane Zn²⁺-transporting four-helix bundle

The design of functional membrane proteins from first principles represents a grand challenge in chemistry and structural biology. Here, we report the design of a membrane-spanning, four-helical bundle that transports first-row transition metal ions Zn(2+) and Co(2+), but not Ca(2+), across membranes. The conduction path was designed to contain two di-metal binding sites that bind with negative cooperativity. X-ray crystallography and solid-state and solution nuclear magnetic resonance indicate that the overall helical bundle is formed from two tightly interacting pairs of helices, which form individual domains that interact weakly along a more dynamic interface. Vesicle flux experiments show that as Zn(2+) ions diffuse down their concentration gradients, protons are antiported. These experiments illustrate the feasibility of designing membrane proteins with predefined structural and dynamic properties.

De novo membrane protein structure prediction

Recent advances in identifying residue-residue contacts from large multiple sequence alignments have enabled impressive gains to be made in the field of protein structure prediction. In this chapter, we discuss these advances and provide a step-by-step guide to applying the latest tools to the de novo modelling of alpha-helical transmembrane proteins. As a practical example, we demonstrate the process of building an accurate 3D model of a G protein-coupled receptor, correctly orientated in the membrane, using only its primary protein sequence.

Evaluating membrane affinity by integrating protein orientations

Energetic interactions of a protein with lipid bilayers determine its propensity to reside in the membrane. Here we seek to evaluate the membrane interactions for EMAPII, a protein found to be released from the cell by unknown mechanisms, as well as several other proteins. Using a knowledge-based coarse-grained membrane potential, we calculate the free energy profiles for these proteins by integrating out the orientation degrees of freedom. Due to the invariance of energy under in-plane rotations about the membrane normal, the orientation space can be reduced to two dimensions and mapped onto the surface of a unit sphere, thus making visualization, sampling and integration more convenient. The integrated free energy profiles determine the relative probabilities along the membrane normal for the proteins regardless of their orientations, and display distinctive characteristics for membrane proteins and water-soluble proteins. The membrane interactions for EMAPII exhibit typical features of a water-soluble protein, with a high energetic barrier to enter or cross the membrane. Our results thus suggest that similar to the non-classical export of FGF1, the release of EMAPII would involve more complicated mechanisms than simple passive diffusion across the membrane.

A Gly-zipper motif mediates homodimerization of the transmembrane domain of the mitochondrial kinase ADCK3

Interactions between α-helices within the hydrophobic environment of lipid bilayers are integral to the folding and function of transmembrane proteins; however, the major forces that mediate these interactions remain debated, and our ability to predict these interactions is still largely untested. We recently demonstrated that the frequent transmembrane association motif GASright, the GxxxG-containing fold of the glycophorin A dimer, is optimal for the formation of extended networks of Cα-H hydrogen bonds, supporting the hypothesis that these bonds are major contributors to association. We also found that optimization of Cα-H hydrogen bonding and interhelical packing is sufficient to computationally predict the structure of known GASright dimers at near atomic level. Here, we demonstrate that this computational method can be used to characterize the structure of a protein not previously known to dimerize, by predicting and validating the transmembrane dimer of ADCK3, a mitochondrial kinase. ADCK3 is involved in the biosynthesis of the redox active lipid, ubiquinone, and human ADCK3 mutations cause a cerebellar ataxia associated with ubiquinone deficiency, but the biochemical functions of ADCK3 remain largely undefined. Our experimental analyses show that the transmembrane helix of ADCK3 oligomerizes, with an interface based on an extended Gly-zipper motif, as predicted by our models. The data provide strong evidence for the hypothesis that optimization of Cα-H hydrogen bonding is an important factor in the association of transmembrane helices. This work also provides a structural foundation for investigating the role of transmembrane association in regulating the biological activity of ADCK3.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

Prediction of linear cationic antimicrobial peptides based on characteristics responsible for their interaction with the membranes

Most available antimicrobial peptides (AMP) prediction methods use common approach for different classes of AMP. Contrary to available approaches, we suggest that a strategy of prediction should be based on the fact that there are several kinds of AMP that vary in mechanisms of action, structure, mode of interaction with membrane, etc. According to our suggestion for each kind of AMP, a particular approach has to be developed in order to get high efficacy. Consequently, in this paper, a particular but the biggest class of AMP, linear cationic antimicrobial peptides (LCAP), has been considered and a newly developed simple method of LCAP prediction described. The aim of this study is the development of a simple method of discrimination of AMP from non-AMP, the efficiency of which will be determined by efficiencies of selected descriptors only and comparison the results of the discrimination procedure with the results obtained by more complicated discriminative methods. As descriptors the physicochemical characteristics responsible for capability of the peptide to interact with an anionic membrane were considered. The following characteristics such as hydrophobicity, amphiphaticity, location of the peptide in relation to membrane, charge density, propensities to disordered structure and aggregation were studied. On the basis of these characteristics, a new simple algorithm of prediction is developed and evaluation of efficacies of the characteristics as descriptors performed. The results show that three descriptors, hydrophobic moment, charge density and location of the peptide along the membranes, can be used as discriminators of LCAPs. For the training set, our method gives the same level of accuracy as more complicated machine learning approaches offered as CAMP database service tools. For the test set accuracy obtained by our method gives even higher value than the one obtained by CAMP prediction tools. The AMP prediction tool based on the considered method is available at http://www.biomedicine.org.ge/dbaasp/.

Structural basis of the proinflammatory signaling complex mediated by TSLP

Thymic stromal lymphopoietin (TSLP), a cytokine produced by epithelial cells at barrier surfaces, is pivotal for the development of widespread chronic inflammatory disorders such as asthma and atopic dermatitis. The structure of the mouse TSLP-mediated signaling complex reveals how TSLP establishes extensive interfaces with its cognate receptor (TSLPR) and the shared interleukin 7 receptor α-chain (IL-7Rα) to evoke membrane-proximal receptor-receptor contacts poised for intracellular signaling. Binding of TSLP to TSLPR is a mechanistic prerequisite for recruitment of IL-7Rα to the high-affinity ternary complex, which we propose is coupled to a structural switch in TSLP at the crossroads of the cytokine-receptor interfaces. Functional interrogation of TSLP-receptor interfaces points to putative interaction hotspots that could be exploited for antagonist design. Finally, we derive the structural rationale for the functional duality of IL-7Rα and establish a consensus for the geometry of ternary complexes mediated by interleukin 2 (IL-2)-family cytokines.

Unique functional properties of conserved arginine residues in the lentivirus lytic peptide domains of the C-terminal tail of HIV-1 gp41

A previous study from our laboratory reported a preferential conservation of arginine relative to lysine in the C-terminal tail (CTT) of HIV-1 envelope (Env). Despite substantial overall sequence variation in the CTT, specific arginines are highly conserved in the lentivirus lytic peptide (LLP) motifs and are scarcely substituted by lysines, in contrast to gp120 and the ectodomain of gp41. However, to date, no explanation has been provided to explain the selective incorporation and conservation of arginines over lysines in these motifs. Herein, we address the functions in virus replication of the most conserved arginines by performing conservative mutations of arginine to lysine in the LLP1 and LLP2 motifs. The presence of lysine in place of arginine in the LLP1 motif resulted in significant impairment of Env expression and consequently virus replication kinetics, Env fusogenicity, and incorporation. By contrast, lysine exchanges in LLP2 only affected the level of Env incorporation and fusogenicity. Our findings demonstrate that the conservative lysine substitutions significantly affect Env functional properties indicating a unique functional role for the highly conserved arginines in the LLP motifs. These results provide for the first time a functional explanation to the preferred incorporation of arginine, relative to lysine, in the CTT of HIV-1 Env. We propose that these arginines may provide unique functions for Env interaction with viral or cellular cofactors that then influence overall Env functional properties.

A frequent, GxxxG-mediated, transmembrane association motif is optimized for the formation of interhelical Cα-H hydrogen bonds

Carbon hydrogen bonds between Cα-H donors and carbonyl acceptors are frequently observed between transmembrane helices (Cα-H···O=C). Networks of these interactions occur often at helix-helix interfaces mediated by GxxxG and similar patterns. Cα-H hydrogen bonds have been hypothesized to be important in membrane protein folding and association, but evidence that they are major determinants of helix association is still lacking. Here we present a comprehensive geometric analysis of homodimeric helices that demonstrates the existence of a single region in conformational space with high propensity for Cα-H···O=C hydrogen bond formation. This region corresponds to the most frequent motif for parallel dimers, GASright, whose best-known example is glycophorin A. The finding suggests a causal link between the high frequency of occurrence of GASright and its propensity for carbon hydrogen bond formation. Investigation of the sequence dependency of the motif determined that Gly residues are required at specific positions where only Gly can act as a donor with its "side chain" Hα. Gly also reduces the steric barrier for non-Gly amino acids at other positions to act as Cα donors, promoting the formation of cooperative hydrogen bonding networks. These findings offer a structural rationale for the occurrence of GxxxG patterns at the GASright interface. The analysis identified the conformational space and the sequence requirement of Cα-H···O=C mediated motifs; we took advantage of these results to develop a structural prediction method. The resulting program, CATM, predicts ab initio the known high-resolution structures of homodimeric GASright motifs at near-atomic level.

Characterization of LdrA (long direct repeat A) protein of Escherichia coli

Reminiscent of apoptotic genes in higher systems, almost all bacteria contain 'suicide' genes, which encode toxins. Toxins inhibit cell growth and may cause cell death. These are coexpressed with their cognate antitoxins from a toxin-antitoxin (TA) operon in normally growing cells. The cellular targets of toxins are diverse. The study of the TA system is important as in addition to their involvement in the normal bacterial physiology, TA systems may also play a role in bacterial pathogenicity. The long direct repeat (Ldr) family of proteins belongs to one such TA system. Here we report characterization of LdrA, the prototypical protein of the Ldr family, and show that it is highly toxic for cell growth. The data suggests that LdrA may exert toxicity by inhibiting ATP synthesis, possibly due to its localization in the cell membrane. This inhibition of energy production leads to the inhibition of biosynthetic reactions in the cell such as DNA replication, transcription and translation, and eventually cell growth.

Structure-activity relationships of the antimicrobial peptide gramicidin S and its analogs: aqueous solubility, self-association, conformation, antimicrobial activity and interaction with model lipid membranes

GS10 [cyclo-(VKLdYPVKLdYP)] is a synthetic analog of the naturally occurring antimicrobial peptide gramicidin (GS) in which the two positively charged ornithine (Orn) residues are replaced by two positively charged lysine (Lys) residues and the two less polar aromatic phenylalanine (Phe) residues are replaced by the more polar tyrosine (Tyr) residues. In this study, we examine the effects of these seemingly conservative modifications to the parent GS molecule on the physical properties of the peptide, and on its interactions with lipid bilayer model and biological membranes, by a variety of biophysical techniques. We show that although GS10 retains the largely β-sheet conformation characteristic of GS, it is less structured in both water and membrane-mimetic solvents. GS10 is also more water soluble and less hydrophobic than GS, as predicted, and also exhibits a reduced tendency for self-association in aqueous solution. Surprisingly, GS10 associates more strongly with zwitterionic and anionic phospholipid bilayer model membranes than does GS, despite its greater water solubility, and the presence of anionic phospholipids and cholesterol (Chol) modestly reduces the association of both GS10 and GS to these model membranes. The strong partitioning of both peptides into lipid bilayers is driven by a large favorable entropy change opposed by a much smaller unfavorable enthalpy change. However, GS10 is also less potent than GS at inducing inverted cubic phases in phospholipid bilayer model membranes and at inhibiting the growth of the cell wall-less bacterium Acholeplasma laidlawii B. These results are discussed in terms of the comparative antibiotic and hemolytic activities of these peptides.

Structural dynamics and topology of phosphorylated phospholamban homopentamer reveal its role in the regulation of calcium transport

Phospholamban (PLN) inhibits the sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA), thereby regulating cardiac diastole. In membranes, PLN assembles into homopentamers that in both the phosphorylated and nonphosphorylated states have been proposed to form ion-selective channels. Here, we determined the structure of the phosphorylated pentamer using a combination of solution and solid-state nuclear magnetic resonance methods. We found that the pinwheel architecture of the homopentamer is preserved upon phosphorylation, with each monomer having an L-shaped conformation. The TM domains form a hydrophobic pore approximately 24 Å long and 2 Å in diameter, which is inconsistent with canonical Ca²⁺-selective channels. Phosphorylation, however, enhances the conformational dynamics of the cytoplasmic region of PLN, causing partial unwinding of the amphipathic helix. We propose that PLN oligomers act as storage for active monomers, keeping SERCA function within a physiological window.

Determination of structural topology of a membrane protein in lipid bilayers using polarization optimized experiments (POE) for static and MAS solid state NMR spectroscopy

The low sensitivity inherent to both the static and magic angle spinning techniques of solid-state NMR (ssNMR) spectroscopy has thus far limited the routine application of multidimensional experiments to determine the structure of membrane proteins in lipid bilayers. Here, we demonstrate the advantage of using a recently developed class of experiments, polarization optimized experiments, for both static and MAS spectroscopy to achieve higher sensitivity and substantial time-savings for 2D and 3D experiments. We used sarcolipin, a single pass membrane protein, reconstituted in oriented bicelles (for oriented ssNMR) and multilamellar vesicles (for MAS ssNMR) as a benchmark. The restraints derived by these experiments are then combined into a hybrid energy function to allow simultaneous determination of structure and topology. The resulting structural ensemble converged to a helical conformation with a backbone RMSD ~0.44 Å, a tilt angle of 24° ± 1°, and an azimuthal angle of 55° ± 6°. This work represents a crucial first step toward obtaining high-resolution structures of large membrane proteins using combined multidimensional oriented solid-state NMR and magic angle spinning solid-state NMR.

Membrane protein orientation and refinement using a knowledge-based statistical potential

BACKGROUND

Recent increases in the number of deposited membrane protein crystal structures necessitate the use of automated computational tools to position them within the lipid bilayer. Identifying the correct orientation allows us to study the complex relationship between sequence, structure and the lipid environment, which is otherwise challenging to investigate using experimental techniques due to the difficulty in crystallising membrane proteins embedded within intact membranes.

RESULTS

We have developed a knowledge-based membrane potential, calculated by the statistical analysis of transmembrane protein structures, coupled with a combination of genetic and direct search algorithms, and demonstrate its use in positioning proteins in membranes, refinement of membrane protein models and in decoy discrimination.

CONCLUSIONS

Our method is able to quickly and accurately orientate both alpha-helical and beta-barrel membrane proteins within the lipid bilayer, showing closer agreement with experimentally determined values than existing approaches. We also demonstrate both consistent and significant refinement of membrane protein models and the effective discrimination between native and decoy structures. Source code is available under an open source license from http://bioinf.cs.ucl.ac.uk/downloads/memembed/.

Structure-function analysis of MurJ reveals a solvent-exposed cavity containing residues essential for peptidoglycan biogenesis in Escherichia coli

Gram-negative bacteria such as Escherichia coli build a peptidoglycan (PG) cell wall in their periplasm using the precursor known as lipid II. Lipid II is a large amphipathic molecule composed of undecaprenyl diphosphate and a disaccharide-pentapeptide that PG-synthesizing enzymes use to build the PG sacculus. During PG biosynthesis, lipid II is synthesized at the cytoplasmic face of the inner membrane and then flipped across the membrane. This translocation of lipid II must be assisted by flippases thought to shield the disaccharide-pentapeptide as it crosses the hydrophobic core of the membrane. The inner membrane protein MurJ is essential for PG biogenesis and homologous to known and putative flippases of the MOP (multidrug/oligo-saccharidyl-lipid/polysaccharide) exporter superfamily, which includes flippases that translocate undecaprenyl diphosphate-linked oligosaccharides across the cytoplasmic membranes of bacteria. Consequently, MurJ has been proposed to function as the lipid II flippase in E. coli. Here, we present a three-dimensional structural model of MurJ generated by the I-TASSER server that suggests that MurJ contains a solvent-exposed cavity within the plane of the membrane. Using in vivo topological studies, we demonstrate that MurJ has 14 transmembrane domains and validate features of the MurJ structural model, including the presence of a solvent-exposed cavity within its transmembrane region. Furthermore, we present functional studies demonstrating that specific charged residues localized in the central cavity are essential for function. Together, our studies support the structural homology of MurJ to MOP exporter proteins, suggesting that MurJ might function as an essential transporter in PG biosynthesis.

Structural adaptations of proteins to different biological membranes

To gain insight into adaptations of proteins to their membranes, intrinsic hydrophobic thicknesses, distributions of different chemical groups and profiles of hydrogen-bonding capacities (α and β) and the dipolarity/polarizability parameter (π*) were calculated for lipid-facing surfaces of 460 integral α-helical, β-barrel and peripheral proteins from eight types of biomembranes. For comparison, polarity profiles were also calculated for ten artificial lipid bilayers that have been previously studied by neutron and X-ray scattering. Estimated hydrophobic thicknesses are 30-31Å for proteins from endoplasmic reticulum, thylakoid, and various bacterial plasma membranes, but differ for proteins from outer bacterial, inner mitochondrial and eukaryotic plasma membranes (23.9, 28.6 and 33.5Å, respectively). Protein and lipid polarity parameters abruptly change in the lipid carbonyl zone that matches the calculated hydrophobic boundaries. Maxima of positively charged protein groups correspond to the location of lipid phosphates at 20-22Å distances from the membrane center. Locations of Tyr atoms coincide with hydrophobic boundaries, while distributions maxima of Trp rings are shifted by 3-4Å toward the membrane center. Distributions of Trp atoms indicate the presence of two 5-8Å-wide midpolar regions with intermediate π* values within the hydrocarbon core, whose size and symmetry depend on the lipid composition of membrane leaflets. Midpolar regions are especially asymmetric in outer bacterial membranes and cell membranes of mesophilic but not hyperthermophilic archaebacteria, indicating the larger width of the central nonpolar region in the later case. In artificial lipid bilayers, midpolar regions are observed up to the level of acyl chain double bonds.

Charge asymmetry in the proteins of the outer membrane

MOTIVATION

Outer membrane beta-barrels (OMBBs) are the proteins found in the outer membrane of bacteria, mitochondria and chloroplasts. There are thousands of beta-barrels reported in genomic databases with ∼2-3% of the genes in gram-negative bacteria encoding these proteins. These proteins have a wide variety of biological functions including active and passive transport, cell adhesion, catalysis and structural anchoring. Of the non-redundant OMBB structures in the Protein Data Bank, half have been solved during the past 5 years. This influx of information provides new opportunities for understanding the chemistry of these proteins. The distribution of charges in proteins in the outer membrane has implications for how the mechanism of outer membrane protein insertion is understood. Understanding the distribution of charges might also assist in organism selection for the heterologous expression of mitochondrial OMBBs.

RESULTS

We find a strong asymmetry in the charge distribution of these proteins. For the outward-facing residues of the beta-barrel within regions of similar amino acid density for both membrane leaflets, the external side of the outer membrane contains almost three times the number of charged residues as the internal side of the outer membrane. Moreover, the lipid bilayer of the outer membrane is asymmetric, and the overall preference for amino acid types to be in the external leaflet of the membrane correlates roughly with the hydrophobicity of the membrane lipids. This preference is demonstrably related to the difference in lipid composition of the external and internal leaflets of the membrane.

Membrane depth-dependent energetic contribution of the tryptophan side chain to the stability of integral membrane proteins

Lipid solvation provides the primary driving force for the insertion and folding of integral membrane proteins. Although the structure of the lipid bilayer is often simplified as a central hydrophobic core sandwiched between two hydrophilic interfacial regions, the complexity of the liquid-crystalline bilayer structure and the gradient of water molecules across the bilayer fine-tune the energetic contributions of individual amino acid residues to the stability of membrane proteins at different depths of the bilayer. The tryptophan side chain is particularly interesting because despite its widely recognized role in anchoring membrane proteins in lipid bilayers, there is little consensus about its hydrophobicity among various experimentally determined hydrophobicity scales. Here we investigated how lipid-facing tryptophan residues located at different depths in the bilayer contribute to the stability of integral membrane proteins using outer membrane protein A (OmpA) as a model. We replaced all lipid-contacting residues of the first transmembrane β-strand of OmpA with alanines and individually incorporated tryptophans in these positions along the strand. By measuring the thermodynamic stability of these proteins, we found that OmpA is slightly more stable when tryptophans are placed in the center of the bilayer and that it is somewhat destabilized as tryptophans approach the interfacial region. However, this trend may be partially reversed when a moderate concentration of urea rather than water is taken as the reference state. The measured stability profiles are driven by similar profiles of the m-value, a parameter that reflects the shielding of hydrophobic surface area from water. Our results indicate that knowledge of the free energy level of the protein's unfolded reference state is important for quantitatively assessing the stability of membrane proteins, which may explain differences in observed profiles between in vivo and in vitro scales.

Structural organization of FtsB, a transmembrane protein of the bacterial divisome

We report the first structural analysis of an integral membrane protein of the bacterial divisome. FtsB is a single-pass membrane protein with a periplasmic coiled coil. Its heterologous association with its partner FtsL represents an essential event for the recruitment of the late components to the division site. Using a combination of mutagenesis, computational modeling, and X-ray crystallography, we determined that FtsB self-associates, and we investigated its structural organization. We found that the transmembrane domain of FtsB homo-oligomerizes through an evolutionarily conserved interaction interface where a polar residue (Gln 16) plays a critical role through the formation of an interhelical hydrogen bond. The crystal structure of the periplasmic domain, solved as a fusion with Gp7, shows that 30 juxta-membrane amino acids of FtsB form a canonical coiled coil. The presence of conserved Gly residue in the linker region suggests that flexibility between the transmembrane and coiled coil domains is functionally important. We hypothesize that the transmembrane helices of FtsB form a stable dimeric core for its association with FtsL into a higher-order oligomer and that FtsL is required to stabilize the periplasmic domain of FtsB, leading to the formation of a complex that is competent for binding to FtsQ, and to their consequent recruitment to the divisome. The study provides an experimentally validated structural model and identifies point mutations that disrupt association, thereby establishing important groundwork for the functional characterization of FtsB in vivo.

Simultaneous prediction of protein secondary structure and transmembrane spans

Prediction of transmembrane spans and secondary structure from the protein sequence is generally the first step in the structural characterization of (membrane) proteins. Preference of a stretch of amino acids in a protein to form secondary structure and being placed in the membrane are correlated. Nevertheless, current methods predict either secondary structure or individual transmembrane states. We introduce a method that simultaneously predicts the secondary structure and transmembrane spans from the protein sequence. This approach not only eliminates the necessity to create a consensus prediction from possibly contradicting outputs of several predictors but bears the potential to predict conformational switches, i.e., sequence regions that have a high probability to change for example from a coil conformation in solution to an α-helical transmembrane state. An artificial neural network was trained on databases of 177 membrane proteins and 6048 soluble proteins. The output is a 3 × 3 dimensional probability matrix for each residue in the sequence that combines three secondary structure types (helix, strand, coil) and three environment types (membrane core, interface, solution). The prediction accuracies are 70.3% for nine possible states, 73.2% for three-state secondary structure prediction, and 94.8% for three-state transmembrane span prediction. These accuracies are comparable to state-of-the-art predictors of secondary structure (e.g., Psipred) or transmembrane placement (e.g., OCTOPUS). The method is available as web server and for download at www.meilerlab.org.

Structural insights into the common γ-chain family of cytokines and receptors from the interleukin-7 pathway

Over the past 13 years, numerous crystal structures of complexes of the common γ-chain (γ(c)) cytokine receptors and their cytokines have been solved. Even with the remarkable progress in the structural biology of γ(c) receptors and their cytokines or interleukins, there are valuable lessons to be learned from the structural and biophysical studies of interleukin-7 (IL-7) and its α-receptor (IL-7Rα) and comparisons with other γ(c) family members. The structure of the IL-7/IL-7Rα complex teaches that interfaces between the γ(c) interleukins and their receptors can vary in size, polarity, and specificity, and that significant conformational changes might be necessary for complexes of interleukins and their receptors to bind the shared, activating γ(c) receptor. Binding, kinetic, and thermodynamic studies of IL-7 and IL-7Rα show that glycosylation and electrostatics can be important to interactions between interleukins and their receptor, even where the glycans and charged residues are distant from the interface. The structure of the IL-7Rα homodimer is a reminder that often-ignored non-activating complexes likely perform roles just as important to signaling as activating complexes. And last but not least, the structural and biophysical studies help explain and potentially treat the diseases caused by aberrant IL-7 signaling.

Assembly of the transmembrane domain of E. coli PhoQ histidine kinase: implications for signal transduction from molecular simulations

The PhoQP two-component system is a signaling complex essential for bacterial virulence and cationic antimicrobial peptide resistance. PhoQ is the histidine kinase chemoreceptor of this tandem machine and assembles in a homodimer conformation spanning the bacterial inner membrane. Currently, a full understanding of the PhoQ signal transduction is hindered by the lack of a complete atomistic structure. In this study, an atomistic model of the key transmembrane (TM) domain is assembled by using molecular simulations, guided by experimental cross-linking data. The formation of a polar pocket involving Asn202 in the lumen of the tetrameric TM bundle is crucial for the assembly and solvation of the domain. Moreover, a concerted displacement of the TM helices at the periplasmic side is found to modulate a rotation at the cytoplasmic end, supporting the transduction of the chemical signal through a combination of scissoring and rotational movement of the TM helices.

Two-component systems (TCSs) are signaling complexes essential for bacterial survival and virulence. PhoQ is the histidine kinase chemoreceptor of the PhoQ-PhoP tandem machine that detects the concentration of cationic species at the inner membrane of Gram-negative bacteria. A full understanding of the PhoQ signal transduction mechanism is currently hindered by the lack of a complete atomistic structure. Here, by using molecular simulations integrated with cross-linking disulfide scanning data, we present the first structural model of the transmembrane (TM) portion of PhoQ from E. coli. Its structural and dynamic features induce a concerted displacement of the TM helices at the periplasmic side, which modulates a rotation at the cytoplasmic end. This supports the idea that signal transduction is promoted through a combination of scissoring and rotational movements of the TM helices. This complex mechanism is the key to understanding how the chemical stimuli sensed by the periplasmic sensor domain trigger, via the relay of the HAMP domain, the histidine auto-phosphorylation and kinase/phosphatase activity at the cytoplasmic end.

Painting proteins blue: β-(1-azulenyl)-L-alanine as a probe for studying protein-protein interactions

We demonstrated that β-(1-azulenyl)-L-alanine, a fluorescent pseudoisosteric analog of tryptophan, exhibits weak environmental dependence and thus allows for using weak intrinsic quenchers, such as methionines, to monitor protein-protein interactions while not perturbing them.

Analyzing the effects of hydrophobic mismatch on transmembrane α-helices using tryptophan fluorescence spectroscopy

Hydrophobic matching between transmembrane protein segments and the lipid bilayer in which they are embedded is a significant factor in the behavior and orientation of such transmembrane segments. The condition of hydrophobic mismatch occurs when the hydrophobic thickness of a lipid bilayer is significantly different than the length of the membrane spanning segment of a protein, resulting in a mismatch. This mismatch can result in altered function of proteins as well as nonnative structural arrangements including effects on transmembrane α-helix tilt angles, oligomerization state, and/or the formation of non-transmembrane topographies. Here, a fluorescence-based protocol is described for testing model transmembrane α-helices and their sensitivity to hydrophobic mismatch by measuring the propensity of these helices to form non-transmembrane structures. Overall, good hydrophobic matching between the bilayer and transmembrane segments is an important factor that must be considered when designing membrane proteins or peptides.

Prediction and design of outer membrane protein-protein interactions

Protein-protein interactions (PPI) play central roles in biological processes, motivating us to understand the structural basis underlying affinity and specificity. In this chapter, we focus on biochemical and computational design strategies of assessing and detecting PPIs of β-barrel outer membrane proteins (OMPs). A few case studies are presented highlighting biochemical techniques used to dissect the energetics of oligomerization and determine amino acids forming the key interactions of the PPI sites. Current computational strategies for detecting/predicting PPIs are introduced, and examples of computational and rational engineering strategies applied to OMPs are presented.

Solvation models and computational prediction of orientations of peptides and proteins in membranes

Membrane-associated peptides and proteins function in the highly heterogeneous environment of the lipid bilayer whose physico-chemical properties change non-monotonically along the bilayer normal. To simulate insertion of peptides and proteins into membranes and correctly reproduce the energetics of this process, an adequate solvation model and physically realistic representation of the lipid bilayer should be employed. We present a brief overview of the existing solvation models and their application for prediction of binding affinities and orientations of proteins in membranes. Particular emphasis is placed on the recently proposed PPM method, the corresponding web server, and the OPM database that were designed for positioning in membranes of integral and peripheral proteins with known three-dimensional structures.

Knowledge-based potential for positioning membrane-associated structures and assessing residue-specific energetic contributions

The complex hydrophobic and hydrophilic milieus of membrane-associated proteins pose experimental and theoretical challenges to their understanding. Here, we produce a nonredundant database to compute knowledge-based asymmetric cross-membrane potentials from the per-residue distributions of C(β), C(γ) and functional group atoms. We predict transmembrane and peripherally associated regions from genomic sequence and position peptides and protein structures relative to the bilayer (available at http://www.degradolab.org/ez). The pseudo-energy topological landscapes underscore positional stability and functional mechanisms demonstrated here for antimicrobial peptides, transmembrane proteins, and viral fusion proteins. Moreover, experimental effects of point mutations on the relative ratio changes of dual-topology proteins are quantitatively reproduced. The functional group potential and the membrane-exposed residues display the largest energetic changes enabling to detect native-like structures from decoys. Hence, focusing on the uniqueness of membrane-associated proteins and peptides, we quantitatively parameterize their cross-membrane propensity, thus facilitating structural refinement, characterization, prediction, and design.

Structural informatics, modeling, and design with an open-source Molecular Software Library (MSL)

We present the Molecular Software Library (MSL), a C++ library for molecular modeling. MSL is a set of tools that supports a large variety of algorithms for the design, modeling, and analysis of macromolecules. Among the main features supported by the library are methods for applying geometric transformations and alignments, the implementation of a rich set of energy functions, side chain optimization, backbone manipulation, calculation of solvent accessible surface area, and other tools. MSL has a number of unique features, such as the ability of storing alternative atomic coordinates (for modeling) and multiple amino acid identities at the same backbone position (for design). It has a straightforward mechanism for extending its energy functions and can work with any type of molecules. Although the code base is large, MSL was created with ease of developing in mind. It allows the rapid implementation of simple tasks while fully supporting the creation of complex applications. Some of the potentialities of the software are demonstrated here with examples that show how to program complex and essential modeling tasks with few lines of code. MSL is an ongoing and evolving project, with new features and improvements being introduced regularly, but it is mature and suitable for production and has been used in numerous protein modeling and design projects. MSL is open-source software, freely downloadable at http://msl-libraries.org. We propose it as a common platform for the development of new molecular algorithms and to promote the distribution, sharing, and reutilization of computational methods.

Protein Structure in Membrane Domains

Of great interest to the academic and pharmaceutical research communities, helical transmembrane proteins are characterized by their ability to dissolve and fold in lipid bilayers—properties conferred by polypeptide spans termed transmembrane domains (TMDs). The apolar nature of TMDs necessitates the use of membrane-mimetic solvents for many structure and folding studies. This review examines the relationship between TMD structure and solvent environment, focusing on principles elucidated largely in membrane-mimetic environments with single-TMD protein and peptide models. Following a brief description of TMD sequence and conformational characteristics gleaned from the structural database, we present an overview of the conceptual models used to study folding in vitro. The impact of sequence and solvent context on the incorporation of TMDs into membranes, and its role in measurements of TMD self-assembly strengths, is then described. We conclude with a discussion of the nonspecific effects of membrane components on TMD stability.

Computational design of membrane proteins

Membrane proteins are involved in a wide variety of cellular processes, and are typically part of the first interaction a cell has with extracellular molecules. As a result, these proteins comprise a majority of known drug targets. Membrane proteins are among the most difficult proteins to obtain and characterize, and a structure-based understanding of their properties can be difficult to elucidate. Notwithstanding, the design of membrane proteins can provide stringent tests of our understanding of these crucial biological systems, as well as introduce novel or targeted functionalities. Computational design methods have been particularly helpful in addressing these issues, and this review discusses recent studies that tailor membrane proteins to display specific structures or functions and examines how redesigned membrane proteins are being used to facilitate structural and functional studies.

Biosynthesis and spectroscopic characterization of 2-TM fragments encompassing the sequence of a human GPCR, the Y4 receptor

This paper presents a divide-and-conquer approach towards obtaining solution structures of G protein-coupled receptors. The human Y4 receptor was dissected into two to three transmembrane helix fragments, which were individually studied by solution NMR. We systematically compared various biosynthetic routes for the expression of the fragments in Escherichia coli and discuss purification strategies. In particular, we have compared the production of transmembrane (TM) fragments as inclusion bodies by using the ΔTrp leader sequence, with membrane-directed expression by using Mistic as the fusion partner, and developed methods for enzymatic cleavage. In addition, direct expression of two-TM fragments into inclusion bodies is a successful route in some cases. With the exception of TM13, we could produce all fragments in isotope-labeled form in quantities sufficient for NMR studies. Almost complete backbone resonance assignment was obtained for the first two helices, as well as for helices 5 and 7, and a high degree was obtained for TM6, while conformational exchange processes resulted in the disappearance of many signals from TM4. In addition, complete assignments were obtained for all residues of the N-terminal domain, as well as the extracellular and cytosolic loops (with the exception of an undecapeptide segment in the second extracellular loop, EC2) and for the complete cytosolic C-terminal tail. In total, backbone resonances of 78 % of all residues were assigned for the Y4 receptor. Predictions of secondary structure based on backbone chemical shifts indicate that most residues from the TM regions adopt helical conformations, with exception of those around polar residues or prolines. However, the domain boundaries differ slightly from those predicted for homology models. We suggest that the obtained chemical shifts might be useful in assigning the full-length receptor.

Modulation of function in a minimalist heme-binding membrane protein

De novo designed heme-binding proteins have been used successfully to recapitulate features of natural hemoproteins. This approach has now been extended to membrane-soluble model proteins. Our group designed a functional hemoprotein, ME1, by engineering a bishistidine binding site into a natural membrane protein, glycophorin A (Cordova et al. in J Am Chem Soc 129:512-518, 2007). ME1 binds iron(III) protoporphyrin IX with submicromolar affinity, has a redox potential of -128 mV, and displays peroxidase activity. Here, we show the effect of aromatic residues in modulating the redox potential in the context of a membrane-soluble model system. We designed aromatic interactions with the heme through a single-point mutant, G25F, in which a phenylalanine is designed to dock against the porphyrin ring. This mutation results in roughly tenfold tighter binding to iron(III) protoporphyrin IX (K(d,app) = 6.5 × 10(-8) M), and lowers the redox potential of the cofactor to -172 mV. This work demonstrates that specific design features aimed at controlling the properties of bound cofactors can be introduced in a minimalist membrane hemoprotein model. The ability to modulate the redox potential of cofactors embedded in artificial membrane proteins is crucial for the design of electron transfer chains across membranes in functional photosynthetic devices.

Structural reorganization of the interleukin-7 signaling complex

We report here an unliganded receptor structure in the common gamma-chain (γ(c)) family of receptors and cytokines. The crystal structure of the unliganded form of the interleukin-7 alpha receptor (IL-7Rα) extracellular domain (ECD) at 2.15 Å resolution reveals a homodimer forming an "X" geometry looking down onto the cell surface with the C termini of the two chains separated by 110 Å and the dimer interface comprising residues critical for IL-7 binding. Further biophysical studies indicate a weak association of the IL-7Rα ECDs but a stronger association between the γ(c)/IL-7Rα ECDs, similar to previous studies of the full-length receptors on CD4(+) T cells. Based on these and previous results, we propose a molecular mechanism detailing the progression from the inactive IL-7Rα homodimer and IL-7Rα-γ(c) heterodimer to the active IL-7-IL-7Rα-γ(c) ternary complex whereby the two receptors undergo at least a 90° rotation away from the cell surface, moving the C termini of IL-7Rα and γ(c) from a distance of 110 Å to less than 30 Å at the cell surface. This molecular mechanism can be used to explain recently discovered IL-7- and γ(c)-independent gain-of-function mutations in IL-7Rα from B- and T-cell acute lymphoblastic leukemia patients. The mechanism may also be applicable to other γ(c) receptors that form inactive homodimers and heterodimers independent of their cytokines.

Neutralizing epitopes in the membrane-proximal external region of HIV-1 gp41 are influenced by the transmembrane domain and the plasma membrane

Failure to elicit broadly neutralizing (bNt) antibodies (Abs) against the membrane-proximal external region of HIV-1 gp41 (MPER) reflects the difficulty of mimicking its neutralization-competent structure (NCS). Here, we analyzed MPER antigenicity in the context of the plasma membrane and identified a role for the gp41 transmembrane domain (TM) in exposing the epitopes of three bNt monoclonal Abs (MAbs) (2F5, 4E10, and Z13e1). We transiently expressed DNA constructs encoding gp41 ectodomain fragments fused to either the TM of the platelet-derived growth factor receptor (PDGFR) or the gp41 TM and cytoplasmic tail domain (CT). Constructs encoding the MPER tethered to the gp41 TM followed by a 27-residue CT fragment (MPER-TM1) produced optimal MAb binding. Critical binding residues for the three Nt MAbs were identified using a panel of 24 MPER-TM1 mutants bearing single amino acid substitutions in the MPER; many were previously shown to affect MAb-mediated viral neutralization. Moreover, non-Nt mutants of MAbs 2F5 and 4E10 exhibited a reduction in binding to MPER-TM1 and yet maintained binding to synthetic MPER peptides, indicating that MPER-TM1 better approximates the MPER NCS than peptides. Replacement of the gp41 TM and CT of MPER-TM1 with the PDGFR TM reduced binding by MAb 4E10, but not 2F5, indicating that the gp41 TM plays a pivotal role in orienting the 4E10 epitope, and more globally, in affecting MPER exposure.

A knowledge-based potential highlights unique features of membrane α-helical and β-barrel protein insertion and folding

Outer membrane β-barrel proteins differ from α-helical inner membrane proteins in lipid environment, secondary structure, and the proposed processes of folding and insertion. It is reasonable to expect that outer membrane proteins may contain primary sequence information specific for their folding and insertion behavior. In previous work, a depth-dependent insertion potential, E(z) , was derived for α-helical inner membrane proteins. We have generated an equivalent potential for TM β-barrel proteins. The similarities and differences between these two potentials provide insight into unique aspects of the folding and insertion of β-barrel membrane proteins. This potential can predict orientation within the membrane and identify functional residues involved in intermolecular interactions.

Computational studies of membrane proteins: models and predictions for biological understanding

We discuss recent progresses in computational studies of membrane proteins based on physical models with parameters derived from bioinformatics analysis. We describe computational identification of membrane proteins and prediction of their topology from sequence, discovery of sequence and spatial motifs, and implications of these discoveries. The detection of evolutionary signal for understanding the substitution pattern of residues in the TM segments and for sequence alignment is also discussed. We further discuss empirical potential functions for energetics of inserting residues in the TM domain, for interactions between TM helices or strands, and their applications in predicting lipid-facing surfaces of the TM domain. Recent progresses in structure predictions of membrane proteins are also reviewed, with further discussions on calculation of ensemble properties such as melting temperature based on simplified state space model. Additional topics include prediction of oligomerization state of membrane proteins, identification of the interfaces for protein-protein interactions, and design of membrane proteins. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Efficient molecular mechanics simulations of the folding, orientation, and assembly of peptides in lipid bilayers using an implicit atomic solvation model

Membrane proteins comprise a significant fraction of the proteomes of sequenced organisms and are the targets of approximately half of marketed drugs. However, in spite of their prevalence and biomedical importance, relatively few experimental structures are available due to technical challenges. Computational simulations can potentially address this deficit by providing structural models of membrane proteins. Solvation within the spatially heterogeneous membrane/solvent environment provides a major component of the energetics driving protein folding and association within the membrane. We have developed an implicit solvation model for membranes that is both computationally efficient and accurate enough to enable molecular mechanics predictions for the folding and association of peptides within the membrane. We derived the new atomic solvation model parameters using an unbiased fitting procedure to experimental data and have applied it to diverse problems in order to test its accuracy and to gain insight into membrane protein folding. First, we predicted the positions and orientations of peptides and complexes within the lipid bilayer and compared the simulation results with solid-state NMR structures. Additionally, we performed folding simulations for a series of host-guest peptides with varying propensities to form alpha helices in a hydrophobic environment and compared the structures with experimental measurements. We were also able to successfully predict the structures of amphipathic peptides as well as the structures for dimeric complexes of short hexapeptides that have experimentally characterized propensities to form beta sheets within the membrane. Finally, we compared calculated relative transfer energies with data from experiments measuring the effects of mutations on the free energies of translocon-mediated insertion of proteins into lipid bilayers and of combined folding and membrane insertion of a beta barrel protein.

Oligomerization state of photosynthetic core complexes is correlated with the dimerization affinity of a transmembrane helix

In the Rhodobacter (Rba.) species of photosynthetic purple bacteria, a single transmembrane α-helix, PufX, is found within the core complex, an essential photosynthetic macromolecular assembly that performs the absorption and the initial processing of light energy. Despite its structural simplicity, many unresolved questions surround PufX, the most important of which is its location within the photosynthetic core complex. One proposed placement of PufX is at the center of a core complex dimer, where two PufX helices associate in the membrane and form a homodimer. Inability for PufX of certain Rba. species to form a homodimer is thought to lead to monomeric core complexes. In the present study, we employ a combination of computational and experimental techniques to test the hypothesized homodimerization of PufX. We carry out a systematic investigation to measure the dimerization affinity of PufX from four Rba. species, Rba. blasticus , Rba. capsulatus , Rba. sphaeroides , and Rba. veldkampii , using a molecular dynamics-based free-energy method, as well as experimental TOXCAT assays. We found that the four PufX helices have substantially different dimerization affinities. Both computational and experimental techniques demonstrate that species with dimeric core complexes have PufX that can potentially form a homodimer, whereas the one species with monomeric core complexes has a PufX with little to no dimerization propensity. Our analysis of the helix-helix interface revealed a number of positions that may be important for PufX dimerization and the formation of a hydrogen-bond network between these GxxxG-containing helices. Our results suggest that the different oligomerization states of core complexes in various Rba. species can be attributed, among other factors, to the different propensity of its PufX helix to homodimerize.

OPM database and PPM web server: resources for positioning of proteins in membranes

The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topology and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approximately 350 organisms that represent approximately 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielectric constant and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calculating spatial positions in membranes of newly determined proteins structures or theoretical models.

Computational design of membrane proteins

This article reviews the recent successes of computational protein design techniques applied to integral membrane proteins. This emerging area is still handicapped by significant difficulties in the experimental characterization of the stability and structure of the designed proteins. Nevertheless, by focusing on oligomeric complexes of single-span transmembrane (TM) peptides with detectable activity, the computational design of membrane proteins has already produced very exciting results. The 'take-home message' is that optimization of van der Waals packing and hydrogen bonding (both 'canonical' and weak Cα-H⋯O bonds) can produce functional structures of remarkable stability and specificity in the membrane.

Anisotropic solvent model of the lipid bilayer. 2. Energetics of insertion of small molecules, peptides, and proteins in membranes

A new computational approach to calculating binding energies and spatial positions of small molecules, peptides, and proteins in the lipid bilayer has been developed. The method combines an anisotropic solvent representation of the lipid bilayer and universal solvation model, which predicts transfer energies of molecules from water to an arbitrary medium with defined polarity properties. The universal solvation model accounts for hydrophobic, van der Waals, hydrogen-bonding, and electrostatic solute-solvent interactions. The lipid bilayer is represented as a fluid anisotropic environment described by profiles of dielectric constant (ε), solvatochromic dipolarity parameter (π*), and hydrogen bonding acidity and basicity parameters (α and β). The polarity profiles were calculated using published distributions of quasi-molecular segments of lipids determined by neutron and X-ray scattering for DOPC bilayer and spin-labeling data that define concentration of water in the lipid acyl chain region. The model also accounts for the preferential solvation of charges and polar groups by water and includes the effect of the hydrophobic mismatch for transmembrane proteins. The method was tested on calculations of binding energies and preferential positions in membranes for small-molecules, peptides and peripheral membrane proteins that have been experimentally studied. The new theoretical approach was implemented in a new version (2.0) of our PPM program and applied for the large-scale calculations of spatial positions in membranes of more than 1000 peripheral and integral proteins. The results of calculations are deposited in the updated OPM database ( http://opm.phar.umich.edu ).

The integrity of the periplasmic domain of the VirA sensor kinase is critical for optimal coordination of the virulence signal response in Agrobacterium tumefaciens

The plant pathogen Agrobacterium tumefaciens responds to three main signals at the plant-bacterium interface: phenolics, such as acetosyringone (AS), monosaccharides, and acidic pH (∼5.5). These signals are transduced via the chromosomally encoded sugar binding protein ChvE and the Ti plasmid-encoded VirA/VirG two-component regulatory system, resulting in the transcriptional activation of the Ti plasmid virulence genes. Here, we present genetic and physical evidence that the periplasmic domain of VirA dimerizes independently of other parts of the protein, and we examine the effects of several engineered mutations in the periplasmic and transmembrane regions of VirA on vir-inducing capacity as indicated by AS sensitivity and maximal level of vir-inducing activity at saturating AS levels. The data indicate that helix-breaking mutations throughout the periplasmic domain of VirA or mutations that reposition the second transmembrane domain (TM2) of VirA relieve the periplasmic domain's repressive effects on the maximal activity of this kinase in response to phenolics, effects normally relieved only when ChvE, sugars, and low pH are also present. Such relief, however, does not sensitize VirA to low concentrations of phenolics, the other major effect of the ChvE-sugar and low pH signals. We further demonstrate that amino acid residues in a small Trg-like motif in the periplasmic domain of VirA are crucial for transmission of the ChvE-sugar signal to the cytoplasmic domain. These experiments provide evidence that small perturbations in the periplasmic domain of VirA can uncouple sugar-mediated changes in AS sensitivity from the sugar-mediated effects on maximal activity.

NMR analysis of the alphaIIb beta3 cytoplasmic interaction suggests a mechanism for integrin regulation

The integrin αIIbβ3 is a transmembrane (TM) heterodimeric adhesion receptor that exists in equilibrium between resting and active ligand binding conformations. In resting αIIbβ3, the TM and cytoplasmic domains of αIIb and β3 form a heterodimer that constrains αIIbβ3 in its resting conformation. To study the structure and dynamics of the cytoplasmic domain heterodimer, we prepared a disulfide-stabilized complex consisting of portions of the TM domains and the full cytoplasmic domains. NMR and hydrogen-deuterium exchange of this complex in micelles showed that the αIIb cytoplasmic domain is largely disordered, but it interacts with and influences the conformation of the β3 cytoplasmic domain. The β3 cytoplasmic domain consists of a stable proximal helix contiguous with the TM helix and two distal amphiphilic helices. To confirm the NMR structure in a membrane-like environment, we studied the β3 cytoplasmic domain tethered to phospholipid bilayers. Hydrogen-deuterium exchange mass spectrometry, as well as circular dichroism spectroscopy, demonstrated that the β3 cytoplasmic domain becomes more ordered and helical under these conditions, consistent with our NMR results. Further, these experiments suggest that the two distal helices associate with lipid bilayers but undergo fluctuations that would allow rapid binding of cytoplasmic proteins regulating integrin activation, such as talin and kindlin-3. Thus, these results provide a framework for understanding the kinetics and thermodynamics of protein interactions involving integrin cytoplasmic domains and suggest that such interactions act in a concerted fashion to influence integrin stalk separation and exposure of extracellular ligand binding sites.