We describe a new approach to calculate the binding of flexible peptides and unfolded proteins to multicomponent lipid membranes. The method is based on the transfer matrix formalism of statistical mechanics recently described as a systematic tool to study DNA-protein-drug binding in gene regulation. Using the energies of interaction of the individual polymer segments with different membrane lipid species and the scaling corrections due to polymer looping, we calculate polymer adsorption characteristics and the degree of sequestration of specific membrane lipids. The method is applied to the effector domain of the MARCKS (myristoylated alanine rich C kinase substrate) protein known to be involved in signal transduction through membrane binding. The calculated binding constants of the MARCKS(151-175) peptide and a series of related peptides to mixed PC/PS/PIP2 membranes are in satisfactory agreement with in vitro experiments. Copyright (c) 2008 European Peptide Society and John Wiley & Sons, Ltd.

Amphipathic alpha-helical peptides are often ascribed an ability to induce curvature stress in lipid membranes. This may lead directly to a bending deformation of the host membrane, or it may promote the formation of defects that involve highly curved lipid layers present in membrane pores, fusion intermediates, and solubilized peptide-micelle complexes. The driving force is the same in all cases: peptides induce a spontaneous curvature in the host lipid layer, the sign of which depends sensitively on the peptide’s structural properties. We provide a quantitative account for this observation on the basis of a molecular-level method. To this end, we consider a lipid membrane with peptides interfacially adsorbed onto one leaflet at high peptide-to-lipid ratio. The peptides are modeled generically as rigid cylinders that interact with the host membrane through a perturbation of the conformational properties of the lipid chains. Through the use of a molecular-level chain packing theory, we calculate the elastic properties, that is, the spontaneous curvature and bending stiffness, of the peptide-decorated lipid membrane as a function of the peptide’s insertion depth. We find a positive spontaneous curvature (preferred bending of the membrane away from the peptide) for small penetration depths of the peptide. At a penetration depth roughly equal to half-insertion into the hydrocarbon core, the spontaneous curvature changes sign, implying negative spontaneous curvature (preferred bending of the membrane toward the peptide) for large penetration depths. Despite thinning of the membrane upon peptide insertion, we find an increase in the bending stiffness. We discuss these findings in terms of how the peptide induces elastic stress.

We present a theory of the dependence on sequence of the three-dimensional size of large single-stranded (ss) RNA molecules. The work is motivated by the fact that the genomes of many viruses are large ssRNA molecules-often several thousand nucleotides long-and that these RNAs are spontaneously packaged into small rigid protein shells. We argue that there has been evolutionary pressure for the genome to have overall spatial properties-including an appropriate radius of gyration, R-g-that facilitate this assembly process. For an arbitrary RNA sequence, we introduce the (thermal) average maximum ladder distance (< MLD >) and use it as a measure of the ‘‘extendedness’’ of the RNA secondary structure. The < MLD > values of viral ssRNAs that package into capsids of fixed size are shown to be consistently smaller than those for randomly permuted sequences of the same length and base composition, and also smaller than those of natural ssRNAs that are not under evolutionary pressure to have a compact native form. By mapping these secondary structures onto a linear polymer model and by using < MLD > as a measure of effective contour length, we predict the R-g values of viral ssRNAs are smaller than those of nonviral sequences. More generally, we predict the average < MLD > values of large nonviral ssRNAs scale as N-0.67 +/- 0.01, where N is the number of nucleotides, and that their R-g values vary as < MLD >(0.5) in an ideal solvent, and hence as N-0.34. An alternative analysis, which explicitly includes all branches, is introduced and shown to yield consistent results.

