Structure of fully hydrated fluid phase DMPC and DLPC lipid bilayers using X-ray scattering from oriented multilamellar arrays and from unilamellar vesicles

Quantitative structures of the fully hydrated fluid phases of dimyristoylphosphatidylcholine (DMPC) and dilauroylphosphatidylcholine (DLPC) were obtained at 30 degrees C. Data for the relative form factors F(q(z)) for DMPC were obtained using a combination of four methods. 1), Volumetric data provided F(0). 2), Diffuse x-ray scattering from oriented stacks of bilayers provided relative form factors |F(q(z))| for high q(z), 0.22 < q(z) < 0.8 A(-1). 3), X-ray scattering from extruded unilamellar vesicles with diameter 600 A provided |F(q(z))| for low q(z), 0.1 < q(z) < 0.3 A(-1). 4), Previous measurements using a liquid crystallographic x-ray method provided |F(2 pi h/D)| for h = 1 and 2 for a range of nearly fully hydrated D-spacings. The data from method 4 overlap and validate the new unilamellar vesicles data for DMPC, so method 4 is not required for DLPC or future studies. We used hybrid electron density models to obtain structural results from these form factors. Comparison of the model electron density profiles with that of gel phase DMPC provides areas per lipid A, 60.6 +/- 0.5 A(2) for DMPC and 63.2 +/- 0.5 A(2) for DLPC. Constraints on the model provided by volume measurements and component volumes obtained from simulations put the electron density profiles rho(z) and the corresponding form factors F(q(z)) on absolute scales. Various thicknesses, such as the hydrophobic thickness and the steric thickness, are obtained and compared to literature values.     

Experimental validation of molecular dynamics simulations of lipid bilayers: a new approach

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.     

Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains

Quantitative structures are obtained at 30 degrees C for the fully hydrated fluid phases of palmitoyloleoylphosphatidylcholine (POPC), with a double bond on the sn-2 hydrocarbon chain, and for dierucoylphosphatidylcholine (di22:1PC), with a double bond on each hydrocarbon chain. The form factors F(qz) for both lipids are obtained using a combination of three methods. (1) Volumetric measurements provide F(0). (2) X-ray scattering from extruded unilamellar vesicles provides /F(qz)/ for low q(z). (3) Diffuse X-ray scattering from oriented stacks of bilayers provides /F(qz)/ for high q(z). Also, data using method (2) are added to our recent data for dioleoylphosphatidylcholine (DOPC) using methods (1) and (3); the new DOPC data agree very well with the recent data and with (4) our older data obtained using a liquid crystallographic X-ray method. We used hybrid electron density models to obtain structural results from these form factors. The result for area per lipid (A) for DOPC 72.4 +/- 0.5 A(2) agrees well with our earlier publications, and we find A = 69.3 +/- 0.5 A2 for di22:1PC and A = 68.3 +/- 1.5 A2 for POPC. We obtain the values for five different average thicknesses: hydrophobic, steric, head-head, phosphate-phosphate and Luzzati. Comparison of the results for these three lipids and for our recent dimyristoylphosphatidylcholine (DMPC) determination provides quantitative measures of the effect of unsaturation on bilayer structure. Our results suggest that lipids with one monounsaturated chain have quantitative bilayer structures closer to lipids with two monounsaturated chains than to lipids with two completely saturated chains.     

Structure of magainin and alamethicin in model membranes studied by x-ray reflectivity

We have investigated the structure of lipid bilayers containing varied molar ratios of different lipids and the antimicrobial peptides magainin and alamethicin. For this structural study, we have used x-ray reflectivity on highly aligned solid-supported multilamellar lipid membranes. The reflectivity curves have been analyzed by semi-kinematical reflectivity theory modeling the bilayer density profile rho(z). Model simulations of the reflectivity curves cover a large range of vertical momentum transfer q(z), and yield excellent agreement between data and theory. The structural changes observed as a function of the molar peptide/lipid concentration P/L are discussed in a comparative way.     

