Search

The cationic lipid 1,2-dioleoyl-3-trimethylammonium propane (DOTAP) is one of the original synthetic cationic lipids used for the liposomal transfection of oligonucleotides in gene therapy. The key structural element of DOTAP is its quaternary ammonium headgroup that is responsible for interactions with both nucleic acids and target cell membranes. Because these interactions are fundamental to the design of a major class of transfection lipids, it is important to understand the structure of DOTAP and how it interacts with halide counterions. Here, we use x-ray and neutron diffraction techniques to examine the structure of DOTAP and how chloride (Cl-) and iodide (I-) counterions alter the hydration properties of the DOTAP headgroup. A problem of particular interest is the poor solubility of DOTAP/I- in water solutions. Our results show that the poor solubility results from very tight binding of the I- counterion to the headgroup and the consequent expulsion of water. The structural principles we report here are important for assessing the suitability of DOTAP and its quaternary ammonium derivatives for transfection., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Biophysical Society. All rights reserved.)

Autism spectrum disorder (ASD) is phenotypically and genetically heterogeneous. We present a CRISPR gene editing strategy to insert a protein tag and premature termination sites creating an induced pluripotent stem cell (iPSC) knockout resource for functional studies of ten ASD-relevant genes (AFF2/FMR2, ANOS1, ASTN2, ATRX, CACNA1C, CHD8, DLGAP2, KCNQ2, SCN2A, TENM1). Neurogenin 2 (NGN2)-directed induction of iPSCs allowed production of excitatory neurons, and mutant proteins were not detectable. RNA sequencing revealed convergence of several neuronal networks. Using both patch-clamp and multi-electrode array approaches, the electrophysiological deficits measured were distinct for different mutations. However, they culminated in a consistent reduction in synaptic activity, including reduced spontaneous excitatory postsynaptic current frequencies in AFF2/FMR2-, ASTN2-, ATRX-, KCNQ2-, and SCN2A-null neurons. Despite ASD susceptibility genes belonging to different gene ontologies, isogenic stem cell resources can reveal common functional phenotypes, such as reduced functional connectivity., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Copy-number variations (CNVs) are strong risk factors for neurodevelopmental and psychiatric disorders. The 15q13.3 microdeletion syndrome region contains up to ten genes and is associated with numerous conditions, including autism spectrum disorder (ASD), epilepsy, schizophrenia, and intellectual disability; however, the mechanisms underlying the pathogenesis of 15q13.3 microdeletion syndrome remain unknown. We combined whole-genome sequencing, human brain gene expression (proteome and transcriptome), and a mouse model with a syntenic heterozygous deletion (Df(h15q13)/+ mice) and determined that the microdeletion results in abnormal development of cortical dendritic spines and dendrite outgrowth. Analysis of large-scale genomic, transcriptomic, and proteomic data identified OTUD7A as a critical gene for brain function. OTUD7A was found to localize to dendritic and spine compartments in cortical neurons, and its reduced levels in Df(h15q13)/+ cortical neurons contributed to the dendritic spine and dendrite outgrowth deficits. Our results reveal OTUD7A as a major regulatory gene for 15q13.3 microdeletion syndrome phenotypes that contribute to the disease mechanism through abnormal cortical neuron morphological development., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

The YidC/Oxa1/Alb3 family of membrane proteins function to insert proteins into membranes in bacteria, mitochondria, and chloroplasts. Recent X-ray structures of YidC from Bacillus halodurans and Escherichia coli revealed a hydrophilic groove that is accessible from the lipid bilayer and the cytoplasm. Here, we explore the water accessibility within the conserved core region of the E. coli YidC using in vivo cysteine alkylation scanning and molecular dynamics (MD) simulations of YidC in POPE/POPG membranes. As expected from the structure, YidC possesses an aqueous membrane cavity localized to the membrane inner leaflet. Both the scanning data and the MD simulations show that the lipid-exposed transmembrane helices 3, 4, and 5 are short, leading to membrane thinning around YidC. Close examination of the MD data reveals previously unrecognized structural features that are likely important for protein stability and function., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

The development of neural connectivity is essential for brain function, and disruption of this process is associated with autism spectrum disorders (ASDs). DIX domain containing 1 (DIXDC1) has previously been implicated in neurodevelopmental disorders, but its role in postnatal brain function remains unknown. Using a knockout mouse model, we determined that DIXDC1 is a regulator of excitatory neuron dendrite development and synapse function in the cortex. We discovered that MARK1, previously linked to ASDs, phosphorylates DIXDC1 to regulate dendrite and spine development through modulation of the cytoskeletal network in an isoform-specific manner. Finally, rare missense variants in DIXDC1 were identified in ASD patient cohorts via genetic sequencing. Interestingly, the variants inhibit DIXDC1 isoform 1 phosphorylation, causing impairment to dendrite and spine growth. These data reveal that DIXDC1 is a regulator of cortical dendrite and synaptic development and provide mechanistic insight into morphological defects associated with neurodevelopmental disorders., (Copyright © 2016 The Author(s). Published by Elsevier Inc. All rights reserved.)

