The mitochondria-targeted peptide SS-31 binds lipid bilayers and modulates surface electrostatics as a key component of its mechanism of action.

Mitochondrial dysfunction underlies many heritable diseases, acquired pathologies, and aging-related declines in health. Szeto-Schiller (SS) peptides comprise a class of amphipathic tetrapeptides that are efficacious toward a wide array of mitochondrial disorders and are believed to target mitochondrial membranes because they are enriched in the anionic phospholipid cardiolipin (CL). However, little is known regarding how SS peptides interact with or alter the physical properties of lipid bilayers. In this study, using biophysical and computational approaches, we have analyzed the interactions of the lead compound SS-31 (elamipretide) with model and mitochondrial membranes. Our results show that this polybasic peptide partitions into the membrane interfacial region with an affinity and a lipid binding density that are directly related to surface charge. We found that SS-31 binding does not destabilize lamellar bilayers even at the highest binding concentrations; however, it did cause saturable alterations in lipid packing. Most notably, SS-31 modulated the surface electrostatics of both model and mitochondrial membranes. We propose nonexclusive mechanisms by which the tuning of surface charge could underpin the mitoprotective properties of SS-31, including alteration of the distribution of ions and basic proteins at the interface, and/or modulation of bilayer physical properties. As a proof of concept, we show that SS-31 alters divalent cation (calcium) distribution within the interfacial region and reduces the energetic burden of calcium stress in mitochondria. The mechanistic details of SS-31 revealed in this study will help inform the development of future compound variants with enhanced efficacy and bioavailability.

Molecular Dynamics Analysis of Cardiolipin and Monolysocardiolipin on Bilayer Properties.

The mitochondrial lipid cardiolipin (CL) contributes to the spatial protein organization and morphological character of the inner mitochondrial membrane. Monolysocardiolipin (MLCL), an intermediate species in the CL remodeling pathway, is enriched in the multisystem disease Barth syndrome. Despite the medical relevance of MLCL, a detailed molecular description that elucidates the structural and dynamic differences between CL and MLCL has not been conducted. To this end, we performed comparative atomistic molecular dynamics studies on bilayers consisting of pure CL or MLCL to elucidate similarities and differences in their molecular and bulk bilayer properties. We describe differential headgroup dynamics and hydrogen bonding patterns between the CL variants and show an increased cohesiveness of MLCL's solvent interfacial region, which may have implications for protein interactions. Finally, using the coarse-grained Martini model, we show that substitution of MLCL for CL in bilayers mimicking mitochondrial composition induces drastic differences in bilayer mechanical properties and curvature-dependent partitioning behavior. Together, the results of this work reveal differences between CL and MLCL at the molecular and mesoscopic levels that may underpin the pathomechanisms of defects in cardiolipin remodeling.

Understanding Conformational Dynamics of Complex Lipid Mixtures Relevant to Biology.

This is a perspective article entitled "Frontiers in computational biophysics: understanding conformational dynamics of complex lipid mixtures relevant to biology" which is following a CECAM meeting with the same name.

Buckling Under Pressure: Curvature-Based Lipid Segregation and Stability Modulation in Cardiolipin-Containing Bilayers.

Mitochondrial metabolic function is affected by the morphology and protein organization of the mitochondrial inner membrane. Cardiolipin (CL) is a unique tetra-acyl lipid that is involved in the maintenance of the highly curved shape of the mitochondrial inner membrane as well as spatial organization of the proteins necessary for respiration and oxidative phosphorylation. Cardiolipin has been suggested to self-organize into lipid domains due to its inverted conical molecular geometry, though the driving forces for this organization are not fully understood. In this work, we use coarse-grained molecular dynamics simulations to study the mechanical properties and lipid dynamics in heterogeneous bilayers both with and without CL, as a function of membrane curvature. We find that incorporation of CL increases bilayer deformability and that CL becomes highly enriched in regions of high negative curvature. We further show that another mitochondrial inverted conical lipid, phosphatidylethanolamine (PE), does not partition or increase the deformability of the membrane in a significant manner. Therefore, CL appears to possess some unique characteristics that cannot be inferred simply from molecular geometry considerations.

Curvature sensing lipid dynamics in a mitochondrial inner membrane model.