Reparameterization of all-atom dipalmitoylphosphatidylcholine lipid parameters enables simulation of fluid bilayers at zero tension

Molecular dynamics simulations of dipalmitoylphosphatidylcholine (DPPC) lipid bilayers using the CHARMM27 force field in the tensionless isothermal-isobaric (NPT) ensemble give highly ordered, gel-like bilayers with an area per lipid of approximately 48 A(2). To obtain fluid (L(alpha)) phase properties of DPPC bilayers represented by the CHARMM energy function in this ensemble, we reparameterized the atomic partial charges in the lipid headgroup and upper parts of the acyl chains. The new charges were determined from the electron structure using both the Mulliken method and the restricted electrostatic potential fitting method. We tested the derived charges in molecular dynamics simulations of a fully hydrated DPPC bilayer. Only the simulation with the new restricted electrostatic potential charges shows significant improvements compared with simulations using the original CHARMM27 force field resulting in an area per lipid of 60.4 +/- 0.1 A(2). Compared to the 48 A(2), the new value of 60.4 A(2) is in fair agreement with the experimental value of 64 A(2). In addition, the simulated order parameter profile and electron density profile are in satisfactory agreement with experimental data. Thus, the biologically more interesting fluid phase of DPPC bilayers can now be simulated in all-atom simulations in the NPT ensemble by employing our modified CHARMM27 force field.     

Closer look at structure of fully hydrated fluid phase DPPC bilayers

X-ray data are presented for the benchmark dipalmitoylphosphatidylcholine lipid bilayer in the most biologically relevant state in which the bilayers are fully hydrated and in the fluid (liquid-crystalline) phase. Form factors F(q(z)) are obtained from a combination of two sample preparations, oriented stacks of bilayers for q(z) extending to 0.85 A(-1) and unilamellar vesicles for smaller q(z). Modeling obtains the electron density profile and values for the area per molecule, for the locations of the component groups, and for the different types of thicknesses of the bilayer, such as the hydrocarbon thickness and the steric thickness.     

Effects of the eukaryotic pore-forming cytolysin Equinatoxin II on lipid membranes and the role of sphingomyelin

Equinatoxin II (EqtII), a protein toxin from the sea anemone Actinia equina, readily creates pores in sphingomyelin-containing lipid membranes. The perturbation by EqtII of model lipid membranes composed of dimyristoylphosphatidycholine and sphingomyelin (10 mol %) was investigated using wideline phosphorus-31 and deuterium NMR. The preferential interaction between EqtII (0.1 and 0.4 mol %) and the individual bilayer lipids was studied by (31)P magic angle spinning NMR, and toxin-induced changes in bilayer morphology were examined by freeze-fracture electron microscopy. Both NMR and EM showed the formation of an additional lipid phase in sphingomyelin-containing mixed lipid multilamellar suspensions with 0.4 mol % EqtII. The new toxin-induced phase consisted of small unilamellar vesicles 20-40 nm in diameter. Deuterium NMR showed that the new lipid phase contains both dimyristoylphosphatidycholine and sphingomyelin. Solid-state (31)P NMR showed an increase in spin-lattice and a decrease in spin-spin relaxation times in mixed-lipid model membranes in the presence of EqtII, consistent with an increase in the intensity of low frequency motions. The (2)H and (31)P spectral intensity distributions confirmed a change in lipid mobility and showed the creation of an isotropic lipid phase, which was identified as the small vesicle structures visible by electron microscopy in the EqtII-lipid suspensions. The toxin appears to enhance slow motions in the membrane lipids and destabilize the membrane. This effect was greatly enhanced in sphingomyelin-containing mixed lipid membranes compared with pure phosphatidylcholine bilayers, suggesting a preferential interaction between the toxin and bilayer sphingomyelin.     

Molecular dynamics simulations of palmitoyloleoylphosphatidylglycerol bilayers