The distribution of peptide conformations in the membrane interface is central to partitioning energetics. Molecular-dynamics simulations enable characterization of in-membrane structural dynamics. Here, we describe melittin partitioning into dioleoylphosphatidylcholine lipids using CHARMM and OPLS force fields. Although the OPLS simulation failed to reproduce experimental results, the CHARMM simulation reported was consistent with experiments. The CHARMM simulation showed melittin to be represented by a narrow distribution of folding states in the membrane interface., (Copyright © 2013 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Lactose permease of Escherichia coli (LacY) catalyzes symport of a galactopyranoside and an H⁺ via an alternating access mechanism. The transition from an inward- to an outward-facing conformation of LacY involves sugar-release followed by deprotonation. Because the transition depends intimately upon the dynamics of LacY in a bilayer environment, molecular dynamics (MD) simulations may be the only means of following the accompanying structural changes in atomic detail. Here, we describe MD simulations of wild-type apo LacY in phosphatidylethanolamine (POPE) lipids that features two protonation states of the critical Glu325. While the protonated system displays configurational stability, deprotonation of Glu325 causes significant structural rearrangements that bring into proximity side chains important for H⁺ translocation and sugar binding and closes the internal cavity. Moreover, protonated LacY in phosphatidylcholine (DMPC) lipids shows that the observed dynamics are lipid-dependent. Together, the simulations describe early dynamics of the inward-to-outward transition of LacY that agree well with experimental data., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

Voltage-dependent ion channels open and close in response to changes in membrane electrical potential due to the motion of their voltage-sensing domains (VSDs). VSD charge displacements within the membrane electric field are observed in electrophysiology experiments as gating currents preceding ionic conduction. The elementary charge motions that give rise to the gating current cannot be observed directly, but appear as discrete current pulses that generate fluctuations in gating current measurements. Here we report direct observation of gating-charge displacements in an atomistic molecular dynamics simulation of the isolated VSD from the KvAP channel in a hydrated lipid bilayer on the timescale (10-μs) expected for elementary gating charge transitions. The results reveal that gating-charge displacements are associated with the water-catalyzed rearrangement of salt bridges between the S4 arginines and a set of conserved acidic side chains on the S1-S3 transmembrane segments in the hydrated interior of the VSD., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

A central feature of the lipid raft concept is the formation of cholesterol-rich lipid domains. The introduction of relatively rigid cholesterol molecules into fluid liquid-disordered (L(d)) phospholipid bilayers can produce liquid-ordered (L(o)) mixtures in which the rigidity of cholesterol causes partial ordering of the flexible hydrocarbon acyl chains of the phospholipids. Several lines of evidence support this concept, but direct structural information about L(o) membranes is lacking. Here we present the structure of L(o) membranes formed from cholesterol and dioleoylphosphatidylcholine (DOPC). Specific deuteration of the DOPC acyl-chain methyl groups and neutron diffraction measurements reveal an extraordinary disorder of the acyl chains of neat L(d) DOPC bilayers. The disorder is so great that >20% of the methyl groups are in intimate contact with water in the bilayer interface. The ordering of the DOPC acyl chains by cholesterol leads to retraction of the methyl groups away from the interface. Molecular dynamics simulations based on experimental systems reveal asymmetric transbilayer distributions of the methyl groups associated with each bilayer leaflet., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The SecY/Sec61 translocon complex, located in the endoplasmic reticulum membrane of eukaryotes (Sec61) or the plasma membrane of prokaryotes (SecY), mediates the transmembrane secretion or insertion of nascent proteins. Mutations that permit the secretion of nascent proteins with defective signal sequences (Prl-phenotype), or interfere with the transmembrane orientation of newly synthesized protein segments, can affect protein topogenesis. The crystallographic structure of SecYEbeta from Methanococcus jannaschii revealed widespread distribution of mutations causing topogenesis defects, but not their molecular mechanisms. Based upon prolonged molecular dynamics simulations of wild-type M. jannaschii SecYEbeta and an extensive sequence-conservation analysis, we show that the closed state of the translocon is stabilized by hydrogen-bonding interactions of numerous highly conserved amino acids. Perturbations induced by mutation at various locations are rapidly relayed to the plug segment that seals the wild-type closed-state translocon, leading to displacement and increased hydration of the plug., (Copyright 2010 Elsevier Ltd. All rights reserved.)

Voltage-sensing domains (VSDs) of voltage-gated potassium (Kv) channels undergo a series of conformational changes upon membrane depolarization, from a down state when the channel is at rest to an up state, all of which lead to the opening of the channel pore. The crystal structures reported to date reveal the pore in an open state and the VSDs in an up state. To gain insights into the structure of the down state, we used a set of experiment-based restraints to generate a model of the down state of the KvAP VSD using molecular-dynamics simulations of the VSD in a lipid bilayer in excess water. The equilibrated VSD configuration is consistent with the biotin-avidin accessibility and internal salt-bridge data used to generate it, and with additional biotin-avidin accessibility data. In the model, both the S3b and S4 segments are displaced approximately 10 A toward the intracellular side with respect to the up-state configuration, but they do not move as a rigid body. Arginine side chains that carry the majority of the gating charge also make large excursions between the up and down states. In both states, arginines interact with water and participate in salt bridges with acidic residues and lipid phosphate groups. An important feature that emerges from the down-state model is that the N-terminal half of the S4 segment adopts a 3(10)-helical conformation, which appears to be necessary to satisfy a complex salt-bridge network., ((c) 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Cell membrane lipids and proteins are heterogeneously distributed in the membrane plane. In recent years, much attention has been paid to the heterogeneous distribution of the lipid components, particularly the formation of cholesterol-rich domains that are thought to be important in signaling processes. This has led to renewed interest in the phase diagrams of complex lipid mixtures, such as three-component mixtures containing phospholipids and cholesterol. We report here a novel fluorescent probe (NBD-R595) that is useful for exploring the phase behaviors of one-, two-, and three-component large unilamellar vesicles. In one-component fluid-phase membranes, the probe has the expected spectral characteristic of monomeric 7-nitrobenzo-2-oxa-1,3-diazol, with a fluorescence maximum of 540 nm when excited at 470 nm. But below the gel-to-liquid crystalline phase transition temperature, an additional emission peak appears at approximately 610 nm, because of Förster resonance energy transfer from NBD-R595 monomers to NBD-R595 Jelley aggregates of limited size formed by the association of 7-nitrobenzo-2-oxa-1,3-diazol moieties. This may be the first report of Förster resonance energy transfer from a single fluorophore in two different physical states. In a test of the probe, we found NBD-R595 to be remarkably sensitive to the molar composition of large unilamellar vesicles formed from cholesterol, distearoylphosphatidylcholine, and dioleoylphosphatidylcholine.

Sphingosine and sphingosine 1-phosphate (S1P) are sphingolipid metabolites that act as signaling messengers to activate or inhibit multiple downstream targets to regulate cell growth, differentiation, and apoptosis. The amphiphilic nature of these compounds leads to aggregation above their critical micelle concentrations (CMCs), which may be important for understanding lysosomal glycosphingolipid storage disorders. We investigated the aggregation of sphingosine and S1P over a comprehensive, physiologically relevant range of pH values, ionic strengths, and lipid concentrations by means of dynamic light scattering, titration, and NMR spectroscopy. The results resolve discrepancies in literature reports of CMC and pK(a) values. At physiological pH, the nominal CMCs of sphingosine and S1P are 0.99 +/- 0.12 microM (pH 7.4) and 14.35 +/- 0.08 microM (pH 7.2), respectively. We find that pH strongly affects the aggregation behavior of sphingosine by changing the ionic and hydrogen-bonding states; the nominal critical aggregation concentrations of protonated and deprotonated sphingosine are 1.71 +/- 0.24 microM and 0.70 +/- 0.02 microM, respectively. NMR measurements revealed that the NH3+-NH2 transition of sphingosine occurs at pH 6.6, and that there is a structural shift in sphingosine aggregates caused by a transition in the predominant hydrogen-bonding network from intramolecular to intermolecular that occurs between pH 6.7 and 9.9.

Intramembrane proteases, which cleave transmembrane (TM) helices, participate in numerous biological processes encompassing all branches of life. Several crystallographic structures of Escherichia coli GlpG rhomboid protease have been determined. In order to understand GlpG dynamics and lipid interactions in a native-like environment, we have examined the molecular dynamics of wild-type and mutant GlpG in different membrane environments. The irregular shape and small hydrophobic thickness of the protein cause significant bilayer deformations that may be important for substrate entry into the active site. Hydrogen-bond interactions with lipids are paramount in protein orientation and dynamics. Mutations in the unusual L1 loop cause changes in protein dynamics and protein orientation that are relayed to the His-Ser catalytic dyad. Similarly, mutations in TM5 change the dynamics and structure of the L1 loop. These results imply that the L1 loop has an important regulatory role in proteolysis.

The innate immune systems of humans and other animals are activated by lipopolysaccharides (LPS), which are glucosamine-based phospholipids that form the outer leaflet of the outer membranes of Gram-negative bacteria. Activation involves interactions of LPS with the innate immunity-receptor comprised of toll-like receptor 4 in complex with so-called MD-2 protein and accessory proteins, such as CD14 and LPS binding protein. The Lipid Metabolites and Pathways Strategy (LIPID MAPS) Consortium has isolated in large amounts a nearly homogeneous LPS, Kdo(2)-Lipid A, and demonstrated that it activates macrophages via toll-like receptor 4. The active form of LPS, monomer or aggregate, is controversial. We have therefore examined the aggregation behavior and other physical properties of Kdo(2)-Lipid A. Differential scanning calorimetry of Kdo(2)-Lipid A suspensions revealed a gel-to-liquid crystalline phase transition at 36.4 degrees C (T(m)). The nominal critical aggregation concentration, determined by dynamic light scattering, was found to be 41.2 +/- 1.6 nM below the T(m) (25 degrees C), but only 8.1 +/- 0.3 nM above the T(m) (37 degrees C). The specific molecular volume of Kdo(2)-Lipid A, obtained by densitometry measurements was found to be 3159 +/- 71 A(3) at 25 degrees C, from which the number of molecules in each aggregate was estimated to be 5.8 x 10(5). The aggregation behavior of Kdo(2)-Lipid A in the presence of lipoprotein-deficient serum suggests that Re LPS monomers and multimers are the active units for the immune system in the CD14-dependent and -independent pathways, respectively.

As a first step toward understanding the principles of the targeting of C2 domains to membranes, we have carried out a molecular dynamics simulation of the C2 domain of cytosolic phospholipase A2 (cPLA2-C2) in a 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer at constant pressure and temperature (NPT, 300 K and 1 atm). Using the high-resolution crystal structure of cPLA2-C2 as a starting point, we embedded two copies of the C2 domain into a pre-equilibrated membrane at the depth and orientation previously defined by electron paramagnetic resonance (EPR). Noting that in the membrane-bound state the three calcium binding loops are complexed to two calcium ions, we initially restrained the calcium ions at the membrane depth determined by EPR. But the depth and orientation of the domains remained within EPR experimental errors when the restraints were later removed. We find that the thermally disordered, chemically heterogeneous interfacial zones of phosphatidylcholine bilayers allow local lipid remodeling to produce a nearly perfect match to the shape and polarity of the C2 domain, thereby enabling the C2 domain to assemble and optimize its own lipid docking site. The result is a cuplike docking site with a hydrophobic bottom and hydrophilic rim. Contrary to expectations, we did not find direct interactions between the protein-bound calcium ions and lipid headgroups, which were sterically excluded from the calcium binding cleft. Rather, the lipid phosphate groups provided outer-sphere calcium coordination through intervening water molecules. These results show that the combined use of high-resolution protein structures, EPR measurements, and molecular dynamics simulations provides a general approach for analyzing the molecular interactions between membrane-docked proteins and lipid bilayers.

Voltage-sensor (VS) domains cause the pore of voltage-gated ion channels to open and close in response to changes in transmembrane potential. Recent experimental studies suggest that VS domains are independent structural units. This independence is revealed dramatically by a voltage-dependent proton-selective channel (Hv), which has a sequence homologous to the VS domains of voltage-gated potassium channels (Kv). Here we show by means of molecular dynamics simulations that the isolated open-state VS domain of the KvAP channel in a lipid membrane has a configuration consistent with a water channel, which we propose as a common feature underlying the conductance of protons, and perhaps other cations, through VS domains.

We have recently shown that current molecular dynamics (MD) atomic force fields are not yet able to produce lipid bilayer structures that agree with experimentally-determined structures within experimental errors. Because of the many advantages offered by experimentally validated simulations, we have developed a novel restraint method for membrane MD simulations that uses experimental diffraction data. The restraints, introduced into the MD force field, act upon specified groups of atoms to restrain their mean positions and widths to values determined experimentally. The method was first tested using a simple liquid argon system, and then applied to a neat dioleoylphosphatidylcholine (DOPC) bilayer at 66% relative humidity and to the same bilayer containing the peptide melittin. Application of experiment-based restraints to the transbilayer double-bond and water distributions of neat DOPC bilayers led to distributions that agreed with the experimental values. Based upon the experimental structure, the restraints improved the simulated structure in some regions while introducing larger differences in others, as might be expected from imperfect force fields. For the DOPC-melittin system, the experimental transbilayer distribution of melittin was used as a restraint. The addition of the peptide caused perturbations of the simulated bilayer structure, but which were larger than observed experimentally. The melittin distribution of the simulation could be fit accurately to a Gaussian with parameters close to the observed ones, indicating that the restraints can be used to produce an ensemble of membrane-bound peptide conformations that are consistent with experiments. Such ensembles pave the way for understanding peptide-bilayer interactions at the atomic level.

A novel protocol has been developed for comparing the structural properties of lipid bilayers determined by simulation with those determined by diffraction experiments, which makes it possible to test critically the ability of molecular dynamics simulations to reproduce experimental data. This model-independent method consists of analyzing data from molecular dynamics bilayer simulations in the same way as experimental data by determining the structure factors of the system and, via Fourier reconstruction, the overall transbilayer scattering-density profiles. Multi-nanosecond molecular dynamics simulations of a dioleoylphosphatidylcholine bilayer at 66% RH (5.4 waters/lipid) were performed in the constant pressure and temperature ensemble using the united-atom GROMACS and the all-atom CHARMM22/27 force fields with the GROMACS and NAMD software packages, respectively. The quality of the simulated bilayer structures was evaluated by comparing simulation with experimental results for bilayer thickness, area/lipid, individual molecular-component distributions, continuous and discrete structure factors, and overall scattering-density profiles. Neither the GROMACS nor the CHARMM22/27 simulations reproduced experimental data within experimental error. The widths of the simulated terminal methyl distributions showed a particularly strong disagreement with the experimentally observed distributions. A comparison of the older CHARMM22 with the newer CHARMM27 force fields shows that significant progress is being made in the development of atomic force fields for describing lipid bilayer systems empirically.

Melittin is arguably the most widely studied amphipathic, membrane-lytic alpha-helical peptide. Although several lines of evidence suggest an interfacial membrane location at low concentrations, melittin's exact position and depth of penetration into the hydrocarbon core are unknown. Furthermore, the structural basis for its lytic action remains largely a matter of conjecture. Using a novel x-ray absolute-scale refinement method, we have now determined the location, orientation, and likely conformation of monomeric melittin in oriented phosphocholine lipid multilayers. Its helical axis is aligned parallel to the bilayer plane at the depth of the glycerol groups, but its average conformation differs from the crystallographic structure. As observed earlier for another amphipathic alpha-helical peptide, the lipid perturbations induced by melittin are remarkably modest. Small bilayer perturbations thus appear to be a general feature of amphipathic helices at low concentrations. In contrast, a dimeric form of melittin causes larger structural perturbations under otherwise identical conditions. These results provide direct structural evidence that self-association of amphipathic helices may be the crucial initial step toward membrane lysis.

Changes in the structure of the hydrocarbon core (HC) of fluid lipid bilayers can reveal how bilayers respond to the partitioning of peptides and other solutes (Jacobs, R. E., and S. H. White. 1989. Biochemistry. 28:3421-3437). The structure of the HC of dioleoylphosphocholine (DOPC) bilayers can be determined from the transbilayer distribution of the double-bonds (Wiener, M. C., and S. H. White. 1992. Biophys. J. 61:434-447). This distribution, representing the time-averaged projection of the double-bond positions onto the bilayer normal (z), can be obtained by means of neutron diffraction and double-bond specific deuteration (Wiener, M. C., G. I. King, and S. H. White. 1991. Biophys. J. 60:568-576). For fully resolved bilayer profiles, a close approximation of the distribution could be obtained by x-ray diffraction and isomorphous bromine labeling at the double-bonds of the DOPC sn-2 acyl chain (Wiener, M. C., and S. H. White. 1991. Biochemistry. 30:6997-7008). We have modified the bromine-labeling approach in a manner that permits determination of the distribution in under-resolved bilayer profiles observed at high water contents. We used this new method to determine the transbilayer distribution of the double-bond bromine labels of DOPC over a hydration range of 5.4 to 16 waters per lipid, which reveals how the HC structure changes with hydration. We found that the transbilayer distributions of the bromines can be described by a pair of Gaussians of 1/e half-width A(Br) located at z = +Z(Br) relative to the bilayer center. For hydrations from 5.4 waters up to 9.4 waters per lipid, Z(Br) decreases from 7.97 +/- 0.27 A to 6.59 +/- 0.15 A, while A(Br) increased from 4.62 +/- 0.62 A to 5.92 +/- 0.37 A, consistent with the expected hydration-induced decrease in HC thickness and increase in area per lipid. After the phosphocholine hydration shell was filled at approximately 12 waters per lipid, we observed a shift in Z(Br) to approximately 7.3 A, indicative of a distinct structural change upon completion of the hydration shell. For hydrations of 12-16 waters per lipid, the bromine distribution remains constant at Z(Br) = 7.33 +/- 0.25 A and A(Br) = 5.35 +/- 0.5 A. The absolute-scale structure factors obtained in the experiments provided an opportunity to test the so-called fluid-minus method of structure-factor scaling. We found that the method is quite satisfactory for determining the phases of structure factors, but not their absolute values.

Many toxins and antimicrobial peptides permeabilize membrane vesicles by forming multimeric pores. Determination of the size of such pores is an important first step for understanding their structure and the mechanism of their self-assembly. We report a simple method for sizing pores in vesicles based on the differential release of co-encapsulated fluorescently labeled dextran markers of two different sizes. The method was tested using the bee venom peptide melittin, which was found to form pores of 25-30 A diameter in palmitoyloleoylphosphatidylcholine (POPC) vesicles at a lipid-to-peptide ratio of 50. This result is consistent with observations on melittin pore formation in erythrocytes (Katsu, T., C. Ninomiya, M. Kuroko, H. Kobayashi, T. Hirota, and Y. Fujita 1988. Action mechanism of amphipathic peptides gramicidin S and melittin on erythrocyte membrane Biochim. Biophys. Acta. 939:57-63).

Tryptophan, proline, and basic amino acids have all been implicated as being important in the assembly and structure of membrane proteins. Indolicidin, an antimicrobial 13-residue peptide-amide isolated from the cytoplasmic granules of bovine neutrophils, is highly enriched in these amino acids: five tryptophans, three prolines, three basic residues, and no acidic residues. Consistent with the likely importance of these amino acids in membrane protein assembly, indolicidin is known to be highly membrane-active and is believed to act by disruption of cell membranes. We have, therefore, examined the interactions of native indolicidin with large unilamellar vesicles (LUV) formed from palmitoyloleoylphosphatidylcholine (POPC), and palmitoyloleoylphosphatidylglycerol (POPG), in order to use it as a model system for studying membrane protein insertion and for evaluating the relative contributions of hydrophobic and electrostatic forces in peptide-bilayer interactions. Equilibrium dialysis measurements indicate that indolicidin binds strongly, but reversibly, to both neutral POPC and anionic POPG vesicles with free energies of transfer of -8.8 +/- 0.2 and -11.5 +/- 0.4 kcal/mol, respectively. The extremely strong partitioning into POPG vesicles necessitated the development of a new equilibrium dialysis method that is described in detail. Tryptophan fluorescence measurements show that indolicidin is located in the bilayer interface and that indole fluorescence is affected by the type of lipid used to form the LUVs. Circular dichroism (CD) measurements reveal unordered conformations in aqueous and bulk organic solutions and a somewhat more ordered, but not alpha-helical, conformation in SDS micelles and lipid bilayers. Fluorescence requenching measurements (Ladokhin et al. 1995. Biophys. J. 69:1964-1971) on vesicles loaded with the fluorophore/quencher pair 8-aminonapthalene-1,3,6 trisulfonic acid (ANTS)/p-xylene-bis-pyridinium bromide (DPX), show that indolicidin induces membrane permeabilization. For anionic POPG, leakage is graded with a high preference for the release of cationic DPX over anionic ANTS. For neutral POPC vesicles no such preference is observed. Leakage induction is more effective with POPG vesicles than with POPC vesicles, as judged by three quantitative measures that are developed in the Appendix.

Agents such as antimicrobial peptides and toxins can permeabilize membrane vesicles to cause leakage of entrapped contents in either a graded or an all-or-none fashion. Determination of which mode of leakage is induced is an important step in understanding the molecular mechanism of membrane permeabilization. Wimley et al. (1994, Protein Sci. 3:1362-1378) have developed a fluorescence method for distinguishing the two modes that makes use of the dye/quencher pair 8-aminonapthalene-1,3,6 trisulfonic acid (ANTS)/p-xylene-bis-pyridinium bromide (DPX) without the usual need for the physical separation of vesicles from released contents. Their "requenching" method establishes the mode of release through the fluorescence changes that occur when DPX is added externally to a solution of vesicles that have released some fraction of their contents. However, the requenching method as originally stated ignored the possibility of preferential release of dye or quencher. Here we extend the theory of the method to take into account preferential release and the effects of graded leakage. The ratio of the rates of release of the cationic quencher DPX and anionic dye 8-aminonapthalene-1,3,6 trisulfonic acid can be estimated by means of the theory. For graded leakage, we show that the release of the markers does not coincide with the fluorescence changes observed in the standard leakage assay. This is true for self-quenching dyes as well and means that 1) the amount of released material will be overestimated and 2) the kinetics will be nonexponential and have artificially high apparent rates. We show how the extended requenching analysis allows the results of leakage experiments to be corrected for artifacts that result from graded and preferential leakage. Experimental evidence is presented for the existence of peptide-induced preferential graded leakage of DPX from both neutral and anionic vesicles.

