Membrane transporter dimerization driven by differential lipid solvation energetics of dissociated and associated states.

Over two-thirds of integral membrane proteins of known structure assemble into oligomers. Yet, the forces that drive the association of these proteins remain to be delineated, as the lipid bilayer is a solvent environment that is both structurally and chemically complex. In this study, we reveal how the lipid solvent defines the dimerization equilibrium of the CLC-ec1 Cl-/H+ antiporter. Integrating experimental and computational approaches, we show that monomers associate to avoid a thinned-membrane defect formed by hydrophobic mismatch at their exposed dimerization interfaces. In this defect, lipids are strongly tilted and less densely packed than in the bulk, with a larger degree of entanglement between opposing leaflets and greater water penetration into the bilayer interior. Dimerization restores the membrane to a near-native state and therefore, appears to be driven by the larger free-energy cost of lipid solvation of the dissociated protomers. Supporting this theory, we demonstrate that addition of short-chain lipids strongly shifts the dimerization equilibrium toward the monomeric state, and show that the cause of this effect is that these lipids preferentially solvate the defect. Importantly, we show that this shift requires only minimal quantities of short-chain lipids, with no measurable impact on either the macroscopic physical state of the membrane or the protein's biological function. Based on these observations, we posit that free-energy differentials for local lipid solvation define membrane-protein association equilibria. With this, we argue that preferential lipid solvation is a plausible cellular mechanism for lipid regulation of oligomerization processes, as it can occur at low concentrations and does not require global changes in membrane properties.

A cell's outer membrane is made of molecules called lipids, which band together to form a flexible thin film, just two molecules thick. This membrane is dotted with proteins that transport materials in to and out of cells. Most of these membrane proteins join with other proteins to form structures known as oligomers. Except, how membrane-bound proteins assemble into oligomers ­ the physical forces driving these molecules to take shape ­ remains unclear. This is partly because the structural, physical and chemical properties of fat-like lipid membranes are radically different to the cell's watery interior. Consequently, the conditions under which membrane oligomers form are distinct from those surrounding proteins inside cells. Membrane proteins are also more difficult to study and characterize than water-soluble proteins inside the cell, and yet many therapeutic drugs such as antibiotics specifically target membrane proteins. Overall, our understanding of how the unique properties of lipid membranes affect the formation of protein structures embedded within, is lacking and warrants further investigation. Now, Chadda, Bernhardt et al. focused on one membrane protein, known as CLC, which tends to exist in pairs ­ or dimers. To understand why these proteins form dimers (a process called dimerization) Chadda, Bernhardt et al. first used computer simulations, and then validated the findings in experimental tests. These complementary approaches demonstrated that the main reason CLC proteins 'dimerize' lies in their interaction with the lipid membrane, and not the attraction of one protein to its partner. When CLC proteins are on their own, they deform the surrounding membrane and create structural defects that put the membrane under strain. But when two CLC proteins join as a dimer, this membrane strain disappears ­ making dimerization the more stable and energetically favorable option. Chadda, Bernhardt et al. also showed that with the addition of a few certain lipids, specifically smaller lipids, cell membranes become more tolerant of protein-induced structural changes. This might explain how cells could use various lipids to fine-tune the activity of membrane proteins by controlling how oligomers form. However, the theory needs to be examined further. Altogether, this work has provided fundamental insights into the physical forces shaping membrane-bound proteins, relevant to researchers studying cell biology and pharmacology alike.

MALTA: a calculator for estimating the coverage with shRNA, CRISPR, and cDNA libraries.

Genetic screens using shRNA, CRISPR, or cDNA libraries rely on adequately transferring the library into cells for further assay. These libraries can have many different elements and each element can be present at different copy numbers within a given pooled library. Calculating how many recipient cells are needed to adequately sample all or most of the different elements within a library is important, especially if one wants to compare the outcomes of different genetic screens that rely on accurately reproducing the starting population of library-containing cells. Here we present a simple application that starts with a list of library elements and their abundance and calculates the minimum sampling number to achieve full transfer of the library to an acceptor cell population to a user-specified level of probability. Users can adjust several input parameters including designating a subpopulation over which the calculation is made. Finally, the program performs a series of Monte Carlo simulations of a user-specified number of picks to produce an empirically determined distribution of each library element.