Phosphatidylglycerol (PG) is one of the important components of biological membranes, but there is a paucity of experimental data to test the accuracy of molecular dynamics (MD) simulations. This work consists of testing the accuracy of the CHARMM36 (C36) lipid force field on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG) lipid bilayers. MD simulations of POPG lipid bilayers are compared to recently available X-ray and neutron scattering and deuterium NMR measurements. Overall, the C36 lipid force field accurately represents the X-ray and neutron form factors, bilayer and hydrocarbon thicknesses and chain deuterium order parameters. The surface area per lipid from MD simulations with C36 (67.7 ± 0.2 Å²) is in excellent agreement with the experimentally determined value of 66.0 ± 1.3 Å². C36 outperforms the lipid force field developed by Berger et al. [15] and suggests that past studies with this force field may result in lateral areas that are too small. Moreover, our studies give some insight into the structural model used in experiments and suggest that the functional form for the head group may not be Gaussian-like. Based on our simulations, the POPG lipid in the C36 lipid force field is well parameterised and can be used for other PG lipids and membrane models with mixed lipids.     

Validation of all-atom phosphatidylcholine lipid force fields in the tensionless NPT ensemble

A recently defined charge set, to be used in conjunction with the all-atom CHARMM27r force field, has been validated for a series of phosphatidylcholine lipids. The work of Sonne et al. successfully replicated experimental bulk membrane behaviour for dipalmitoylphosphatidylcholine (DPPC) under the isothermal-isobaric (NPT) ensemble. Previous studies using the defined CHARMM27r charge set have resulted in lateral membrane contraction when used in the tensionless NPT ensemble, forcing the lipids to adopt a more ordered conformation than predicted experimentally. The current study has extended the newly defined charge set to 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC) and 1-palmitoyl-2-docosahexaenoyl-sn-glycero-3-phosphatidylcholine (PDPC). Molecular dynamics simulations were run for each of the lipids (including DPPC) using both the CHARMM27r charge set and the newly defined modified charge set. In all three cases a significant improvement was seen in both bulk membrane properties and individual atomistic effects. Membrane width, area per lipid and the depth of water penetration were all seen to converge to experimental values. Deuterium order parameters generated with the new charge set showed increased disorder across the width of the bilayer and reflected both results from experiment and similar simulations run with united atom models. These newly validated models can now find use in mixed biological simulations under the tensionless ensemble without concern for lateral contraction.     

Radiation hybrid comparative mapping between human chromosome 17 and porcine chromosome 12 demonstrates conservation of gene order

A comparative study of human chromosome 17 (HSA17) and pig chromosome 12 (SSC12) was conducted using both somatic cell hybrid panel (SCHP) and radiation hybrid (RH) panel analysis. Sequences from an expressed sequence tag (EST) project in pig reproduction were examined and six genes and ESTs originally believed to map to HSA17 were selected for this study. The genes/ESTs were TATA box binding protein-associated factor (TAF2N/RBP56), alpha-2-plasmin inhibitor (SERPINF2/PLI), H3 histone family 3B (H3F3B), aminopeptidase puromycin sensitive (NPEPPS), an expressed sequence tag (ESTMI015) and P311 protein (P311). The SCHP analysis mapped five genes/ESTs (TAF2N, H3F3B, SERPINF2, NPEPPS and ESTMI015) to SSC12q11-q15 and SSC12p11-p15 with 100% concordance, and assigned P311 to SSC2 (1/2q24)-q29 with 100% concordance. Radiation hybrid analysis of all six genes confirmed the SCHP mapping results, with average retention frequency of 25%. Recent human sequence data demonstrated that P311 is actually located on HSA5q. As HSA5q and SSC2q show conserved syntenic regions predicted from bi-directional painting, our P311 mapping data is consistent with these results. An expanded comparative SSC12 RH map integrating the five new type I markers and 23 previously mapped loci was established using a LOD score threshold of 4.8. The gene order of the five genes/ESTs on the SSC12 framework RH map (H3F3B-ESTMI015-NPEPPS-TAF2N-SERPINF2) is identical to the HSA17 GB4 map but with inversion of the map as conventionally drawn.     

Effects of MHC and gender on lupus-like autoimmunity in Nba2 congenic mice