We present in this paper the complete structure of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in the L alpha phase (66% RH, 23 degrees C) obtained by the joint refinement of neutron and x-ray lamellar diffraction data. The structural details obtained have previously required a large number of neutron diffraction experiments, using numerous specifically-deuterated phospholipid isomorphs (Büldt et al., 1978. Nature (Lond.). 271:182-184). The joint-refinement approach minimizes specific deuteration by utilizing independent neutron and x-ray data sets. The method yields a quasimolecular structure consisting of a series of multiatomic fragments that are each represented by one or several Gaussian distributions whose positions and widths can be determined to within 0.06 to 0.52 A exclusive of the methylene region. The image of DOPC at 66% RH (5.36 +/- 0.08 waters per lipid) is consistent with many aspects of bilayer structure previously determined by structural and spectroscopic studies. The most striking feature of the structure is the large amount of transbilayer thermal motion suggested by the widths and overlaps of the Gaussian envelopes of the quasimolecular fragments. We discuss the "dynamic bilayer thickness" which describes the minimum effective thickness of the hydrocarbon permeability barrier in terms of the thermal motion of the water. A gradient of thermal motion exists that increases in either direction away from the glycerol backbone which is the most constrained portion of the bilayer. The steric interactions between headgroups of apposed bilayers, expected at the hydration level of our experiments, are clearly revealed. A useful consequence of the quasimolecular structure is that average boundaries within bilayers calculated using composition and volumetric data and ad hoc assumptions can be related to the positions of the principal structural groups. Several measures of "bilayer thickness" in common use can be identified as the positions of the cholines for Luzzati's d1 (Luzzati and Husson. 1962. J. Cell Biol. 12:207-219) and the glycerols for Small's dL (Small. 1967. J. Lipid Res. 8:551-556). We do not know if these relations will be true at other hydrations or for other lipids. Of particular interest is the fact that the position of the carbonyl groups marks the average hydrocarbon/headgroup boundary. It must be emphasized, however, that this region of the bilayer must be generally characterized as one of tumultuous chemical heterogeneity because of the thermal motion of the bilayer.

