Computer modelling studies of the bilayer/water interface.

This review summarises high resolution studies on the interface of lamellar lipid bilayers composed of the most typical lipid molecules which constitute the lipid matrix of biomembranes. The presented results were obtained predominantly by computer modelling methods. Whenever possible, the results were compared with experimental results obtained for similar systems. The first and main section of the review is concerned with the bilayer-water interface and is divided into four subsections. The first describes the simplest case, where the interface consists only of lipid head groups and water molecules and focuses on interactions between the lipid heads and water molecules; the second describes the interface containing also mono- and divalent ions and concentrates on lipid-ion interactions; the third describes direct inter-lipid interactions. These three subsections are followed by a discussion on the network of direct and indirect inter-lipid interactions at the bilayer interface. The second section summarises recent computer simulation studies on the interactions of antibacterial membrane active compounds with various models of the bacterial outer membrane. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.

Communication: consistent picture of lateral subdiffusion in lipid bilayers: molecular dynamics simulation and exact results.

This communication presents a molecular dynamics simulation study of a bilayer consisting of 128 dioleoyl-sn-glycero-3-phosphocholine molecules, which focusses on the center-of-mass diffusion of the lipid molecules parallel to the membrane plane. The analysis of the simulation results is performed within the framework of the generalized Langevin equation and leads to a consistent picture of subdiffusion. The mean square displacement of the lipid molecules evolves as ∝ t(α), with α between 0.5 and 0.6, and the fractional diffusion coefficient is close to the experimental value for a similar system obtained by fluorescence correlation spectroscopy. We show that the long-time tails of the lateral velocity autocorrelation function and the associated memory function agree well with exact results which have been recently derived by asymptotic analysis [G. Kneller, J. Chem. Phys. 134, 224106 (2011)]. In this context, we define characteristic time scales for these two quantities.

Structural and functional characterization of SplA, an exclusively specific protease of Staphylococcus aureus.

Staphylococcus aureus is a dangerous human pathogen whose antibiotic resistance is steadily increasing and no efficient vaccine is as yet available. This serious threat drives extensive studies on staphylococcal physiology and pathogenicity pathways, especially virulence factors. Spl (serine protease-like) proteins encoded by an operon containing up to six genes are a good example of poorly characterized secreted proteins probably involved in virulence. In the present study, we describe an efficient heterologous expression system for SplA and detailed biochemical and structural characterization of the recombinant SplA protease. The enzyme shares a significant sequence homology to V8 protease and epidermolytic toxins which are well documented staphylococcal virulence factors. SplA has a very narrow substrate specificity apparently imposed by the precise recognition of three amino acid residues positioned N-terminal to the hydrolysed peptide bond. To explain determinants of this extended specificity we resolve the crystal structure of SplA and define the consensus model of substrate binding. Furthermore we demonstrate that artificial N-terminal elongation of mature SplA mimicking a naturally present signal peptide abolishes enzymatic activity. The probable physiological role of the process is discussed. Of interest, even though precise N-terminal trimming is a common regulatory mechanism among S1 family enzymes, the crystal structure of SplA reveals novel significantly different mechanistic details.

Electron paramagnetic resonance spectroscopy reveals alterations in the redox state of endogenous copper and iron complexes in photodynamic stress-induced ischemic mouse liver.

Divalent copper and iron cations have been acknowledged for their catalytic roles in physiological processes critical for homeostasis maintenance. Being redox-active, these metals act as cofactors in the enzymatic reactions of electron transfer. However, under pathophysiological conditions, owing to their high redox potentials, they may exacerbate stress-induced injury. This could be particularly hazardous to the liver - the main body reservoir of these two metals. Surprisingly, the involvement of Cu and Fe in liver pathology still remains poorly understood. Hypoxic stress in the tissue may act as a stimulus that mobilizes these ions from their hepatic stores, aggravating the systemic injury. Since ischemia poses a serious complication in liver surgery (e.g. transplantation) we aimed to reveal the status of Cu and Fe via spectroscopic analysis of mouse ischemic liver tissue. Herein, we establish a novel non-surgical model of focal liver ischemia, achieved by applying light locally when a photosensitizer is administered systemically. Photodynamic treatment results in clear-cut areas of the ischemic hepatic tissue, as confirmed by ultrasound scans, mean velocity measurements, 3D modelling of vasculature and (immuno)histological analysis. For reference, we assessed the samples collected from the animals which developed transient systemic endotoxemic stress induced by a non-lethal dose of lipopolysaccharide. The electron paramagnetic resonance (EPR) spectra recorded in situ in the liver samples reveal a dramatic increase in the level of Cu adducts solely in the ischemic tissues. In contrast, other typical free radical components of the liver EPR spectra, such as reduced Riske clusters are not detected; these differences are not followed by changes in the blood EPR spectra. Taken together, our results suggest that local ischemic stress affects paramagnetic species containing redox-active metals. Moreover, because in our model hepatic vascular flow is impaired, these effects are only local (confined to the liver) and are not propagated systemically.