The lupus-like disease that develops in hybrids of NZB and NZW mice is genetically complex, involving both MHC- and non-MHC-encoded genes. Studies in this model have indicated that the H2d/z MHC type, compared with H2d/d or H2z/z, is critical for disease development. C57BL/6 (B6) mice (H2b/b) congenic for NZB autoimmunity 2 (Nba2), a NZB-derived susceptibility locus on distal chromosome 1, produce autoantibodies to nuclear Ags, but do not develop kidney disease. Crossing B6.Nba2 to NZW results in H2b/z F1 offspring that develop severe lupus nephritis. Despite the importance of H2z in past studies, we found no enhancement of autoantibody production or nephritis in H2b/z vs H2b/b B6.Nba2 mice, and inheritance of H2z/z markedly suppressed autoantibody production. (B6.Nba2 x NZW)F1 mice, compared with MHC-matched B6.Nba2 mice, produced higher levels of IgG autoantibodies to chromatin, but not to dsDNA. Although progressive renal damage with proteinuria only occurred in F1 mice, kidneys of some B6.Nba2 mice showed similar extensive IgG and C3 deposition. We also studied male and female B6.Nba2 and F1 mice with different MHC combinations to determine whether increased susceptibility to lupus among females was also expressed within the context of the Nba2 locus. Regardless of MHC or the presence of NZW genes, females produced higher levels of antinuclear autoantibodies, and female F1 mice developed severe proteinuria with higher frequencies. Together, these studies help to clarify particular genetic and sex-specific influences on the pathogenesis of lupus nephritis.     

Effects of force fields on the conformational and dynamic properties of amyloid β(1‐40) dimer explored by replica exchange molecular dynamics simulations

The conformational space and structural ensembles of amyloid beta (Aβ) peptides and their oligomers in solution are inherently disordered and proven to be challenging to study. Optimum force field selection for molecular dynamics (MD) simulations and the biophysical relevance of results are still unknown. We compared the conformational space of the Aβ(1‐40) dimers by 300 ns replica exchange MD simulations at physiological temperature (310 K) using: the AMBER‐ff99sb‐ILDN, AMBER‐ff99sb*‐ILDN, AMBER‐ff99sb‐NMR, and CHARMM22* force fields. Statistical comparisons of simulation results to experimental data and previously published simulations utilizing the CHARMM22* and CHARMM36 force fields were performed. All force fields yield sampled ensembles of conformations with collision cross sectional areas for the dimer that are statistically significantly larger than experimental results. All force fields, with the exception of AMBER‐ff99sb‐ILDN (8.8 ± 6.4%) and CHARMM36 (2.7 ± 4.2%), tend to overestimate the α‐helical content compared to experimental CD (5.3 ± 5.2%). Using the AMBER‐ff99sb‐NMR force field resulted in the greatest degree of variance (41.3 ± 12.9%). Except for the AMBER‐ff99sb‐NMR force field, the others tended to under estimate the expected amount of β‐sheet and over estimate the amount of turn/bend/random coil conformations. All force fields, with the exception AMBER‐ff99sb‐NMR, reproduce a theoretically expected β‐sheet‐turn‐β‐sheet conformational motif, however, only the CHARMM22* and CHARMM36 force fields yield results compatible with collapse of the central and C‐terminal hydrophobic cores from residues 17‐21 and 30‐36. Although analyses of essential subspace sampling showed only minor variations between force fields, secondary structures of lowest energy conformers are different.     

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.     

Identification of a putative gamma-aminobutyric acid (GABA) receptor subunit rho2 cDNA and colocalization of the genes encoding rho2 (GABRR2) and rho1 (GABRR1) to human chromosome 6q14-q21 and mouse chromosome 4.