We continue in this paper the presentation of theoretical and experimental methods for the joint refinement of neutron and x-ray lamellar diffraction data for the analysis of fluid (L alpha phase) bilayer structure (Wiener, M. C., and S. H. White. 1991 a, b, c. Biophys. J. 59:162-173 and 174-185; Biochemistry. 30:6997-7008; Wiener, M. C., G. I. King, and S. H. White. Biophys. J. 60: 568-576). We show how to obtain the distribution and packing of the terminal methyls in the interior of a fluid dioleoylphosphatidylcholine bilayer (66% RH) by combining x-ray and neutron scattering-length transbilayer profiles with no a priori assumptions about the functional form of the distribution. We find that the methyls can be represented by a Gaussian function with 1/e-halfwidth of 2.95 +/- 0.28 A situated at the bilayer center. There is substantial mixing of the methyls and methylenes in the bilayer center. The Gaussian representation of the methyl distribution is narrower and has a different shape than predicted by several simulations of fluid bilayers (Gruen, D. W. R., and E. H. B. de Lacey. 1984. Surfactants in Solution, Vol. 1. Plenum Publishing Corp., New York. 279-306; de Loof, H., et al. 1991. Biochemistry. 30:2099-2133) but this may be due to the smaller area/lipid of our experiments and the presence of the double-bonds. Determination of the absolute specific volume of DOPC and an analysis of bulk alkane volumetric data over a range of hydrostatic pressures lead to estimates of methylene and methyl volumes at the bilayer center of 27 +/- 1 A3 and 57.2 +/- 3.6 A3, respectively. This result provides direct confirmation of the common assumption that the molecular packing of methyl and methylene groups in bilayers is the same as in bulk liquid alkanes.