Is the tilt of the lipid head group correlated with the number of intermolecular interactions at the bilayer interface?

Lipid and water molecules comprising the bilayer form an integral entity owing to only weak physical interactions. At the bilayer interface, these interactions chiefly involve hydrogen bonding and charge pairing. Lipid head groups make hydrogen bonds (H-bonds) predominantly with water, whereas interlipid H-bonds and charge pairs are less numerous. Both interlipid H-bonding and charge pairing depend on the distance and relative orientation of the interacting head groups. In this computational paper, correlations are analysed between the orientation of the lipid head group and the number of interlipid interactions at the interface of a bilayer made of galactolipids, forming direct interlipid H-bonds, and of phosphatidylcholines forming interlipid charge pairs. The correlations are not strong, however, in both bilayers they show a similar trend.

Stretching of buckled filaments by thermal fluctuations.

We study the buckling instability of filaments or elastic rods in two spatial dimensions in the presence of thermal fluctuations. We present an analytical solution based on a renormalization-like procedure where we integrate out short wavelength fluctuations in order to obtain an effective theory governing the buckling instability. We calculate the resulting shift of the critical force by fluctuation effects and the average projected filament length parallel to the force direction as a function of the applied force and of the contour length of the filament. We find that, in the buckled state, thermal fluctuations lead to an increase in the mean projected length of the filament in the force direction. As a function of the contour length, the mean projected length exhibits a cusp at the buckling instability, which becomes rounded by thermal fluctuations. Our analytic results are confirmed by Monte Carlo simulations.

Assessing gastric toxicity of xanthone derivatives of anti-inflammatory activity using simulation and experimental approaches.

Xanthones are tricyclic compounds of natural or synthetic origin exhibiting a broad spectrum of therapeutic activities. Three synthetic xanthone derivatives (KS1, KS2, and KS3) with properties typical for nonsteroidal anti-inflammatory drugs (NSAID) were objects of the presented model study. NSAIDs are in common use however; several of them exhibit gastric toxicity predominantly resulting from their direct interactions with the outermost lipid layer of the gastric mucosa that impair its hydrophobic barrier property. Among the studied xanthones, gastric toxicity of only KS2 has been determined in previous pharmacological studies, and it is low. In this study, carried out using X-ray diffraction and computer simulation, a palmitoyloleoylphosphatidylcholine-cholesterol bilayer (POPC-Chol) was used as a model of a hydrophobic layer of lipids protecting gastric mucosa as POPC and Chol are the main lipids in human mucus. X-ray diffraction data were used to validate the computer model. The aim of the study was to assess potential gastric toxicity of the xanthones by analysing their atomic level interactions with lipids, ions, and water in the lipid bilayer and their effect on the bilayer physicochemical properties. The results show that xanthones have small effect on the bilayer properties except for its rigidity whereas their interactions with water, ions, and lipids depend on their protonation state and for a given state, are similar for all the xanthones. As gastric toxicity of KS2 is low, based on MD simulations one can predict that toxicity of KS1 and KS3 is also low.

A computer model of a polyunsaturated monogalactolipid bilayer.

1,2-di-O-acyl-3-O-ß-D-galactopyranosyl-sn-glycerol (MGDG) is the main lipid component of thylakoid membranes of higher plants and algae. This monogalactolipid is thought of as a non-bilayer lipid but actually it can form both lamellar and nonlamellar phases. In this study, molecular dynamics (MD) simulations of the fully hydrated di-18:3 MGDG bilayer in the lamellar phase were carried out at 310 and 295 K for 200 and 450 ns, respectively, using the GROMACS 4 software package and OPLS-AA force field. At both temperatures, the lamellar phase of the systems was stable. The pure di-18:3 MGDG bilayer is the first step towards creating a computer model of the lipid matrix of the thylakoid membrane and the main aim of this study was to validate the computer model of di-18:3 MGDG in the bilayer and also to assess the properties of the bilayer. However, only a few of the predicted properties could be compared with those derived experimentally and in other MD simulations because of insufficient amount of such data. Thus, direct validation of the MGDG bilayer proved difficult. Therefore, in the validation process also an indirect approach was used, in which a computer model of the 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayer simulated at the same temperatures using the same force field as the MGDG bilayer was assessed. Successful validation of the DOPC bilayer parameterized in the OPLS-AA force field and similar properties of the MGDG molecules in the pure 18:3 MGDG and in binary 18:3 MGDG-PC bilayers indicate that the computer model of the MDGD molecule is faithful and the MGDG bilayer is representative on the time scales covered in these MD simulations.