Informatic Analysis of Sequence Data from Batch Yeast 2-Hybrid Screens.

We have adapted the yeast 2-hybrid assay to simultaneously uncover dozens of transient and static protein interactions within a single screen utilizing high-throughput short-read DNA sequencing. The resulting sequence datasets can not only track what genes in a population that are enriched during selection for positive yeast 2-hybrid interactions, but also give detailed information about the relevant subdomains of proteins sufficient for interaction. Here, we describe a full suite of stand-alone software programs that allow non-experts to perform all the bioinformatics and statistical steps to process and analyze DNA sequence fastq files from a batch yeast 2-hybrid assay. The processing steps covered by these software include: 1) mapping and counting sequence reads corresponding to each candidate protein encoded within a yeast 2-hybrid prey library; 2) a statistical analysis program that evaluates the enrichment profiles; and 3) tools to examine the translational frame and position within the coding region of each enriched plasmid that encodes the interacting proteins of interest.

Cellular encoding of Cy dyes for single-molecule imaging.

A general method is described for the site-specific genetic encoding of cyanine dyes as non-canonical amino acids (Cy-ncAAs) into proteins. The approach relies on an improved technique for nonsense suppression with in vitro misacylated orthogonal tRNA. The data show that Cy-ncAAs (based on Cy3 and Cy5) are tolerated by the eukaryotic ribosome in cell-free and whole-cell environments and can be incorporated into soluble and membrane proteins. In the context of the Xenopus laevis oocyte expression system, this technique yields ion channels with encoded Cy-ncAAs that are trafficked to the plasma membrane where they display robust function and distinct fluorescent signals as detected by TIRF microscopy. This is the first demonstration of an encoded cyanine dye as a ncAA in a eukaryotic expression system and opens the door for the analysis of proteins with single-molecule resolution in a cellular environment.

DEEPN as an Approach for Batch Processing of Yeast 2-Hybrid Interactions.

We adapted the yeast 2-hybrid assay to simultaneously uncover multiple transient protein interactions within a single screen by using a strategy termed DEEPN (dynamic enrichment for evaluation of protein networks). This approach incorporates high-throughput DNA sequencing and computation to follow competition among a plasmid population encoding interacting partners. To demonstrate the capacity of DEEPN, we identify a wide range of ubiquitin-binding proteins, including interactors that we verify biochemically. To demonstrate the specificity of DEEPN, we show that DEEPN allows simultaneous comparison of candidate interactors across multiple bait proteins, allowing differential interactions to be identified. This feature was used to identify interactors that distinguish between GTP- and GDP-bound conformations of Rab5.

The dimerization equilibrium of a ClC Cl(-)/H(+) antiporter in lipid bilayers.

Interactions between membrane protein interfaces in lipid bilayers play an important role in membrane protein folding but quantification of the strength of these interactions has been challenging. Studying dimerization of ClC-type transporters offers a new approach to the problem, as individual subunits adopt a stable and functionally verifiable fold that constrains the system to two states - monomer or dimer. Here, we use single-molecule photobleaching analysis to measure the probability of ClC-ec1 subunit capture into liposomes during extrusion of large, multilamellar membranes. The capture statistics describe a monomer to dimer transition that is dependent on the subunit/lipid mole fraction density and follows an equilibrium dimerization isotherm. This allows for the measurement of the free energy of ClC-ec1 dimerization in lipid bilayers, revealing that it is one of the strongest membrane protein complexes measured so far, and introduces it as new type of dimerization model to investigate the physical forces that drive membrane protein association in membranes.

Optimization of an elastic network augmented coarse grained model to study CCMV capsid deformation.