We described in two previous papers a method for the joint refinement of the structure of fluid bilayers using neutron and x-ray diffraction data (Wiener, M. C., and S. H. White 1991a, b. Biophys. J. 59: 162-173 and 174-185). An essential part of the method is the appropriate scaling of the diffraction data. Here we describe the scaling of the neutron data and the determination of the transbilayer distribution of double bonds in liquid-crystalline (L alpha phase) phospholipid bilayers of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC). The distribution was determined by neutron diffraction of oriented multilayers (66% RH) of DOPC specifically deuterated at the 9- and 10-position of both acyl chains. The double-bond distribution is described accurately by a pair of Gaussian functions each located at a position Zcc = 7.88 +/- 0.09 A from the bilayer center with 1/e-halfwidths of Acc = 4.29 +/- 0.16 A. Previously, we determined the transbilayer distribution of bromine atoms in a specifically halogenated lipid, 1-oleoyl-2-9,10-dibromostearoyl-sn-glycero-3-phosphocholine (OBPC), and showed it to be an isomorphous replacement for DOPC (Wiener, M. C., and S. H. White, 1991c. Biochemistry. In press). A comparison of the double-bond and bromine profiles indicates that the positions of the centers of the deuterated double bond and the brominated methylene Gaussian distributions are equal within experimental error and that each label undergoes similar average thermal motions with respect to the bilayer normal. The observation that the average position of a label on both acyl chains (the deuterated double bonds) is similar to the average position of a label on the 2-chain alone (the brominated methylenes) indicates that the maximum separation along the bilayer normal between the double bonds of the acyl chains is 1 A or less. The fully-resolved transbilayer water distribution, previously determined at lower resolution (Jacobs, R. E., and S. H. White. 1989. Biochemistry. 28:3421-3437), was obtained from the analysis of neutron diffraction data of DOPC hydrated with a D20/H20 mixture. The water distribution is described accurately by a pair of Gaussian functions each located at a position Zw = 22.51 +/- 0.77 A from the bilayer center with 1/e-half widths of Aw = 4.63 +/- 0.48A. We present the relative absolute neutron and x-ray structure factors of DOPC at 66% RH that will be used to solve the complete structure of DOPC which will be presented in a later paper of this series.