Properties of water hydrating the galactolipid and phospholipid bilayers: a molecular dynamics simulation study.

Molecular dynamics simulations of 1,2-di-O-acyl-3-O-ß-D-galactopyranosyl-sn-glycerol (MGDG) and 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) bilayers were carried out to compare the effect of the lipid head group's chemical structure on the dynamics and orientational order of the water molecules hydrating the bilayer. The effect of the bilayers on the diffusion of water is strong for the neighbouring water molecules i.e., those located not further than 4 Å from any bilayer atom. This is because the neighbouring water molecules are predominantly hydrogen bonded to the lipid oxygen atoms and their mobility is limited to a confined spatial volume. The choline group of DOPC and the galactose group of MGDG affect water diffusion less than the polar groups located deeper in the bilayer interface, and similarly. The latter is an unexpected result since interactions of water with these groups have a vastly different origin. The least affected by the bilayer lipids is the lateral diffusion of unbound water in the bilayer plane (x,y-plane) - it is because the diffusion is not confined by the periodic boundary conditions, whereas that perpendicular to the plane is. Interactions of water molecules with lipid groups also enforce certain orientations of water dipole moments. The profile of an average water orientation along the bilayer normal for the MGDG bilayer differs from that for the DOPC bilayer. In the DOPC bilayer, the ordering effect of the lipid head groups extends further into the water phase than in the MGDG bilayer, whereas inside the bilayer/water interface, ordering of the water dipoles in the MGDG bilayer is higher. It is possible that differences in the profiles of an average water orientation across the bilayer in the DOPC and MGDG bilayers are responsible for differences in the lateral pressure profiles of these bilayers.

Unifying autocatalytic and zeroth-order branching models for growing actin networks.

The directed polymerization of actin networks is an essential element of many biological processes, including cell migration. Different theoretical models considering the interplay between the underlying processes of polymerization, capping, and branching have resulted in conflicting predictions. One of the main reasons for this discrepancy is the assumption of a branching reaction that is either first order (autocatalytic) or zeroth order in the number of existing filaments. Here we introduce a unifying framework from which the two established scenarios emerge as limiting cases for low and high filament numbers. A smooth transition between the two cases is found at intermediate conditions. We also derive a threshold for the capping rate above which autocatalytic growth is predicted at sufficiently low filament number. Below the threshold, zeroth-order characteristics are predicted to dominate the dynamics of the network for all accessible filament numbers. Together, these mechanisms allow cells to grow stable actin networks over a large range of different conditions.

Orientation of lutein in a lipid bilayer - revisited.

Lutein is present in the human retina and lens, where it plays a protective role. As lutein is associated with the lipid matrix of biomembranes, the role depends on its membrane location. Experimental studies predicted two orientations of lutein in a phosphatidylcholine (PC) bilayer: vertical and horizontal. Using a molecular dynamics simulation, we observed, in two different PC bilayers, both orientations of lutein, and in each bilayer, a single change from vertical to horizontal orientation or vice versa. Both orientations were stabilized by hydrogen bonding of lutein OH groups with mainly carbonyl but also phosphate oxygen atoms of PC.

Temperature- and hydration-dependent protein dynamics in photosystem II of green plants studied by quasielastic neutron scattering.

Protein dynamics in hydrated and vacuum-dried photosystem II (PS II) membrane fragments from spinach has been investigated by quasielastic neutron scattering (QENS) in the temperature range between 5 and 300 K. Three distinct temperature ranges can be clearly distinguished by active type(s) of protein dynamics: (A) At low temperatures (T < 120 K), the protein dynamics of both dry and hydrated PS II is characterized by harmonic vibrational motions. (B) In the intermediate temperature range (120 < T < 240 K), the total mean square displacement total slightly deviates from the predicted linear behavior. The QENS data indicate that this deviation, which is virtually independent of the extent of hydration, is due to a partial onset of diffusive protein motions. (C) At temperatures above 240 K, the protein flexibility drastically changes because of the onset of diffusive (large-amplitude) protein motions. This dynamical transition is clearly hydration-dependent since it is strongly suppressed in dry PS II. The thermally activated onset of protein flexibility as monitored by QENS is found to be strictly correlated with the temperature-dependent increase of the electron transport efficiency from Q(A)(-) to QB (Garbers et al. (1998) Biochemistry 37, 11399-11404). Analogously, the freezing of protein mobility by dehydration in dry PS II appears to be responsible for the blockage of Q(A)(-) reoxidation by Q(B) at hydration values lower than 45% r.h. (Kaminskaya et al. (2003) Biochemistry 42, 8119-8132). Similar effects were observed for reactions of the water-oxidizing complex as outlined in the Discussion section.