The cellular cytoskeleton is a complex dynamical network that constantly remodels as cells divide and move. This reorganization process occurs not only at the cell membrane, but also in the cell interior (bulk). During locomotion, regulated actin assembly near the plasma membrane produces lamellipodia and filopodia. Therefore, most in vitro experiments explore phenomena taking place in the vicinity of a surface. To understand how the molecular machinery of a cell self-organizes in a more general way, we studied bulk polymerization of actin in the presence of actin-related protein 2/3 complex and a nucleation promoting factor as a model for actin assembly in the cell interior separate from membranes. Bulk polymerization of actin in the presence of the verprolin homology, cofilin homology, and acidic region, domain of Wiskott-Aldrich syndrome protein, and actin-related protein 2/3 complex results in spontaneous formation of diffuse aster-like structures. In the presence of fascin these asters transition into stars with bundles of actin filaments growing from the surface, similar to star-like structures recently observed in vivo. The transition from asters to stars depends on the ratio [fascin]/[G actin]. The polarity of the actin filaments during the transition is preserved, as in the transition from lamellipodia to filopodia. Capping protein inhibits star formation. Based on these experiments and kinetic Monte Carlo simulations, we propose a model for the spontaneous self-assembly of asters and their transition into stars. This mechanism may apply to the transition from lamellipodia to filopodia in vivo.

Peripheral proteins can trigger the formation of domains in mixed. uid- like lipid membranes. We analyze the mechanism underlying this process for proteins that bind electrostatically onto a. at two- component membrane, composed of charged and neutral lipid species. Of particular interest are membranes in whichthe hydrocarbon lipid tails tend to segregate owing to nonideal chain mixing,but the ( protein- free) lipid membrane is nevertheless stable due to the electrostatic repulsion between the charged lipid headgroups. The adsorption of charged, say basic, proteins onto a membrane containing anionic lipids induces local lipid demixing, whereby charged lipids migrate toward ( or away from) the adsorption site, so as to minimize the electrostatic binding free energy. Apart from reducing lipid headgroup repulsion, this process creates a gradient in lipid composition around the adsorption zone, and hence a line energy whose magnitude depends on the protein’s size and charge and the extent of lipid chain nonideality. Above a certain critical lipid nonideality, the line energy islarge enough to induce domain formation, i. e., protein aggregation and, concomitantly, macroscopic lipid phase separation. We quantitatively analyze the thermodynamic stability of the dressed membrane based on nonlinear Poisson- Boltzmann theory, accounting for both the microscopic characteristics of the proteins and lipid composition modulations at and around the adsorption zone. Spinodal surfaces and critical points of the dressed membranes are calculated for several different model proteins of spherical and disklike shapes. Among the models studied we. nd the most substantial protein- induced membrane destabilization for disk- like proteins whose charges are concentrated in the membrane- facing surface. If additional charges reside on the side faces of the proteins, direct protein- protein repulsion diminishes considerably the propensity fordomain formation. Generally, a highly charged. at face of a macroion appears most ef. cient in inducing large compositional gradients, hence a large and unfavorable line energy and consequently lateral macroion aggregation and, concomitantly, macroscopic lipid phase separation.

Fluid membranes containing charged lipids enhance binding of oppositely charged proteins by mobilizing these lipids into the interaction zone, overcoming the concomitant entropic losses due to lipid segregation and lower conformational freedom upon macromolecule adsorption. We study this energetic-entropic interplay using Monte Carlo simulations and theory. Our model system consists of a flexible cationic polyelectrolyte, interacting, via Debye-Huckel and short-ranged repulsive potentials, with membranes containing neutral lipids, 1% tetravalent, and 10% ( or 1%) monovalent anionic lipids. Adsorption onto a fluid membrane is invariably stronger than to an equally charged frozen or uniform membrane. Although monovalent lipids may suffice for binding rigid macromolecules, polyvalent counter-lipids ( e. g., phosphatidylinositol 4,5 bisphosphate), whose entropy loss upon localization is negligible, are crucial for binding flexible macromolecules, which lose conformational entropy upon adsorption. Extending Rosenbluth’s Monte Carlo scheme we directly simulate polymer adsorption on fluid membranes. Yet, we argue that similar information could be derived from a biased superposition of quenched membrane simulations. Using a simple cell model we account for surface concentration effects, and show that the average adsorption probabilities on annealed and quenched membranes coincide at vanishing surface concentrations. We discuss the relevance of our model to the electrostatic-switch mechanism of, e. g., the myristoylated alanine-rich C kinase substrate protein.