Cardiolipin (CL) is a negatively charged four acyl chain lipid, associated with energy production in bacterial and mitochondrial membranes. Due to the shape of CL, negative curvatures of aggregates are favorable if the charges in the head group can be reduced. The phase polymorphism of CL, and of associated derivatives with 2, 3, 4, or 5 chains, has been determined previously and offers a model system in which micellar, lamellar, and inverse hexagonal phases can be observed. We present an extension to a previously established coarse-grained molecular dynamics model with the aim of reproducing the different CL phases with two adjustable parameters: the number of acyl chains and the effective head group charge. With molecular dynamics simulations of large lipid systems, we observed transitions between different phases on the nanosecond to microsecond time scale. Charge screening by high salt or low pH was successfully modeled by a reduction of phosphate charge, which led to the adoption of aggregates with more negative curvature. Although specific ion binding at the interface and other atomistic details are sacrificed in the coarse-grained model, we found that it captures general phase features over a large range of aggregate geometries.

Cardiolipin is a key lipid component in the inner mitochondrial membrane, where the lipid is involved in energy production, cristae structure, and mechanisms in the apoptotic pathway. In this article we used molecular dynamics computer simulations to investigate cardiolipin and its effect on the structure of lipid bilayers. Three cardiolipin/POPC bilayers with different lipid compositions were simulated: 100, 9.2, and 0% cardiolipin. We found strong association of sodium counterions to the carbonyl groups of both lipid types, leaving in the case of 9.2% cardiolipin virtually no ions in the aqueous compartment. Although binding occurred primarily at the carbonyl position, there was a preference to bind to the carbonyl groups of cardiolipin. Ion binding and the small headgroup of cardiolipin gave a strong ordering of the hydrocarbon chains. We found significant effects in the water dipole orientation and water dipole potential which can compensate for the electrostatic repulsion that otherwise should force charged lipids apart. Several parameters relevant for the molecular structure of cardiolipin were calculated and compared with results from analyses of coarse-grained simulations and available X-ray structural data.

Lipid shape and charge are connected with the physical properties and the biological function of membranes. Cardiolipin, a double phospholipid with four chains and the potential of changing its charge with pH, is crucially connected with mitochondrial inner membrane shape, and recent experiments suggest that local pH changes allow highly curved local geometries. Here, we use a coarse-grained molecular dynamics model to investigate the mechanical properties of cardiolipin bilayers, systematically varying the headgroup charge and the composition in mixtures with zwitterionic 1,2-dioleoyl-glycero-3-phosphatidylcholine (DOPC) or 1,2-dioleoyl-glycero-3-phosphatidylethanolamine (DOPE). Low cardiolipin charge, corresponding to low pH, was found to induce bending moduli on the order of kBT and curved microdomains. On the length scale investigated, in contrast to continuum theoretical models, we found the area modulus and bending modulus to be inversely correlated for mixtures of cardiolipin and DOPC/DOPE, explainable by changes in the effective headgroup volume.

Lipid membranes composed of 1,2-di-(9Z,12Z,15Z)-octadecatrienoyl-3-O-β-D-galactosyl-sn-glycerol or monogalactosyldiacylglycerol (MGDG) and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) were studied by means of molecular dynamics (MD) computer simulations. Three lipid compositions were considered: 0%, 20%, and 45% MGDG (by mole) denoted as MG-0, MG-20, and MG-45, respectively. The article is focused on the calculation of NMR dipolar interactions, which were confronted with previously reported experimental couplings. Dynamical processes and orientational distributions relevant for the averaging of dipolar interactions were evaluated. Furthermore, several parameters important for characterization of the bilayer structure, molecular organization, and dynamics were investigated. In general, only a minor change in DMPC properties was observed upon the increased MGDG/DMPC ratio, whereas properties related to MGDG undergo a more pronounced change. This effect was ascribed to the fact that DMPC is a bilayer (L(α)) forming lipid, whereas MGDG prefers a reverse hexagonal (H(II)) arrangement.

Cardiolipin is a key lipid component in many biological membranes. Proton conduction and proton-lipid interactions on the membrane surface are thought to be central to mitochondrial energy production. However, details on the cardiolipin headgroup structure are lacking and the protonation state of this lipid at physiological pH is not fully established. Here we present ab initio DFT calculations of the cardiolipin (CL) headgroup and its 2'-deoxy derivative (dCL), with the aim of establishing a connection between structure and acid-base equilibrium in CL. Furthermore, we investigate the effects of solvation on the molecular conformations. In our model, both CL and dCL showed a significant gap between the two pK(a) values, with pK(a2) above the physiological range, and intramolecular hydrogen bonds were found to play a central role in the conformations of both molecules. This behavior was also observed experimentally in CL. Structures derived from the DFT calculations were compared with those obtained experimentally, collected for CL in the Protein Data Bank, and conformations from previous as well as new molecular dynamics simulations of cardiolipin bilayers. Transition states for proton transfer in CL were investigated, and we estimate that protons can exchange between phosphate groups with an approximate 4-5 kcal/mol barrier. Computed NMR and IR spectral properties were found to be in reasonable agreement with experimental results available in the literature.