No Arabic abstract
This paper has been withdrawn.
Many biological systems fold thin sheets of lipid membrane into complex three-dimensional structures. This microscopic origami is often mediated by the adsorption and self-assembly of proteins on a membrane. As a model system to study adsorption-mediated interactions, we study the collective behavior of micrometric particles adhered to a lipid vesicle. We estimate the colloidal interactions using a maximum likelihood analysis of particle trajectories. When the particles are highly wrapped by a tense membrane, we observe strong long-range attractions with a typical binding energy of 150 $k_B T$ and significant forces extending a few microns.
Typical biomolecular systems such as cellular membranes, DNA, and protein complexes are highly charged. Thus, efficient and accurate treatment of electrostatic interactions is of great importance in computational modelling of such systems. We have employed the GROMACS simulation package to perform extensive benchmarking of different commonly used electrostatic schemes on a range of computer architectures (Pentium-4, IBM Power 4, and Apple/IBM G5) for single processor and parallel performance up to 8 nodes - we have also tested the scalability on four different networks, namely Infiniband, GigaBit Ethernet, Fast Ethernet, and nearly uniform memory architecture, i.e., communication between CPUs is possible by directly reading from or writing to other CPUs local memory. It turns out that the particle-mesh Ewald method (PME) performs surprisingly well and offers competitive performance unless parallel runs on PC hardware with older network infrastructure are needed. Lipid bilayers of sizes 128, 512 and 2048 lipid molecules were used as the test systems representing typical cases encountered in biomolecular simulations. Our results enable an accurate prediction of computational speed on most current computing systems, both for serial and parallel runs. These results should be helpful in, for example, choosing the most suitable configuration for a small departmental computer cluster.
We use a coarse-grained molecular dynamics simulation to investigate the interaction between neutral or charged nanoparticles (NPs) and a vesicle consisting of neutral and negatively charged lipids. We focus on the interaction strengths of hydrophilic and hydrophobic attraction and electrostatic interactions between a lipid molecule and an NP. A neutral NP passes through the lipid membrane when the hydrophobic interaction is sufficiently strong. As the valence of the positively charged NP increases, the membrane permeation speed of the NP is increased compared with the neutral NP and charged lipids are accumulated around the charged NP. A charged NP with a high valence passes through the lipid membrane via a transient channel formed by charged lipids or transport-like endocytosis. These permeation processes can be classified based on analyses of the density correlation function. When the non-electrostatic interaction parameters are large enough, a negatively charged NP can be adsorbed on the membrane and a neutral lipid-rich region is formed directly below the NP. The NP is spontaneously incorporated into the vesicle under various conditions and the incorporation is mediated by the membrane curvature. We reveal how the NPs behavior depends on the NP valence, size, and the non-electrostatic interaction parameters.
We present two models for electron-driven uphill proton transport across lipid membranes, with the electron energy converted to the proton gradient via the electrostatic interaction. In the first model, associated with the cytochrome c oxidase complex in the inner mitochondria membranes, the electrostatic coupling to the site occupied by an electron lowers the energy level of the proton-binding site, making the proton transfer possible. In the second model, roughly describing the redox loop in a nitrate respiration of E. coli bacteria, an electron displaces a proton from the negative side of the membrane to a shuttle, which subsequently diffuses across the membrane and unloads the proton to its positive side. We show that both models can be described by the same approach, which can be significantly simplified if the system is separated into several clusters, with strong Coulomb interaction inside each cluster and weak transfer couplings between them. We derive and solve the equations of motion for the electron and proton creation/annihilation operators, taking into account the appropriate Coulomb terms, tunnel couplings, and the interaction with the environment. For the second model, these equations of motion are solved jointly with a Langevin-type equation for the shuttle position. We obtain expressions for the electron and proton currents and determine their dependence on the electron and proton voltage build-ups, on-site charging energies, reorganization energies, temperature, and other system parameters. We show that the quantum yield in our models can be up to 100% and the power-conversion efficiency can reach 35%.
We report a high energy-resolution neutron backscattering study to investigate slow motions on nanosecond time scales in highly oriented solid supported phospholipid bilayers of the model system DMPC -d54 (deuterated 1,2-dimyristoyl-sn-glycero-3-phoshatidylcholine), hydrated with heavy water. This technique allows to discriminate the onset of mobility at different length scales for the different molecular components, as e.g.@ the lipid acyl-chains and the hydration water in between the membrane stacks, respectively, and provides a benchmark test regarding the feasibility of neutron backscattering investigations on these sample systems. We discuss freezing of the lipid acyl-chains, as observed by this technique, and observe a second freezing transition which we attribute to the hydration water.