This is the first in a series of papers concerned with methods for the determination of the structures of fluid phospholipid bilayers in the liquid-crystalline (L alpha) phase. The basic approach is the joint refinement of quasimolecular models (King and White, 1986. Biophys. J. 49:1047-1054) using x-ray and neutron diffraction data. We present here (a) the rationale for quasimolecular models, (b) the nature of the resolution problem for thermally disordered bilayers, and (c) an analysis of the resolution of experiments in which Gaussian functions are used to describe the distribution of submolecular components. We show that multilamellar liquid-crystalline bilayers are best described by the convolution of a perfect lattice function with a thermally disordered bilayer unit cell. Lamellar diffraction measurements on such a system generally yield only 5-10 orders of diffraction data from which transbilayer profiles of the unit cell can be constructed. The canonical resolution of these transbilayer profiles, defined as the Bragg spacing divided by the index of the highest recorded diffraction order, is typically 5-10 A. Using simple model calculations, we show that the canonical resolution is a measure of the widths of the distributions of constituents of the unit cell rather than a measure of the spatial separation of the distributions. The widths provide a measure of the thermal motion of the bilayer constituents which can be described by Gaussian functions. The equilibrium positions of the centers of the distributions can be determined with a precision of 0.1-0.5 A based upon typical experimental errors.