Screening of a genomic DNA library with a portion of the cDNA encoding the gamma-aminobutyric acid (GABA) receptor subunit rho1 identified two distinct clones. DNA sequencing revealed that one clone contained a single exon from the rho1 gene (GABBR1) while the second clone encompassed an exon with 96% identity to the rho1 gene. Screening of a human retina cDNA library with oligonucleotides specific for the exon in the second clone identified a 3-kb cDNA with an open reading frame of 1395 bp. The predicted amino acid sequence of this cDNA demonstrates 30 to 38% similarity to alpha, beta, gamma, and delta GABA receptor subunits and 74% similarity to the GABA rho1 subunit suggesting that the newly isolated cDNA encodes a new member of the rho subunit family, tentatively named GABA rho2. Polymerase chain reaction (PCR) amplification of rho1 and rho2 gene sequences from DNA of three somatic cell hybrid panels maps both genes to human chromosome 6, bands q14 to q21. Tight linkage was also demonstrated between restriction fragment length variants (RFLVs) from each rho gene and the Tsha locus on mouse chromosome 4, which is homologous to the CGA locus on human chromosome 6q12-q21. These two lines of evidence confirm that GABRR1 and newly identified GABRR2 map to the same region on human chromosome 6. This close physical association and high degree of sequence similarity raises the possibility that one rho gene arose from the other by duplication.     

Structure of the fully hydrated gel phase of dipalmitoylphosphatidylcholine.

X-ray diffraction intensities for lamellar repeats, h = 1 to 7, and wide-angle x-ray scattering were measured for the gel phase of fully hydrated dipalmitoylphosphatidylcholine. A hybrid model, which represents the electron density along the lamellar repeat direction as a continuous function composed of constant strips and superimposed Gaussians, is defined. The data were used to determine the best parameters in hybrid models and also in the older strip models. The most successful results were obtained when the density of the methylene region was constrained to the value obtained from the wide-angle scattering. Further analysis utilized the lipid volume obtained from absolute specific volume measurements. Together with the fundamental relations derived in the previous paper, the electron density modeling yielded the headgroup volume (340 +/- 10 A3) and the methylene volume (25.3 +/- 0.2A3). The results were in agreement whether the hybrid model or the strip model was used and whether our data or the data of Torbet and Wilkins were used. Additional structural results, such as the area (45.9 +/- 2.0 A2) and the number of waters of hydration per lipid (10.6 +/- 2.0), required one additional piece of information, which we took to be the tilt angle theta, which is 30 +/- 3 degrees from other experiments in the literature. Absolute electron density profiles, which clearly indicate two features in the headgroup region, are presented. The analysis yielded an accurate value of F(0), which contributes to the continuous scattering transform F(X), which is also given.     

The assembly of ionic currents in a thalamic neuron. I. The three-dimensional model

We have previously discussed qualitative models for bursting and thalamic neurons that were obtained by modifying a simple two-dimensional model for repetitive firing. In this paper we report the results of making a similar sequence of modifications to a more elaborate six-dimensional model of repetitive firing which is based on the Hodgkin-Huxley equations. To do this we first reduce the six-dimensional model to a two-dimensional model that resembles our original two-dimensional qualitative model. This is achieved by defining a new variable, which we call q. We then add a subthreshold inward current and a subthreshold outward current having a variable, z, that changes slowly. This gives a three-dimensional (v,q,z) model of the Hodgkin-Huxley type, which we refer to as the z-model. Depending on the choice of parameter values this model resembles our previous models of bursting and thalamic neurons. At each stage in the development of these models we return to the corresponding seven-dimensional model to confirm that we can obtain similar solutions by using the complete system of equations. The analysis of the three-dimensional model involves a state diagram and a stability diagram. The state diagram shows the projection of the phase path from v,q,z space into the v,z plane, together with the projections of the curves z = 0 and v = q = 0. The stability of the points on the curve v = q = 0, which we call the v, q nullcurve, is determined by the stability diagram. Taken together the state and stability diagrams show how to assemble the ionic currents to produce a given firing pattern.     

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari–Szabo squared generalized order parameter (O²) values of amide N? H bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O² values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari–Szabo generalized order parameters fall below an R² of 0.8.     

Molecular Dynamics and NMR Spectroscopy Studies of E.?coli Lipopolysaccharide Structure and Dynamics