Membrane curvature is essential for many cellular structures and processes, and factors such as leaflet asymmetry, lipid composition, and proteins all play important roles. Cardiolipin is the signature lipid of mitochondrial membranes and is essential for maintaining the highly curved shapes of the inner mitochondrial membrane (IMM) and the spatial arrangement of membrane proteins. In this study, we investigate the partitioning behavior of various lipids present in the IMM using coarse-grained molecular dynamics simulations. This study explores curved bilayer systems containing phosphatidylcholine (PC), phosphatidylethanolamine (PE), and cardiolipin (CDL) in binary and ternary component mixtures. Curvature properties such as mean and Gaussian curvatures, as well as the distribution of lipids into the various curved regions of the cristae models, are quantified. Overall, this work represents an advance beyond previous studies on lipid curvature sensing by simulating these systems in a geometry that has the morphological features and scales of curvature consistent with regions of the IMM. We find that CDL has a stronger preference for accumulating in regions of negative curvature than PE lipids, in agreement with previous results. Furthermore, we find lipid partitioning propensity is dominated by sensitivity to mean curvature, while there is a weaker correlation with Gaussian curvature.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus.

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 µs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ1), a construct where the 8 C-terminal residues of γ are fused to γ1 (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

A computational analysis of the interaction of lattice and intramolecular vibrational modes in crystalline alpha-RDX.

The vibrational spectrum of a computer model of crystalline RDX was studied using a 216-molecule periodic supercell, allowing for intra- and intermolecular degrees of freedom using the force field by Boyd et al. [J. Chem. Phys. 124, 104508 (2006)]. The normal modes were analyzed with regard to their activity involving molecule center-of-mass translations and rotations, as well as 15 intramolecular degrees of freedom, including bond stretches, bend and dihedral angle variations, and out-of-plane motions of the nitro groups. We correlate center-of-mass motions with the occupation of internal degrees of freedom for all of the normal modes in the model with particular attention to correlations between nitro rotations and lattice modes. Transfer of lattice energy to internal degrees of freedom can occur through doorway modes and is significant for the initiation of detonation. Several clusters of potential doorway modes are found which involve significant lattice motion as well as nitro rotations. Such groups of modes have been found in the ranges between 186 and 220 and between 420 and 434 cm(-1). Symmetry properties and details of the involved molecular motions are described.

BUMPy: A Model-Independent Tool for Constructing Lipid Bilayers of Varying Curvature and Composition.

Molecular dynamics is a powerful tool to investigate atomistic and mesoscopic phenomena in lipid bilayer systems. These studies have progressed with the advent of increased computational power, and efforts are now increasingly being directed toward investigating the role of curvature and bilayer morphology, as these are critical features of biological processes. Computational studies of lipid bilayers benefit from tools that can create starting configurations for molecular dynamics simulations, but the majority of such tools are restricted to generating flat bilayers. Generating curved bilayer configurations comes with practical complications and potential ramifications on physical properties in the simulated system if the bilayer is initiated in a high-strain state. We present a new tool for creating curved lipid bilayers that combines flexibility of shape, force field, model resolution, and bilayer composition. A key aspect of our approach is the use of the monolayer pivotal plane location to accurately estimate interleaflet area differences in a curved bilayer. Our tool is named BUMPy (Building Unique Membranes in Python), is written in Python, is fast, and has a simple command line interface.

Stability of Norwalk Virus Capsid Protein Interfaces Evaluated by in Silico Nanoindentation.

Norwalk virus causes severe gastroenteritis for which there is currently no specific anti-viral therapy. A stage of the infection process is uncoating of the protein capsid to expose the viral genome and allow for viral replication. A mechanical characterization of the Norwalk virus may provide important information relating to the mechanism of uncoating. The mechanical strength of the Norwalk virus has previously been investigated using atomic force microscopy (AFM) nanoindentation experiments. Those experiments cannot resolve specific molecular interactions, and therefore, we have employed a molecular modeling approach to gain insights into the potential uncoating mechanism of the Norwalk capsid. In this study, we perform simulated nanoindentation using a coarse-grained structure-based model, which provides an estimate of the spring constant in good agreement with the experimentally determined value. We further analyze the fracture mechanisms and determine weak interfaces in the capsid structure, which are potential sites to inhibit uncoating by stabilization of these weak interfaces. We conclude by identifying potential target sites at the junction of a weak protein-protein interface.

The observation of the first vibrational overtone of dihydrogen in the luminescence of zeolites at low temperatures.

Dehydrated zeolites NaA, NaY, and barium-exchanged NaY luminesce when irradiated with the 1064 nm laser light of an FT-Raman instrument. When hydrogen is adsorbed in the zeolite, the luminescence is altered in several ways. Most remarkable is the appearance of "absorptions" in the positions of the first vibrational overtones of H(2) and HD. Although these features are in the expected positions with reasonable band profiles for overtone absorptions, the large extinction of the luminescent intensity requires a more efficient mechanism than simple reabsorption of emitted photons. In addition to the appearance of holes in the luminescence spectrum, other luminescent features are substantially quenched by the presence of hydrogen or, in one case, augmented.