This is the second of two papers describing a method for the joint refinement of the structure of fluid bilayers using x-ray and neutron diffraction data. We showed in the first paper (Wiener, M. C., and S. H. White. 1990. Biophys. J. 59:162-173) that fluid bilayers generally consist of a nearly perfect lattice of thermally disordered unit cells and that the canonical resolution d/hmax is a measure of the widths of quasimolecular components represented by simple Gaussian functions. The thermal disorder makes possible a "composition space" representation in which the quasimolecular Gaussian distributions describe the number or probability of occupancy per unit length across the width of the bilayer of each component. This representation permits the joint refinement of neutron and x-ray lamellar diffraction data by means of a single quasimolecular structure that is fit simultaneously to both diffraction data sets. Scaling of each component by the appropriate neutron or x-ray scattering length maps the composition space profile to the appropriate scattering length space for comparison to experimental data. Other extensive properties, such as mass, can also be obtained by an appropriate scaling of the refined composition space structure. Based upon simple bilayer models involving crystal and liquid crystal structural information, we estimate that a fluid bilayer with hmax observed diffraction orders will be accurately represented by a structure with approximately hmax quasimolecular components. Strategies for assignment of quasimolecular components are demonstrated through detailed parsing of a phospholipid molecule based upon the one-dimensional projection of the crystal structure of dimyristoylphosphatidylcholine. Finally, we discuss in detail the number of experimental variables required for the composition space joint refinement. We find fluid bilayer structures to be marginally determined by the experimental data. The analysis of errors, which takes on particular importance under these circumstances, is also discussed.

We consider in this paper the statistical distribution of hydrophobic residues along the length of protein chains. For this purpose we used a binary hydrophobicity scale which assigns hydrophobic residues a value of one and non-hydrophobes a value of zero. The resulting binary sequences are tested for randomness using the standard run test. For the majority of the 5,247 proteins examined, the distribution of hydrophobic residues along a sequence cannot be distinguished from that expected for a random distribution. This suggests that (a) functional proteins may have originated from random sequences, (b) the folding of proteins into compact structures may be much more permissive with less sequence specificity than previously thought, and (c) the clusters of hydrophobic residues along chains which are revealed by hydrophobicity plots are a natural consequence of a random distribution and can be conveniently described by binomial statistics.

Stable myogenic cell lines have been derived at a high frequency by transfection of a cloned multipotential mouse embryo cell line, C3H 10T1/2, with cloned human DNA linked to a selectable neomycin resistance gene. The myogenic phenotype remains linked to neomycin resistance during secondary transfections. Although proliferative in growth conditions, these cell lines maintain the ability to differentiate and express muscle-specific proteins. We conclude that there is a simple genetic basis for myogenic determination and that a single gene, myd, converts 10T1/2 cells to a myoblast lineage. Southern blot analysis demonstrates nonidentity of myd and the MyoD1 gene. Northern blot analysis shows that myd-transfected myogenic lineages express MyoD1 mRNA while parental 10T1/2 cells do not. These results suggest that a dependent regulatory gene pathway mediates myogenic determination and differentiation.

The formation of planar bilayer membranes from lipid monolayers as described by Montal and Mueller (Proc. Natl. Acad. Sci. 1972. 69:3561) is analyzed. Bilayers absolutely free of alkane solvents or other nonpolar hydrocarbons can be formed on polytetrafluoroethylene (PTFE) (e.g. Teflon) septa only if certain boundary conditions are satisfied. Measurements have been made of the contact angles between monolayer-coated water and PTFE in the presence and absence of alkane solvents. The measurement suggest that the boundary conditions for formation of stable bilayers can be satisfied only when a nonpolar solvent is present. We conclude that the bilayer must be surrounded by a torus of alkane solvent, petroleum jelly, or silicone grease depending upon the details of technique used to form the bilayer. The non-polar solvent used in the formation of the bilayer may or may not be present in the bilayer depending upon the water solubility and size of the solvent molecule relative to the size of the alkyl chain of the lipid. Detailed sketches describing the formation of bilayers from monolayers are presented.

A simple technique for forming "black" lipid bilayer membranes containing negligible amounts of alkyl solvent is described. The membranes are formed by the method of Mueller et al (Circulation. 1962. 26:1167.) from glyceryl monooleate (GMO) dispersed in squalene. The squalene forms an annulus to satisfy the boundary conditions of the planar bilayer but does not appear to dissolve noticeably in the bilayer itself. The specific geometric capacitance (Cg) of the membranes at 20 degrees C formed by this technique is 0.7771 +/- 0.0048 muF/cm2. Theoretical estimates of Cg for solvent-free bilayers range from 0.75 to 0.81 muF/cm2. Alkane-free GMO bilayers formed from n-octadecane by the solvent freeze-out method of White (Biochim. Biophys. Acta. 1974. 356:8) have values of Cg = 0.7903 +/- 0.0013 muF/cm2 at 20.5 degrees C. The agreement between the various values of Cg strongly suggests that the bilayers are free of squalene. DC potentials applied to the bilayers have no detectable effect on the value of Cg, as expected for solvent-free films. The ability to form bilayers essentially free of the solvent used in the forming solution makes it possible to determine the area per molecule of the surface active lipid in the bilayer. The area per molecule of GMO at 20 degrees C is estimated to be 37.9 +/- 0.2 A2.