Lipopolysaccharide (LPS), a component of Gram-negative bacterial outer membranes, comprises three regions: lipid A, core oligosaccharide, and O-antigen polysaccharide. Using the CHARMM36 lipid and carbohydrate force fields, we have constructed a model of an Escherichia coli R1 (core) O6 (antigen) LPS molecule. Several all-atom bilayers are built and simulated with lipid A only (LIPA) and varying lengths of 0 (LPS0), 5 (LPS5), and 10 (LPS10) O6 antigen repeating units; a single unit of O6 antigen contains five sugar residues. From ¹H,¹H-NOESY experiments, cross-relaxation rates are obtained from an O-antigen polysaccharide sample. Although some experimental deviations are due to spin-diffusion, the remaining effective proton-proton distances show generally very good agreement between NMR experiments and molecular dynamics simulations. The simulation results show that increasing the LPS molecular length has an impact on LPS structure and dynamics and also on LPS bilayer properties. Terminal residues in a LPS bilayer are more flexible and extended along the membrane normal. As the core and O-antigen are added, per-lipid area increases and lipid bilayer order decreases. In addition, results from mixed LPS0/5 and LPS0/10 bilayer simulations show that the LPS O-antigen conformations at a higher concentration of LPS5 and LPS10 are more orthogonal to the membrane and less flexible. The O-antigen concentration of mixed LPS bilayers does not have a significant effect on per-lipid area and hydrophobic thickness. Analysis of ion and water penetration shows that water molecules can penetrate inside the inner core region, and hydration is critical to maintain the integrity of the bilayer structure.     

Application of a SCC-DFTB QM/MM approach to the investigation of the catalytic mechanism of fatty acid amide hydrolase

Self-consistent charge density functional tight binding (SCC-DFTB) is a promising method for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of enzyme-catalyzed reactions. The acylation reaction of fatty acid amide hydrolase (FAAH), a promising drug target, was investigated by applying a SCC-DFTB/CHARMM27 scheme. Calculated potential energy barriers resulted in reasonable agreement with experiments for oleamide (OA) and oleoylmethyl ester (OME) substrates, outperforming previous calculations performed at the PM3/CHARMM22 level. Furthermore, the experimental preference of FAAH in hydrolyzing OA faster than OME was adequately reproduced by calculations. All these findings indicate that the SCC-DFTB/CHARMM27 approach can be successfully applied to mechanistic investigations of FAAH-catalyzed reactions.     

Lipid Configurations from Molecular Dynamics Simulations

The extent to which current force fields faithfully reproduce conformational properties of lipids in bilayer membranes, and whether these reflect the structural principles established for phospholipids in bilayer crystals, are central to biomembrane simulations. We determine the distribution of dihedral angles in palmitoyl-oleoyl phosphatidylcholine from molecular dynamics simulations of hydrated fluid bilayer membranes. We compare results from the widely used lipid force field of Berger et al. with those from the most recent C36 release of the CHARMM force field for lipids. Only the CHARMM force field produces the chain inequivalence with sn-1 as leading chain that is characteristic of glycerolipid packing in fluid bilayers. The exposure and high partial charge of the backbone carbonyls in Berger lipids leads to artifactual binding of Na⁺ ions reported in the literature. Both force fields predict coupled, near-symmetrical distributions of headgroup dihedral angles, which is compatible with models of interconverting mirror-image conformations used originally to interpret NMR order parameters. The Berger force field produces rotamer populations that correspond to the headgroup conformation found in a phosphatidylcholine lipid bilayer crystal, whereas CHARMM36 rotamer populations are closer to the more relaxed crystal conformations of phosphatidylethanolamine and glycerophosphocholine. CHARMM36 alone predicts the correct relative signs of the time-average headgroup order parameters, and reasonably reproduces the full range of NMR data from the phosphate diester to the choline methyls. There is strong motivation to seek further experimental criteria for verifying predicted conformational distributions in the choline headgroup, including the ³¹P chemical shift anisotropy and ¹⁴N and CD₃ NMR quadrupole splittings.