The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.

Although an X-ray crystal structure of lactose permease (LacY) has been presented with bound galactopyranoside, neither the sugar nor the residues ligating the sugar can be identified with precision at ~3.5 Å. Therefore, additional evidence is important for identifying side chains likely to be involved in binding. On the basis of a clue from site-directed alkylation suggesting that Asn272, Gly268, and Val264 on one face of helix VIII might participate in galactoside binding, molecular dynamics simulations were conducted initially. The simulations indicate that Asn272 (helix VIII) is sufficiently close to the galactopyranosyl ring of a docked lactose analogue to play an important role in binding, the backbone at Gly268 may be involved, and Val264 does not interact with the bound sugar. When the three side chains are subjected to site-directed mutagenesis, with the sole exception of mutant Asn272 → Gln, various other replacements for Asn272 either markedly decrease affinity for the substrate (i.e., high KD) or abolish binding altogether. However, mutant Gly268 → Ala exhibits a moderate 8-fold decrease in affinity, and binding by mutant Val264 → Ala is affected only minimally. Thus, Asn272 and possibly Gly268 may comprise additional components of the galactoside-binding site in LacY.

Camelid nanobodies (Nbs) raised against the outward-facing conformer of a double-Trp mutant of the lactose permease of Escherichia coli (LacY) stabilize the permease in outward-facing conformations. Isothermal titration calorimetry is applied herein to dissect the binding thermodynamics of two Nbs, one that markedly improves access to the sugar-binding site and another that dramatically increases the affinity for galactoside. The findings presented here show that both enthalpy and entropy contribute favorably to binding of the Nbs to wild-type (WT) LacY and that binding of Nb to double-Trp mutant G46W/G262W is driven by a greater enthalpy at an entropic penalty. Thermodynamic analyses support the interpretation that WT LacY is stabilized in outward-facing conformations like the double-Trp mutant with closure of the water-filled cytoplasmic cavity through conformational selection. The LacY conformational transition required for ligand binding is reflected by a favorable entropy increase. Molecular dynamics simulations further suggest that the entropy increase likely stems from release of immobilized water molecules primarily from the cytoplasmic cavity upon closure.

On the periplasmic side of LacY, two conserved Gly-Gly pairs in helices II and XI (Gly46 and Gly370, respectively) and helices V and VIII (Gly159 and Gly262, respectively) allow close packing of each helix pair in the outward (periplasmic)-closed conformation. Previous studies demonstrate that replacing one Gly residue in each Gly-Gly pair with Trp leads to opening of the periplasmic cavity with abrogation of transport activity, but an increased rate of galactoside binding. To further investigate the role of the Gly-Gly pairs, 11 double-replacement mutants were constructed for each pair at positions 46 (helix II) and 262 (helix VIII). Replacement with Ala or Ser results in decreased but significant transport activity, while replacements with Thr, Val, Leu, Asn, Gln, Tyr, Trp, Glu, or Lys exhibit very little or no transport. Remarkably, however, the double mutants bind galactoside with affinities 10-20-fold higher than that of the pseudo-WT or WT LacY. Moreover, site-directed alkylation of a periplasmic Cys replacement indicates that the periplasmic cavity becomes readily accessible in the double-replacement mutants. Molecular dynamics simulations with the WT and double-Leu mutant in the inward-open/outward-closed conformation provide support for this interpretation.

Heavy metals in cells are typically regulated by PIB-type ATPases. The first structure of the class, a Cu(+)-ATPase from Legionella pneumophila (LpCopA), outlined a copper transport pathway across the membrane, which was inferred to be occluded. Here we show by molecular dynamics simulations that extracellular water solvated the transmembrane (TM) domain, results indicative of a Cu(+)-release pathway. Furthermore, a new LpCopA crystal structure determined at 2.8-Å resolution, trapped in the preceding E2P state, delineated the same passage, and site-directed-mutagenesis activity assays support a functional role for the conduit. The structural similarities between the TM domains of the two conformations suggest that Cu(+)-ATPases couple dephosphorylation and ion extrusion differently than do the well-characterized PII-type ATPases. The ion pathway explains why certain Menkes' and Wilson's disease mutations impair protein function and points to a site for inhibitors targeting pathogens.

ClC-1 protein channels facilitate rapid passage of chloride ions across cellular membranes, thereby orchestrating skeletal muscle excitability. Malfunction of ClC-1 is associated with myotonia congenita, a disease impairing muscle relaxation. Here, we present the cryo-electron microscopy (cryo-EM) structure of human ClC-1, uncovering an architecture reminiscent of that of bovine ClC-K and CLC transporters. The chloride conducting pathway exhibits distinct features, including a central glutamate residue ("fast gate") known to confer voltage-dependence (a mechanistic feature not present in ClC-K), linked to a somewhat rearranged central tyrosine and a narrower aperture of the pore toward the extracellular vestibule. These characteristics agree with the lower chloride flux of ClC-1 compared with ClC-K and enable us to propose a model for chloride passage in voltage-dependent CLC channels. Comparison of structures derived from protein studied in different experimental conditions supports the notion that pH and adenine nucleotides regulate ClC-1 through interactions between the so-called cystathionine-β-synthase (CBS) domains and the intracellular vestibule ("slow gating"). The structure also provides a framework for analysis of mutations causing myotonia congenita and reveals a striking correlation between mutated residues and the phenotypic effect on voltage gating, opening avenues for rational design of therapies against ClC-1-related diseases.