The major protective coat of most viruses is a highly symmetric protein capsid that forms spontaneously from many copies of identical proteins. Structural and mechanical properties of such capsids, as well as their self-assembly process, have been studied experimentally and theoretically, including modeling efforts by computer simulations on various scales. Atomistic models include specific details of local protein binding but are limited in system size and accessible time, while coarse grained (CG) models do get access to longer time and length scales but often lack the specific local interactions. Multi-scale models aim at bridging this gap by systematically connecting different levels of resolution. Here, a CG model for CCMV (Cowpea Chlorotic Mottle Virus), a virus with an icosahedral shell of 180 identical protein monomers, is developed, where parameters are derived from atomistic simulations of capsid protein dimers in aqueous solution. In particular, a new method is introduced to combine the MARTINI CG model with a supportive elastic network based on structural fluctuations of individual monomers. In the parametrization process, both network connectivity and strength are optimized. This elastic-network optimized CG model, which solely relies on atomistic data of small units (dimers), is able to correctly predict inter-protein conformational flexibility and properties of larger capsid fragments of 20 and more subunits. Furthermore, it is shown that this CG model reproduces experimental (Atomic Force Microscopy) indentation measurements of the entire viral capsid. Thus it is shown that one obvious goal for hierarchical modeling, namely predicting mechanical properties of larger protein complexes from models that are carefully parametrized on elastic properties of smaller units, is achievable.

Molecular dynamics simulation of the unfolding of individual bacteriorhodopsin helices in sodium dodecyl sulfate micelles.

We report molecular dynamics simulations of the trends in the changes in secondary structure of the seven individual helices of bacteriorhodopsin when inserted into sodium dodecyl sulfate (SDS) micelles, and their dependence on the amino acid sequence. The results indicate that the partitioning of the helices in the micelles and their stability are dependent on the hydrophobicity of the transmembrane segments. Helices A, B, and E are stable and retain their initial secondary structure throughout the 100 ns simulation time. In contrast, helices C, D, F, and G show structural perturbations within the first 10 ns. The instabilities are localized near charged residues within the transmembrane segments. The overall structural instability of the helix is correlated with its partitioning to the surface of the micelle and its interaction with polar groups there. The in silico experiments were performed to complement the in vitro experiments that examined the partial denaturation of bacteriorhodopsin in SDS described in the preceding article (DOI 10.1021/bi201769z ). The simulations are consistent with the trends revealed by the experimental results but strongly underestimate the extent of helix to extended coil transformation. The reason may be either that the sampling time was not sufficiently long or, more interestingly, that interhelix residue interactions play a role in the unfolding of the helices.

Secondary and tertiary structure of bacteriorhodopsin in the SDS denatured state.

We characterized the structure of partially unfolded bacteriorhodopsin in sodium dodecyl sulfate (SDS) micelles and compared it with its in vitro refolded structure after reconstitution with dimyristoylphosphatidylcholine/3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (DMPC/CHAPS). Intrahelical and interhelical distances were mapped in the protein using strategically located spin-label pairs at helical ends, assayed by pulsed electron paramagnetic resonance spectroscopy (double electron-electron spin resonance, DEER). We find that in SDS the intrahelical end-to-end distances exhibit broad distributions, suggesting a heterogeneous ensemble of conformations with differing secondary structures. Nevertheless, a majority of the denatured population retains end-to-end distances similar to those in the native state. In contrast, the observed greatly increased interhelical distances, in addition to their very broad distributions, suggest that in the SDS micelles very little of the native tertiary structure remains.

Structural changes in bacteriorhodopsin during in vitro refolding from a partially denatured state.

We report on the formation of the secondary and tertiary structure of bacteriorhodopsin during its in vitro refolding from an SDS-denatured state. We used the mobility of single spin labels in seven samples, attached at various locations to six of the seven helical segments to engineered cysteine residues, to follow coil-to-helix formation. Distance measurements obtained by spin dipolar quenching in six samples labeled at either the cytoplasmic or extracellular ends of pairs of helices revealed the time dependence of the recovery of the transmembrane helical bundle. The secondary structure in the majority of the helical segments refolds with a time constant of <100-140 ms. Recovery of the tertiary structure is achieved by sequential association of the helices and occurs in at least three distinct steps with time constants of 1), well below 1 s; 2), 3-4 s; and 3), 60-130 s (the latter depending on the helical pair). The slowest of these processes occurs in concert with recovery of the retinal chromophore.