Neutron diffraction methods provide information about the distribution of matter in biological and model membrane systems. The information is derived from plots (profiles) of scattering length density along an axis normal to the membrane plane. Without the use of specific deuteration, the generally low resolution of the profiles limits their interpretation in terms of specific chemical constituents (e.g., lipid headgroup, lipid hydrocarbon, protein, and water). A fundamental and useful structural assignment to make is the boundary between the headgroup and hydrocarbon regions of bilayers. We demonstrate here that strip-function model representations of neutron scattering length density profiles of bilayers are sufficient to determine accurately the position of the headgroup-hydrocarbon boundary. The resulting hydrocarbon thickness of the bilayer is useful for determining the area per lipid molecule and consequently the molecular packing arrangements of the membrane constituents. We analyze data obtained from dioleoylphosphatidylcholine (DOPC) bilayers at 66% RH using standard Fourier profile analyses and from DOPC deuterated specifically at the C-2 carbon of the acyl chains using difference Fourier analysis. We demonstrate that strip-function models accurately define the positions of the C-2 carbons and thus the hydrocarbon thickness (dhc) of the bilayer. We then show, using quasi-molecular models, that the strip-model analysis probably provides an accurate measure of dhc because of the exceptionally high scattering length density difference between the carbonyl and methylene groups.

We describe in this paper measurements of the partial specific volumes of lipid (upsilon L) and water (upsilon W) in mechanical mixtures of egg yolk lecithin and water over a range of hydrations (5-55 wt% water) that includes the region at which excess water appears. upsilon L and upsilon W are found to be very nearly 1 cm3/g over the entire range. The water activities of the mixtures were also determined and found to be the same as for lipid deposited as oriented multilayers on solid substrates.

Temperature-dependent structural changes in planar bilayer membranes formed from glycerol monooleate (GMO) dispersed in various n-alkane solvents (C12-C17) have been studied using precise measurements of specific geometric capacitance (Cg). Cg generally increases as temperature (T) decreases. A change in the slope of Cg(T) occurs between 15 and 18 degrees C for all solvent systems examined. Measurements of the interfacial tension (gamma) of the bulk GMO-alkane dispersions against 0.1 M NaCl show that gamma generally decreases with decreasing temperature. The data can be fitted with two straight lines of different slope which intersect on the average at 17 degrees C. Pagano et al. (1973, Science (Wash. D.C.). 181:557) have shown using calorimetry that GMO has a phase transition at about 15 degrees C. Thus, the changes in Cg and gamma with temperature are likely to result from a GMO phase transition. A second structural change is observed to occur between 5 and 10 degrees C which has not been detected calorimetrically. Calculations of Cg based on various estimates of the hydrocarbon dielectric coefficient (epsilon-b) and/or hydrocarbon thickness (delta-b) leads to models for the structure of the bilayer above and below the phase transition temperature.

A method is described for measuring the specific capacitance (C(m)) of lipid bilayer membranes with an estimated experimental error of only 1%. The gross capacitance was measured with an AC Wheatstone bridge and a photographic technique was used to determine the area of thin membrane. The results of measurements on oxidized cholesterol-decane membranes formed in 1 x 10(-2) M KCl show that C(m) depends upon temperature, voltage, time, and the age of the bulk membrane solutions. For a freshly thinned membrane (from 5 week old solution), C(m) increases exponentially from an initial value of 0.432 +/-0.021 (SD) muF/cm(2) with a time constant of approximately 15 min. A 100 mv potential applied across the membrane for 10-20 min prior to making measurements eliminated this time dependence and produced final-state membranes. C(m) of final-state membranes depends upon applied voltage (V(a)) and obeys the equation C(m) = C(0) + betaV(a) (2) where V(a) approximately V(DC) + V(rms) (AC). C(0) and beta depend upon temperature; C(0) decreases linearly with temperature while beta increases linearly. At 20 degrees C, C(0) = 0.559 +/-0.01 (SD) muF/cm(2) and beta = 0.0123 +/-0.0036 (SD) (muF/cm(2))/(mv(2)) and at 34 degrees C, C(0) = 0.472 +/-0.01 and beta = 0.0382 +/-0.0039. These variations in C(m) are interpreted as resulting from thickness changes. The possibility that they result from diffuse layer and/or membrane dielectric phenomena is discussed and found to be unlikely. The results are discussed in terms of membrane stability by constructing hypothetical potential energy vs. thickness curves.

The characteristics and behavior of the torus (annulus) surrounding planar lipid bilayer membranes formed across a cylindrical aperture are analyzed using equations for the shape and volume of the annulus derived by the methods of variational calculus. The analysis leads to the following results: (a) Design criteria for the aperture can be established. (b) The transition region between thin film and thick annulus can be defined quantitatively and its effect on the measurement of specific capacitance determined. (c) At fixed annulus volume the diameter of the thin membrane is a function of the thin film-annulus contact angle. This suggests a new method for examining changes in free energy of the thin film, and explains why the area of thin film increases reversibly when potentials are present across the film. (d) In the absence of buoyant forces, the equations for the shape and volume of the annulus consist of incomplete elliptic integrals of the first and second kinds; however, the shape of the annulus in the transition region can be described with good accuracy by an approximate equation of greater simplicity.